Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/boundaryops/normal_derivative.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 `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