From: fred Date: Fri, 6 Dec 1996 22:41:30 +0000 (+0000) Subject: flower-1.0.12 X-Git-Tag: release/1.5.59~6642 X-Git-Url: https://git.donarmstrong.com/?a=commitdiff_plain;h=b65a088f3b9c1048ee3dcf237b956c1220252025;p=lilypond.git flower-1.0.12 --- 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; } - - - -