Changeset 8 for libDC.cpp

Timestamp:

01/23/08 11:32:31 (16 years ago)

Author:

smidl

Message:

Kalmany funkci, PF nefunkci

Files:

• libDC.cpp (modified) (3 diffs)

libDC.cpp

@@ -115 +115 @@
 }
 
-void ldmat::mult_qform( const mat &C, bool trans ) {
+void ldmat::mult_sym( const mat &C, bool trans ) {
 
 //TODO better
 
-        it_assert_debug( C.cols()==L.cols(), "ldmat::mult_qform wrong input argument" );
+        it_assert_debug( C.cols()==L.cols(), "ldmat::mult_sym wrong input argument" );
         mat Ct=C;
 
@@ -128 +128 @@
                 Ct = C.transpose();
                 Ct *= this->to_mat();
+                Ct *= C;
+        }
+
+        ldmat Lnew=ldmat( Ct );
+        L = Lnew.L;
+        D = Lnew.D;
+}
+
+void ldmat::mult_sym( const mat &C, ldmat &U, bool trans ) {
+
+//TODO better
+
+//TODO input test
+
+        mat Ct=C;
+
+        if ( trans==false ) { // return C*this*C'
+                Ct *= U.to_mat();
+                Ct *= C.transpose();
+        } else {        // return C'*this*C
+                Ct = C.transpose();
+                Ct *= U.to_mat();
                 Ct *= C;
         }
@@ -232 +254 @@
                 kD = 1.0;
                 *kr = 0.0;
-                if ( *Df < -threshold ) it_warning( "Problem in dydr: subraction of dyad results in negative definitness. Likely mistake in calling function." );
+                if ( *Df < -threshold ) {
+                it_warning( "Problem in dydr: subraction of dyad results in negative definitness. Likely mistake in calling function." );}
                 *Df = 0.0;
         }