Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/boundaryops/boundary_restriction.jl @ 1652:65b2d2c72fbc feature/sbp_operators/laplace_curvilinear

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add boundary restriction operator for mapped grid
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 26 Jun 2024 12:54:29 +0200
parents 08f06bfacd5c
children 1f42944d4a72
comparison
equal deleted inserted replaced
1651:707fc9761c2b 1652:65b2d2c72fbc
23 function boundary_restriction(g::EquidistantGrid, stencil_set::StencilSet, boundary) 23 function boundary_restriction(g::EquidistantGrid, stencil_set::StencilSet, boundary)
24 closure_stencil = parse_stencil(stencil_set["e"]["closure"]) 24 closure_stencil = parse_stencil(stencil_set["e"]["closure"])
25 converted_stencil = convert(Stencil{eltype(g)}, closure_stencil) 25 converted_stencil = convert(Stencil{eltype(g)}, closure_stencil)
26 return BoundaryOperator(g, converted_stencil, boundary) 26 return BoundaryOperator(g, converted_stencil, boundary)
27 end 27 end
28
29 function boundary_restriction(g::MappedGrid, stencil_set::StencilSet, boundary)
30 return boundary_restriction(logicalgrid(g), stencil_set, boundary)
31 end