From b65a088f3b9c1048ee3dcf237b956c1220252025 Mon Sep 17 00:00:00 2001
From: fred <fred>
Date: Fri, 6 Dec 1996 22:41:30 +0000
Subject: [PATCH] flower-1.0.12

---
 flower/choleski.cc | 10 ++--------
 1 file changed, 2 insertions(+), 8 deletions(-)

diff --git a/flower/choleski.cc b/flower/choleski.cc
index 82740a27c1..4202a9b1e5 100644
--- 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;
 }
-
-
-
-
-- 
2.39.5