Changeset 18 for bdm/math

Timestamp:

02/01/08 10:49:29 (16 years ago)

Author:

smidl

Message:

update

Location:

bdm/math

Files:

• Makefile (modified) (1 diff)
• libDC.cpp (modified) (3 diffs)

bdm/math/Makefile

@@ -2 +2 @@
 
 clean:
-        rm *.o
+        $(RM) *.o

bdm/math/libDC.cpp

@@ -31 +31 @@
 
         dim = V.cols();
+        mat F;
+        vec D0;
         it_assert_debug( dim == V.rows(),"ldmat::ldmat matrix V is not square!" );
-
         
+        //decompose V in cholesky
+        D0 = ones(dim);
+        
+
+        using std::cout;
+        cout<<V;
+        chol(V,F);
+        // L and D will be allocated by ldform()
+        this->ldform(F,D0);     
 }
 
@@ -194 +204 @@
 
         it_assert_debug( A.cols()==dim,"ldmat::ldform A is not compatible" );
-        it_assert_debug( D.length()==A.rows(),"ldmat::ldform Vector D must have the length as row count of A" );
+        it_assert_debug( D0.length()==A.rows(),"ldmat::ldform Vector D must have the length as row count of A" );
 
         L=concat_vertical( zeros( n,n ), diag( sqrt( D0 ) )*A );
@@ -202 +212 @@
         
         vec w=zeros( n );
-        vec v;
+        vec v=zeros(n);
         double sum, beta, pom;