[Toon-members] TooN Cholesky.h

[Ethan Eade]

CVSROOT:        /cvsroot/toon
Module name:    TooN
Changes by:     Ethan Eade <ethaneade>  06/07/16 19:43:16

Modified files:
        .              : Cholesky.h 

Log message:
        Another silly bug.

CVSWeb URLs:
http://cvs.savannah.gnu.org/viewcvs/TooN/Cholesky.h?cvsroot=toon&r1=1.8&r2=1.9

Patches:
Index: Cholesky.h
===================================================================
RCS file: /cvsroot/toon/TooN/Cholesky.h,v
retrieving revision 1.8
retrieving revision 1.9
diff -u -b -r1.8 -r1.9
--- Cholesky.h  14 Jul 2006 18:45:37 -0000      1.8
+++ Cholesky.h  16 Jul 2006 19:43:16 -0000      1.9
@@ -82,8 +82,9 @@
        {
            for (int col = 0; col<S; col++) {
                Vector<S> t,x;
-               for (int i=col; i<S; i++) {
-                   double psum = 1;
+               t[col] = invdiag[i];
+               for (int i=col+1; i<S; i++) {
+                   double psum = 0;
                    for (int j=col; j<i; j++)
                        psum -= L[i][j]*t[j];               
                    t[i] = invdiag[i] * psum;