Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/d2.jl @ 637:4a81812150f4 feature/volume_and_boundary_operators

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Change qudrature closure from tuple of reals to tuple of Stencils. Also remove parametrization of stencil width in D2 since this was illformed for the 2nd order case.
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Sun, 03 Jan 2021 18:15:14 +0100
parents b05f542b2e8f
children f13d45c10f55
comparison
equal deleted inserted replaced
636:a1dfaf305f41 637:4a81812150f4
3 @enum Parity begin 3 @enum Parity begin
4 odd = -1 4 odd = -1
5 even = 1 5 even = 1
6 end 6 end
7 7
8 struct D2{T,N,M,K} <: ConstantStencilOperator 8 struct D2{T,M} <: ConstantStencilOperator
9 quadratureClosure::NTuple{M,T} 9 innerStencil::Stencil{T}
10 innerStencil::Stencil{T,N} 10 closureStencils::NTuple{M,Stencil{T}}
11 closureStencils::NTuple{M,Stencil{T,K}} 11 eClosure::Stencil{T}
12 eClosure::Stencil{T,M} 12 dClosure::Stencil{T}
13 dClosure::Stencil{T,M} 13 quadratureClosure::NTuple{M,Stencil{T}}
14 parity::Parity 14 parity::Parity
15 end 15 end
16 16
17 function closuresize(D::D2)::Int 17 closuresize(D::D2{T,M}) where {T,M} = M
18 return length(D.quadratureClosure)
19 end