Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/boundaryops/normal_derivative.jl @ 975:5be8e25c81b3 feature/tensormapping_application_promotion

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge default
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 15 Mar 2022 07:37:11 +0100
parents 38d1752a9aff
children 7bf3121c6864 1ba8a398af9c
comparison
equal deleted inserted replaced
957:86889fc5b63f 975:5be8e25c81b3
1 """ 1 """
2 normal_derivative(grid::EquidistantGrid, closure_stencil::Stencil, boundary::CartesianBoundary) 2 normal_derivative(grid, closure_stencil::Stencil, boundary)
3 normal_derivative(grid::EquidistantGrid{1}, closure_stencil::Stencil, region::Region)
4 3
5 Creates the normal derivative boundary operator `d` as a `TensorMapping` 4 Creates the normal derivative boundary operator `d` as a `TensorMapping`
6 5
7 `d` is the normal derivative of a grid function at the boundary specified by `boundary` or `region` using some `closure_stencil`. 6 `d` computes the normal derivative of a grid function on `boundary` a `Stencil` `closure_stencil`.
8 `d'` is the prolongation of the normal derivative of a grid function to the whole grid using the same `closure_stencil`. 7 `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 8 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. 9 a `BoundaryOperator`.
10
11 See also: [`boundary_operator`](@ref).
11 """ 12 """
12 function normal_derivative(grid::EquidistantGrid, closure_stencil, boundary::CartesianBoundary) 13 function normal_derivative(grid, closure_stencil::Stencil, boundary)
13 direction = dim(boundary) 14 direction = dim(boundary)
14 h_inv = inverse_spacing(grid)[direction] 15 h_inv = inverse_spacing(grid)[direction]
15 return SbpOperators.boundary_operator(grid, scale(closure_stencil,h_inv), boundary) 16 return SbpOperators.boundary_operator(grid, scale(closure_stencil,h_inv), boundary)
16 end 17 end
17 normal_derivative(grid::EquidistantGrid{1}, closure_stencil, region::Region) = normal_derivative(grid, closure_stencil, CartesianBoundary{1,typeof(region)}()) 18
18 export normal_derivative 19 """
20 normal_derivative(grid, stencil_set, boundary)
21
22 Creates a `normal_derivative` operator on `grid` given a parsed TOML
23 `stencil_set`.
24 """
25 normal_derivative(grid, stencil_set, boundary) = normal_derivative(grid, parse_stencil(stencil_set["d1"]["closure"]), boundary)