Changeset 884

Timestamp:

03/29/10 23:01:29 (14 years ago)

Author:

smidl

Message:

Kalman correction

Files:

• library/bdm/estim/kalman.cpp (modified) (1 diff)

library/bdm/estim/kalman.cpp

@@ -263 +263 @@
         //add diagonals for rgr
         int active_x = y_block_size;
+        int active_Bcol = 0;
         for ( int r = 0; r < urv.length(); r++ ) {
                 if (u_block_sizes(r)>0){
                         diagonal_part ( this->A, active_x + urv.size ( r ), active_x, u_block_sizes ( r ) - urv.size ( r ) );
-                        this->B.set_submatrix ( active_x, 0, eye ( urv.size ( r ) ) );
+                        this->B.set_submatrix ( active_x, active_Bcol, eye ( urv.size ( r ) ) );
+                        active_Bcol+=u_block_sizes(r);
                 }
                 active_x += u_block_sizes ( r );