Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/volumeops/volume_operator.jl @ 1356:49d03d1169ef feature/variable_derivatives

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Remove code designed to help inference because it seems to have been solved elsewhere
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 08 Feb 2023 21:32:48 +0100
parents 102ebdaf7c11
children 4684c7f1c4cb
comparison
equal deleted inserted replaced
1355:102ebdaf7c11 1356:49d03d1169ef
30 function LazyTensors.apply(op::VolumeOperator, v::AbstractVector, i::Index{Upper}) 30 function LazyTensors.apply(op::VolumeOperator, v::AbstractVector, i::Index{Upper})
31 return @inbounds Int(op.parity)*apply_stencil_backwards(op.closure_stencils[op.size[1]-Int(i)+1], v, Int(i)) 31 return @inbounds Int(op.parity)*apply_stencil_backwards(op.closure_stencils[op.size[1]-Int(i)+1], v, Int(i))
32 end 32 end
33 33
34 function LazyTensors.apply(op::VolumeOperator, v::AbstractVector, i) 34 function LazyTensors.apply(op::VolumeOperator, v::AbstractVector, i)
35 return LazyTensors.apply_with_region(op, v, closure_size(op), op.size[1], i) 35 r = getregion(i, closure_size(op), op.size[1])
36 return LazyTensors.apply(op, v, Index(i, r))
36 end 37 end
37 # TODO: Move this to LazyTensors when we have the region communication down. 38 # TODO: Move this to LazyTensors when we have the region communication down.