Show
Ignore:
Timestamp:
08/19/09 16:54:24 (15 years ago)
Author:
vbarta
Message:

using own error macros (basically copied from IT++, but never aborting)

Files:
1 modified

Legend:

Unmodified
Added
Removed
  • library/bdm/math/chmat.cpp

    r495 r565  
    1616        qr ( Z, R ); 
    1717        Ch = R ( 0, Ch.rows() - 1, 0, Ch.cols() - 1 ); 
    18 }; 
     18} 
     19 
    1920mat chmat::to_mat() const { 
    2021        mat F = Ch.T() * Ch; 
    2122        return F; 
    22 }; 
     23} 
     24 
    2325void chmat::mult_sym ( const mat &C ) { 
    24         it_assert_debug ( C.rows() == dim, "Wrong dimension of U" ); 
     26        bdm_assert_debug ( C.rows() == dim, "Wrong dimension of U" ); 
    2527        if ( !qr ( Ch*C.T(), Ch ) ) { 
    26                 it_warning ( "QR unstable in chmat mult_sym" ); 
     28                bdm_warning ( "QR unstable in chmat mult_sym" ); 
    2729        } 
    28 }; 
     30} 
     31 
    2932void chmat::mult_sym ( const mat &C , chmat &U ) const { 
    30         it_assert_debug ( C.rows() == U.dim, "Wrong dimension of U" ); 
     33        bdm_assert_debug ( C.rows() == U.dim, "Wrong dimension of U" ); 
    3134        if ( !qr ( Ch*C.T(), U.Ch ) ) { 
    32                 it_warning ( "QR unstable in chmat mult_sym" ); 
     35                bdm_warning ( "QR unstable in chmat mult_sym" ); 
    3336        } 
    34 }; 
     37} 
     38 
    3539void chmat::mult_sym_t ( const mat &C ) { 
    36         it_assert_debug ( C.cols() == dim, "Wrong dimension of U" ); 
     40        bdm_assert_debug ( C.cols() == dim, "Wrong dimension of U" ); 
    3741        if ( !qr ( Ch*C, Ch ) ) { 
    38                 it_warning ( "QR unstable in chmat mult_sym" ); 
     42                bdm_warning ( "QR unstable in chmat mult_sym" ); 
    3943        } 
    40 }; 
     44} 
     45 
    4146void chmat::mult_sym_t ( const mat &C, chmat &U ) const { 
    42         it_assert_debug ( C.cols() == U.dim, "Wrong dimension of U" ); 
     47        bdm_assert_debug ( C.cols() == U.dim, "Wrong dimension of U" ); 
    4348        if ( !qr ( Ch*C, U.Ch ) ) { 
    44                 it_warning ( "QR unstable in chmat mult_sym" ); 
     49                bdm_warning ( "QR unstable in chmat mult_sym" ); 
    4550        } 
    4651}