Changeset 1341
- Timestamp:
- Jul 29, 2004, 12:44:39 PM (22 years ago)
- Location:
- trunk/psLib/src
- Files:
-
- 10 edited
-
dataManip/psFunctions.c (modified) (11 diffs)
-
dataManip/psFunctions.h (modified) (9 diffs)
-
dataManip/psStats.c (modified) (1 diff)
-
image/psImageStats.c (modified) (4 diffs)
-
imageops/psImageStats.c (modified) (4 diffs)
-
math/psPolynomial.c (modified) (11 diffs)
-
math/psPolynomial.h (modified) (9 diffs)
-
math/psSpline.c (modified) (11 diffs)
-
math/psSpline.h (modified) (9 diffs)
-
math/psStats.c (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psLib/src/dataManip/psFunctions.c
r1281 r1341 7 7 * polynomials. It also contains a Gaussian functions. 8 8 * 9 * @version $Revision: 1.1 6$ $Name: not supported by cvs2svn $10 * @date $Date: 2004-07-2 2 23:45:16$9 * @version $Revision: 1.17 $ $Name: not supported by cvs2svn $ 10 * @date $Date: 2004-07-29 22:44:34 $ 11 11 * 12 12 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 344 344 *****************************************************************************/ 345 345 float 346 ps EvalPolynomial1D(float x,346 psPolynomial1DEval(float x, 347 347 const psPolynomial1D *myPoly) 348 348 { … … 352 352 353 353 if (NULL == myPoly) { 354 psAbort(__func__, "ps EvalPolynomial1D(): myPoly is NULL\n");354 psAbort(__func__, "psPolynomial1DEval(): myPoly is NULL\n"); 355 355 } 356 356 … … 362 362 363 363 if (NULL == myPoly->coeff) { 364 psAbort(__func__, "ps EvalPolynomial1D(): myPoly->coeff is NULL\n");364 psAbort(__func__, "psPolynomial1DEval(): myPoly->coeff is NULL\n"); 365 365 } 366 366 … … 375 375 376 376 float 377 ps EvalPolynomial2D(float x,377 psPolynomial2DEval(float x, 378 378 float y, 379 379 const psPolynomial2D *myPoly) … … 398 398 399 399 float 400 ps EvalPolynomial3D(float x,400 psPolynomial3DEval(float x, 401 401 float y, 402 402 float z, … … 428 428 429 429 float 430 ps EvalPolynomial4D(float w,430 psPolynomial4DEval(float w, 431 431 float x, 432 432 float y, … … 678 678 *****************************************************************************/ 679 679 double 680 psD EvalPolynomial1D(double x,680 psDPolynomial1DEval(double x, 681 681 const psDPolynomial1D *myPoly) 682 682 { … … 700 700 701 701 double 702 psD EvalPolynomial2D(double x,702 psDPolynomial2DEval(double x, 703 703 double y, 704 704 const psDPolynomial2D *myPoly) … … 723 723 724 724 double 725 psD EvalPolynomial3D(double x,725 psDPolynomial3DEval(double x, 726 726 double y, 727 727 double z, … … 753 753 754 754 double 755 psD EvalPolynomial4D(double w,755 psDPolynomial4DEval(double w, 756 756 double x, 757 757 double y, -
trunk/psLib/src/dataManip/psFunctions.h
r1298 r1341 12 12 * @author George Gusciora, MHPCC 13 13 * 14 * @version $Revision: 1. 9$ $Name: not supported by cvs2svn $15 * @date $Date: 2004-07-2 7 19:51:55$14 * @version $Revision: 1.10 $ $Name: not supported by cvs2svn $ 15 * @date $Date: 2004-07-29 22:44:34 $ 16 16 * 17 17 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 106 106 /** Evaluate 1D polynomial */ 107 107 float 108 ps EvalPolynomial1D(float x, ///< Value at which to evaluate108 psPolynomial1DEval(float x, ///< Value at which to evaluate 109 109 const psPolynomial1D *myPoly ///< Coefficients for the polynomial 110 110 ); … … 112 112 /** Evaluate 2D polynomial */ 113 113 float 114 ps EvalPolynomial2D(float x, ///< Value x at which to evaluate114 psPolynomial2DEval(float x, ///< Value x at which to evaluate 115 115 float y, ///< Value y at which to evaluate 116 116 const psPolynomial2D *myPoly ///< Coefficients for the polynomial … … 119 119 /** Evaluate 3D polynomial */ 120 120 float 121 ps EvalPolynomial3D(float x, ///< Value x at which to evaluate121 psPolynomial3DEval(float x, ///< Value x at which to evaluate 122 122 float y, ///< Value y at which to evaluate 123 123 float z, ///< Value z at which to evaluate … … 127 127 /** Evaluate 4D polynomial */ 128 128 float 129 ps EvalPolynomial4D(float w, ///< Value w at which to evaluate129 psPolynomial4DEval(float w, ///< Value w at which to evaluate 130 130 float x, ///< Value x at which to evaluate 131 131 float y, ///< Value y at which to evaluate … … 193 193 /** Evaluate 1D polynomial (double precision) */ 194 194 double 195 psD EvalPolynomial1D(double x, ///< Value at which to evaluate195 psDPolynomial1DEval(double x, ///< Value at which to evaluate 196 196 const psDPolynomial1D *myPoly ///< Coefficients for the polynomial 197 197 ); … … 199 199 /** Evaluate 2D polynomial (double precision) */ 200 200 double 201 psD EvalPolynomial2D(double x, ///< Value x at which to evaluate201 psDPolynomial2DEval(double x, ///< Value x at which to evaluate 202 202 double y, ///< Value y at which to evaluate 203 203 const psDPolynomial2D *myPoly ///< Coefficients for the polynomial … … 206 206 /** Evaluate 3D polynomial (double precision) */ 207 207 double 208 psD EvalPolynomial3D(double x, ///< Value x at which to evaluate208 psDPolynomial3DEval(double x, ///< Value x at which to evaluate 209 209 double y, ///< Value y at which to evaluate 210 210 double z, ///< Value z at which to evaluate … … 214 214 /** Evaluate 4D polynomial (double precision) */ 215 215 double 216 psD EvalPolynomial4D(double w, ///< Value w at which to evaluate216 psDPolynomial4DEval(double w, ///< Value w at which to evaluate 217 217 double x, ///< Value x at which to evaluate 218 218 double y, ///< Value y at which to evaluate -
trunk/psLib/src/dataManip/psStats.c
r1336 r1341 1049 1049 oldMidpoint = midpoint; 1050 1050 1051 f = ps EvalPolynomial1D(midpoint, myPoly);1051 f = psPolynomial1DEval(midpoint, myPoly); 1052 1052 // printf("p_ps1DPolyMedian() iteration %d. f(%f) is %f\n", numIterations, midpoint, f); 1053 1053 if (fabs(f - getThisValue) <= FLT_EPSILON) { -
trunk/psLib/src/image/psImageStats.c
r1323 r1341 200 200 sums[i][j]: This will contain the sum of 201 201 input->data.F32[x][y] * 202 ps EvalPolynomial1D((float) x, chebPolys[i]) *203 ps EvalPolynomial1D((float) y, chebPolys[j]);202 psPolynomial1DEval((float) x, chebPolys[i]) * 203 psPolynomial1DEval((float) y, chebPolys[j]); 204 204 over all pixels (x,y) in the image. 205 205 *****************************************************************************/ … … 242 242 tmp = 0.0; 243 243 for (x=0;x<input->numRows;x++) { 244 tmp+= ps EvalPolynomial1D(rScalingFactors[x], chebPolys[i]) *245 ps EvalPolynomial1D(rScalingFactors[x], chebPolys[j]);244 tmp+= psPolynomial1DEval(rScalingFactors[x], chebPolys[i]) * 245 psPolynomial1DEval(rScalingFactors[x], chebPolys[j]); 246 246 247 247 } … … 257 257 for (y=0;y<input->numCols;y++) { 258 258 sums[i][j]+= input->data.F32[x][y] * 259 ps EvalPolynomial1D(rScalingFactors[x], chebPolys[i]) *260 ps EvalPolynomial1D(cScalingFactors[y], chebPolys[j]);259 psPolynomial1DEval(rScalingFactors[x], chebPolys[i]) * 260 psPolynomial1DEval(cScalingFactors[y], chebPolys[j]); 261 261 } 262 262 } … … 343 343 for (i=0;i<coeffs->nX;i++) { 344 344 for (j=0;j<coeffs->nY;j++) { 345 polySum+= ps EvalPolynomial1D(rScalingFactors[x], chebPolys[i]) *346 ps EvalPolynomial1D(cScalingFactors[y], chebPolys[j]) *345 polySum+= psPolynomial1DEval(rScalingFactors[x], chebPolys[i]) * 346 psPolynomial1DEval(cScalingFactors[y], chebPolys[j]) * 347 347 coeffs->coeff[i][j]; 348 348 -
trunk/psLib/src/imageops/psImageStats.c
r1323 r1341 200 200 sums[i][j]: This will contain the sum of 201 201 input->data.F32[x][y] * 202 ps EvalPolynomial1D((float) x, chebPolys[i]) *203 ps EvalPolynomial1D((float) y, chebPolys[j]);202 psPolynomial1DEval((float) x, chebPolys[i]) * 203 psPolynomial1DEval((float) y, chebPolys[j]); 204 204 over all pixels (x,y) in the image. 205 205 *****************************************************************************/ … … 242 242 tmp = 0.0; 243 243 for (x=0;x<input->numRows;x++) { 244 tmp+= ps EvalPolynomial1D(rScalingFactors[x], chebPolys[i]) *245 ps EvalPolynomial1D(rScalingFactors[x], chebPolys[j]);244 tmp+= psPolynomial1DEval(rScalingFactors[x], chebPolys[i]) * 245 psPolynomial1DEval(rScalingFactors[x], chebPolys[j]); 246 246 247 247 } … … 257 257 for (y=0;y<input->numCols;y++) { 258 258 sums[i][j]+= input->data.F32[x][y] * 259 ps EvalPolynomial1D(rScalingFactors[x], chebPolys[i]) *260 ps EvalPolynomial1D(cScalingFactors[y], chebPolys[j]);259 psPolynomial1DEval(rScalingFactors[x], chebPolys[i]) * 260 psPolynomial1DEval(cScalingFactors[y], chebPolys[j]); 261 261 } 262 262 } … … 343 343 for (i=0;i<coeffs->nX;i++) { 344 344 for (j=0;j<coeffs->nY;j++) { 345 polySum+= ps EvalPolynomial1D(rScalingFactors[x], chebPolys[i]) *346 ps EvalPolynomial1D(cScalingFactors[y], chebPolys[j]) *345 polySum+= psPolynomial1DEval(rScalingFactors[x], chebPolys[i]) * 346 psPolynomial1DEval(cScalingFactors[y], chebPolys[j]) * 347 347 coeffs->coeff[i][j]; 348 348 -
trunk/psLib/src/math/psPolynomial.c
r1281 r1341 7 7 * polynomials. It also contains a Gaussian functions. 8 8 * 9 * @version $Revision: 1.1 6$ $Name: not supported by cvs2svn $10 * @date $Date: 2004-07-2 2 23:45:16$9 * @version $Revision: 1.17 $ $Name: not supported by cvs2svn $ 10 * @date $Date: 2004-07-29 22:44:34 $ 11 11 * 12 12 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 344 344 *****************************************************************************/ 345 345 float 346 ps EvalPolynomial1D(float x,346 psPolynomial1DEval(float x, 347 347 const psPolynomial1D *myPoly) 348 348 { … … 352 352 353 353 if (NULL == myPoly) { 354 psAbort(__func__, "ps EvalPolynomial1D(): myPoly is NULL\n");354 psAbort(__func__, "psPolynomial1DEval(): myPoly is NULL\n"); 355 355 } 356 356 … … 362 362 363 363 if (NULL == myPoly->coeff) { 364 psAbort(__func__, "ps EvalPolynomial1D(): myPoly->coeff is NULL\n");364 psAbort(__func__, "psPolynomial1DEval(): myPoly->coeff is NULL\n"); 365 365 } 366 366 … … 375 375 376 376 float 377 ps EvalPolynomial2D(float x,377 psPolynomial2DEval(float x, 378 378 float y, 379 379 const psPolynomial2D *myPoly) … … 398 398 399 399 float 400 ps EvalPolynomial3D(float x,400 psPolynomial3DEval(float x, 401 401 float y, 402 402 float z, … … 428 428 429 429 float 430 ps EvalPolynomial4D(float w,430 psPolynomial4DEval(float w, 431 431 float x, 432 432 float y, … … 678 678 *****************************************************************************/ 679 679 double 680 psD EvalPolynomial1D(double x,680 psDPolynomial1DEval(double x, 681 681 const psDPolynomial1D *myPoly) 682 682 { … … 700 700 701 701 double 702 psD EvalPolynomial2D(double x,702 psDPolynomial2DEval(double x, 703 703 double y, 704 704 const psDPolynomial2D *myPoly) … … 723 723 724 724 double 725 psD EvalPolynomial3D(double x,725 psDPolynomial3DEval(double x, 726 726 double y, 727 727 double z, … … 753 753 754 754 double 755 psD EvalPolynomial4D(double w,755 psDPolynomial4DEval(double w, 756 756 double x, 757 757 double y, -
trunk/psLib/src/math/psPolynomial.h
r1298 r1341 12 12 * @author George Gusciora, MHPCC 13 13 * 14 * @version $Revision: 1. 9$ $Name: not supported by cvs2svn $15 * @date $Date: 2004-07-2 7 19:51:55$14 * @version $Revision: 1.10 $ $Name: not supported by cvs2svn $ 15 * @date $Date: 2004-07-29 22:44:34 $ 16 16 * 17 17 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 106 106 /** Evaluate 1D polynomial */ 107 107 float 108 ps EvalPolynomial1D(float x, ///< Value at which to evaluate108 psPolynomial1DEval(float x, ///< Value at which to evaluate 109 109 const psPolynomial1D *myPoly ///< Coefficients for the polynomial 110 110 ); … … 112 112 /** Evaluate 2D polynomial */ 113 113 float 114 ps EvalPolynomial2D(float x, ///< Value x at which to evaluate114 psPolynomial2DEval(float x, ///< Value x at which to evaluate 115 115 float y, ///< Value y at which to evaluate 116 116 const psPolynomial2D *myPoly ///< Coefficients for the polynomial … … 119 119 /** Evaluate 3D polynomial */ 120 120 float 121 ps EvalPolynomial3D(float x, ///< Value x at which to evaluate121 psPolynomial3DEval(float x, ///< Value x at which to evaluate 122 122 float y, ///< Value y at which to evaluate 123 123 float z, ///< Value z at which to evaluate … … 127 127 /** Evaluate 4D polynomial */ 128 128 float 129 ps EvalPolynomial4D(float w, ///< Value w at which to evaluate129 psPolynomial4DEval(float w, ///< Value w at which to evaluate 130 130 float x, ///< Value x at which to evaluate 131 131 float y, ///< Value y at which to evaluate … … 193 193 /** Evaluate 1D polynomial (double precision) */ 194 194 double 195 psD EvalPolynomial1D(double x, ///< Value at which to evaluate195 psDPolynomial1DEval(double x, ///< Value at which to evaluate 196 196 const psDPolynomial1D *myPoly ///< Coefficients for the polynomial 197 197 ); … … 199 199 /** Evaluate 2D polynomial (double precision) */ 200 200 double 201 psD EvalPolynomial2D(double x, ///< Value x at which to evaluate201 psDPolynomial2DEval(double x, ///< Value x at which to evaluate 202 202 double y, ///< Value y at which to evaluate 203 203 const psDPolynomial2D *myPoly ///< Coefficients for the polynomial … … 206 206 /** Evaluate 3D polynomial (double precision) */ 207 207 double 208 psD EvalPolynomial3D(double x, ///< Value x at which to evaluate208 psDPolynomial3DEval(double x, ///< Value x at which to evaluate 209 209 double y, ///< Value y at which to evaluate 210 210 double z, ///< Value z at which to evaluate … … 214 214 /** Evaluate 4D polynomial (double precision) */ 215 215 double 216 psD EvalPolynomial4D(double w, ///< Value w at which to evaluate216 psDPolynomial4DEval(double w, ///< Value w at which to evaluate 217 217 double x, ///< Value x at which to evaluate 218 218 double y, ///< Value y at which to evaluate -
trunk/psLib/src/math/psSpline.c
r1281 r1341 7 7 * polynomials. It also contains a Gaussian functions. 8 8 * 9 * @version $Revision: 1.1 6$ $Name: not supported by cvs2svn $10 * @date $Date: 2004-07-2 2 23:45:16$9 * @version $Revision: 1.17 $ $Name: not supported by cvs2svn $ 10 * @date $Date: 2004-07-29 22:44:34 $ 11 11 * 12 12 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 344 344 *****************************************************************************/ 345 345 float 346 ps EvalPolynomial1D(float x,346 psPolynomial1DEval(float x, 347 347 const psPolynomial1D *myPoly) 348 348 { … … 352 352 353 353 if (NULL == myPoly) { 354 psAbort(__func__, "ps EvalPolynomial1D(): myPoly is NULL\n");354 psAbort(__func__, "psPolynomial1DEval(): myPoly is NULL\n"); 355 355 } 356 356 … … 362 362 363 363 if (NULL == myPoly->coeff) { 364 psAbort(__func__, "ps EvalPolynomial1D(): myPoly->coeff is NULL\n");364 psAbort(__func__, "psPolynomial1DEval(): myPoly->coeff is NULL\n"); 365 365 } 366 366 … … 375 375 376 376 float 377 ps EvalPolynomial2D(float x,377 psPolynomial2DEval(float x, 378 378 float y, 379 379 const psPolynomial2D *myPoly) … … 398 398 399 399 float 400 ps EvalPolynomial3D(float x,400 psPolynomial3DEval(float x, 401 401 float y, 402 402 float z, … … 428 428 429 429 float 430 ps EvalPolynomial4D(float w,430 psPolynomial4DEval(float w, 431 431 float x, 432 432 float y, … … 678 678 *****************************************************************************/ 679 679 double 680 psD EvalPolynomial1D(double x,680 psDPolynomial1DEval(double x, 681 681 const psDPolynomial1D *myPoly) 682 682 { … … 700 700 701 701 double 702 psD EvalPolynomial2D(double x,702 psDPolynomial2DEval(double x, 703 703 double y, 704 704 const psDPolynomial2D *myPoly) … … 723 723 724 724 double 725 psD EvalPolynomial3D(double x,725 psDPolynomial3DEval(double x, 726 726 double y, 727 727 double z, … … 753 753 754 754 double 755 psD EvalPolynomial4D(double w,755 psDPolynomial4DEval(double w, 756 756 double x, 757 757 double y, -
trunk/psLib/src/math/psSpline.h
r1298 r1341 12 12 * @author George Gusciora, MHPCC 13 13 * 14 * @version $Revision: 1. 9$ $Name: not supported by cvs2svn $15 * @date $Date: 2004-07-2 7 19:51:55$14 * @version $Revision: 1.10 $ $Name: not supported by cvs2svn $ 15 * @date $Date: 2004-07-29 22:44:34 $ 16 16 * 17 17 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 106 106 /** Evaluate 1D polynomial */ 107 107 float 108 ps EvalPolynomial1D(float x, ///< Value at which to evaluate108 psPolynomial1DEval(float x, ///< Value at which to evaluate 109 109 const psPolynomial1D *myPoly ///< Coefficients for the polynomial 110 110 ); … … 112 112 /** Evaluate 2D polynomial */ 113 113 float 114 ps EvalPolynomial2D(float x, ///< Value x at which to evaluate114 psPolynomial2DEval(float x, ///< Value x at which to evaluate 115 115 float y, ///< Value y at which to evaluate 116 116 const psPolynomial2D *myPoly ///< Coefficients for the polynomial … … 119 119 /** Evaluate 3D polynomial */ 120 120 float 121 ps EvalPolynomial3D(float x, ///< Value x at which to evaluate121 psPolynomial3DEval(float x, ///< Value x at which to evaluate 122 122 float y, ///< Value y at which to evaluate 123 123 float z, ///< Value z at which to evaluate … … 127 127 /** Evaluate 4D polynomial */ 128 128 float 129 ps EvalPolynomial4D(float w, ///< Value w at which to evaluate129 psPolynomial4DEval(float w, ///< Value w at which to evaluate 130 130 float x, ///< Value x at which to evaluate 131 131 float y, ///< Value y at which to evaluate … … 193 193 /** Evaluate 1D polynomial (double precision) */ 194 194 double 195 psD EvalPolynomial1D(double x, ///< Value at which to evaluate195 psDPolynomial1DEval(double x, ///< Value at which to evaluate 196 196 const psDPolynomial1D *myPoly ///< Coefficients for the polynomial 197 197 ); … … 199 199 /** Evaluate 2D polynomial (double precision) */ 200 200 double 201 psD EvalPolynomial2D(double x, ///< Value x at which to evaluate201 psDPolynomial2DEval(double x, ///< Value x at which to evaluate 202 202 double y, ///< Value y at which to evaluate 203 203 const psDPolynomial2D *myPoly ///< Coefficients for the polynomial … … 206 206 /** Evaluate 3D polynomial (double precision) */ 207 207 double 208 psD EvalPolynomial3D(double x, ///< Value x at which to evaluate208 psDPolynomial3DEval(double x, ///< Value x at which to evaluate 209 209 double y, ///< Value y at which to evaluate 210 210 double z, ///< Value z at which to evaluate … … 214 214 /** Evaluate 4D polynomial (double precision) */ 215 215 double 216 psD EvalPolynomial4D(double w, ///< Value w at which to evaluate216 psDPolynomial4DEval(double w, ///< Value w at which to evaluate 217 217 double x, ///< Value x at which to evaluate 218 218 double y, ///< Value y at which to evaluate -
trunk/psLib/src/math/psStats.c
r1336 r1341 1049 1049 oldMidpoint = midpoint; 1050 1050 1051 f = ps EvalPolynomial1D(midpoint, myPoly);1051 f = psPolynomial1DEval(midpoint, myPoly); 1052 1052 // printf("p_ps1DPolyMedian() iteration %d. f(%f) is %f\n", numIterations, midpoint, f); 1053 1053 if (fabs(f - getThisValue) <= FLT_EPSILON) {
Note:
See TracChangeset
for help on using the changeset viewer.
