Changeset 76 for bdm/math/chmat.h

Timestamp:

04/18/08 14:01:20 (16 years ago)

Author:

smidl

Message:

oprava chol + invqform

Files:

• bdm/math/chmat.h (modified) (3 diffs)

bdm/math/chmat.h

@@ -21 +21 @@
 /*! \brief Symmetric matrix stored in square root decomposition using upper cholesky
 
-This matrix represent $A=Ch Ch'$ where only the upper triangle is stored;
+This matrix represent \f$A=Ch' Ch\f$ where only the upper triangle \f$Ch\f$ is stored;
 
 */
@@ -39 +39 @@
         vec sqrt_mult ( const vec &v ) const;
         double qform ( const vec &v ) const;
+        double invqform ( const vec &v ) const;
         void clear();
         void add ( const chmat &A2, double w=1.0 ) {};
         //!Inversion in the same form, i.e. cholesky
-        void inv ( chmat &Inv ) const   { ( Inv.Ch ) = itpp::inv ( Ch );}; //ltuinv is wrong?
+        void inv ( chmat &Inv ) const   { ( Inv.Ch ) = itpp::inv ( Ch ).T();}; //Fixme: can be more efficient
         ;
 //      void inv ( mat &Inv );
@@ -52 +53 @@
         //! Default constructor
         chmat ( const vec &v) : sqmat ( v.length() ),Ch ( diag(sqrt(v)) ) {};
+        //! Copy constructor
+        chmat ( const chmat &Ch0) : sqmat ( Ch0.dim),Ch ( Ch0.dim,Ch0.dim ) {Ch=Ch0.Ch;};
         //! Default constructor (m3k:cholform)
         chmat ( const mat &M ) : sqmat ( M.rows() ),Ch ( M.rows(),M.cols() ) {