Changeset 76 for bdm/math/chmat.cpp

Timestamp:

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

Author:

smidl

Message:

oprava chol + invqform

Files:

• bdm/math/chmat.cpp (modified) (2 diffs)

bdm/math/chmat.cpp

@@ -11 +11 @@
         mat Z;
         mat R;
-        mat Q;
-        Z = concat_vertical ( Ch,mat(w*v) );
-         qr ( Z,Q,R );
+        mat V(1,v.length());
+        V.set_row(0,v*w);
+        Z = concat_vertical ( Ch,V );
+        qr ( Z,R );
         Ch = R ( 0, Ch.rows()-1, 0, Ch.cols()-1 );
 };
-inline mat chmat::to_mat(){mat F=Ch.T()*Ch;return F;};
+inline mat chmat::to_mat() {mat F=Ch.T() *Ch;return F;};
 void chmat::mult_sym ( const mat &C ) {
-it_error("not implemented");
+        it_error ( "not implemented" );
 };
-void chmat::mult_sym ( const mat &C , chmat &U) const {
-it_error("not implemented");
+void chmat::mult_sym ( const mat &C , chmat &U ) const {
+        it_error ( "not implemented" );
 };
 void chmat::mult_sym_t ( const mat &C ) {
-it_error("not implemented");
+        it_error ( "not implemented" );
 };
 void chmat::mult_sym_t ( const mat &C, chmat &U ) const {
-it_error("not implemented");
+        it_error ( "not implemented" );
 };
 double chmat::logdet() const {
         double ldet=0.0; int i;
         //sum of logs of (possibly negative!) diagonal entries
-        for (i=0;i<Ch.rows();i++ ) {ldet+=log ( std::fabs(Ch ( i,i )) );}
+        for ( i=0;i<Ch.rows();i++ ) {ldet+=log ( std::fabs ( Ch ( i,i ) ) );}
         return 2*ldet; //compensate for Ch being sqrt()
 };
@@ -38 +39 @@
 inline vec chmat::sqrt_mult ( const vec &v ) const {vec pom; pom = Ch*v; return pom;};
 inline double chmat::qform ( const vec &v ) const {vec pom; pom = Ch*v; return pom*pom;};
+inline double chmat::invqform ( const vec &v ) const {
+        vec pom(v.length());
+        forward_substitution(Ch.T(),v,pom);
+        return pom*pom;
+};
 inline void chmat::clear() {Ch.clear();};