view src/SbpOperators/d2.jl @ 638:08b2c7a2d063 feature/volume_and_boundary_operators

Implement the Quadrature operator as a VolumeOperator. Make DiagonalQuadrature a special case of the general Quadrature operator. Update tests.
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Mon, 04 Jan 2021 09:32:11 +0100
parents 4a81812150f4
children f13d45c10f55
line wrap: on
line source

export D2, closuresize

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

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

closuresize(D::D2{T,M}) where {T,M} = M