Mercurial > repos > public > sbplib_julia
comparison src/SbpOperators/boundaryops/normal_derivative.jl @ 651:67639b1c99ea
Merged feature/volume_and_boundary_operators
author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
---|---|
date | Wed, 20 Jan 2021 17:52:55 +0100 |
parents | bfc893d03cbf |
children | 1accc3e051d0 |
comparison
equal
deleted
inserted
replaced
615:52749b687a67 | 651:67639b1c99ea |
---|---|
1 """ | |
2 NormalDerivative(grid::EquidistantGrid, closure_stencil::Stencil, boundary::CartesianBoundary) | |
3 NormalDerivative(grid::EquidistantGrid{1}, closure_stencil::Stencil, region::Region) | |
4 | |
5 Creates the normal derivative boundary operator `d` as a `TensorMapping` | |
6 | |
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`. | |
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. | |
11 """ | |
12 function NormalDerivative(grid::EquidistantGrid, closure_stencil::Stencil, boundary::CartesianBoundary) | |
13 direction = dim(boundary) | |
14 h_inv = inverse_spacing(grid)[direction] | |
15 return SbpOperators.boundary_operator(grid, scale(closure_stencil,h_inv), boundary) | |
16 end | |
17 NormalDerivative(grid::EquidistantGrid{1}, closure_stencil::Stencil, region::Region) = NormalDerivative(grid, closure_stencil, CartesianBoundary{1,typeof(region)}()) | |
18 export NormalDerivative |