Changeset 8245 for trunk/psLib/src/math/psSpline.c
- Timestamp:
- Aug 8, 2006, 4:26:44 PM (20 years ago)
- File:
-
- 1 edited
-
trunk/psLib/src/math/psSpline.c (modified) (13 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psLib/src/math/psSpline.c
r8232 r8245 6 6 * This file contains the routines that allocate, free, and evaluate splines. 7 7 * 8 * @version $Revision: 1.15 1$ $Name: not supported by cvs2svn $9 * @date $Date: 2006-08-0 8 23:32:23$8 * @version $Revision: 1.152 $ $Name: not supported by cvs2svn $ 9 * @date $Date: 2006-08-09 02:26:44 $ 10 10 * 11 11 * Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii … … 119 119 const psVector* y) ///< Coordinates 120 120 { 121 psTrace( __func__, 4, "---- %s() begin ----\n", __func__);121 psTrace("psLib.math", 4, "---- %s() begin ----\n", __func__); 122 122 if (psTraceGetLevel(__func__) >= 6) { 123 123 p_psVectorPrint(1, (psVector *) x, "x"); … … 143 143 u[i] = ((6.0 * u[i] / (X[i+1] - X[i-1])) - (sig * u[i-1])) / p; 144 144 145 psTrace( __func__, 6, "X[%d] is %f\n", i, X[i]);146 psTrace( __func__, 6, "Y[%d] is %f\n", i, Y[i]);147 psTrace( __func__, 6, "u[%d] is %f\n", i, u[i]);145 psTrace("psLib.math", 6, "X[%d] is %f\n", i, X[i]); 146 psTrace("psLib.math", 6, "Y[%d] is %f\n", i, Y[i]); 147 psTrace("psLib.math", 6, "u[%d] is %f\n", i, u[i]); 148 148 } 149 149 … … 154 154 for (psS32 k=(n-2);k>=0;k--) { 155 155 derivs2[k] = derivs2[k] * derivs2[k+1] + u[k]; 156 psTrace( __func__, 6, "derivs2[%d] is %f\n", k, derivs2[k]);156 psTrace("psLib.math", 6, "derivs2[%d] is %f\n", k, derivs2[k]); 157 157 } 158 158 psFree(u); 159 psTrace( __func__, 4, "---- %s() end ----\n", __func__);159 psTrace("psLib.math", 4, "---- %s() end ----\n", __func__); 160 160 return(derivs2); 161 161 } … … 210 210 const psVector* y) ///< Coordinates. 211 211 { 212 psTrace( __func__, 3, "---- %s() begin ----\n", __func__);212 psTrace("psLib.math", 3, "---- %s() begin ----\n", __func__); 213 213 PS_ASSERT_VECTOR_NON_NULL(y, NULL); 214 214 PS_ASSERT_VECTOR_TYPE_F32_OR_F64(y, NULL); 215 215 PS_ASSERT_INT_LARGER_THAN_OR_EQUAL(y->n, 2, NULL); 216 216 psS32 numSplines = (y->n)-1; 217 psTrace( __func__, 5, "numSplines is %d\n", numSplines);217 psTrace("psLib.math", 5, "numSplines is %d\n", numSplines); 218 218 219 219 // … … 281 281 i, i+1, xPtr->data.F32[i], xPtr->data.F32[i+1]); 282 282 } 283 psTrace( __func__, 6, "x data (%f - %f) (%f)\n", xPtr->data.F32[i], xPtr->data.F32[i+1], H);283 psTrace("psLib.math", 6, "x data (%f - %f) (%f)\n", xPtr->data.F32[i], xPtr->data.F32[i+1], H); 284 284 // 285 285 // ******** Calculate 0-order term ******** … … 334 334 spline->spline[i]->coeff[3]+= spline->p_psDeriv2[i+1] / (6.0 * H); 335 335 336 psTrace( __func__, 6, "(spline->spline[%u])->coeff[0] is %f\n", i, spline->spline[i]->coeff[0]);337 psTrace( __func__, 6, "(spline->spline[%u])->coeff[1] is %f\n", i, spline->spline[i]->coeff[1]);338 psTrace( __func__, 6, "(spline->spline[%u])->coeff[2] is %f\n", i, spline->spline[i]->coeff[2]);339 psTrace( __func__, 6, "(spline->spline[%u])->coeff[3] is %f\n", i, spline->spline[i]->coeff[3]);336 psTrace("psLib.math", 6, "(spline->spline[%u])->coeff[0] is %f\n", i, spline->spline[i]->coeff[0]); 337 psTrace("psLib.math", 6, "(spline->spline[%u])->coeff[1] is %f\n", i, spline->spline[i]->coeff[1]); 338 psTrace("psLib.math", 6, "(spline->spline[%u])->coeff[2] is %f\n", i, spline->spline[i]->coeff[2]); 339 psTrace("psLib.math", 6, "(spline->spline[%u])->coeff[3] is %f\n", i, spline->spline[i]->coeff[3]); 340 340 } 341 341 … … 343 343 psFree(yPtr); 344 344 } 345 psTrace( __func__, 3, "---- %s() end ----\n", __func__);345 psTrace("psLib.math", 3, "---- %s() end ----\n", __func__); 346 346 return(spline); 347 347 } … … 365 365 float x) 366 366 { 367 psTrace( __func__, 3, "---- %s() begin ----\n", __func__);367 psTrace("psLib.math", 3, "---- %s() begin ----\n", __func__); 368 368 PS_ASSERT_PTR_NON_NULL(spline, NAN); 369 369 PS_ASSERT_INT_NONNEGATIVE(spline->n, NAN); … … 379 379 380 380 psS32 binNum = (x < spline->knots->data.F32[0]) ? 0 : n-1; 381 psTrace( __func__, 3, "---- %s() end ----\n", __func__);381 psTrace("psLib.math", 3, "---- %s() end ----\n", __func__); 382 382 return(psPolynomial1DEval(spline->spline[binNum], x)); 383 383 } … … 392 392 } 393 393 394 psTrace( __func__, 3, "---- %s() end ----\n", __func__);394 psTrace("psLib.math", 3, "---- %s() end ----\n", __func__); 395 395 return(psPolynomial1DEval(spline->spline[binNum], x)); 396 396 } … … 403 403 const psVector *x) 404 404 { 405 psTrace( __func__, 3, "---- %s() begin ----\n", __func__);405 psTrace("psLib.math", 3, "---- %s() begin ----\n", __func__); 406 406 PS_ASSERT_PTR_NON_NULL(spline, NULL); 407 407 PS_ASSERT_VECTOR_TYPE(spline->knots, PS_TYPE_F32, NULL); … … 426 426 } 427 427 428 psTrace( __func__, 3, "---- %s() end ----\n", __func__);428 psTrace("psLib.math", 3, "---- %s() end ----\n", __func__); 429 429 return(tmpVector); 430 430 }
Note:
See TracChangeset
for help on using the changeset viewer.
