Changeset 565 for library/bdm/math/square_mat.cpp
- Timestamp:
- 08/19/09 16:54:24 (15 years ago)
- Files:
-
- 1 modified
Legend:
- Unmodified
- Added
- Removed
-
library/bdm/math/square_mat.cpp
r536 r565 35 35 }; 36 36 fsqmat::fsqmat ( const mat &M0 ) : sqmat ( M0.cols() ) { 37 it_assert_debug ( ( M0.cols() == M0.rows() ), "M0 must be square" );37 bdm_assert_debug ( ( M0.cols() == M0.rows() ), "M0 must be square" ); 38 38 M = M0; 39 39 }; … … 65 65 ldmat::ldmat ( const mat &V ) : sqmat ( V.cols() ) { 66 66 67 it_assert_debug ( dim == V.rows(), "ldmat::ldmat matrix V is not square!" );67 bdm_assert_debug ( dim == V.rows(), "ldmat::ldmat matrix V is not square!" ); 68 68 69 69 // L and D will be allocated by ldform() … … 81 81 double *rraw = r._data(); 82 82 83 it_assert_debug ( v.length() == dim, "LD::ldupdt vector v is not compatible with this ld." );83 bdm_assert_debug ( v.length() == dim, "LD::ldupdt vector v is not compatible with this ld." ); 84 84 85 85 for ( int i = dim - 1; i >= 0; i-- ) { … … 123 123 int dim = D.length(); 124 124 125 it_assert_debug ( ld2.D.length() == dim, "LD.add() incompatible sizes of LDs;" );125 bdm_assert_debug ( ld2.D.length() == dim, "LD.add() incompatible sizes of LDs" ); 126 126 127 127 //Fixme can be done more efficiently either via dydr or ldform … … 231 231 int mn = ( m < n ) ? m : n ; 232 232 233 // it_assert_debug( A.cols()==dim,"ldmat::ldform A is not compatible" ); 234 it_assert_debug ( D0.length() == A.rows(), "ldmat::ldform Vector D must have the length as row count of A" ); 233 bdm_assert_debug ( D0.length() == A.rows(), "ldmat::ldform Vector D must have the length as row count of A" ); 235 234 236 235 L = concat_vertical ( zeros ( n, n ), diag ( sqrt ( D0 ) ) * A ); … … 389 388 *kr = 0.0; 390 389 if ( *Df < -threshold ) { 391 it_warning ( "Problem in dydr: subraction of dyad results in negative definitness. Likely mistake in calling function." );390 bdm_warning ( "Problem in dydr: subraction of dyad results in negative definitness. Likely mistake in calling function." ); 392 391 } 393 392 *Df = 0.0;