Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/boundaryops/normal_derivative.jl @ 780:3b29b2ff1f0e operator_storage_array_of_table

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge in refactor/sbp_operators_method_signatures
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 19 Jul 2021 08:47:33 +0200
parents 5eb1edef8a7b
children 86776d06b883
comparison
equal deleted inserted replaced
772:bea2feebbeca 780:3b29b2ff1f0e
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