Fix bug in HLLC 3D.

pkestene · Feb 18, 2024 · a7e369d · a7e369d
1 parent 866aab9
commit a7e369d
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/shared/RiemannSolvers.h b/src/shared/RiemannSolvers.h
@@ -474,7 +474,7 @@ void riemann_hllc(const HydroState& qleft,
   real_t ecinr = HALF_F*rr*ur*ur;
   ecinr += HALF_F*rr*qright[IV]*qright[IV];
   if (std::is_same<HydroState,HydroState3d>::value)
-    ecinl += HALF_F*rr*qright[IW]*qright[IW];
+    ecinr += HALF_F*rr*qright[IW]*qright[IW];
 
   real_t etotr = pr*entho+ecinr;
   real_t ptotr = pr;