Changeset 455

Timestamp:

07/30/09 16:44:28 (16 years ago)

Author:

smidl

Message:

test mult_sym_t(.,.)

Location:

library

Files:

• bdm/math/chmat.cpp (modified) (1 diff)
• tests/square_mat_test.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_assert_debug(C.cols()==dim, "Wrong dimension of U");
+        it_assert_debug(C.rows()==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 {
-        it_assert_debug(C.cols()==U.dim, "Wrong dimension of U");
+        it_assert_debug(C.rows()==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_assert_debug(C.rows()==dim, "Wrong dimension of U");
+        it_assert_debug(C.cols()==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_assert_debug(C.rows()==U.dim, "Wrong dimension of U");
+        it_assert_debug(C.cols()==U.dim, "Wrong dimension of U");
         if(!qr(Ch*C, U.Ch)) {it_warning("QR unstable in chmat mult_sym");}
 };

library/tests/square_mat_test.cpp

@@ -84 +84 @@
     CHECK_CLOSE(res, sqmat2.to_mat(), epsilon);
 
+    mat C = randu(sz, sz-1);
+         TMatrix CAC(sz-1);
+    sqmat.mult_sym_t(C,CAC);
+    res = (C.T() * A) * C;
+    CHECK_CLOSE(res, CAC.to_mat(), epsilon);
+
     sqmat2 = sqmat;
     sqmat2.mult_sym_t(B);