release: 0.0.2

[lilypond.git] / choleski.cc

diff --git a/choleski.cc b/choleski.cc
index db3fe82d60e370eb7bd7a6e06b749dcb776f4a6b..2d8c68aa66f157e4a289d6f997a3bf5e41a2bea8 100644 (file)
--- a/choleski.cc
+++ b/choleski.cc
@@ -18,6 +18,8 @@ Choleski_decomposition::solve(Vector rhs)const
        assert(D(i));
        y(i) /= D(i);
     }
+
+    // backward subst
     Vector x(n);
     for (int i=n-1; i >= 0; i--) {
        Real sum(0.0);
@@ -28,6 +30,10 @@ Choleski_decomposition::solve(Vector rhs)const
     return x;
 }
 
+/*
+  Standard matrix algorithm.
+  */
+
 Choleski_decomposition::Choleski_decomposition(Matrix P)
     : L(P.dim()), D(P.dim())
 {
@@ -49,9 +55,9 @@ Choleski_decomposition::Choleski_decomposition(Matrix P)
        D(k) = d;
     }
 
-    #ifdef NDEBUG
+#ifdef NDEBUG
     assert((original()-P).norm() < EPS);
-    #endif
+#endif
 }
      
 Matrix
@@ -79,7 +85,11 @@ Choleski_decomposition::inverse() const
     Matrix I1(n), I2(original());
     I1.unit();
     assert((I1-original()*invm).norm() < EPS);
-    #endif
+#endif
     
     return invm;
 }
+
+
+
+