probability_and_matlab_extensions.h

#ifndef PROBABILITY_AND_MATLAB_EXTENSIONS_H
#define PROBABILITY_AND_MATLAB_EXTENSIONS_H

#include "matlab.h"
#include "probability.h"


void plotDensity(const Density& density,const char* densityName = NULL,const char* x1_name="X",const char* x2_name="Y"){

    MatlabHandle matlab;

    matlab.pass(x1_name,"x1_name");
    matlab.pass(x2_name,"x2_name");
    matlab.pass(densityName,"densityName");

    matlab.pass<float>(density.getExpectedN(),"expectedN");
    matlab.pass<float>(density.getExpectedM(),"expectedM");
    matlab.pass<float>(density.getExpectedMN(),"expectedMN");
    matlab.pass<float>(density.getExpectedNN(),"expectedNN");
    matlab.pass<float>(density.getExpectedMM(),"expectedMM");

    matlab.pass<int>(density.getN(),"N");
    matlab.pass<int>(density.getM(),"M");
    matlab.pass<int>(density.getnCenter(),"nCenter");
    matlab.pass<int>(density.getmCenter(),"mCenter");
    matlab.pass<float>(density.getTotalWeight(),"totalWeight");

    matlab.pass<float>(density.getColPtr(),density.getM(),"Mdensity");
    matlab.pass<float>(density.getRowPtr(),density.getN(),"Ndensity");
    matlab.pass<float>(density.getMatrixPtr(),density.getM(),density.getN(),"MNdensity");

    matlab.pass<float>(density.getExpectedNgivenM_ptr(),density.getM(),"expectedNgivenM");
    matlab.pass<float>(density.getExpectedNNgivenM_ptr(),density.getM(),"expectedNNgivenM");
    matlab.pass<float>(density.getExpectedMgivenN_ptr(),density.getN(),"expectedMgivenN");
    matlab.pass<float>(density.getExpectedMMgivenN_ptr(),density.getN(),"expectedMMgivenN");

    matlab.loadScriptandSend("displayDensity.m");

}

#endif