Changeset 75 for bdm/math/libDC.cpp

Timestamp:

04/18/08 14:00:29 (16 years ago)

Author:

smidl

Message:

oprava operaci na ld

Files:

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

bdm/math/libDC.cpp

@@ -162 +162 @@
         double x = 0.0, sum;
         int i,j;
+        vec s(v.length());
+        vec S=L*v;
 
         for ( i=0; i<D.length(); i++ ) { //rows of L
                 sum = 0.0;
-                for ( j=i; j<=i; j++ ){sum+=L( i,j )*v( i );}
+                for ( j=0; j<=i; j++ ){sum+=L( i,j )*v( j );}
                 x +=D( i )*sum*sum;
+                s(i)=sum;
+        };
+        return x;
+}
+
+double ldmat::invqform( const vec &v ) const {
+        double x = 0.0;
+        int i;
+        vec pom(v.length());
+        
+        backward_substitution(L.T(),v,pom);
+        
+        for ( i=0; i<D.length(); i++ ) { //rows of L
+                x +=pom(i)*pom(i)/D(i);
         };
         return x;
@@ -298 +314 @@
                         j = i + k; //change in .m 1+1=2, here 0+0+1=1
                         s = L( j, i );
-                        for ( m = i + 1; m < ( j - 1 ); m++ ) {
+                        for ( m = i + 1; m < ( j ); m++ ) {
                                 s += L( m, i ) * Il( j, m );
                         }