Mercurial > repos > public > sbplib_julia
view stencil.jl @ 133:c4cc2a1a68b9 stencil_index
Close branch
author | Jonatan Werpers <jonatan@werpers.com> |
---|---|
date | Thu, 21 Feb 2019 16:20:55 +0100 |
parents | 84b1ad5a3755 |
children |
line wrap: on
line source
struct Stencil{T<:Real,N} range::Tuple{Int,Int} weights::NTuple{N,T} end function flip(s::Stencil) range = (-s.range[2], -s.range[1]) return Stencil(range, reverse(s.weights)) end # Provides index into the Stencil based on offset for the root element function Base.getindex(s::Stencil, i::Int) @boundscheck if i < s.range[1] || s.range[2] < i return eltype(s.weights)(0) end return s.weights[1 + i - s.range[1]] end Base.@propagate_inbounds function apply(s::Stencil, v::AbstractVector, i::Int) w = zero(eltype(v)) for j ∈ s.range[1]:s.range[2] @inbounds weight = s[j] w += weight*v[i+j] end return w end