Mercurial > repos > public > sbplib_julia

diff src/SbpOperators/boundaryops/boundary_restriction.jl @ 1099:05a25a5063bb refactor/sbpoperators/inflation

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Try to remove volume_operator and boundary_operator methods
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 21 Mar 2022 12:51:39 +0100
parents 1ba8a398af9c
children 157a78959e5d
line wrap: on
line diff
--- a/src/SbpOperators/boundaryops/boundary_restriction.jl	Mon Mar 21 10:04:15 2022 +0100
+++ b/src/SbpOperators/boundaryops/boundary_restriction.jl	Mon Mar 21 12:51:39 2022 +0100
@@ -15,7 +15,9 @@
 """
 function boundary_restriction(grid, closure_stencil::Stencil, boundary)
     converted_stencil = convert(Stencil{eltype(grid)}, closure_stencil)
-    return SbpOperators.boundary_operator(grid, converted_stencil, boundary)
+
+    op = BoundaryOperator(restrict(grid, dim(boundary)), converted_stencil, region(boundary))
+    return LazyTensors.inflate(op, size(grid), dim(boundary))
 end
 
 """