Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/volumeops/quadratures/inverse_quadrature.jl @ 675:1ce3a104afc8 feature/boundary_quads

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge in default
author Jonatan Werpers <jonatan@werpers.com>
date Sun, 07 Feb 2021 21:28:53 +0100
parents e14627e79a54
children 0bec3c4e78c0
comparison
equal deleted inserted replaced
670:538ccbaeb1f8 675:1ce3a104afc8
30 Creates the inverse of the diagonal quadrature operator defined by the inner stencil 30 Creates the inverse of the diagonal quadrature operator defined by the inner stencil
31 1/h and a set of 1-element closure stencils in `closure_stencils`. Note that 31 1/h and a set of 1-element closure stencils in `closure_stencils`. Note that
32 the closure stencils are those of the quadrature operator (and not the inverse). 32 the closure stencils are those of the quadrature operator (and not the inverse).
33 """ 33 """
34 function InverseDiagonalQuadrature(grid::EquidistantGrid, closure_stencils::NTuple{M,Stencil{T,1}}) where {T,M} 34 function InverseDiagonalQuadrature(grid::EquidistantGrid, closure_stencils::NTuple{M,Stencil{T,1}}) where {T,M}
35 inv_inner_stencil = Stencil(Tuple{T}(1),center=1) 35 inv_inner_stencil = Stencil(one(T), center=1)
36 inv_closure_stencils = reciprocal_stencil.(closure_stencils) 36 inv_closure_stencils = reciprocal_stencil.(closure_stencils)
37 return InverseQuadrature(grid, inv_inner_stencil, inv_closure_stencils) 37 return InverseQuadrature(grid, inv_inner_stencil, inv_closure_stencils)
38 end 38 end
39 export InverseDiagonalQuadrature 39 export InverseDiagonalQuadrature
40 40