Changeset 701

Timestamp:

11/04/09 22:55:04 (15 years ago)

Author:

smidl

Message:

improving estimator

Files:

• applications/bdmtoolbox/mex/estimator.cpp (modified) (2 diffs)
• applications/bdmtoolbox/tutorial/userguide/pdfds_example.m (modified) (3 diffs)
• library/bdm/bdmroot.h (modified) (1 diff)

applications/bdmtoolbox/mex/estimator.cpp

@@ -148 +148 @@
                         Ename="Est"+num2str ( i );
                 }
-
+                if (Es(i)->posterior()._log_level()<1){
+                        const_cast<epdf&>(Es (i) ->posterior()).set_log_level (1);
+                }
+                
                 Es ( i )->log_register ( *L,Ename ); // estimate
         }
@@ -182 +185 @@
 
                 for ( int i=0; i<Es.length(); i++ ) {
-                        Es ( i )->bayes ( Dls ( i )->pushdown ( dt ), Dlsc(i) ->pushdown(dt) );         // update estimates
+                        if (tK + Es ( i )->_rvc().mint() > 0 ) {
+                                Es ( i )->bayes ( Dls ( i )->pushdown ( dt ), Dlsc(i) ->pushdown(dt) );         // update estimates
+                        }
                         Es ( i )->log_write ();
                 }

applications/bdmtoolbox/tutorial/userguide/pdfds_example.m

@@ -9 +9 @@
 fy.A     = [0.5, -0.9];
 fy.const = 0;
-fy.R     = 1e-1;
+fy.R     = 1e-2;
 
 % create f(u_t| )
@@ -18 +18 @@
 
 % create DS
-DS.class = 'MpdfDS';
-DS.mpdf.class  = 'mprod';
-DS.mpdf.mpdfs  = {fy, fu};
+DS.class = 'PdfDS';
+DS.pdf.class  = 'mprod';
+DS.pdf.pdfs  = {fy, fu};
 DS.init_rv = RVtimes([y,y,y], [-1,-2,-3]);
 DS.init_values = [0.1, 0.2, 0.3];
@@ -34 +34 @@
 true_theta=[fy.A fy.const]; 
 true_R=[fy.R]; 
-save mpdfds_results Data drv y u true_theta true_R
+save pdfds_results Data drv y u true_theta true_R

library/bdm/bdmroot.h

@@ -85 +85 @@
         virtual void validate() {
         }
+        //! access function
+        int _log_level() const {return log_level;}
         
 };