diff --git a/euler/src/eval.rs b/euler/src/eval.rs index 46bc5a3..74f6140 100644 --- a/euler/src/eval.rs +++ b/euler/src/eval.rs @@ -91,7 +91,7 @@ impl<'a, D: Dimension, BP: EvaluatorPressure> Evaluator *e = p / (GAMMA - 1.0) + rho * (u*u + v*v) / 2.0; }); - azip!((rho in &rho, rhou in &mut rhou) *rhou /= rho); - azip!((rho in &rho, rhov in &mut rhov) *rhov /= rho); + azip!((rho in &rho, rhou in &mut rhou) *rhou *= rho); + azip!((rho in &rho, rhov in &mut rhov) *rhov *= rho); } }