example_mathmatrix_kldivergence.cpp

Programa para el testeo de las funciones.

Programa para el testeo de las funciones.Para compilar o código example_mathmatrix_kldivergence.cpp:

g++ -static -o example_mathmatrix_kldivergence example_mathmatrix_kldivergence.cpp -lpdsramm 

Para executar o programa:

./example_mathmatrix_kldivergence

Retornando por consola:

Pds::KLDivergence(0.0,0.0):0
Pds::KLDivergence(0.0,1.0):inf
Pds::KLDivergence(1.0,0.0):inf
Pds::KLDivergence(1.0,1.0):0

example_mathmatrix_kldivergence.png

example_mathmatrix_kldivergence_cases.png

Código example_mathmatrix_kldivergence.cpp:

 
#include <Pds/Ra>
    
int main(void)
{
    unsigned int N=128;
    Pds::Vector p=Pds::LinSpace(0.0,1,N); 
    Pds::Vector q=Pds::LinSpace(0.0,1,N);
    Pds::Matrix P,Q;
    
    Pds::Meshgrid (p,q,P,Q);
    Pds::Matrix Z=Pds::KLDivergence(P,Q);
    
    Pds::Octave::XLabel="P";
    Pds::Octave::YLabel="Q";
    Pds::Octave::Plot::Surf(P,Q,Z,"testando.m","example_mathmatrix_kldivergence.png");
    
    std::cout<<"Pds::KLDivergence(0.0,0.0):"<<Pds::KLDivergence(0.0,0.0)<<"\n";
    std::cout<<"Pds::KLDivergence(0.0,1.0):"<<Pds::KLDivergence(0.0,1.0)<<"\n";
    std::cout<<"Pds::KLDivergence(1.0,0.0):"<<Pds::KLDivergence(1.0,0.0)<<"\n";
    std::cout<<"Pds::KLDivergence(1.0,1.0):"<<Pds::KLDivergence(1.0,1.0)<<"\n";
    
    
    Pds::Octave::XLabel="q";
    Pds::Octave::YLabel="KL";
    Pds::Octave::Legend={"p=0","p=1"};
    Pds::Octave::Plot::CurveXYXY(   q,Z.GetRowAsColVector(0),
                                    q,Z.GetRowAsColVector(N-1),
                                    "testando.m",
                                    "example_mathmatrix_kldivergence_cases.png");
    return 0;
}