Changeset 700 for library/bdm/estim/particles.cpp

Timestamp:

11/04/09 22:54:58 (15 years ago)

Author:

smidl

Message:

Making tutorial/userguide example work again (changes of mpdf and bayes)

Files:

• library/bdm/estim/particles.cpp (modified) (2 diffs)

library/bdm/estim/particles.cpp

@@ -78 +78 @@
         int n=pf->__w().length();
         vec &lls = pf->_lls();
+        Array<vec> &samples=pf->__samples();
         
         // generate samples - time step
@@ -84 +85 @@
         #pragma parallel for
         for ( i = 0; i < n; i++ ) {
-                BMs(i) -> bayes(this2bm.pushdown(yt), this2bm.get_cond(yt,cond));
+                vec bm_cond(BMs(i)->dimensionc());
+                this2bm.fill_cond(yt,cond, bm_cond);
+                pf2bm.filldown(samples(i), bm_cond);
+                BMs(i) -> bayes(this2bm.pushdown(yt), bm_cond);
                 lls ( i ) += BMs(i)->_ll();
         }