We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent b9a8bd9 commit b8632baCopy full SHA for b8632ba
2 files changed
src/extrapolation/constant.jl
@@ -7,15 +7,9 @@ ConstantExtrapolation{T,ITP,IT,GT}(::Type{T}, N, itp::ITP, ::Type{IT}, ::Type{GT
7
extrapolate{T,N,IT,GT}(itp::AbstractInterpolation{T,N,IT,GT}, ::Type{Flat}) =
8
ConstantExtrapolation(T,N,itp,IT,GT)
9
10
-function extrap_prep{T,ITP,IT}(etp::Type{ConstantExtrapolation{T,1,ITP,IT,OnGrid}}, x)
11
- :(x = clamp(x, 1, size(etp,1)))
12
-end
13
-function extrap_prep{T,ITP,IT}(etp::Type{ConstantExtrapolation{T,1,ITP,IT,OnCell}}, x)
14
- :(x = clamp(x, .5, size(etp,1)+.5))
15
16
function extrap_prep{T,N,ITP,IT}(etp::Type{ConstantExtrapolation{T,N,ITP,IT,OnGrid}}, xs...)
17
- :(@nexprs $N d->(xs[d] = clamp(xs[d], 1, size(etp,d))))
+ :(@nexprs $N d->(x_d = clamp(xs[d], one($xs[d]), size(etp,d))))
18
end
19
function extrap_prep{T,N,ITP,IT}(etp::Type{ConstantExtrapolation{T,N,ITP,IT,OnCell}}, xs...)
20
- :(@nexprs $N d->(xs[d] = clamp(xs[d], .5, size(etp,d)+.5)))
+ :(@nexprs $N d->(x_d = clamp(xs[d], 1//2, size(etp,d)+1//2)))
21
src/extrapolation/periodic.jl
@@ -8,5 +8,5 @@ extrapolate{T,N,IT,GT}(itp::AbstractInterpolation{T,N,IT,GT}, ::Type{Periodic})
PeriodicExtrapolation(T,N,itp,IT,GT)
function extrap_prep{T,N,ITP,IT,GT}(etp::Type{PeriodicExtrapolation{T,N,ITP,IT,GT}}, xs...)
- :(@nexprs $N d->(; x_d = mod(xs[d]-one($xs[d]), size(etp.itp,d)-one($xs[d]))+one($xs[d])))
+ :(@nexprs $N d->(x_d = mod(xs[d]-one($xs[d]), size(etp.itp,d)-one($xs[d]))+one($xs[d])))
0 commit comments