Mercurial > repos > public > sbplib_julia
diff src/SbpOperators/laplace/secondderivative.jl @ 543:1a53eb83ed24 refactor/tensor_index_coupling
Remove some unecessary Index types
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Fri, 27 Nov 2020 12:02:47 +0100 |
| parents | 011ca1639153 |
| children | 9330338d6ab5 |
line wrap: on
line diff
--- a/src/SbpOperators/laplace/secondderivative.jl Fri Nov 27 11:27:37 2020 +0100 +++ b/src/SbpOperators/laplace/secondderivative.jl Fri Nov 27 12:02:47 2020 +0100 @@ -20,10 +20,6 @@ LazyTensors.range_size(D2::SecondDerivative) = D2.size LazyTensors.domain_size(D2::SecondDerivative) = D2.size -#TODO: The 1D tensor mappings should not have to dispatch on 1D tuples if we write LazyTensor.apply for vararg right?!?! -# Currently have to index the Tuple{Index} in each method in order to call the stencil methods which is ugly. -# I thought I::Vararg{Index,R} fell back to just Index for R = 1 - # Apply for different regions Lower/Interior/Upper or Unknown region function LazyTensors.apply(D2::SecondDerivative{T}, v::AbstractVector{T}, I::Index{Lower}) where T return @inbounds D2.h_inv*D2.h_inv*apply_stencil(D2.closureStencils[Int(I)], v, Int(I))