define fundamental derivative G

AMReX-Astro · Aug 2, 2023 · a38ff14 · a38ff14
1 parent 249e109
commit a38ff14
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/EOS/gamma_law/actual_eos.H b/EOS/gamma_law/actual_eos.H
@@ -265,7 +265,7 @@ void actual_eos (I input, T& state)
 
             // sound speed
             state.cs = std::sqrt(eos_gamma * state.p * rhoinv);
-            state.d2pdr2_s = state.cs * state.cs * (std::pow(state.cs, 2) / state.p - rhoinv);
+            state.G = 0.5 * (1.0 + (rho / state.p) * state.cs * state.cs);
         }
     }
 

diff --git a/interfaces/eos_type.H b/interfaces/eos_type.H
@@ -27,7 +27,7 @@ struct eos_t:eos_base_t {
     amrex::Real dsdr{};
     amrex::Real dpde{};
     amrex::Real dpdr_e{};
-    amrex::Real d2pdr2_s{};
+    amrex::Real G{}; // fundmental derivative (Thompson 1971)
 
     amrex::Real cv{};
     amrex::Real cp{};
@@ -146,7 +146,7 @@ struct chem_eos_t:eos_base_t {
     amrex::Real dpdr_e{};
     amrex::Real dedT{};
     amrex::Real dedr{};
-    amrex::Real d2pdr2_s{};
+    amrex::Real G{};
 
     amrex::Real mu{};
     amrex::Real cs{};