Changeset 1890

Timestamp:

03/26/08 11:42:25 (8 months ago)

Author:

peter

Message:

source/punch_do.cpp:

Bug-fix for PR61 - division by zero in the PUNCH HYDROGEN IONIZATION command.

Files:

• branches/c07.02_branch/source/punch_do.cpp (modified) (1 diff)

branches/c07.02_branch/source/punch_do.cpp

@@ -1692 +1692 @@
                                           stage/MAX2(1e-37,fout), 
                                           /* simple H+ */
-                                          iso.gamnc[ipH_LIKE][ipHYDROGEN][ipH1s]*dense.xIonDense[ipHYDROGEN][0]/(dense.eden* dense.xIonDense[ipHYDROGEN][1]) ,
+                                          iso.gamnc[ipH_LIKE][ipHYDROGEN][ipH1s]*dense.xIonDense[ipHYDROGEN][0]/SDIV(dense.eden*dense.xIonDense[ipHYDROGEN][1]),
                                           secondaries.csupra[ipHYDROGEN][0]);