Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/volumeops/quadratures/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 2a95beb9ef1d e14627e79a54
children 43cf58c69f91
comparison
equal deleted inserted replaced
670:538ccbaeb1f8 675:1ce3a104afc8
24 return H 24 return H
25 end 25 end
26 export quadrature 26 export quadrature
27 27
28 function quadrature(grid::EquidistantGrid, closure_stencils::NTuple{M,Stencil{T}}) where {M,T} 28 function quadrature(grid::EquidistantGrid, closure_stencils::NTuple{M,Stencil{T}}) where {M,T}
29 inner_stencil = Stencil(Tuple{T}(1),center=1) 29 inner_stencil = CenteredStencil(one(T))
30 return quadrature(grid, inner_stencil, closure_stencils) 30 return quadrature(grid, inner_stencil, closure_stencils)
31 end 31 end
32 32
33 """ 33 """
34 boundary_quadrature(grid::EquidistantGrid, inner_stencil, closure_stencils, id::CartesianBoundary) 34 boundary_quadrature(grid::EquidistantGrid, inner_stencil, closure_stencils, id::CartesianBoundary)
50 function boundary_quadrature(grid::EquidistantGrid{1}, inner_stencil::Stencil{T}, closure_stencils::NTuple{M,Stencil{T}}, id::CartesianBoundary{1}) where {M,T} 50 function boundary_quadrature(grid::EquidistantGrid{1}, inner_stencil::Stencil{T}, closure_stencils::NTuple{M,Stencil{T}}, id::CartesianBoundary{1}) where {M,T}
51 return IdentityMapping{T}() 51 return IdentityMapping{T}()
52 end 52 end
53 53
54 function boundary_quadrature(grid::EquidistantGrid, closure_stencils::NTuple{M,Stencil{T}}, id::CartesianBoundary) where {M,T} 54 function boundary_quadrature(grid::EquidistantGrid, closure_stencils::NTuple{M,Stencil{T}}, id::CartesianBoundary) where {M,T}
55 inner_stencil = Stencil(Tuple{T}(1),center=1) 55 inner_stencil = CenteredStencil(one(T))
56 return boundary_quadrature(grid,inner_stencil,closure_stencils,id) 56 return boundary_quadrature(grid,inner_stencil,closure_stencils,id)
57 end 57 end
58 58
59 """ 59 """
60 orthogonal_grid(grid,dim) 60 orthogonal_grid(grid,dim)