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

Timestamp:

07/28/09 15:31:57 (15 years ago)

Author:

smidl

Message:

missing functions in chmat - square_matrix test now passes for chmat

Files:

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

library/bdm/math/chmat.cpp

@@ -17 +17 @@
 mat chmat::to_mat() const {mat F=Ch.T() *Ch;return F;};
 void chmat::mult_sym ( const mat &C ) {
-        it_error ( "not implemented" );
+        it_assert_debug(C.cols()==dim, "Wrong dimension of U");
+        if(!qr(Ch*C.T(), Ch)) {it_warning("QR unstable in chmat mult_sym");}
 };
 void chmat::mult_sym ( const mat &C , chmat &U ) const {
-        mat Z=C*Ch;
-        U.Ch= chol(Z*Z.T());
+        it_assert_debug(C.cols()==U.dim, "Wrong dimension of U");
+        if(!qr(Ch*C.T(), U.Ch)) {it_warning("QR unstable in chmat mult_sym");}
 };
 void chmat::mult_sym_t ( const mat &C ) {
-        it_error ( "not implemented" );
+        it_assert_debug(C.rows()==dim, "Wrong dimension of U");
+        if(!qr(Ch*C, Ch)) {it_warning("QR unstable in chmat mult_sym");}
 };
 void chmat::mult_sym_t ( const mat &C, chmat &U ) const {
-        it_error ( "not implemented" );
+        it_assert_debug(C.rows()==U.dim, "Wrong dimension of U");
+        if(!qr(Ch*C, U.Ch)) {it_warning("QR unstable in chmat mult_sym");}
 };
 double chmat::logdet() const {