Index: trunk/psLib/src/dataManip/psFunctions.c =================================================================== --- trunk/psLib/src/dataManip/psFunctions.c (revision 2431) +++ trunk/psLib/src/dataManip/psFunctions.c (revision 2433) @@ -7,6 +7,6 @@ * polynomials. It also contains a Gaussian functions. * - * @version $Revision: 1.68 $ $Name: not supported by cvs2svn $ - * @date $Date: 2004-11-24 20:34:13 $ + * @version $Revision: 1.69 $ $Name: not supported by cvs2svn $ + * @date $Date: 2004-11-24 21:59:43 $ * * Copyright 2004 Maui High Performance Computing Center, University of Hawaii @@ -399,6 +399,6 @@ if (myPoly->mask[loop_x][loop_y] == 0) { polySum += myPoly->coeff[loop_x][loop_y] * - psPolynomial1DEval(x, chebPolys[loop_x]) * - psPolynomial1DEval(y, chebPolys[loop_y]); + psPolynomial1DEval(chebPolys[loop_x], x) * + psPolynomial1DEval(chebPolys[loop_y], y); } } @@ -465,7 +465,7 @@ if (myPoly->mask[loop_x][loop_y][loop_z] == 0) { polySum += myPoly->coeff[loop_x][loop_y][loop_z] * - psPolynomial1DEval(x, chebPolys[loop_x]) * - psPolynomial1DEval(y, chebPolys[loop_y]) * - psPolynomial1DEval(z, chebPolys[loop_z]); + psPolynomial1DEval(chebPolys[loop_x], x) * + psPolynomial1DEval(chebPolys[loop_y], y) * + psPolynomial1DEval(chebPolys[loop_z], z); } } @@ -545,8 +545,8 @@ if (myPoly->mask[loop_w][loop_x][loop_y][loop_z] == 0) { polySum += myPoly->coeff[loop_w][loop_x][loop_y][loop_z] * - psPolynomial1DEval(w, chebPolys[loop_w]) * - psPolynomial1DEval(x, chebPolys[loop_x]) * - psPolynomial1DEval(y, chebPolys[loop_y]) * - psPolynomial1DEval(z, chebPolys[loop_z]); + psPolynomial1DEval(chebPolys[loop_w], w) * + psPolynomial1DEval(chebPolys[loop_x], x) * + psPolynomial1DEval(chebPolys[loop_y], y) * + psPolynomial1DEval(chebPolys[loop_z], z); } } @@ -651,6 +651,6 @@ if (myPoly->mask[loop_x][loop_y] == 0) { polySum += myPoly->coeff[loop_x][loop_y] * - psPolynomial1DEval(x, chebPolys[loop_x]) * - psPolynomial1DEval(y, chebPolys[loop_y]); + psPolynomial1DEval(chebPolys[loop_x], x) * + psPolynomial1DEval(chebPolys[loop_y], y); } } @@ -718,7 +718,7 @@ if (myPoly->mask[loop_x][loop_y][loop_z] == 0) { polySum += myPoly->coeff[loop_x][loop_y][loop_z] * - psPolynomial1DEval(x, chebPolys[loop_x]) * - psPolynomial1DEval(y, chebPolys[loop_y]) * - psPolynomial1DEval(z, chebPolys[loop_z]); + psPolynomial1DEval(chebPolys[loop_x], x) * + psPolynomial1DEval(chebPolys[loop_y], y) * + psPolynomial1DEval(chebPolys[loop_z], z); } } @@ -798,8 +798,8 @@ if (myPoly->mask[loop_w][loop_x][loop_y][loop_z] == 0) { polySum += myPoly->coeff[loop_w][loop_x][loop_y][loop_z] * - psPolynomial1DEval(w, chebPolys[loop_w]) * - psPolynomial1DEval(x, chebPolys[loop_x]) * - psPolynomial1DEval(y, chebPolys[loop_y]) * - psPolynomial1DEval(z, chebPolys[loop_z]); + psPolynomial1DEval(chebPolys[loop_w], w) * + psPolynomial1DEval(chebPolys[loop_x], x) * + psPolynomial1DEval(chebPolys[loop_y], y) * + psPolynomial1DEval(chebPolys[loop_z], z); } } @@ -1160,5 +1160,5 @@ } -float psPolynomial1DEval(float x, const psPolynomial1D* myPoly) +float psPolynomial1DEval(const psPolynomial1D* myPoly, float x) { PS_POLY_CHECK_NULL(myPoly, NAN); @@ -1176,6 +1176,6 @@ } -psVector *psPolynomial1DEvalVector(const psVector *x, - const psPolynomial1D *myPoly) +psVector *psPolynomial1DEvalVector(const psPolynomial1D *myPoly, + const psVector *x) { PS_POLY_CHECK_NULL(myPoly, NULL); @@ -1189,5 +1189,5 @@ tmp = psVectorAlloc(x->n, PS_TYPE_F32); for (i=0;in;i++) { - tmp->data.F32[i] = psPolynomial1DEval(x->data.F32[i], myPoly); + tmp->data.F32[i] = psPolynomial1DEval(myPoly, x->data.F32[i]); } @@ -1198,5 +1198,5 @@ } -float psPolynomial2DEval(float x, float y, const psPolynomial2D* myPoly) +float psPolynomial2DEval(const psPolynomial2D* myPoly, float x, float y) { PS_POLY_CHECK_NULL(myPoly, NAN); @@ -1214,7 +1214,8 @@ } -psVector *psPolynomial2DEvalVector(const psVector *x, - const psVector *y, - const psPolynomial2D *myPoly) +psVector *psPolynomial2DEvalVector(const psPolynomial2D *myPoly, + const psVector *x, + const psVector *y) + { PS_POLY_CHECK_NULL(myPoly, NULL); @@ -1237,7 +1238,8 @@ tmp = psVectorAlloc(vecLen, PS_TYPE_F32); for (i=0;idata.F32[i] = psPolynomial2DEval(myX->data.F32[i], - myY->data.F32[i], - myPoly); + tmp->data.F32[i] = psPolynomial2DEval(myPoly, + myX->data.F32[i], + myY->data.F32[i]); + } @@ -1251,5 +1253,5 @@ } -float psPolynomial3DEval(float x, float y, float z, const psPolynomial3D* myPoly) +float psPolynomial3DEval(const psPolynomial3D* myPoly, float x, float y, float z) { PS_POLY_CHECK_NULL(myPoly, NAN); @@ -1267,8 +1269,9 @@ } -psVector *psPolynomial3DEvalVector(const psVector *x, +psVector *psPolynomial3DEvalVector(const psPolynomial3D *myPoly, + const psVector *x, const psVector *y, - const psVector *z, - const psPolynomial3D *myPoly) + const psVector *z) + { PS_POLY_CHECK_NULL(myPoly, NULL); @@ -1297,8 +1300,9 @@ tmp = psVectorAlloc(vecLen, PS_TYPE_F32); for (i=0;idata.F32[i] = psPolynomial3DEval(myX->data.F32[i], + tmp->data.F32[i] = psPolynomial3DEval(myPoly, + myX->data.F32[i], myY->data.F32[i], - myZ->data.F32[i], - myPoly); + myZ->data.F32[i]); + } @@ -1315,5 +1319,5 @@ } -float psPolynomial4DEval(float w, float x, float y, float z, const psPolynomial4D* myPoly) +float psPolynomial4DEval(const psPolynomial4D* myPoly, float w, float x, float y, float z) { PS_POLY_CHECK_NULL(myPoly, NAN); @@ -1331,9 +1335,9 @@ } -psVector *psPolynomial4DEvalVector(const psVector *w, +psVector *psPolynomial4DEvalVector(const psPolynomial4D *myPoly, + const psVector *w, const psVector *x, const psVector *y, - const psVector *z, - const psPolynomial4D *myPoly) + const psVector *z) { PS_POLY_CHECK_NULL(myPoly, NULL); @@ -1368,9 +1372,10 @@ tmp = psVectorAlloc(vecLen, PS_TYPE_F32); for (i=0;idata.F32[i] = psPolynomial4DEval(myW->data.F32[i], + tmp->data.F32[i] = psPolynomial4DEval(myPoly, + myW->data.F32[i], myX->data.F32[i], myY->data.F32[i], - myZ->data.F32[i], - myPoly); + myZ->data.F32[i]); + } @@ -1555,5 +1560,5 @@ -double psDPolynomial1DEval(double x, const psDPolynomial1D* myPoly) +double psDPolynomial1DEval(const psDPolynomial1D* myPoly, double x) { PS_POLY_CHECK_NULL(myPoly, NAN); @@ -1571,6 +1576,7 @@ } -psVector *psDPolynomial1DEvalVector(const psVector *x, - const psDPolynomial1D *myPoly) +psVector *psDPolynomial1DEvalVector(const psDPolynomial1D *myPoly, + const psVector *x) + { PS_POLY_CHECK_NULL(myPoly, NULL); @@ -1584,5 +1590,6 @@ tmp = psVectorAlloc(x->n, PS_TYPE_F64); for (i=0;in;i++) { - tmp->data.F64[i] = psDPolynomial1DEval(x->data.F64[i], myPoly); + tmp->data.F64[i] = psDPolynomial1DEval(myPoly, + x->data.F64[i]); } @@ -1594,5 +1601,7 @@ -double psDPolynomial2DEval(double x, double y, const psDPolynomial2D* myPoly) +double psDPolynomial2DEval(const psDPolynomial2D* myPoly, + double x, + double y) { PS_POLY_CHECK_NULL(myPoly, NAN); @@ -1610,7 +1619,7 @@ } -psVector *psDPolynomial2DEvalVector(const psVector *x, - const psVector *y, - const psDPolynomial2D *myPoly) +psVector *psDPolynomial2DEvalVector(const psDPolynomial2D *myPoly, + const psVector *x, + const psVector *y) { PS_POLY_CHECK_NULL(myPoly, NULL); @@ -1633,7 +1642,8 @@ tmp = psVectorAlloc(vecLen, PS_TYPE_F64); for (i=0;idata.F64[i] = psDPolynomial2DEval(myX->data.F64[i], - myY->data.F64[i], - myPoly); + tmp->data.F64[i] = psDPolynomial2DEval(myPoly, + myX->data.F64[i], + myY->data.F64[i]); + } @@ -1648,5 +1658,8 @@ -double psDPolynomial3DEval(double x, double y, double z, const psDPolynomial3D* myPoly) +double psDPolynomial3DEval(const psDPolynomial3D* myPoly, + double x, + double y, + double z) { PS_POLY_CHECK_NULL(myPoly, NAN); @@ -1664,8 +1677,9 @@ } -psVector *psDPolynomial3DEvalVector(const psVector *x, +psVector *psDPolynomial3DEvalVector(const psDPolynomial3D *myPoly, + const psVector *x, const psVector *y, - const psVector *z, - const psDPolynomial3D *myPoly) + const psVector *z) + { PS_POLY_CHECK_NULL(myPoly, NULL); @@ -1694,8 +1708,9 @@ tmp = psVectorAlloc(vecLen, PS_TYPE_F64); for (i=0;idata.F64[i] = psDPolynomial3DEval(myX->data.F64[i], + tmp->data.F64[i] = psDPolynomial3DEval(myPoly, + myX->data.F64[i], myY->data.F64[i], - myZ->data.F64[i], - myPoly); + myZ->data.F64[i]); + } @@ -1712,5 +1727,9 @@ } -double psDPolynomial4DEval(double w, double x, double y, double z, const psDPolynomial4D* myPoly) +double psDPolynomial4DEval(const psDPolynomial4D* myPoly, + double w, + double x, + double y, + double z) { PS_POLY_CHECK_NULL(myPoly, NAN); @@ -1728,9 +1747,9 @@ } -psVector *psDPolynomial4DEvalVector(const psVector *w, +psVector *psDPolynomial4DEvalVector(const psDPolynomial4D *myPoly, + const psVector *w, const psVector *x, const psVector *y, - const psVector *z, - const psDPolynomial4D *myPoly) + const psVector *z) { PS_POLY_CHECK_NULL(myPoly, NULL); @@ -1765,9 +1784,10 @@ tmp = psVectorAlloc(vecLen, PS_TYPE_F64); for (i=0;idata.F64[i] = psDPolynomial4DEval(myW->data.F64[i], + tmp->data.F64[i] = psDPolynomial4DEval(myPoly, + myW->data.F64[i], myX->data.F64[i], myY->data.F64[i], - myZ->data.F64[i], - myPoly); + myZ->data.F64[i]); + } @@ -2110,11 +2130,14 @@ if (x < (spline->domains)[0]) { - return(psPolynomial1DEval(x, spline->spline[0])); + return(psPolynomial1DEval(spline->spline[0], + x)); } else if (x > (spline->domains)[n-1]) { - return(psPolynomial1DEval(x, spline->spline[n-1])); - } - } - - return(psPolynomial1DEval(x, spline->spline[binNum])); + return(psPolynomial1DEval(spline->spline[n-1], + x)); + } + } + + return(psPolynomial1DEval(spline->spline[binNum], + x)); }