Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/boundaryops/boundary_restriction.jl @ 875:067a322e4f73 laplace_benchmarks

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge with feature/laplace_opset
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 27 Jan 2022 10:55:08 +0100
parents 86776d06b883
children 6b47a9ee1632
comparison
equal deleted inserted replaced
874:7e9ebd572deb 875:067a322e4f73
7 `e` is the restriction of a grid function to the boundary specified by `boundary` or `region` using some `closure_stencil`. 7 `e` is the restriction of a grid function to the boundary specified by `boundary` or `region` using some `closure_stencil`.
8 `e'` is the prolongation of a grid function on the boundary to the whole grid using the same `closure_stencil`. 8 `e'` is the prolongation of a grid function on the boundary to the whole grid using the same `closure_stencil`.
9 On a one-dimensional `grid`, `e` is a `BoundaryOperator`. On a multi-dimensional `grid`, `e` is the inflation of 9 On a one-dimensional `grid`, `e` is a `BoundaryOperator`. On a multi-dimensional `grid`, `e` is the inflation of
10 a `BoundaryOperator`. Also see the documentation of `SbpOperators.boundary_operator(...)` for more details. 10 a `BoundaryOperator`. Also see the documentation of `SbpOperators.boundary_operator(...)` for more details.
11 """ 11 """
12 boundary_restriction(grid::EquidistantGrid, closure_stencil::Stencil, boundary::CartesianBoundary) = SbpOperators.boundary_operator(grid, closure_stencil, boundary) 12 function boundary_restriction(grid::EquidistantGrid, closure_stencil, boundary::CartesianBoundary)
13 boundary_restriction(grid::EquidistantGrid{1}, closure_stencil::Stencil, region::Region) = boundary_restriction(grid, closure_stencil, CartesianBoundary{1,typeof(region)}()) 13 converted_stencil = convert(Stencil{eltype(grid)}, closure_stencil)
14 14 return SbpOperators.boundary_operator(grid, converted_stencil, boundary)
15 export boundary_restriction 15 end
16 boundary_restriction(grid::EquidistantGrid{1}, closure_stencil, region::Region) = boundary_restriction(grid, closure_stencil, CartesianBoundary{1,typeof(region)}())