Context Navigation

bm_bayes.cpp

Timestamp:

05/27/10 13:08:31 (14 years ago)

Author:

smidl

Message:

new mexes + cleanup

Files:

r978	r1002
14	14	#include <estim/kalman.h>
15	15	#include <estim/particles.h>
	16	#include <estim/mixtures.h>
16	17	#include <mex/mex_BM.h>
17	18
18	19	void mexFunction ( int n_output, mxArray output[], int n_input, const mxArray input[] ) {
19	20	// Check the number of inputs and output arguments
20		if ( n_input<1 ) mexErrMsgTxt ( "Usage:\n"
21		"bm_struct=bm_bayes(BM_struct, yt, cond)\n"
22		" bm_struct = struct('class','BM_offspring',...); % description of Bayesian Model (BM)\n"
23		" yt = [...]; % vector of observed data\n"
24		" cond = [...]; % vector of conditioning variables\n"
25		"output:\n"
26		" bm_struct = struct('class','BM_offsprin',...); % description of BM afret update." );
	21	if ( n_input<1 ) mexErrMsgTxt (
	22	"Usage:\n"
	23	"[bm_struct,logevidence]=bm_bayes(BM_struct, yt, cond)\n"
	24	" BM_struct = struct('class','BM_offspring',...); % description of Bayesian Model (BM)\n"
	25	" yt = [...]; % vector of observed data\n"
	26	" cond = [...]; % vector of conditioning variables\n"
	27	"output:\n"
	28	" bm_struct = struct('class','BM_offsprin',...); % description of BM afret update.\n"
	29	" logevidence = 0.0; % logarithm of normalizing constant f(yt\|cond)"
	30	);
27	31
28	32	RV::clear_all();
…	…
50	54	UI::save(bm.get(), Out);
51	55	output[0]=UImxArray::create_mxArray(Out);
	56	if (n_output>1){
	57	mxArray *le=mxCreateDoubleScalar(mxREAL);
	58	double2mxArray(bm->_ll(), le);
	59	output[1]=le;
	60	}
52	61	} else {
53	62	mexErrMsgTxt ( "Given object is not BM" );