Changeset 744 for library/bdm/estim/kalman.cpp
- Timestamp:
- 11/28/09 14:38:48 (15 years ago)
- Files:
-
- 1 modified
Legend:
- Unmodified
- Added
- Removed
-
library/bdm/estim/kalman.cpp
r739 r744 264 264 int active_x = y_block_size; 265 265 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 } 268 270 active_x += u_block_sizes ( r ); 269 271 }