Changeset 437 for library/bdm/math/chmat.h

Timestamp:

07/28/09 15:31:57 (15 years ago)

Author:

smidl

Message:

missing functions in chmat - square_matrix test now passes for chmat

Files:

• library/bdm/math/chmat.h (modified) (2 diffs)

library/bdm/math/chmat.h

@@ -40 +40 @@
         void clear();
         //! add another chmat \c A2 with weight \c w.
-        void add ( const chmat &A2, double w=1.0 ) {};
+        void add ( const chmat &A2, double w=1.0 ) {
+                it_assert_debug(dim==A2.dim,"Matrices of unequal dimension");
+                mat pre=concat_vertical(Ch,sqrt(w)*A2.Ch); 
+                mat post=zeros(pre.rows(),pre.cols()); 
+                if(!qr(pre,post)) {it_warning("Unstable QR in chmat add");}
+                Ch=post(0,dim-1,0,dim-1);
+        };
         //!Inversion in the same form, i.e. cholesky
         void inv ( chmat &Inv ) const   { ( Inv.Ch ) = itpp::inv ( Ch ).T();}; //Fixme: can be more efficient
@@ -80 +86 @@
         //! Operators
         chmat& operator += ( const chmat &A2 );
-        chmat& operator -= ( const chmat &A2 ) ;
-        chmat& operator * ( const chmat &A2 );
-        chmat& operator * ( const double &d ){Ch*d; return *this;};
+        chmat& operator -= ( const chmat &A2 );
+        chmat& operator * ( const double &d ){Ch*sqrt(d); return *this;};
         chmat& operator = ( const chmat &A2 ){Ch=A2.Ch;dim=A2.dim;return *this;}
+        chmat& operator *= ( double x ){Ch*=sqrt(x); return *this;};
 };