diff options
-rw-r--r-- | Examples/Advect/src/advect_initial_shock.cc | 2 | ||||
-rw-r--r-- | Examples/Burgers/src/burgers_initial_shock.cc | 2 | ||||
-rw-r--r-- | Examples/Euler/src/euler_initial_shock.cc | 8 | ||||
-rw-r--r-- | Examples/EulerAuto/src/eulerauto_initial_shock.cc | 8 |
4 files changed, 10 insertions, 10 deletions
diff --git a/Examples/Advect/src/advect_initial_shock.cc b/Examples/Advect/src/advect_initial_shock.cc index c599296..5c9706e 100644 --- a/Examples/Advect/src/advect_initial_shock.cc +++ b/Examples/Advect/src/advect_initial_shock.cc @@ -120,7 +120,7 @@ static void advect_initial_shock_Body(cGH const * restrict const cctkGH, int con CCTK_REAL v3L = 0; - CCTK_REAL rhoL = ToReal(amp)*UnitStep(-0.5 + xL); + CCTK_REAL rhoL = StepFunction(-0.5 + xL)*ToReal(amp); /* Copy local copies back to grid functions */ rho[index] = rhoL; diff --git a/Examples/Burgers/src/burgers_initial_shock.cc b/Examples/Burgers/src/burgers_initial_shock.cc index a0c5b6e..3baf23f 100644 --- a/Examples/Burgers/src/burgers_initial_shock.cc +++ b/Examples/Burgers/src/burgers_initial_shock.cc @@ -93,7 +93,7 @@ static void burgers_initial_shock_Body(cGH const * restrict const cctkGH, int co /* Precompute derivatives */ /* Calculate temporaries and grid functions */ - CCTK_REAL uL = Piecewise(List(List(ToReal(uL0),xL < 0.5)),0. + + CCTK_REAL uL = ToReal(uL0) + StepFunction(-0.5 + xL)*(-ToReal(uL0) + ToReal(uR0)); /* Copy local copies back to grid functions */ 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 */ 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 */ |