Changeset 738 for library/bdm/math/chmat.cpp

Timestamp:

11/25/09 12:46:08 (15 years ago)

Author:

mido

Message:

a few moves of code from h to cpp, however, only part of the whole library is done

Files:

• library/bdm/math/chmat.cpp (modified) (1 diff)

library/bdm/math/chmat.cpp

@@ -5 +5 @@
 
 namespace bdm {
+
+void chmat::add ( const chmat &A2, double w ) {
+        bdm_assert_debug ( dim == A2.dim, "Matrices of unequal dimension" );
+        mat pre = concat_vertical ( Ch, sqrt ( w ) * A2.Ch );
+        mat post = zeros ( pre.rows(), pre.cols() );
+        if ( !qr ( pre, post ) ) {
+                bdm_warning ( "Unstable QR in chmat add" );
+        }
+        Ch = post ( 0, dim - 1, 0, dim - 1 );
+};
 
 void chmat::opupdt ( const vec &v, double w ) {