Mercurial > repos > public > sbplib_julia
changeset 1458:b411493597f3 bugfix/sbp_operators/stencil_return_type
Update op for backwards apply too
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Sun, 26 Nov 2023 23:07:21 +0100 |
| parents | 6e9ac7925f20 |
| children | 450110ed5d18 |
| files | src/SbpOperators/stencil.jl |
| diffstat | 1 files changed, 1 insertions(+), 1 deletions(-) [+] |
line wrap: on
line diff
--- a/src/SbpOperators/stencil.jl Sun Nov 26 21:14:44 2023 +0100 +++ b/src/SbpOperators/stencil.jl Sun Nov 26 23:07:21 2023 +0100 @@ -79,7 +79,7 @@ end Base.@propagate_inbounds @inline function apply_stencil_backwards(s::Stencil, v::AbstractVector, i::Int) - T = Base.promote_op(*, eltype(s), eltype(v)) + T = Base.promote_op((a₁,v₁,a₂,v₂)->a₁*v₁ + a₂*v₂, eltype(s), eltype(v), eltype(s), eltype(v)) w = zero(T) @simd for k ∈ length(s):-1:1 w += s.weights[k]*v[i - s.range[k]]