diff options
author | rhaas <rhaas@8e189c6b-2ab8-4400-aa02-70a9cfce18b9> | 2014-03-13 03:02:09 +0000 |
---|---|---|
committer | rhaas <rhaas@8e189c6b-2ab8-4400-aa02-70a9cfce18b9> | 2014-03-13 03:02:09 +0000 |
commit | ae4ca16a00a0f48fd60a7431edab6f3250cbde69 (patch) | |
tree | 94bf0e62b4825e3862e2a71a89e61359bd88ddd7 | |
parent | 93d7a40a8f849bdd8a39a8e2e941c94795eb1f50 (diff) |
EOS_Omni nuc_eos_cxx: ensure cs2 is 0 or positive
From: Christian Ott <cott@tapir.caltech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEOS/EOS_Omni/trunk@105 8e189c6b-2ab8-4400-aa02-70a9cfce18b9
-rw-r--r-- | src/nuc_eos_cxx/nuc_eos_press_cs2.cc | 10 |
1 files changed, 8 insertions, 2 deletions
diff --git a/src/nuc_eos_cxx/nuc_eos_press_cs2.cc b/src/nuc_eos_cxx/nuc_eos_press_cs2.cc index 35de771..a741143 100644 --- a/src/nuc_eos_cxx/nuc_eos_press_cs2.cc +++ b/src/nuc_eos_cxx/nuc_eos_press_cs2.cc @@ -63,7 +63,10 @@ void CCTK_FNAME(nuc_eos_m_kt1_press_eps_cs2)(const int *restrict n_in, 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 @@ void CCTK_FNAME(nuc_eos_m_kt0_press_cs2)(const int *restrict n_in, 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 } |