From 227ea812f7f17c04762944e71860c3ff696250d6 Mon Sep 17 00:00:00 2001 From: Ian Hinder Date: Thu, 15 Dec 2011 19:03:44 +0100 Subject: Regenerate hydro examples --- Examples/Euler/src/euler_initial_shock.cc | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) (limited to 'Examples/Euler/src') diff --git a/Examples/Euler/src/euler_initial_shock.cc b/Examples/Euler/src/euler_initial_shock.cc index b8fe537..d20ce9a 100644 --- a/Examples/Euler/src/euler_initial_shock.cc +++ b/Examples/Euler/src/euler_initial_shock.cc @@ -114,17 +114,17 @@ static void euler_initial_shock_Body(cGH const * restrict const cctkGH, int cons /* Precompute derivatives */ /* Calculate temporaries and grid functions */ - CCTK_REAL rhoL = Piecewise(List(List(ToReal(rhoL0),xL < 0.5)),0. + - ToReal(rhoR0)); + CCTK_REAL rhoL = ToReal(rhoL0) + StepFunction(-0.5 + + xL)*(-ToReal(rhoL0) + ToReal(rhoR0)); - CCTK_REAL v1L = Piecewise(List(List(ToReal(vL0),xL < 0.5)),0. + + CCTK_REAL v1L = ToReal(vL0) + StepFunction(-0.5 + xL)*(-ToReal(vL0) + ToReal(vR0)); CCTK_REAL v2L = 0; CCTK_REAL v3L = 0; - CCTK_REAL pL = Piecewise(List(List(ToReal(pL0),xL < 0.5)),0. + + CCTK_REAL pL = ToReal(pL0) + StepFunction(-0.5 + xL)*(-ToReal(pL0) + ToReal(pR0)); /* Copy local copies back to grid functions */ -- cgit v1.2.3