Change UnitStep to StepFunction in older hydro examples

UnitStep is converted to Piecewise in Mathematica 8, so we use our own StepFunction instead.

4 files changed, 8 insertions, 8 deletions

diff --git a/Examples/Advect.m b/Examples/Advect.m
index 2bbb334..0463686 100644
--- a/Examples/Advect.m
+++ b/Examples/Advect.m
@@ -79,7 +79,7 @@ initialShockCalc =
     v1 -> v0,
     v2 -> 0,
     v3 -> 0,
-    rho -> amp UnitStep[x-0.5]
+    rho -> amp StepFunction[x-0.5]
   }
 };
 
diff --git a/Examples/Burgers.m b/Examples/Burgers.m
index 71d7930..3d89261 100644
--- a/Examples/Burgers.m
+++ b/Examples/Burgers.m
@@ -75,7 +75,7 @@ initialShockCalc =
   ConditionalOnKeyword -> {"initial_data", "shock"},
   Equations ->
   {
-    u -> uR0 UnitStep[x-0.5] + uL0 (1-UnitStep[x-0.5])
+    u -> uR0 StepFunction[x-0.5] + uL0 (1-StepFunction[x-0.5])
   }
 };
 
diff --git a/Examples/Euler.m b/Examples/Euler.m
index e3c22ad..b2716cd 100644
--- a/Examples/Euler.m
+++ b/Examples/Euler.m
@@ -88,11 +88,11 @@ initialShockCalc =
   ConditionalOnKeyword -> {"initial_data", "shock"},
   Equations ->
   {
-    rho -> rhoR0 UnitStep[x-0.5] + rhoL0 (1-UnitStep[x-0.5]),
-    v[1] -> vR0 UnitStep[x-0.5] + vL0 (1-UnitStep[x-0.5]),
+    rho -> rhoR0 StepFunction[x-0.5] + rhoL0 (1-StepFunction[x-0.5]),
+    v[1] -> vR0 StepFunction[x-0.5] + vL0 (1-StepFunction[x-0.5]),
     v[2] -> 0,
     v[3] -> 0,
-    p -> pR0 UnitStep[x-0.5] + pL0 (1-UnitStep[x-0.5])
+    p -> pR0 StepFunction[x-0.5] + pL0 (1-StepFunction[x-0.5])
   }
 };
 
diff --git a/Examples/EulerAuto.m b/Examples/EulerAuto.m
index e021d36..8eb69de 100644
--- a/Examples/EulerAuto.m
+++ b/Examples/EulerAuto.m
@@ -36,11 +36,11 @@ initialShockCalc =
   ConditionalOnKeyword -> {"initial_data", "shock"},
   Equations ->
   {
-    rho -> rhoR0 UnitStep[x-0.5] + rhoL0 (1-UnitStep[x-0.5]),
-    v[1] -> vR0 UnitStep[x-0.5] + vL0 (1-UnitStep[x-0.5]),
+    rho -> rhoR0 StepFunction[x-0.5] + rhoL0 (1-StepFunction[x-0.5]),
+    v[1] -> vR0 StepFunction[x-0.5] + vL0 (1-StepFunction[x-0.5]),
     v[2] -> 0,
     v[3] -> 0,
-    p -> pR0 UnitStep[x-0.5] + pL0 (1-UnitStep[x-0.5])
+    p -> pR0 StepFunction[x-0.5] + pL0 (1-StepFunction[x-0.5])
   }
 };