view src/SbpOperators/boundary_conditions/sat.jl @ 1949:34dcad2f2cd6 feature/grids/parameter_spaces

Close before merge
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 11 Feb 2025 08:52:45 +0100
parents 77f192b05b1d
children
line wrap: on
line source

"""
    sat_tensors(op, grid, bc::BoundaryCondition; kwargs...)

The penalty tensor and boundary operator used to construct a
simultaneous-approximation-term for imposing `bc` related to `op`.

For `penalty_tensor, L  = sat_tensors(...)` then `SAT(u,g) =
penalty_tensor*(L*u - g)`  where `g` is the boundary data.
"""
function sat_tensors end


"""
    sat(op, grid, bc::BoundaryCondition; kwargs...)

Simultaneous-Approximation-Term for a general `bc` to `op`. Returns a function
`SAT(u,g)` weakly imposing `bc` when added to `op*u`.

Internally `sat_tensors(op, grid, bc; ...)` is called to construct the
necessary parts for the SAT.
"""
function sat(op, grid, bc::BoundaryCondition; kwargs...)
    penalty_tensor, L = sat_tensors(op, grid, bc; kwargs...)
    return SAT(u, g) = penalty_tensor*(L*u - g)
end