Mercurial > repos > public > sbplib_julia

changeset 380:81053b1992b6

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Remove parity field from Secondderivative
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 30 Sep 2020 21:35:10 +0200
parents de4746d6d126
children dacbcba33d7d
files src/SbpOperators/laplace/secondderivative.jl
diffstat 1 files changed, 2 insertions(+), 3 deletions(-) [+]
line wrap: on
line diff
--- a/src/SbpOperators/laplace/secondderivative.jl	Wed Sep 30 21:30:29 2020 +0200
+++ b/src/SbpOperators/laplace/secondderivative.jl	Wed Sep 30 21:35:10 2020 +0200
@@ -8,14 +8,13 @@
     h_inv::T # The grid spacing could be included in the stencil already. Preferable?
     innerStencil::Stencil{T,N}
     closureStencils::NTuple{M,Stencil{T,K}}
-    parity::Parity
     size::NTuple{1,Int}
 end
 export SecondDerivative
 
 function SecondDerivative(grid::EquidistantGrid{1}, innerStencil, closureStencils)
     h_inv = grid.inverse_spacing[1]
-    return SecondDerivative(h_inv, innerStencil, closureStencils, even, size(grid))
+    return SecondDerivative(h_inv, innerStencil, closureStencils, size(grid))
 end
 
 LazyTensors.range_size(D2::SecondDerivative) = D2.size
@@ -36,7 +35,7 @@
 
 function LazyTensors.apply(D2::SecondDerivative{T}, v::AbstractVector{T}, I::Index{Upper}) where T
     N = length(v) # TODO: Use domain_size here instead? N = domain_size(D2,size(v))
-    return @inbounds D2.h_inv*D2.h_inv*Int(D2.parity)*apply_stencil_backwards(D2.closureStencils[N-Int(I)+1], v, Int(I))
+    return @inbounds D2.h_inv*D2.h_inv*apply_stencil_backwards(D2.closureStencils[N-Int(I)+1], v, Int(I))
 end
 
 function LazyTensors.apply(D2::SecondDerivative{T}, v::AbstractVector{T}, index::Index{Unknown}) where T