User: rhaas Date: 2014/03/12 10:02 PM
Modified: /trunk/src/nuc_eos_cxx/ nuc_eos_press_cs2.cc
Log: EOS_Omni nuc_eos_cxx: ensure cs2 is 0 or positive
From: Christian Ott cott@tapir.caltech.edu
File Changes:
Directory: /trunk/src/nuc_eos_cxx/ ==================================
File [modified]: nuc_eos_press_cs2.cc Delta lines: +8 -2 =================================================================== --- trunk/src/nuc_eos_cxx/nuc_eos_press_cs2.cc 2014-03-13 03:02:05 UTC (rev 104) +++ trunk/src/nuc_eos_cxx/nuc_eos_press_cs2.cc 2014-03-13 03:02:09 UTC (rev 105) @@ -63,7 +63,10 @@ prs[i] = exp(prs[i]); eps[i] = exp(eps[i]) - energy_shift; #if HAVEGR - cs2[i] = rho[i] * cs2[i] / (rho[i] + rho[i] * eps[i] + prs[i]); + cs2[i] = MAX(rho[i] * cs2[i] / (rho[i] + rho[i] * eps[i] + prs[i]), + 0.0e0); +#else + cs2[i] = MAX(cs2[i], 0.0e0); #endif }
@@ -175,7 +178,10 @@ for(int i=0;i<n;i++) { prs[i] = exp(prs[i]); #if HAVEGR - cs2[i] = rho[i] * cs2[i] / (rho[i] + rho[i]*eps[i] + prs[i]); + cs2[i] = MAX(rho[i] * cs2[i] / (rho[i] + rho[i] * eps[i] + prs[i]), + 0.0e0); +#else + cs2[i] = MAX(cs2[i], 0.0e0); #endif }