mex_datasource.h

#ifndef MXDS_H
#define MXDS_H


#include "../bdm/bdmerror.h"
#include "../bdm/base/datasources.h"
#include "mex_parser.h"

namespace bdm {
class mxArrayDS : public MemDS {
        public:
        mxArrayDS ():MemDS() {};

        void from_setting ( const Setting &set ) {
                Data = mxArray2mat ( mexGetVariable ( "base", set["varname"] ) );
/*              UI::get ( rowid, set, "rids" , UI::compulsory );
                bdm_assert_debug ( max ( rowid ) <= Data.rows(), "MemDS rowid is too high for given Dat." );

                UI::get ( delays, set, "tds", UI::compulsory );
                time = max ( delays );
                bdm_assert_debug ( time < Data.cols(), "MemDS delays are too high." );
*/
                //set MemDS
                rowid = linspace(0,Data.rows()-1);
                delays = zeros_i(Data.rows());
                dtsize=rowid.length();
                utsize=0;
                
                shared_ptr<RV> r = UI::build<RV> ( set, "rv", UI::optional );
                RV ru = RV();
                if (r){
                        set_rvs ( *r, ru );
                } else {
                        RV def((const char*)set["varname"],Data.rows());
                        set_rvs(def, ru);
                }
        }


        // TODO dodelat void to_setting( Setting &set ) const;
};

UIREGISTER ( mxArrayDS );
SHAREDPTR ( mxArrayDS );

class mexDS : public DS {
        protected:
                string step_name;
                string input_name;
                vec dt;
                vec ut;
        public:
        mexDS ():DS() {};

        void from_setting ( const Setting &set ) {
                UI::get(step_name, set, "step_name", UI::compulsory);
                UI::get(input_name, set, "input_name", UI::compulsory);
                
                shared_ptr<RV> ry = UI::build<RV> ( set, "rv_out", UI::compulsory );
                shared_ptr<RV> ru = UI::build<RV> ( set, "rv_in", UI::compulsory);

                dtsize=ry->_dsize();
                utsize=ru->_dsize();

                set_drv(*ry, *ru);
        }
        
        void step() {
                mxArray* tmp;
                mxArray* tmp2;
                // write inputs to variable input_name
                mxArray* mxinp= mexGetVariable ( "global", input_name.c_str()) ;
                vec2mxArray(ut, mxinp);
                // call function step_name
                mexCallMATLAB(1, &tmp, 0, (mxArray **) &tmp2, step_name.c_str());
                // save its results
                dt=mxArray2vec(tmp);
        }
        void write(vec &ut0){ ut=ut0;}
        void getdata(vec &dt_out){dt_out = dt;  }


        // TODO dodelat void to_setting( Setting &set ) const;
};

UIREGISTER ( mexDS );
SHAREDPTR ( mexDS );

}
#endif //MXDS_H

---