flower-1.0.12

diff --git a/flower/choleski.cc b/flower/choleski.cc
index 82740a27c18addcbeca4278ac401dca765c2b773..4202a9b1e5eed0a3b64349f7fcb1954ea5a5b073 100644 (file)
--- a/flower/choleski.cc
+++ b/flower/choleski.cc
@@ -15,13 +15,11 @@ Choleski_decomposition::solve(Vector rhs)const
            sum += y(j) * L(i,j);
        y(i) = (rhs(i) - sum)/L(i,i);
     }
-    for (int i=0; i < n; i++) {
-       assert(D(i));
+    for (int i=0; i < n; i++)
        y(i) /= D(i);
-    }
 
     // backward subst
-    Vector x(n);
+    Vector &x(rhs);            // using input as return val.
     for (int i=n-1; i >= 0; i--) {
        Real sum(0.0);
        for (int j=i+1; j < n; j++)
@@ -91,7 +89,3 @@ Choleski_decomposition::inverse() const
     
     return invm;
 }
-
-
-
-