Changeset 4901 for trunk/psphot/src/psPolynomials.c

Timestamp:

Aug 29, 2005, 5:44:55 PM (21 years ago)

Author:

eugene

Message:

various psphot updates

File:

• trunk/psphot/src/psPolynomials.c (modified) (7 diffs)

trunk/psphot/src/psPolynomials.c

@@ -42 +42 @@
 }
 
+static void psPolynomial1DFree (psPolynomial1D *poly) {
+
+  if (poly == NULL) return;
+  psFree (poly->coeff);
+  psFree (poly->coeffErr);
+  psFree (poly->mask);
+  return;
+}
+
 // XXX EAM : use Nterm = Norder + 1 definition  
 // XXX EAM : should we provide both order and nterms in struct?
@@ -54 +63 @@
 
     newPoly = (psPolynomial1D* ) psAlloc(sizeof(psPolynomial1D));
-    // p_psMemSetDeallocator(newPoly, (psFreeFcn) polynomial1DFree);
-    // XXX EAM : me, being lazy
-
     newPoly->type = type;
     newPoly->n = nOrder;
@@ -67 +73 @@
         newPoly->mask[i] = 0;
     }
+    p_psMemSetDeallocator(newPoly, (psFreeFcn) psPolynomial1DFree);
     return(newPoly);
 }
@@ -185 +192 @@
             myPoly->coeff[k] = coeffs->data.F64[k];
         }
+        psFree (ALUD);
+        psFree (outPerm);
+        psFree (coeffs);
     }
 
@@ -262 +272 @@
 }
 
+static void psPolynomial2DFree(psPolynomial2D *poly) {
+
+  if (poly == NULL) return;
+  for (int i = 0; i < poly->nX + 1; i++) {
+    psFree (poly->coeff[i]);
+    psFree (poly->coeffErr[i]);
+    psFree (poly->mask[i]);
+  }
+  psFree (poly->coeff);
+  psFree (poly->coeffErr);
+  psFree (poly->mask);
+  return;
+}
+
 // XXX EAM : use Nterm = Norder + 1 definition  
 // the user requests a polynomial of order Norder
@@ -277 +301 @@
 
     newPoly = (psPolynomial2D* ) psAlloc(sizeof(psPolynomial2D));
-    // p_psMemSetDeallocator(newPoly, (psFreeFcn) polynomial2DFree);
-    // XXX EAM : me, being lazy
-
     newPoly->type = type;
     newPoly->nX = nXorder;
@@ -299 +320 @@
         }
     }
+    p_psMemSetDeallocator(newPoly, (psFreeFcn) psPolynomial2DFree);
     return(newPoly);
 }