Mercurial > repos > public > sbplib_julia

comparison Notes.md @ 1137:6757cc9ba22e feature/boundary_conditions

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Update notes
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Tue, 11 Oct 2022 18:16:53 +0200
parents 396278072f18
children ea2e8254820a
comparison
equal deleted inserted replaced
1136:470a70a6c1e6 1137:6757cc9ba22e
1 # Notes 1 # Notes
2
3 ## Boundary Conditions and SATs
4
5 Types for boundary conditions:
6
7 * abstract type `BoundaryDataType`
8 * abstract type `BoundaryCondition{T<:BoundaryDataType}`
9 * concrete types `ConstantBoundaryData <: BoundaryDataType` and similar
10 * concrete types `NeumannCondition{BDT<:BoundaryDataType} <: BoundaryCondition{BDT}` and similar
11 The concrete `BoundaryDataType` subtypes are "thin types" wrapping the boundary data, and are used to indicate how the boundary data should be used in e.g. sat routines. The concrete `BoundaryCondition{BDT}` subtypes are used for assembling the tensors used to construct e.g. a SAT.
12
13 SAT methods:
14 There are multiple options for what the SAT methods could return.
15 * (Current) a function which returns a `LazyTensorApplication`, e.g. `f = sat(grid,op,bc)`. The the resulting `LazyTensorApplication` can then be added to scheme i.e. `scheme = op*u + f(u)`. This is how one typically would write the SAT in the litterature. Depdending on the type of data in the BC, e.g. time-depdendent etc one can return f(u,t).
16 * `LazyTensor`s `closure, penalty = sat(grid,op,bc)` like in the matlab version. Probably the most general one. Up to the user to make use of the returned `LazyTensor`s. One can for example then easily include the closures to the operator and have eg. `D = (op + closure)*u`.
17 * A `LazyTensor` for closure, and a `LazyArray` for `penalty*data`. Mix of the above.
18
19 It is not clear if any of these are preferable as it currently stands.
2 20
3 ## Reading operators 21 ## Reading operators
4 22
5 Jonatan's suggestion is to add methods to `Laplace`, `SecondDerivative` and 23 Jonatan's suggestion is to add methods to `Laplace`, `SecondDerivative` and
6 similar functions that take in a filename from which to read stencils. These 24 similar functions that take in a filename from which to read stencils. These