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

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

@@ -315 +315 @@
     float z0 = 0.0;
     float z1 = pow((1.0 / limit), (1.0 / ALPHA));
-    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[PM_PAR_I0], flux, axes.major, PAR[PM_PAR_SXX], PAR[PM_PAR_SXY], PAR[PM_PAR_SYY]);
     if (PAR[PM_PAR_7] < 0.0) z1 *= 2.0;