Changeset 31031 for branches/eam_branches/ipp-20110213/psModules/src/objects/models/pmModel_SERSIC.c

Timestamp:

Mar 24, 2011, 9:22:53 AM (15 years ago)

Author:

eugene

Message:

better information on the radius assertion failure

File:

• branches/eam_branches/ipp-20110213/psModules/src/objects/models/pmModel_SERSIC.c (modified) (1 diff)

branches/eam_branches/ipp-20110213/psModules/src/objects/models/pmModel_SERSIC.c

@@ -346 +346 @@
     psF64 radius = axes.major * sqrt (2.0) * pow(zn, 0.5 / PAR[PM_PAR_7]);
 
-    // fprintf (stderr, "sersic model %f %f, n %f, radius: %f, zn: %f, f/Io: %f, major: %f\n", PAR[PM_PAR_XPOS], PAR[PM_PAR_YPOS], PAR[PM_PAR_7], radius, zn, flux/PAR[PM_PAR_I0], axes.major);
-
-    psAssert (isfinite(radius), "fix this code: z should not be nan for %f", PAR[PM_PAR_7]);
+    psAssert (isfinite(radius), "fix this code: radius should not be nan for Io = %f, flux = %f, major = %f (%f, %f, %f), par 7 = %f", 
+              PAR[PM_PAR_I0], flux, axes.major, PAR[PM_PAR_SXX], PAR[PM_PAR_SXY], PAR[PM_PAR_SYY], PAR[PM_PAR_7]);
     return (radius);
 }