example_ceo_model.cpp

 
#include <Pds/Ra>
#include <Pds/It>
 
int main(void)
{
    unsigned int L=100000;
    Pds::Vector p_exp;
    Pds::Vector Us(L);  
    
    // Dados do canal
    double ps=0.5;
    Pds::Vector p={0.1,0.2,0.3,0.4};
    
    // Source Us
    Us.FillRandC(ps);
    
    // Input channel U
    Pds::Matrix U=Pds::ChannelModel::Bsc(Us,p);
    
    // Second fusion method: NaiveBayes 
    Pds::Vector Unb=Pds::Ceo::Binary::Fusion::NaiveBayes(ps,p,U);
    
    // Results
    p.T().Print("p.T():\t");
    std::cout<<"BERnb:\t"<<Us.Xor(Unb).Mean()<<"\n";
    
    return 0;
}