Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/boundary_conditions/sat.jl @ 1635:b62770cec7d0 update/julia_1.10.3

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, 25 Jun 2024 15:32:19 +0200
parents 77f192b05b1d
children
comparison
equal deleted inserted replaced
1593:7f5e4d2a5112 1635:b62770cec7d0
1 """
2 sat_tensors(op, grid, bc::BoundaryCondition; kwargs...)
3
4 The penalty tensor and boundary operator used to construct a
5 simultaneous-approximation-term for imposing `bc` related to `op`.
6
7 For `penalty_tensor, L = sat_tensors(...)` then `SAT(u,g) =
8 penalty_tensor*(L*u - g)` where `g` is the boundary data.
9 """
10 function sat_tensors end
11
12
13 """
14 sat(op, grid, bc::BoundaryCondition; kwargs...)
15
16 Simultaneous-Approximation-Term for a general `bc` to `op`. Returns a function
17 `SAT(u,g)` weakly imposing `bc` when added to `op*u`.
18
19 Internally `sat_tensors(op, grid, bc; ...)` is called to construct the
20 necessary parts for the SAT.
21 """
22 function sat(op, grid, bc::BoundaryCondition; kwargs...)
23 penalty_tensor, L = sat_tensors(op, grid, bc; kwargs...)
24 return SAT(u, g) = penalty_tensor*(L*u - g)
25 end