Show
Ignore:
Timestamp:
01/14/10 21:37:57 (14 years ago)
Author:
smidl
Message:

arena experiment + numerical fixes

Files:
1 modified

Legend:

Unmodified
Added
Removed
  • library/bdm/mpdm/arx_agent.h

    r776 r787  
    6565                } 
    6666                void adapt(const vec &glob_dt){                  
    67                         dlDt.filldown(glob_dt, dt); 
     67                        dlDt.filldown(glob_dt, dt);  
     68                        dlDt.store_data(glob_dt); // do not ever use dlDt again!  
     69                         
    6870                        lqg_arx->adapt(dt); 
    6971                        // remove old predictors 
     
    8082                                merger->merge(); 
    8183                                enorm<chmat> joint_pred;  
    82                                 joint_pred.set_parameters(merger->mean(), merger->covariance()); 
     84                                mat Cov=merger->covariance(); 
     85                                if (sumsum(Cov)==0.0){ 
     86                                        bdm_error("merging failed"); 
     87                                } 
     88                                joint_pred.set_parameters(merger->mean(), Cov); 
    8389                                joint_pred.set_rv(merger->_rv()); 
    8490                                 
     
    96102                        ut=lqg_arx->ctrlaction(dt); 
    97103                        dlU.filldown(ut,glob_ut); 
    98                         // save dat for next time 
    99                         dlDt.store_data(dt); 
    100104                } 
    101105                virtual void ds_register(const DS &ds){