Changeset 1889

Show
Ignore:
Timestamp:
03/26/08 11:34:00 (3 months ago)
Author:
peter
Message:

source/punch_do.cpp:

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

Files:

Legend:

Unmodified
Added
Removed
Modified
Copied
Moved

  • trunk/source/punch_do.cpp

    r1852 r1889  
    18501850                                                  stage/MAX2(1e-37,fout),  
    18511851                                                  /* simple H+ */ 
    1852                                                   iso.gamnc[ipH_LIKE][ipHYDROGEN][ipH1s]*dense.xIonDense[ipHYDROGEN][0]/(dense.eden* dense.xIonDense[ipHYDROGEN][1]) , 
     1852                                                  safe_div( iso.gamnc[ipH_LIKE][ipHYDROGEN][ipH1s]*dense.xIonDense[ipHYDROGEN][0], dense.eden*dense.xIonDense[ipHYDROGEN][1] ), 
    18531853                                                  secondaries.csupra[ipHYDROGEN][0]); 
    18541854