Changeset 374 for library/mex/mexBM.cpp

Timestamp:

06/11/09 17:15:06 (15 years ago)

Author:

miro

Message:

mexBM works. An example added (it ignores any underlying theory).
config2mxstruct probably doesn't handle vectors and matrices, to be fixed

Files:

• library/mex/mexBM.cpp (modified) (6 diffs)

library/mex/mexBM.cpp

@@ -1 +1 @@
 #include <itpp/itmex.h>
 #include <estim/arx.h>
+#include "config2mxstruct.h"
 
 using namespace bdm;
@@ -12 +13 @@
                 void from_setting(const Setting &S)  {
                         name = (const char *) S["name"];
-                        string fname = name+"_new";
-                        mexCallMATLAB(1, &data, 0, 0, fname.c_str());
+                        UImxConfig conf(S);
+                        data = mxDuplicateArray(conf.mxconfig);
+                        //mexCallMATLAB(1, &data, 0, 0, (name+"_new").c_str());
                         //TODO (future...):
                         //mxArray * init_data = setting2mxarray S["init_data"];
@@ -35 +37 @@
                 mexBM() {}
 
+                mxArray *get_data() {
+                        //mexCallMATLAB(0, NULL, 1, &data, "dump");
+                        return mxDuplicateArray(data);
+                }
+
                 void from_setting(const Setting &S)  {
                         name = (const char *) S["name"];
-                        string fname = name+"_new";
-                        mexCallMATLAB(1, &data, 0, 0, (name+"_new").c_str());
+                        UImxConfig conf(S);
+                        data = mxDuplicateArray(conf.mxconfig);
+                        //string fname = name+"_new";
+                        //mexCallMATLAB(1, &data, 0, 0, (name+"_new").c_str());
                         //the following works as long as the posterior is the
                         //only member object there could be a structure of 
@@ -48 +57 @@
                 void bayes(const vec &dt)  {
                 //void bayes()  {
-                        string fname = name+"_bayes";
-                        mexCallMATLAB(1, &data, 1, &data, fname.c_str());
+                        mxArray *tmp, *old;
+                        mxArray *in[2];
+                        in[0] = data;
+                        in[1] = mxCreateDoubleMatrix(dt.size(), 1, mxREAL);
+                        vec2mxArray(dt, in[1]);
+                        mexCallMATLAB(1, &tmp, 2, in, (name+"_bayes").c_str());
+                        old = data;
+                        data = mxDuplicateArray(tmp);
+                        if (old) mxDestroyArray(old);
+                        if (tmp) mxDestroyArray(tmp);
+                        //mexCallMATLAB(0, NULL, 1, &data, "dump");
                 } 
                 const mexEpdf& posterior() const  {
@@ -72 +90 @@
         config.readFile(filename.c_str());
         mexBM mb;
-        mb.from_setting(config.getRoot());
-        vec a = "1.0 2.0 3.0";
-        mb.bayes(a);
+        Setting &root = config.getRoot();
+        mb.from_setting(root);
+        vec dt = "18.0";
+        mb.bayes(dt);
+        if(n_output>0) output[0] = mb.get_data();
         
 
@@ -80 +100 @@
 
     // Create output vectors
-        output[0] = mxCreateDoubleMatrix(1,1, mxREAL);
+        //output[0] = mxCreateDoubleMatrix(1,1, mxREAL);
 
     // Convert the IT++ format to Matlab format for output