Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/volumeops/quadratures/inverse_quadrature.jl @ 674:621460cf8279

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge feature/stencil_convenience.
author Jonatan Werpers <jonatan@werpers.com>
date Sun, 07 Feb 2021 21:16:40 +0100
parents e14627e79a54
children 0bec3c4e78c0
comparison
equal deleted inserted replaced
664:7d7c1d636de3 674:621460cf8279
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