Skip to content

Commit 967ce42

Browse files
Fix missing SSPRK interpolation function
Not sure how this got dropped in the move.
1 parent b03a2a3 commit 967ce42

File tree

3 files changed

+12
-4
lines changed

3 files changed

+12
-4
lines changed

lib/OrdinaryDiffEqSSPRK/src/OrdinaryDiffEqSSPRK.jl

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -29,7 +29,7 @@ include("addsteps.jl")
2929
include("functions.jl")
3030

3131
export SSPRK53_2N2, SSPRK22, SSPRK53, SSPRK63, SSPRK83, SSPRK43, SSPRK432, SSPRKMSVS32,
32-
SSPRK54, SSPRK53_2N1, SSPRK104, SSPRK932, SSPRKMSVS43, SSPRK73, SSPRK53_H,
32+
SSPRK54, SSPRK53_2N1, SSPRK104, SSPRK932, SSPRKMSVS43, SSPRK73, SSPRK53_H,
3333
SSPRK33, SHLDDRK_2N, KYKSSPRK42, SHLDDRK52
3434

3535
end

lib/OrdinaryDiffEqSSPRK/src/algorithms.jl

Lines changed: 0 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -1,6 +1,3 @@
1-
using OrdinaryDiffEq: explicit_rk_docstring, trivial_limiter!
2-
using Static: False
3-
41
@doc explicit_rk_docstring(
52
"A third-order, five-stage explicit strong stability preserving (SSP) low-storage method.
63
Fixed timestep only.",

lib/OrdinaryDiffEqSSPRK/src/interpolants.jl

Lines changed: 11 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -123,6 +123,17 @@ end
123123
c10diff2invdt2+k[1][idxs] * b10diff2invdt
124124
end
125125

126+
@muladd function _ode_interpolant!(out, Θ, dt, y₀, y₁, k,
127+
cache::Union{SSPRK22ConstantCache, SSPRK22Cache,
128+
SSPRK33ConstantCache, SSPRK33Cache,
129+
SSPRK43ConstantCache, SSPRK43Cache,
130+
SSPRK432ConstantCache, SSPRK432Cache},
131+
idxs, T::Type{Val{2}}, differential_vars::Nothing)
132+
@ssprkpre2
133+
@inbounds @.. broadcast=false out=(y₁[idxs] - y₀[idxs]) * c10diff2invdt2 + k[1][idxs] * b10diff2invdt
134+
out
135+
end
136+
126137
@muladd function _ode_interpolant!(out, Θ, dt, y₀, y₁, k,
127138
cache::Union{SSPRK22ConstantCache, SSPRK22Cache,
128139
SSPRK33ConstantCache, SSPRK33Cache,

0 commit comments

Comments
 (0)