Changeset 455 for library/bdm/math/chmat.cpp
- Timestamp:
- 07/30/09 16:44:28 (15 years ago)
- Files:
-
- 1 modified
Legend:
- Unmodified
- Added
- Removed
-
library/bdm/math/chmat.cpp
r437 r455 17 17 mat chmat::to_mat() const {mat F=Ch.T() *Ch;return F;}; 18 18 void chmat::mult_sym ( const mat &C ) { 19 it_assert_debug(C. cols()==dim, "Wrong dimension of U");19 it_assert_debug(C.rows()==dim, "Wrong dimension of U"); 20 20 if(!qr(Ch*C.T(), Ch)) {it_warning("QR unstable in chmat mult_sym");} 21 21 }; 22 22 void chmat::mult_sym ( const mat &C , chmat &U ) const { 23 it_assert_debug(C. cols()==U.dim, "Wrong dimension of U");23 it_assert_debug(C.rows()==U.dim, "Wrong dimension of U"); 24 24 if(!qr(Ch*C.T(), U.Ch)) {it_warning("QR unstable in chmat mult_sym");} 25 25 }; 26 26 void chmat::mult_sym_t ( const mat &C ) { 27 it_assert_debug(C. rows()==dim, "Wrong dimension of U");27 it_assert_debug(C.cols()==dim, "Wrong dimension of U"); 28 28 if(!qr(Ch*C, Ch)) {it_warning("QR unstable in chmat mult_sym");} 29 29 }; 30 30 void chmat::mult_sym_t ( const mat &C, chmat &U ) const { 31 it_assert_debug(C. rows()==U.dim, "Wrong dimension of U");31 it_assert_debug(C.cols()==U.dim, "Wrong dimension of U"); 32 32 if(!qr(Ch*C, U.Ch)) {it_warning("QR unstable in chmat mult_sym");} 33 33 };