example_ceo_model_entropy_joint.cpp

 
#include <Pds/Ra>
#include <Pds/It>
    
int main(void)
{
    unsigned int M=5;
    unsigned int L=24;
    
    double ps;
    double rho;
    double jh,js;
    
    Pds::Vector Jh(L);
    Pds::Vector Js(L);
    Pds::Vector RhoI(L);
    
    Pds::Vector Rho=Pds::LinSpace(0,0.5,L);
    Pds::Vector p=Pds::Ones(M,1);
    
    ps=0.5;
    for(unsigned int l=0;l<L;l++)
    {
        jh=Pds::Ceo::Binary::Entropy::HjOmegaM(p*Rho[l],ps);
        js=Pds::Ceo::Binary::Entropy::HjsOmegaM(Rho[l],M,ps);
        rho=Pds::Ceo::Binary::Entropy::HjsOmegaMInv(js,M,ps);
        
        Jh.SetRaw(l,jh);
        Js.SetRaw(l,js);
        RhoI.SetRaw(l,rho);
    }
    
    Pds::Octave::XLabel="Rho";
    Pds::Octave::YLabel="Entropy Joint M="+std::to_string(M);
    Pds::Octave::Legend={"Jh","Js"};
    Pds::Octave::Plot::CurveXYXY(Rho,Jh,Rho,Js,"temporal.m","example_ceo_model_entropy_joint.png");
    
    Pds::Octave::XLabel="Rho";
    Pds::Octave::YLabel="Inv\\{H(Rho,M)\\}";
    Pds::Octave::Plot::CurveXY(Rho,RhoI,"temporal.m","example_ceo_model_entropy_joint_rho.png");
    return 0;
}