Changeset 613 for library/bdm/base/datasources.h

Timestamp:

09/13/09 23:14:31 (15 years ago)

Author:

smidl

Message:

documentation

Files:

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

library/bdm/base/datasources.h

@@ -128 +128 @@
                         iepdf=UI::build<epdf> ( set,"epdf",UI::compulsory );
                         bdm_assert(iepdf->isnamed(), "Input epdf must be named, check if RV is given correctly");
-                        dt =  iepdf->sample();
+                        dt =  zeros(iepdf->dimension());
                         dtsize=dt.length();
                         set_drv(iepdf->_rv(),RV());
                         utsize =0;
+                }
+                void validate() {
+                        dt = iepdf->sample();
                 }
 };
@@ -172 +175 @@
                 \code
                 class = "MpdfDS";
-                mpdf = {class="mpdf_offspring", ...}// mpdf to simulate
+                mpdf = {class="mpdf_offspring", ...};  // mpdf to simulate
+                --- optional ---
+                init_rv = {class="RV",names=...};      // define what rv to initialize - typically delayed values!
+                init_values = [...];                   // vector of initial values corresponding to init_rv
                 \endcode
 
+                If init_rv is not given, init_values are set to zero.
                 */
                 void from_setting ( const Setting &set ) {
@@ -188 +195 @@
                         ut2rgr.set_connection(impdf->_rvc(), Urv); 
                         yt2rgr.set_connection(impdf->_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 ( impdf->dimension() );
@@ -196 +216 @@
                         utsize=ut.length();
                         dtsize = ytsize+utsize;
+                        validate();
+                }
+                void validate() {
+                        //taken from sample() - shift of history is not done here
+                        ut2rgr.filldown ( ut,rgr );
+                        yt2rgr.filldown ( yt,rgr );
+                        yt=impdf->samplecond ( rgr );
                 }
 };