Merge pull request #102 from LLNL/bug-fix/GSPH-unitsError

Bug fix/gsph units error
LLNL · Dec 10, 2021 · aa28fb4 · aa28fb4
2 parents dc88a2b + 38241ca
commit aa28fb4
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/GSPH/RiemannSolvers/HLLC.cc b/src/GSPH/RiemannSolvers/HLLC.cc
@@ -155,7 +155,7 @@ interfaceState(const int i,
 
   }else{ // if ci & cj too small punt to normal av
     const auto uij = std::min((vi-vj).dot(rhatij),0.0);
-    Pstar = 0.5 * (uij*uij)/(rhoi+rhoj);
+    Pstar += 0.25 * (rhoi+rhoj) * (uij*uij);
   }
 }// Scalar interface class