view SbpOperators/src/d2.jl @ 273:babc4288e6a6 boundary_conditions

Dispatch apply(Quadrature) etc on region indices and create utility wrapper function dispatching on standard ints
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Mon, 06 Jan 2020 10:48:38 +0100
parents 51e7de109c25
children
line wrap: on
line source

export D2, closuresize, readOperator

@enum Parity begin
    odd = -1
    even = 1
end

struct D2{T,N,M,K} <: ConstantStencilOperator
    quadratureClosure::NTuple{M,T}
    innerStencil::Stencil{T,N}
    closureStencils::NTuple{M,Stencil{T,K}}
    eClosure::Stencil{T,M}
    dClosure::Stencil{T,M}
    parity::Parity
end

function closuresize(D::D2)::Int
    return length(D.quadratureClosure)
end