Changeset 25683
- Timestamp:
- Sep 29, 2009, 10:10:32 PM (17 years ago)
- Location:
- branches/eam_branches/20090715/psModules/src/objects/models
- Files:
-
- 6 edited
-
pmModel_GAUSS.c (modified) (1 diff)
-
pmModel_PGAUSS.c (modified) (1 diff)
-
pmModel_PS1_V1.c (modified) (1 diff)
-
pmModel_QGAUSS.c (modified) (1 diff)
-
pmModel_RGAUSS.c (modified) (1 diff)
-
pmModel_SERSIC.c (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
branches/eam_branches/20090715/psModules/src/objects/models/pmModel_GAUSS.c
r25679 r25683 260 260 psEllipseAxes axes = psEllipseShapeToAxes (shape, 20.0); 261 261 psF64 radius = axes.major * sqrt (2.0 * log(PAR[PM_PAR_I0] / flux)); 262 psAssert (isfinite(radius , "fix this code: radius should not be nan for %f", PAR[PM_PAR_I0]));262 psAssert (isfinite(radius), "fix this code: radius should not be nan for %f", PAR[PM_PAR_I0]); 263 263 264 264 return (radius); -
branches/eam_branches/20090715/psModules/src/objects/models/pmModel_PGAUSS.c
r25679 r25683 287 287 // choose a z value guaranteed to be beyond our limit 288 288 float z0 = pow((1.0 / limit), (1.0 / 3.0)); 289 psAssert (isfinite(z0 , "fix this code: z0 should not be nan for %f", PAR[PM_PAR_I0]));289 psAssert (isfinite(z0), "fix this code: z0 should not be nan for %f", PAR[PM_PAR_I0]); 290 290 float z1 = (1.0 / limit); 291 psAssert (isfinite(z1 , "fix this code: z1 should not be nan for %f", PAR[PM_PAR_I0]));291 psAssert (isfinite(z1), "fix this code: z1 should not be nan for %f", PAR[PM_PAR_I0]); 292 292 z1 = PS_MAX (z0, z1); 293 293 z0 = 0.0; -
branches/eam_branches/20090715/psModules/src/objects/models/pmModel_PS1_V1.c
r25679 r25683 325 325 float z0 = 0.0; 326 326 float z1 = pow((1.0 / limit), (1.0 / ALPHA)); 327 psAssert (isfinite(z1 , "fix this code: z1 should not be nan for %f", PAR[PM_PAR_7]));327 psAssert (isfinite(z1), "fix this code: z1 should not be nan for %f", PAR[PM_PAR_7]); 328 328 if (PAR[PM_PAR_7] < 0.0) z1 *= 2.0; 329 329 -
branches/eam_branches/20090715/psModules/src/objects/models/pmModel_QGAUSS.c
r25679 r25683 321 321 // choose a z value guaranteed to be beyond our limit 322 322 float z0 = pow((1.0 / limit), (1.0 / 2.25)); 323 psAssert (isfinite(z0 , "fix this code: z0 should not be nan for %f", PAR[PM_PAR_7]));323 psAssert (isfinite(z0), "fix this code: z0 should not be nan for %f", PAR[PM_PAR_7]); 324 324 float z1 = (1.0 / limit) / PAR[PM_PAR_7]; 325 psAssert (isfinite(z1 , "fix this code: z1 should not be nan for %f", PAR[PM_PAR_7]));325 psAssert (isfinite(z1), "fix this code: z1 should not be nan for %f", PAR[PM_PAR_7]); 326 326 z1 = PS_MAX (z0, z1); 327 327 z0 = 0.0; -
branches/eam_branches/20090715/psModules/src/objects/models/pmModel_RGAUSS.c
r25679 r25683 314 314 // choose a z value guaranteed to be beyond our limit 315 315 float z0 = pow((1.0 / limit), (1.0 / PAR[PM_PAR_7])); 316 psAssert (isfinite(z0 , "fix this code: z0 should not be nan for %f", PAR[PM_PAR_7]));316 psAssert (isfinite(z0), "fix this code: z0 should not be nan for %f", PAR[PM_PAR_7]); 317 317 float z1 = (1.0 / limit); 318 psAssert (isfinite(z1 , "fix this code: z1 should not be nan for %f", PAR[PM_PAR_7]));318 psAssert (isfinite(z1), "fix this code: z1 should not be nan for %f", PAR[PM_PAR_7]); 319 319 z1 = PS_MAX (z0, z1); 320 320 z0 = 0.0; -
branches/eam_branches/20090715/psModules/src/objects/models/pmModel_SERSIC.c
r25679 r25683 324 324 325 325 psF64 z = pow (-log(limit), (1.0 / PAR[PM_PAR_7])); 326 psAssert (isfinite(z , "fix this code: z should not be nan for %f", PAR[PM_PAR_7]));326 psAssert (isfinite(z), "fix this code: z should not be nan for %f", PAR[PM_PAR_7]); 327 327 328 328 psF64 radius = sigma * sqrt (2.0 * z); 329 psAssert (isfinite(radius , "fix this code: radius should not be nan for %f, %f", PAR[PM_PAR_7], sigma));329 psAssert (isfinite(radius), "fix this code: radius should not be nan for %f, %f", PAR[PM_PAR_7], sigma); 330 330 331 331 if (isnan(radius))
Note:
See TracChangeset
for help on using the changeset viewer.
