Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/boundaryops/normal_derivative.jl @ 769:0158c3fd521c operator_storage_array_of_table

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge in default
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 15 Jul 2021 00:06:16 +0200
parents 1accc3e051d0
children 5eb1edef8a7b
comparison
equal deleted inserted replaced
768:7c87a33963c5 769:0158c3fd521c
1 """ 1 """
2 NormalDerivative(grid::EquidistantGrid, closure_stencil::Stencil, boundary::CartesianBoundary) 2 normal_derivative(grid::EquidistantGrid, closure_stencil::Stencil, boundary::CartesianBoundary)
3 NormalDerivative(grid::EquidistantGrid{1}, closure_stencil::Stencil, region::Region) 3 normal_derivative(grid::EquidistantGrid{1}, closure_stencil::Stencil, region::Region)
4 4
5 Creates the normal derivative boundary operator `d` as a `TensorMapping` 5 Creates the normal derivative boundary operator `d` as a `TensorMapping`
6 6
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 NormalDerivative(grid::EquidistantGrid, closure_stencil::Stencil, boundary::CartesianBoundary) 12 function normal_derivative(grid::EquidistantGrid, closure_stencil::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 NormalDerivative(grid::EquidistantGrid{1}, closure_stencil::Stencil, region::Region) = NormalDerivative(grid, closure_stencil, CartesianBoundary{1,typeof(region)}()) 17 normal_derivative(grid::EquidistantGrid{1}, closure_stencil::Stencil, region::Region) = normal_derivative(grid, closure_stencil, CartesianBoundary{1,typeof(region)}())
18 export NormalDerivative 18 export normal_derivative