Changeset 5812

Timestamp:

Dec 19, 2005, 9:44:27 AM (20 years ago)

Author:

magnier

Message:

fixed coeffErr

File:

• branches/eam_rel9_b1/psLib/src/math/psMinimize.c (modified) (3 diffs)

branches/eam_rel9_b1/psLib/src/math/psMinimize.c

@@ -10 +10 @@
  *  @author EAM, IfA
  *
- *  @version $Revision: 1.147.4.1.2.3 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2005-12-19 00:51:36 $
+ *  @version $Revision: 1.147.4.1.2.4 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-12-19 19:44:27 $
  *
  *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
@@ -2714 +2714 @@
             for (int iz = 0; iz < nZterm; iz++) {
                 int nx = ix+iy*nXterm+iz*nXterm*nYterm;
-                int ny = jx+jy*nXterm+jz*nXterm*nYterm;
                 myPoly->coeff[ix][iy][iz] = B->data.F64[nx];
-                myPoly->coeffErr[ix][iy][iz] = sqrt(A->data.F64[nx][ny]);
+                myPoly->coeffErr[ix][iy][iz] = sqrt(A->data.F64[nx][nx]);
             }
         }
@@ -3176 +3175 @@
                 for (int it = 0; it < nTterm; it++) {
                     int nx = ix+iy*nXterm+iz*nXterm*nYterm+it*nXterm*nYterm*nZterm;
-                    int ny = jx+jy*nXterm+jz*nXterm*nYterm+jt*nXterm*nYterm*nZterm;
                     myPoly->coeff[ix][iy][iz][it] = B->data.F64[nx];
-                    myPoly->coeffErr[ix][iy][iz][it] = sqrt(A->data.F64[nx][ny]);
+                    myPoly->coeffErr[ix][iy][iz][it] = sqrt(A->data.F64[nx][nx]);
                 }
             }