From 92ce7d11645bead38983194dc7c26831e0172966 Mon Sep 17 00:00:00 2001 From: Ian Hinder Date: Thu, 30 Sep 2010 11:40:03 +0100 Subject: Euler.m: Fix to use reconstructed En --- Examples/Euler.m | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) (limited to 'Examples/Euler.m') diff --git a/Examples/Euler.m b/Examples/Euler.m index 27f872b..75f678f 100644 --- a/Examples/Euler.m +++ b/Examples/Euler.m @@ -160,7 +160,7 @@ eulerDenFlux[rho_, v_, p_, i_] := eulerSFlux[rho_, v_, p_, {i_, j_}] := rho v[i] v[j] + p Euc[i,j]; -eulerEnFlux[rho_, v_, p_, i_] := +eulerEnFlux[rho_, v_, p_, En_, i_] := v[i](En + p); zeroRHSCalc[] := @@ -225,10 +225,10 @@ fluxCalc[i_] := ShiftMinus[pRight, i], {i, uj}] + alpha (ShiftMinus[SRight[uj],i] - SLeft[uj])), - EnF -> 1/2 (eulerEnFlux[rhoLeft, vLeft, pLeft, i] + + EnF -> 1/2 (eulerEnFlux[rhoLeft, vLeft, pLeft, EnLeft, i] + eulerEnFlux[ShiftMinus[rhoRight,i], vRightTemp, - ShiftMinus[pRight,i], i] + + ShiftMinus[pRight,i], ShiftMinus[EnRight,i], i] + alpha (ShiftMinus[EnRight,i] - EnLeft)) } }; -- cgit v1.2.3