mulimoen
/
SummationByParts
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Fix bug with u,v->rhou,rhov conversion

This commit is contained in:
Magnus Ulimoen 2021-06-29 21:31:40 +02:00
parent 0878e75b8d
commit f6f70fe79a
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
euler/src/eval.rs
View File

@ -91,7 +91,7 @@ impl<'a, D: Dimension, BP: EvaluatorPressure<D>> Evaluator<D>
*e = p / (GAMMA - 1.0) + rho * (u*u + v*v) / 2.0; *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, rhou in &mut rhou) *rhou *= rho);
azip!((rho in &rho, rhov in &mut rhov) *rhov /= rho); azip!((rho in &rho, rhov in &mut rhov) *rhov *= rho);
} }
} }