Changeset 1020 for trunk/psLib/src/dataManip/psFunctions.c

Timestamp:

Jun 14, 2004, 9:33:09 AM (22 years ago)

Author:

gusciora

Message:

...

File:

• trunk/psLib/src/dataManip/psFunctions.c (modified) (17 diffs)

trunk/psLib/src/dataManip/psFunctions.c

@@ -4 +4 @@
  *         routines.
  *
- *  @author George Gusciora, MHPCC
+ *  This file will hold the functions for allocated, freeing, and evaluating
+ *  polynomials.  It also contains a Gaussian functions.
  *
- *  @version $Revision: 1.11 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2004-06-07 00:32:53 $
+ *  @version $Revision: 1.12 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2004-06-14 19:32:42 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -103 +104 @@
 psPolynomial1D *psPolynomial1DAlloc(int n)
 {
+    int i = 0;
     psPolynomial1D *newPoly = NULL;
 
@@ -110 +112 @@
     newPoly->coeffErr = (float *) psAlloc(n * sizeof(float));
     newPoly->mask     = (char *)  psAlloc(n * sizeof(char));
+    for (i=0;i<n;i++) {
+        newPoly->coeff[i] = 0.0;
+        newPoly->coeffErr[i] = 0.0;
+        newPoly->mask[i] = 0.0;
+    }
 
     return(newPoly);
@@ -117 +124 @@
 {
     int x = 0;
+    int y = 0;
     psPolynomial2D *newPoly = NULL;
 
@@ -131 +139 @@
         newPoly->mask[x]     = (char *)  psAlloc(nY * sizeof(char));
     }
-
+    for (x=0;x<nX;x++) {
+        for (y=0;y<nY;y++) {
+            newPoly->coeff[x][y] = 0.0;
+            newPoly->coeffErr[x][y] = 0.0;
+            newPoly->mask[x][y] = 0;
+        }
+    }
 
     return(newPoly);
@@ -140 +154 @@
     int x = 0;
     int y = 0;
+    int z = 0;
     psPolynomial3D *newPoly = NULL;
 
@@ -160 +175 @@
         }
     }
+    for (x=0;x<nX;x++) {
+        for (y=0;y<nY;y++) {
+            for (z=0;z<nZ;z++) {
+                newPoly->coeff[x][y][z] = 0.0;
+                newPoly->coeffErr[x][y][z] = 0.0;
+                newPoly->mask[x][y][z] = 0;
+            }
+        }
+    }
 
     return(newPoly);
@@ -169 +193 @@
     int x = 0;
     int y = 0;
+    int z = 0;
     psPolynomial4D *newPoly = NULL;
 
@@ -195 +220 @@
         }
     }
+    for (w=0;w<nW;w++) {
+        for (x=0;x<nX;x++) {
+            for (y=0;y<nY;y++) {
+                for (z=0;z<nZ;z++) {
+                    newPoly->coeff[w][x][y][z] = 0.0;
+                    newPoly->coeffErr[w][x][y][z] = 0.0;
+                    newPoly->mask[w][x][y][z] = 0;
+                }
+            }
+        }
+    }
 
     return(newPoly);
@@ -394 +430 @@
 psDPolynomial1D *psDPolynomial1DAlloc(int n)
 {
+    int i = 0;
     psDPolynomial1D *newPoly = NULL;
 
@@ -401 +438 @@
     newPoly->coeffErr = (double *) psAlloc(n * sizeof(double));
     newPoly->mask     = (char *)  psAlloc(n * sizeof(char));
+    for (i=0;i<n;i++) {
+        newPoly->coeff[i] = 0.0;
+        newPoly->coeffErr[i] = 0.0;
+        newPoly->mask[i] = 0.0;
+    }
 
     return(newPoly);
@@ -408 +450 @@
 {
     int x = 0;
+    int y = 0;
     psDPolynomial2D *newPoly = NULL;
 
@@ -422 +465 @@
         newPoly->mask[x]     = (char *)  psAlloc(nY * sizeof(char));
     }
+    for (x=0;x<nX;x++) {
+        for (y=0;y<nY;y++) {
+            newPoly->coeff[x][y] = 0.0;
+            newPoly->coeffErr[x][y] = 0.0;
+            newPoly->mask[x][y] = 0;
+        }
+    }
 
     return(newPoly);
@@ -430 +480 @@
     int x = 0;
     int y = 0;
+    int z = 0;
     psDPolynomial3D *newPoly = NULL;
 
@@ -450 +501 @@
         }
     }
+    for (x=0;x<nX;x++) {
+        for (y=0;y<nY;y++) {
+            for (z=0;z<nZ;z++) {
+                newPoly->coeff[x][y][z] = 0.0;
+                newPoly->coeffErr[x][y][z] = 0.0;
+                newPoly->mask[x][y][z] = 0;
+            }
+        }
+    }
 
     return(newPoly);
@@ -459 +519 @@
     int x = 0;
     int y = 0;
+    int z = 0;
     psDPolynomial4D *newPoly = NULL;
 
@@ -485 +546 @@
         }
     }
+    for (w=0;w<nW;w++) {
+        for (x=0;x<nX;x++) {
+            for (y=0;y<nY;y++) {
+                for (z=0;z<nZ;z++) {
+                    newPoly->coeff[w][x][y][z] = 0.0;
+                    newPoly->coeffErr[w][x][y][z] = 0.0;
+                    newPoly->mask[w][x][y][z] = 0;
+                }
+            }
+        }
+    }
 
     return(newPoly);