Changeset 1440 for trunk/psLib/src/dataManip/psFunctions.c
- Timestamp:
- Aug 9, 2004, 1:34:58 PM (22 years ago)
- File:
-
- 1 edited
-
trunk/psLib/src/dataManip/psFunctions.c (modified) (29 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psLib/src/dataManip/psFunctions.c
r1407 r1440 8 8 * polynomials. It also contains a Gaussian functions. 9 9 * 10 * @version $Revision: 1.2 1$ $Name: not supported by cvs2svn $11 * @date $Date: 2004-08-0 7 00:06:06$10 * @version $Revision: 1.22 $ $Name: not supported by cvs2svn $ 11 * @date $Date: 2004-08-09 23:34:57 $ 12 12 * 13 13 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 49 49 50 50 /*****************************************************************************/ 51 static void polynomial1DFree(psPolynomial1D * myPoly);52 static void polynomial2DFree(psPolynomial2D * myPoly);53 static void polynomial3DFree(psPolynomial3D * myPoly);54 static void polynomial4DFree(psPolynomial4D * myPoly);55 static void dPolynomial1DFree(psDPolynomial1D * myPoly);56 static void dPolynomial2DFree(psDPolynomial2D * myPoly);57 static void dPolynomial3DFree(psDPolynomial3D * myPoly);58 static void dPolynomial4DFree(psDPolynomial4D * myPoly);51 static void polynomial1DFree(psPolynomial1D* myPoly); 52 static void polynomial2DFree(psPolynomial2D* myPoly); 53 static void polynomial3DFree(psPolynomial3D* myPoly); 54 static void polynomial4DFree(psPolynomial4D* myPoly); 55 static void dPolynomial1DFree(psDPolynomial1D* myPoly); 56 static void dPolynomial2DFree(psDPolynomial2D* myPoly); 57 static void dPolynomial3DFree(psDPolynomial3D* myPoly); 58 static void dPolynomial4DFree(psDPolynomial4D* myPoly); 59 59 60 60 /*****************************************************************************/ … … 120 120 NOTE: There is no way to seed the random generator. 121 121 *****************************************************************************/ 122 psVector *psGaussianDev(float mean, float sigma, int Npts)123 { 124 psVector *gauss = NULL;122 psVector* psGaussianDev(float mean, float sigma, int Npts) 123 { 124 psVector* gauss = NULL; 125 125 const gsl_rng_type *T = NULL; 126 126 gsl_rng *r = NULL; … … 144 144 This routine must allocate memory for the polynomial structures. 145 145 *****************************************************************************/ 146 psPolynomial1D *psPolynomial1DAlloc(int n)146 psPolynomial1D* psPolynomial1DAlloc(int n) 147 147 { 148 148 int i = 0; 149 psPolynomial1D *newPoly = NULL;150 151 newPoly = (psPolynomial1D *) psAlloc(sizeof(psPolynomial1D));149 psPolynomial1D* newPoly = NULL; 150 151 newPoly = (psPolynomial1D* ) psAlloc(sizeof(psPolynomial1D)); 152 152 p_psMemSetDeallocator(newPoly, (psFreeFcn) polynomial1DFree); 153 153 newPoly->n = n; … … 164 164 } 165 165 166 psPolynomial2D *psPolynomial2DAlloc(int nX, int nY)167 { 168 int x = 0; 169 int y = 0; 170 psPolynomial2D *newPoly = NULL;171 172 newPoly = (psPolynomial2D *) psAlloc(sizeof(psPolynomial2D));166 psPolynomial2D* psPolynomial2DAlloc(int nX, int nY) 167 { 168 int x = 0; 169 int y = 0; 170 psPolynomial2D* newPoly = NULL; 171 172 newPoly = (psPolynomial2D* ) psAlloc(sizeof(psPolynomial2D)); 173 173 p_psMemSetDeallocator(newPoly, (psFreeFcn) polynomial2DFree); 174 174 newPoly->nX = nX; … … 194 194 } 195 195 196 psPolynomial3D *psPolynomial3DAlloc(int nX, int nY, int nZ)196 psPolynomial3D* psPolynomial3DAlloc(int nX, int nY, int nZ) 197 197 { 198 198 int x = 0; 199 199 int y = 0; 200 200 int z = 0; 201 psPolynomial3D *newPoly = NULL;202 203 newPoly = (psPolynomial3D *) psAlloc(sizeof(psPolynomial3D));201 psPolynomial3D* newPoly = NULL; 202 203 newPoly = (psPolynomial3D* ) psAlloc(sizeof(psPolynomial3D)); 204 204 p_psMemSetDeallocator(newPoly, (psFreeFcn) polynomial3DFree); 205 205 newPoly->nX = nX; … … 233 233 } 234 234 235 psPolynomial4D *psPolynomial4DAlloc(int nW, int nX, int nY, int nZ)235 psPolynomial4D* psPolynomial4DAlloc(int nW, int nX, int nY, int nZ) 236 236 { 237 237 int w = 0; … … 239 239 int y = 0; 240 240 int z = 0; 241 psPolynomial4D *newPoly = NULL;242 243 newPoly = (psPolynomial4D *) psAlloc(sizeof(psPolynomial4D));241 psPolynomial4D* newPoly = NULL; 242 243 newPoly = (psPolynomial4D* ) psAlloc(sizeof(psPolynomial4D)); 244 244 p_psMemSetDeallocator(newPoly, (psFreeFcn) polynomial4DFree); 245 245 newPoly->nW = nW; … … 281 281 } 282 282 283 static void polynomial1DFree(psPolynomial1D * myPoly)283 static void polynomial1DFree(psPolynomial1D* myPoly) 284 284 { 285 285 psFree(myPoly->coeff); … … 288 288 } 289 289 290 static void polynomial2DFree(psPolynomial2D * myPoly)290 static void polynomial2DFree(psPolynomial2D* myPoly) 291 291 { 292 292 int x = 0; … … 302 302 } 303 303 304 static void polynomial3DFree(psPolynomial3D * myPoly)304 static void polynomial3DFree(psPolynomial3D* myPoly) 305 305 { 306 306 int x = 0; … … 323 323 } 324 324 325 static void polynomial4DFree(psPolynomial4D * myPoly)325 static void polynomial4DFree(psPolynomial4D* myPoly) 326 326 { 327 327 int w = 0; … … 353 353 Polynomial coefficients will be accessed in [w][x][y][z] fashion. 354 354 *****************************************************************************/ 355 float psPolynomial1DEval(float x, const psPolynomial1D * myPoly)355 float psPolynomial1DEval(float x, const psPolynomial1D* myPoly) 356 356 { 357 357 int loop_x = 0; … … 380 380 } 381 381 382 float psPolynomial2DEval(float x, float y, const psPolynomial2D * myPoly)382 float psPolynomial2DEval(float x, float y, const psPolynomial2D* myPoly) 383 383 { 384 384 int loop_x = 0; … … 400 400 } 401 401 402 float psPolynomial3DEval(float x, float y, float z, const psPolynomial3D * myPoly)402 float psPolynomial3DEval(float x, float y, float z, const psPolynomial3D* myPoly) 403 403 { 404 404 int loop_x = 0; … … 426 426 } 427 427 428 float psPolynomial4DEval(float w, float x, float y, float z, const psPolynomial4D * myPoly)428 float psPolynomial4DEval(float w, float x, float y, float z, const psPolynomial4D* myPoly) 429 429 { 430 430 int loop_w = 0; … … 458 458 } 459 459 460 psDPolynomial1D *psDPolynomial1DAlloc(int n)460 psDPolynomial1D* psDPolynomial1DAlloc(int n) 461 461 { 462 462 int i = 0; 463 psDPolynomial1D *newPoly = NULL;464 465 newPoly = (psDPolynomial1D *) psAlloc(sizeof(psDPolynomial1D));463 psDPolynomial1D* newPoly = NULL; 464 465 newPoly = (psDPolynomial1D* ) psAlloc(sizeof(psDPolynomial1D)); 466 466 p_psMemSetDeallocator(newPoly, (psFreeFcn) dPolynomial1DFree); 467 467 newPoly->n = n; … … 478 478 } 479 479 480 psDPolynomial2D *psDPolynomial2DAlloc(int nX, int nY)481 { 482 int x = 0; 483 int y = 0; 484 psDPolynomial2D *newPoly = NULL;485 486 newPoly = (psDPolynomial2D *) psAlloc(sizeof(psDPolynomial2D));480 psDPolynomial2D* psDPolynomial2DAlloc(int nX, int nY) 481 { 482 int x = 0; 483 int y = 0; 484 psDPolynomial2D* newPoly = NULL; 485 486 newPoly = (psDPolynomial2D* ) psAlloc(sizeof(psDPolynomial2D)); 487 487 p_psMemSetDeallocator(newPoly, (psFreeFcn) dPolynomial2DFree); 488 488 newPoly->nX = nX; … … 508 508 } 509 509 510 psDPolynomial3D *psDPolynomial3DAlloc(int nX, int nY, int nZ)510 psDPolynomial3D* psDPolynomial3DAlloc(int nX, int nY, int nZ) 511 511 { 512 512 int x = 0; 513 513 int y = 0; 514 514 int z = 0; 515 psDPolynomial3D *newPoly = NULL;516 517 newPoly = (psDPolynomial3D *) psAlloc(sizeof(psDPolynomial3D));515 psDPolynomial3D* newPoly = NULL; 516 517 newPoly = (psDPolynomial3D* ) psAlloc(sizeof(psDPolynomial3D)); 518 518 p_psMemSetDeallocator(newPoly, (psFreeFcn) dPolynomial3DFree); 519 519 newPoly->nX = nX; … … 547 547 } 548 548 549 psDPolynomial4D *psDPolynomial4DAlloc(int nW, int nX, int nY, int nZ)549 psDPolynomial4D* psDPolynomial4DAlloc(int nW, int nX, int nY, int nZ) 550 550 { 551 551 int w = 0; … … 553 553 int y = 0; 554 554 int z = 0; 555 psDPolynomial4D *newPoly = NULL;556 557 newPoly = (psDPolynomial4D *) psAlloc(sizeof(psDPolynomial4D));555 psDPolynomial4D* newPoly = NULL; 556 557 newPoly = (psDPolynomial4D* ) psAlloc(sizeof(psDPolynomial4D)); 558 558 p_psMemSetDeallocator(newPoly, (psFreeFcn) dPolynomial4DFree); 559 559 newPoly->nW = nW; … … 595 595 } 596 596 597 static void dPolynomial1DFree(psDPolynomial1D * myPoly)597 static void dPolynomial1DFree(psDPolynomial1D* myPoly) 598 598 { 599 599 psFree(myPoly->coeff); … … 602 602 } 603 603 604 static void dPolynomial2DFree(psDPolynomial2D * myPoly)604 static void dPolynomial2DFree(psDPolynomial2D* myPoly) 605 605 { 606 606 int x = 0; … … 616 616 } 617 617 618 static void dPolynomial3DFree(psDPolynomial3D * myPoly)618 static void dPolynomial3DFree(psDPolynomial3D* myPoly) 619 619 { 620 620 int x = 0; … … 637 637 } 638 638 639 static void dPolynomial4DFree(psDPolynomial4D * myPoly)639 static void dPolynomial4DFree(psDPolynomial4D* myPoly) 640 640 { 641 641 int w = 0; … … 667 667 Polynomial coefficients will be accessed in [w][x][y][z] fashion. 668 668 *****************************************************************************/ 669 double psDPolynomial1DEval(double x, const psDPolynomial1D * myPoly)669 double psDPolynomial1DEval(double x, const psDPolynomial1D* myPoly) 670 670 { 671 671 int loop_x = 0; … … 686 686 } 687 687 688 double psDPolynomial2DEval(double x, double y, const psDPolynomial2D * myPoly)688 double psDPolynomial2DEval(double x, double y, const psDPolynomial2D* myPoly) 689 689 { 690 690 int loop_x = 0; … … 706 706 } 707 707 708 double psDPolynomial3DEval(double x, double y, double z, const psDPolynomial3D * myPoly)708 double psDPolynomial3DEval(double x, double y, double z, const psDPolynomial3D* myPoly) 709 709 { 710 710 int loop_x = 0; … … 732 732 } 733 733 734 double psDPolynomial4DEval(double w, double x, double y, double z, const psDPolynomial4D * myPoly)734 double psDPolynomial4DEval(double w, double x, double y, double z, const psDPolynomial4D* myPoly) 735 735 { 736 736 int loop_w = 0;
Note:
See TracChangeset
for help on using the changeset viewer.
