allow other real types

JoshuaLampert · Feb 1, 2024 · 529413a · 529413a
1 parent 46e6514
commit 529413a
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/src/equations/bbm_bbm_variable_bathymetry_1d.jl b/src/equations/bbm_bbm_variable_bathymetry_1d.jl
@@ -67,7 +67,7 @@ function initial_condition_manufactured(x, t,
                                         mesh)
     eta = exp(t) * cospi(2 * (x - 2 * t))
     v = exp(t / 2) * sinpi(2 * (x - t / 2))
-    D = 5.0 + 2.0 * cospi(2 * x)
+    D = 5 + 2 * cospi(2 * x)
     return SVector(eta, v, D)
 end
 
@@ -97,7 +97,7 @@ function source_terms_manufactured(q, x, t, equations::BBMBBMVariableEquations1D
            4 * (a4 + 2 * pi * a3) * (16 * sinpi(x)^4 - 26 * sinpi(x)^2 + 7)) * exp(t / 2) /
           3 - exp(t / 2) * a4 / 2 - pi * exp(t / 2) * a3 - pi * exp(t) * a5
 
-    return SVector(dq1, dq2, 0.0)
+    return SVector(dq1, dq2, zero(dq1))
 end
 
 """

diff --git a/src/equations/svaerd_kalisch_1d.jl b/src/equations/svaerd_kalisch_1d.jl
@@ -157,7 +157,7 @@ function source_terms_manufactured(q, x, t, equations::SvaerdKalischEquations1D)
           (a19) * a8 * a3 / 2 + 3 * pi * a18 * exp(t) * a3 * a4 -
           2.5 * a14 * a16 * a15 * a8 * a3 / (a13)
 
-    return SVector(dq1, dq2, 0.0)
+    return SVector(dq1, dq2, zero(dq1))
 end
 
 function create_cache(mesh,