Changeset 1921 for trunk/psLib/src/math/psPolynomial.c

Timestamp:

Sep 28, 2004, 1:27:37 PM (22 years ago)

Author:

gusciora

Message:

Chebyshev polynomila fitting routines.

File:

• trunk/psLib/src/math/psPolynomial.c (modified) (6 diffs)

trunk/psLib/src/math/psPolynomial.c

@@ -7 +7 @@
  *  polynomials.  It also contains a Gaussian functions.
  *
- *  @version $Revision: 1.42 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2004-09-27 23:41:42 $
+ *  @version $Revision: 1.43 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2004-09-28 23:27:37 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -486 +486 @@
 float p_psChebPolynomial1DEval(float x, const psPolynomial1D* myPoly)
 {
+    /*
     psVector *d;
     int n;
@@ -507 +508 @@
     psFree(d);
     return(tmp);
-    /*
-        int n;
-        int i;
-        float tmp;
-        psPolynomial1D **chebPolys = NULL;
-     
-        n = myPoly->n;
-        chebPolys = CreateChebyshevPolys(n);
-     
-        tmp = 0.0;
-        for (i=0;i<myPoly->n;i++) {
-            tmp+= (myPoly->coeff[i] * psPolynomial1DEval(x, chebPolys[i]));
-        }
-        tmp-= (myPoly->coeff[0]/2.0);
-        return(tmp);
     */
+
+    int n;
+    int i;
+    float tmp;
+    psPolynomial1D **chebPolys = NULL;
+
+    n = myPoly->n;
+    chebPolys = CreateChebyshevPolys(n);
+
+    tmp = 0.0;
+    for (i=0;i<myPoly->n;i++) {
+        tmp+= (myPoly->coeff[i] * psPolynomial1DEval(x, chebPolys[i]));
+        //            printf("HMMM: psPolynomial1DEval(%f, chebPolys[%d]) is %f\n", x, i, psPolynomial1DEval(x, chebPolys[i]));
+    }
+    tmp-= (myPoly->coeff[0]/2.0);
+
+
+    return(tmp);
 }
 
@@ -1942 +1946 @@
 
 /*****************************************************************************
-p_psInterpolate1D(): This routine will take as input psVectors domain and
+p_psVectorInterpolate(): This routine will take as input psVectors domain and
 range, and the x value, assumed to lie with the domain vector.  It produces
 as output the LaGrange interpolated value of a polynomial of the specified
@@ -1959 +1963 @@
             "---- p_psVectorInterpolate() begin ----\n");
 
-    if (domain->type.type != range->type.type != x->type.type) {
+    if ((domain->type.type != range->type.type) ||
+            (domain->type.type != x->type.type)) {
         // XXX psError
+        printf("domain->type.type != range->type.type != x->type.type\n");
+        printf("%d %d %d\n", domain->type.type, range->type.type, x->type.type);
     }
     if (domain->n != range->n) {
         // XXX psError
+        printf("domain->n != range->n\n");
     }
     if (order > (domain->n - 1)) {
         // XXX psError: not enough data points for order-order interpolation.
+        printf("not enough data points for order-order interpolation.\n");
     }
 
@@ -1979 +1988 @@
     } else {
         // XXX psError: type not supported
-    }
-
+        printf("XXX psError: type not supported\n");
+    }
+
+    printf("return(NULL)\n");
     psTrace(".psLib.dataManip.psFunctions.p_psVectorInterpolate", 4,
             "---- p_psVectorInterpolate() end ----\n");