Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/constantstenciloperator.jl @ 543:1a53eb83ed24 refactor/tensor_index_coupling

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Remove some unecessary Index types
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 27 Nov 2020 12:02:47 +0100
parents 011ca1639153
children d1929491180b
comparison
equal deleted inserted replaced
542:011ca1639153 543:1a53eb83ed24
60 -h_inv*apply_stencil_backwards(op.dClosure,v,length(v)) 60 -h_inv*apply_stencil_backwards(op.dClosure,v,length(v))
61 end 61 end
62 62
63 export apply_normal_derivative_transpose 63 export apply_normal_derivative_transpose
64 64
65 function apply_normal_derivative(op::ConstantStencilOperator, h_inv::Real, v::Number, i::Index, N::Integer, ::Type{Lower}) 65 function apply_normal_derivative(op::ConstantStencilOperator, h_inv::Real, v::Number, i, N::Integer, ::Type{Lower})
66 @boundscheck if !(0<length(Int(i)) <= N) 66 @boundscheck if !(0<length(Int(i)) <= N)
67 throw(BoundsError()) 67 throw(BoundsError())
68 end 68 end
69 h_inv*op.dClosure[Int(i)-1]*v 69 h_inv*op.dClosure[Int(i)-1]*v
70 end 70 end
71 71
72 function apply_normal_derivative(op::ConstantStencilOperator, h_inv::Real, v::Number, i::Index, N::Integer, ::Type{Upper}) 72 function apply_normal_derivative(op::ConstantStencilOperator, h_inv::Real, v::Number, i, N::Integer, ::Type{Upper})
73 @boundscheck if !(0<length(Int(i)) <= N) 73 @boundscheck if !(0<length(Int(i)) <= N)
74 throw(BoundsError()) 74 throw(BoundsError())
75 end 75 end
76 -h_inv*op.dClosure[N-Int(i)]*v 76 -h_inv*op.dClosure[N-Int(i)]*v
77 end 77 end