diff options
author | Ian Hinder <ian.hinder@aei.mpg.de> | 2011-12-15 19:03:44 +0100 |
---|---|---|
committer | Ian Hinder <ian.hinder@aei.mpg.de> | 2011-12-15 19:03:44 +0100 |
commit | 227ea812f7f17c04762944e71860c3ff696250d6 (patch) | |
tree | a9378b89be164c146bbace1e320b8e6820071606 /Examples/EulerAuto | |
parent | 7bbffe843d3a60b2f3432e492700cf4d1f3ff7eb (diff) |
Regenerate hydro examples
Diffstat (limited to 'Examples/EulerAuto')
-rw-r--r-- | Examples/EulerAuto/src/eulerauto_initial_shock.cc | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/Examples/EulerAuto/src/eulerauto_initial_shock.cc b/Examples/EulerAuto/src/eulerauto_initial_shock.cc index d4bf071..f6893b7 100644 --- a/Examples/EulerAuto/src/eulerauto_initial_shock.cc +++ b/Examples/EulerAuto/src/eulerauto_initial_shock.cc @@ -93,17 +93,17 @@ static void eulerauto_initial_shock_Body(cGH const * restrict const cctkGH, int /* 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 */ |