Changeset 31031 for branches/eam_branches/ipp-20110213/psModules/src/objects/models/pmModel_RGAUSS.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_RGAUSS.c (modified) (1 diff)

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

@@ -305 +305 @@
     // choose a z value guaranteed to be beyond our limit
     float z0 = pow((1.0 / limit), (1.0 / PAR[PM_PAR_7]));
-    psAssert (isfinite(z0), "fix this code: z0 should not be nan for %f", PAR[PM_PAR_7]);
+    psAssert (isfinite(z0), "fix this code: z0 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]);
+
     float z1 = (1.0 / limit);
-    psAssert (isfinite(z1), "fix this code: z1 should not be nan for %f", PAR[PM_PAR_7]);
+    psAssert (isfinite(z1), "fix this code: z1 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]);
+
     z1 = PS_MAX (z0, z1);
     z0 = 0.0;