Legend:
- Unmodified
- Added
- Removed
-
bdm/math/chmat.h
r95 r98 42 42 double invqform ( const vec &v ) const; 43 43 void clear(); 44 //! add another chmat \c A2 with weight \c w. 44 45 void add ( const chmat &A2, double w=1.0 ) {}; 45 46 //!Inversion in the same form, i.e. cholesky -
bdm/math/libDC.cpp
r75 r98 213 213 int mn = (m<n) ? m :n ; 214 214 215 it_assert_debug( A.cols()==dim,"ldmat::ldform A is not compatible" );215 // it_assert_debug( A.cols()==dim,"ldmat::ldform A is not compatible" ); 216 216 it_assert_debug( D0.length()==A.rows(),"ldmat::ldform Vector D must have the length as row count of A" ); 217 217 … … 299 299 L.del_rows(0,m-1); 300 300 D.del(0,m-1); 301 302 dim = L.rows(); 301 303 } 302 304 -
bdm/math/libDC.h
r85 r98 245 245 void setL (const vec &nL){L=nL;} 246 246 247 //! Access functions 248 const vec& _D() const {return D;} 249 //! Access functions 250 const mat& _L() const {return L;} 251 247 252 //! add another ldmat matrix 248 253 ldmat& operator += ( const ldmat &ldA ); … … 254 259 //! print both \c L and \c D 255 260 friend std::ostream &operator<< ( std::ostream &os, const ldmat &sq ); 256 257 261 258 262 protected: