Changeset 283 for bdm/stat/libDS.h

Timestamp:

02/24/09 14:14:01 (15 years ago)

Author:

smidl

Message:

get rid of BMcond + adaptation in doprava/

Files:

• bdm/stat/libDS.h (modified) (5 diffs)

bdm/stat/libDS.h

@@ -28 +28 @@
 */
 class MemDS : public DS {
+        protected:
         //! internal matrix of data
         mat Data;
@@ -45 +46 @@
         void step();
         //!Default constructor
+        MemDS () {};
         MemDS ( mat &Dat, ivec &rowid, ivec &delays );
+};
+
+/*! Read Data Matrix from an IT file
+
+*/
+class FileDS: public MemDS {
+
+public:
+        FileDS ( const string &fname, const string &varname ) :MemDS() {
+                it_file it ( fname );
+                it << Name ( varname ); 
+                it >> Data;
+                time =0;
+                //rowid and delays are ignored
+        }
+        void getdata ( vec &dt ) {
+                it_assert_debug ( dt.length() ==Data.rows(),"" );
+                dt = Data.get_col(time);
+        };
+        void getdata ( vec &dt, const ivec &indeces ){
+                it_assert_debug ( dt.length() ==indeces.length(),"" );
+                vec tmp(indeces.length());
+                tmp = Data.get_col(time);
+                dt = tmp(indeces);
+        };
+        //! returns number of data in the file;
+        int ndat(){return Data.cols();}
 };
 
@@ -96 +125 @@
         { model.set_parameters ( Th0, mu0, sqR0 );};
         //! Set
-        void set_drv(RV &yrv, RV &urv, RV &rrv){
+        void set_drv ( RV &yrv, RV &urv, RV &rrv ) {
                 Rrv = rrv;
                 Urv = urv;
-                dt_size = yrv._dsize()+urv._dsize();
-                
-                RV drv = concat(yrv,urv);
+                dt_size = yrv._dsize() +urv._dsize();
+
+                RV drv = concat ( yrv,urv );
                 Drv = drv;
                 int td = rrv.mint();
-                H.set_size(drv._dsize()*(-td+1));
-                U.set_size(Urv._dsize());
-                for (int i=-1;i>=td;i--){
-                        drv.t(-1);
-                        Drv.add(drv); //shift u1
+                H.set_size ( drv._dsize() * ( -td+1 ) );
+                U.set_size ( Urv._dsize() );
+                for ( int i=-1;i>=td;i-- ) {
+                        drv.t ( -1 );
+                        Drv.add ( drv ); //shift u1
                 }
-                rgrlnk.set_connection(rrv,Drv);
-                
+                rgrlnk.set_connection ( rrv,Drv );
+
                 dtsize = Drv._dsize();
                 utsize = Urv._dsize();
@@ -121 +150 @@
         virtual void log_add ( logger &L ) {
                 //DS::log_add ( L ); too long!!
-                L_dt=L.add ( Drv(0,dt_size),"" );
+                L_dt=L.add ( Drv ( 0,dt_size ),"" );
                 L_ut=L.add ( Urv,"" );
 
@@ -131 +160 @@
         virtual void logit ( logger &L ) {
                 //DS::logit ( L );
-                L.logit( L_dt, H.left(dt_size));
-                L.logit(L_ut, U);
-                
+                L.logit ( L_dt, H.left ( dt_size ) );
+                L.logit ( L_ut, U );
+
                 mat &A =model._A();
                 mat R =model._R();