Mercurial > repos > public > sbplib_julia
diff src/SbpOperators/volumeops/laplace/laplace.jl @ 1396:35840a0681d1 feature/boundary_conditions
Start drafting new implemenentation of boundary conditions
| author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
|---|---|
| date | Wed, 26 Jul 2023 23:11:02 +0200 |
| parents | bdcdbd4ea9cd |
| children | 8d60d045c2a2 |
line wrap: on
line diff
--- a/src/SbpOperators/volumeops/laplace/laplace.jl Wed Jul 26 21:35:50 2023 +0200 +++ b/src/SbpOperators/volumeops/laplace/laplace.jl Wed Jul 26 23:11:02 2023 +0200 @@ -57,8 +57,8 @@ """ sat_tensors(Δ::Laplace, g::TensorGrid, bc::NeumannCondition) -Returns anonymous functions for construction the `LazyTensorApplication`s -recuired in order to impose a Neumann boundary condition. +Returns the LazyTensors required to impose a Neumann condition +SAT = sat_op(d*u - g) See also: [`sat`,`NeumannCondition`](@ref). """ @@ -70,7 +70,6 @@ e = boundary_restriction(g, set, id) d = normal_derivative(g, set, id) - closure(u) = H⁻¹*e'*Hᵧ*d*u - penalty(g) = -H⁻¹*e'*Hᵧ*g - return closure, penalty + sat_tensor = H⁻¹∘e'∘Hᵧ + return sat_tensor, d end