Changeset 679 for library/bdm/estim/particles.cpp
- Timestamp:
- 10/23/09 00:05:25 (15 years ago)
- Files:
-
- 1 modified
Legend:
- Unmodified
- Added
- Removed
-
library/bdm/estim/particles.cpp
r665 r679 46 46 } 47 47 48 void PF::bayes ( const vec &dt ) { 49 vec yt=dt.left(obs->dimension()); 50 vec ut=dt.get(obs->dimension(),dt.length()-1); 48 void PF::bayes( const vec &yt, const vec &cond ) { 49 const vec &ut=cond; //todo 51 50 52 51 int i; … … 74 73 // } 75 74 76 void MPF::bayes ( const vec & dt) {75 void MPF::bayes ( const vec &yt, const vec &cond ) { 77 76 // follows PF::bayes in most places!!! 78 77 int i; … … 85 84 #pragma parallel for 86 85 for ( i = 0; i < n; i++ ) { 87 BMs(i) -> condition(pf->posterior()._sample(i)); 88 BMs(i) -> bayes(dt); 86 BMs(i) -> bayes(this2bm.pushdown(yt), this2bm.get_cond(yt,cond)); 89 87 lls ( i ) += BMs(i)->_ll(); 90 88 }