Mercurial > repos > public > sbplib_julia

diff src/SbpOperators/boundaryops/normal_derivative.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/normal_derivative.jl	Mon Mar 21 10:04:15 2022 +0100
+++ b/src/SbpOperators/boundaryops/normal_derivative.jl	Mon Mar 21 12:51:39 2022 +0100
@@ -13,7 +13,9 @@
 function normal_derivative(grid, closure_stencil::Stencil, boundary)
     direction = dim(boundary)
     h_inv = inverse_spacing(grid)[direction]
-    return SbpOperators.boundary_operator(grid, scale(closure_stencil,h_inv), boundary)
+
+    op = BoundaryOperator(restrict(grid, dim(boundary)), scale(closure_stencil,h_inv), region(boundary))
+    return LazyTensors.inflate(op, size(grid), dim(boundary))
 end
 
 """