Changeset 265 for bdm/stat/libDS.cpp
- Timestamp:
- 02/09/09 23:14:58 (15 years ago)
- Files:
-
- 1 modified
Legend:
- Unmodified
- Added
- Removed
-
bdm/stat/libDS.cpp
r263 r265 44 44 void ArxDS::step() { 45 45 double tmp; 46 //shift history 47 H.shift_right ( 0, Drv.count() +Urv.count() ); 48 49 H.set_subvector ( Drv.count(), U ); // write U after Drv 50 46 51 //get regressor 47 52 rgr = rgrlnk.get_val ( H ); 48 53 // Eval Y 49 Y = model.samplecond ( rgr,tmp);54 H.set_subvector ( 0, model.samplecond ( rgr,tmp ) ); 50 55 51 //shift history52 H.shift_right ( 0, Drv.count() +Urv.count() );53 //fill new data54 H.set_subvector ( Drv.count(),Y );55 56 //Leaving U.length() empty slots - these will be filled by write()57 56 } 58 57 … … 60 59 RV fullrgr ( const RV &drv0, const RV &urv0, const RV &rrv0 ) { 61 60 RV T ( urv0 ); 62 RV pom = concat (drv0, urv0);61 RV pom = concat ( drv0, urv0 ); 63 62 int mint = rrv0.mint(); 64 63 for ( int i=0; i>mint; i-- ) { … … 69 68 } 70 69 71 ArxDS::ArxDS ( RV &drv, RV &urv, RV &rrv ) : 72 DS (drv,urv), Rrv(rrv), Hrv( fullrgr ( drv,urv,rrv )), //RVs73 Y(drv.count()), H(Hrv.count()),rgr ( rrv.count() ), //tmp variables74 rgrlnk (Hrv ,rrv ) ,model ( drv,rrv ) {70 ArxDS::ArxDS ( RV &drv, RV &urv, RV &rrv ) : 71 DS ( drv,urv ), Rrv ( rrv ), Hrv ( concat(drv,fullrgr ( drv,urv,rrv )) ), //RVs 72 H ( Hrv.count() ) ,U ( urv.count() ),rgr ( rrv.count() ), //tmp variables 73 rgrlnk ( rrv, Hrv ) ,model ( drv,rrv ) { 75 74 }