comparison src/SbpOperators/volumeops/quadratures/inverse_quadrature.jl @ 765:fdd48f6ace1c operator_storage_array_of_table

Merge in default
author Jonatan Werpers <jonatan@werpers.com>
date Sun, 07 Feb 2021 21:32:21 +0100
parents e14627e79a54
children 0bec3c4e78c0
comparison
equal deleted inserted replaced
764:d34b515b0ae7 765:fdd48f6ace1c
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