Mercurial > repos > public > sbplib_julia

diff src/SbpOperators/stencil.jl @ 1453:e1222fbb7c4d bugfix/sbp_operators/stencil_return_type

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix broken test using Core.Compiler.return_type. Needs to be evaluated for performance
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 24 Nov 2023 22:42:50 +0100
parents 14cb97284373
children
line wrap: on
line diff
--- a/src/SbpOperators/stencil.jl	Fri Nov 24 21:49:07 2023 +0100
+++ b/src/SbpOperators/stencil.jl	Fri Nov 24 22:42:50 2023 +0100
@@ -69,7 +69,8 @@
 end
 
 Base.@propagate_inbounds @inline function apply_stencil(s::Stencil, v::AbstractVector, i::Int)
-    w = zero(promote_type(eltype(s),eltype(v)))
+    T = Core.Compiler.return_type(*, Tuple{eltype(s),eltype(v)})
+    w = zero(T)
     @simd for k ∈ 1:length(s)
         w += s.weights[k]*v[i + s.range[k]]
     end
@@ -78,7 +79,8 @@
 end
 
 Base.@propagate_inbounds @inline function apply_stencil_backwards(s::Stencil, v::AbstractVector, i::Int)
-    w = zero(promote_type(eltype(s),eltype(v)))
+    T = Core.Compiler.return_type(*, Tuple{eltype(s),eltype(v)})
+    w = zero(T)
     @simd for k ∈ length(s):-1:1
         w += s.weights[k]*v[i - s.range[k]]
     end