Mercurial > repos > public > sbplib_julia

view SbpOperators/src/d2.jl @ 328:9cc5d1498b2d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Refactor 1D diagonal inner product in quadrature.jl to separate file. Write tests for quadratures. Clean up laplace and secondderivative
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 24 Sep 2020 22:31:48 +0200
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