Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/boundaryops/normal_derivative.jl @ 866:1784b1c0af3e feature/laplace_opset

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge with default
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Wed, 19 Jan 2022 14:44:24 +0100
parents 5eb1edef8a7b
children 86776d06b883
comparison
equal deleted inserted replaced
865:545a6c1a0a0e 866:1784b1c0af3e
7 `d` is the normal derivative of a grid function at the boundary specified by `boundary` or `region` using some `closure_stencil`. 7 `d` is the normal derivative of a grid function at the boundary specified by `boundary` or `region` using some `closure_stencil`.
8 `d'` is the prolongation of the normal derivative of a grid function to the whole grid using the same `closure_stencil`. 8 `d'` is the prolongation of the normal derivative of a grid function to the whole grid using the same `closure_stencil`.
9 On a one-dimensional `grid`, `d` is a `BoundaryOperator`. On a multi-dimensional `grid`, `d` is the inflation of 9 On a one-dimensional `grid`, `d` is a `BoundaryOperator`. On a multi-dimensional `grid`, `d` 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 function normal_derivative(grid::EquidistantGrid, closure_stencil::Stencil, boundary::CartesianBoundary) 12 function normal_derivative(grid::EquidistantGrid, closure_stencil, boundary::CartesianBoundary)
13 direction = dim(boundary) 13 direction = dim(boundary)
14 h_inv = inverse_spacing(grid)[direction] 14 h_inv = inverse_spacing(grid)[direction]
15 return SbpOperators.boundary_operator(grid, scale(closure_stencil,h_inv), boundary) 15 return SbpOperators.boundary_operator(grid, scale(closure_stencil,h_inv), boundary)
16 end 16 end
17 normal_derivative(grid::EquidistantGrid{1}, closure_stencil::Stencil, region::Region) = normal_derivative(grid, closure_stencil, CartesianBoundary{1,typeof(region)}()) 17 normal_derivative(grid::EquidistantGrid{1}, closure_stencil, region::Region) = normal_derivative(grid, closure_stencil, CartesianBoundary{1,typeof(region)}())
18 export normal_derivative 18 export normal_derivative