|
Revision 943, 1.4 kB
(checked in by smidl, 14 years ago)
|
|
syntax of epredictor
|
| Line | |
|---|
| 1 | /*! |
|---|
| 2 | \file |
|---|
| 3 | \brief wrapper function for epdf.marginal() |
|---|
| 4 | |
|---|
| 5 | |
|---|
| 6 | */ |
|---|
| 7 | |
|---|
| 8 | #include <stat/emix.h> |
|---|
| 9 | |
|---|
| 10 | using namespace bdm; |
|---|
| 11 | |
|---|
| 12 | #ifdef MEX |
|---|
| 13 | #include <mex/mex_parser.h> |
|---|
| 14 | #include <mex/mex_pdf.h> |
|---|
| 15 | |
|---|
| 16 | void mexFunction ( int n_output, mxArray *output[], int n_input, const mxArray *input[] ) { |
|---|
| 17 | // Check the number of inputs and output arguments |
|---|
| 18 | if ( n_input<1 ) mexErrMsgTxt ( "Usage:\n" |
|---|
| 19 | "vec=pdf_marginal(pdf_struct, RVs)\n" |
|---|
| 20 | " pdf_struct = struct('class','pdf_offspring',...); % description of pdf\n" |
|---|
| 21 | " RVs = RV(...); % description of variables of teh marginal" |
|---|
| 22 | "output:\n" |
|---|
| 23 | " pdf_struct = struct('class','marginal_type',...); % description of marginal pdf." ); |
|---|
| 24 | |
|---|
| 25 | RV::clear_all(); |
|---|
| 26 | //CONFIG |
|---|
| 27 | if (n_input<2){mexErrMsgTxt("RVs not given"); } |
|---|
| 28 | |
|---|
| 29 | UImxArray Cfg(input[0]); |
|---|
| 30 | UImxArray RVs(input[1]); |
|---|
| 31 | //Cfg.writeFile("epdf_marginal.cfg"); |
|---|
| 32 | |
|---|
| 33 | shared_ptr<epdf> ep=UI::build<epdf>(Cfg); |
|---|
| 34 | |
|---|
| 35 | if (ep){ |
|---|
| 36 | //if ( n_output<1 ) mexErrMsgTxt ( "No output - nothing to do!" ); |
|---|
| 37 | shared_ptr<RV> rv=UI::build<RV>(RVs); |
|---|
| 38 | shared_ptr<epdf> marginal = ep->marginal(*rv); |
|---|
| 39 | UImxArray Out; |
|---|
| 40 | UI::save(marginal, Out); |
|---|
| 41 | output[0]=UImxArray::create_mxArray(Out); |
|---|
| 42 | } else { |
|---|
| 43 | mexErrMsgTxt ( "Given object is not epdf" ); |
|---|
| 44 | } |
|---|
| 45 | } |
|---|
| 46 | #endif |
|---|
