Changeset 536 for library/bdm/estim/kalman.h
- Timestamp:
- 08/16/09 18:13:31 (15 years ago)
- Files:
-
- 1 modified
Legend:
- Unmodified
- Added
- Removed
-
library/bdm/estim/kalman.h
r529 r536 125 125 return est; 126 126 } 127 const enorm<sq_T>* _e() const {128 return &est;129 }130 127 //!access function 131 128 mat& __K() { … … 214 211 return E; 215 212 }; 216 const enorm<fsqmat>* _e() const {217 return &E;218 };219 213 const mat _R() { 220 214 return P; … … 277 271 // TODO dodelat void to_setting( Setting &set ) const; 278 272 273 const enorm<chmat>& posterior() {return est;} 279 274 }; 280 275 … … 423 418 424 419 est.set_rv ( RV ( "MM", A ( 0 )->posterior().dimension(), 0 ) ); 425 est.set_parameters ( A ( 0 )-> _e()->mean(), A ( 0 )->_e()->_R() );420 est.set_parameters ( A ( 0 )->posterior().mean(), A ( 0 )->posterior()._R() ); 426 421 } 427 422 void bayes ( const vec &dt ) { … … 439 434 case 1: { 440 435 int mi = max_index ( w ); 441 const enorm<chmat> * st = ( Models ( mi )->_e() );442 est.set_parameters ( st ->mean(), st->_R() );436 const enorm<chmat> &st = Models ( mi )->posterior() ; 437 est.set_parameters ( st.mean(), st._R() ); 443 438 } 444 439 break; … … 451 446 } 452 447 } 453 //all posterior densities are equal => return the first one 454 const enorm<chmat>* _e() const { 455 return &est; 456 } 457 //all posterior densities are equal => return the first one 448 //! posterior density 458 449 const enorm<chmat>& posterior() const { 459 450 return est;