Mercurial > repos > public > sbplib_julia
changeset 1484:8d60d045c2a2 feature/boundary_conditions
Add todo and minor edits to variable naming and docstring
| author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
|---|---|
| date | Mon, 25 Dec 2023 19:25:10 +0100 |
| parents | abc5cc0aec94 |
| children | e96ee7d7ac9c |
| files | src/SbpOperators/volumeops/laplace/laplace.jl |
| diffstat | 1 files changed, 6 insertions(+), 5 deletions(-) [+] |
line wrap: on
line diff
--- a/src/SbpOperators/volumeops/laplace/laplace.jl Mon Dec 25 19:24:08 2023 +0100 +++ b/src/SbpOperators/volumeops/laplace/laplace.jl Mon Dec 25 19:25:10 2023 +0100 @@ -53,12 +53,13 @@ end laplace(g::EquidistantGrid, stencil_set) = second_derivative(g, stencil_set) +# TODO: Add sat_tensor for Diirichlet condition """ -sat_tensors(Δ::Laplace, g::TensorGrid, bc::NeumannCondition) +sat_tensors(Δ::Laplace, g::Grid, bc::NeumannCondition) -Returns the LazyTensors required to impose a Neumann condition -SAT = sat_op(d*u - g) +The operators required to construct the SAT for imposing Neumann condition + See also: [`sat`,`NeumannCondition`](@ref). """ @@ -70,6 +71,6 @@ e = boundary_restriction(g, set, id) d = normal_derivative(g, set, id) - sat_tensor = H⁻¹∘e'∘Hᵧ - return sat_tensor, d + sat_op = H⁻¹∘e'∘Hᵧ + return sat_op, d end