Show
Ignore:
Timestamp:
11/28/09 14:38:48 (15 years ago)
Author:
smidl
Message:

Working unitsteps and controlloop + corresponding fixes

Files:
1 modified

Legend:

Unmodified
Added
Removed
  • library/bdm/estim/kalman.cpp

    r739 r744  
    264264        int active_x = y_block_size; 
    265265        for ( int r = 0; r < urv.length(); r++ ) { 
    266                 diagonal_part ( this->A, active_x + urv.size ( r ), active_x, u_block_sizes ( r ) - urv.size ( r ) ); 
    267                 this->B.set_submatrix ( active_x, 0, eye ( urv.size ( r ) ) ); 
     266                if (u_block_sizes(r)>0){ 
     267                        diagonal_part ( this->A, active_x + urv.size ( r ), active_x, u_block_sizes ( r ) - urv.size ( r ) ); 
     268                        this->B.set_submatrix ( active_x, 0, eye ( urv.size ( r ) ) ); 
     269                } 
    268270                active_x += u_block_sizes ( r ); 
    269271        }