Changeset 1080 for library/bdm/base

Timestamp:

06/10/10 22:13:50 (14 years ago)

Author:

smidl

Message:

fix tests

Files:

• library/bdm/base/datasources.h (modified) (2 diffs)

library/bdm/base/datasources.h

@@ -248 +248 @@
     //! numeric values of regressor
     vec rgr;
-
+        
+        //! RV for initial conditions
+        RV init_rv;
+        //! initial conditions
+        vec init_values;
+                
 public:
     void step();
@@ -279 +284 @@
     */
     void from_setting ( const Setting &set ) {
-        ipdf = UI::build<pdf> ( set, "pdf", UI::compulsory );
-
-        RV Yrv = ipdf->_rv();
-        // get unique rvs form rvc
-        RV rgrv0 = ipdf->_rvc().remove_time();
-        // input is what in not in Yrv
-        Urv = rgrv0.subt ( Yrv );
-        Drv=  concat(Yrv,Urv);
-        // connect input and output to rvc
-        ut2rgr.set_connection ( ipdf->_rvc(), Urv );
-        yt2rgr.set_connection ( ipdf->_rvc(), Yrv );
-
-        //set history - if given
-        shared_ptr<RV> rv_ini = UI::build<RV> ( set, "init_rv", UI::optional );
-        if ( rv_ini ) { // check if
-            vec val;
-            UI::get ( val, set, "init_values", UI::optional );
-            if ( val.length() != rv_ini->_dsize() ) {
-                bdm_error ( "init_rv and init_values fields have incompatible sizes" );
-            } else {
-                ut2rgr.set_history ( *rv_ini, val );
-                yt2rgr.set_history ( *rv_ini, val );
-            }
-        }
-
-        yt = zeros ( ipdf->dimension() );
-        rgr = zeros ( ipdf->dimensionc() );
-        ut = zeros ( Urv._dsize() );
-
-        utsize = ut.length();
-        dtsize = yt.length() + utsize;
+                ipdf = UI::build<pdf> ( set, "pdf", UI::compulsory );
+                UI::get ( init_rv, set, "init_rv", UI::optional);
+                UI::get ( init_values, set, "init_values", UI::optional);
     }
 
     void validate() {
-        DS::validate();
+                RV Yrv = ipdf->_rv();
+                // get unique rvs form rvc
+                RV rgrv0 = ipdf->_rvc().remove_time();
+                // input is what in not in Yrv
+                Urv = rgrv0.subt ( Yrv );
+                Drv=  concat(Yrv,Urv);
+                // connect input and output to rvc
+                ut2rgr.set_connection ( ipdf->_rvc(), Urv );
+                yt2rgr.set_connection ( ipdf->_rvc(), Yrv );
+                
+                if ( init_rv._dsize()>0 ) { // check if
+                        if ( init_values.length() != init_rv._dsize() ) {
+                                bdm_error ( "init_rv and init_values fields have incompatible sizes" );
+                        } else {
+                                ut2rgr.set_history ( init_rv, init_values );
+                                yt2rgr.set_history ( init_rv, init_values );
+                        }
+                }
+                
+                yt = zeros ( ipdf->dimension() );
+                rgr = zeros ( ipdf->dimensionc() );
+                ut = zeros ( Urv._dsize() );
+                
+                utsize = ut.length();
+                dtsize = yt.length() + utsize;
+                
+                DS::validate();
 
         //taken from sample() - shift of history is not done here