Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/volumeops/laplace/laplace.jl @ 1040:7fc8df5157a7 refactor/lazy_tensors

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, 22 Mar 2022 14:23:55 +0100
parents 1ba8a398af9c 99d1f5651d0b
children b4ee47f2aafb 6530fceef37c
comparison
equal deleted inserted replaced
1037:9e76bf19904c 1040:7fc8df5157a7
1 """ 1 """
2 Laplace{T, Dim, TM} <: LazyTensor{T, Dim, Dim} 2 Laplace{T, Dim, TM} <: LazyTensor{T, Dim, Dim}
3 3
4 Implements the Laplace operator, approximating ∑d²/xᵢ² , i = 1,...,`Dim` as a 4 Implements the Laplace operator, approximating ∑d²/xᵢ² , i = 1,...,`Dim` as a
5 `LazyTensor`. Additionally `Laplace` stores the stencil set (parsed from TOML) 5 `LazyTensor`. Additionally `Laplace` stores the `StencilSet`
6 used to construct the `LazyTensor`. 6 used to construct the `LazyTensor `.
7 """ 7 """
8 struct Laplace{T, Dim, TM<:LazyTensor{T, Dim, Dim}} <: LazyTensor{T, Dim, Dim} 8 struct Laplace{T, Dim, TM<:LazyTensor{T, Dim, Dim}} <: LazyTensor{T, Dim, Dim}
9 D::TM # Difference operator 9 D::TM # Difference operator
10 stencil_set # Stencil set of the operator 10 stencil_set::StencilSet # Stencil set of the operator
11 end 11 end
12 12
13 """ 13 """
14 Laplace(grid::Equidistant, stencil_set) 14 Laplace(grid::Equidistant, stencil_set)
15 15
16 Creates the `Laplace` operator `Δ` on `grid` given a parsed TOML 16 Creates the `Laplace` operator `Δ` on `grid` given a `stencil_set`.
17 `stencil_set`. See also [`laplace`](@ref). 17
18 See also [`laplace`](@ref).
18 """ 19 """
19 function Laplace(grid::EquidistantGrid, stencil_set) 20 function Laplace(grid::EquidistantGrid, stencil_set::StencilSet)
20 inner_stencil = parse_stencil(stencil_set["D2"]["inner_stencil"]) 21 inner_stencil = parse_stencil(stencil_set["D2"]["inner_stencil"])
21 closure_stencils = parse_stencil.(stencil_set["D2"]["closure_stencils"]) 22 closure_stencils = parse_stencil.(stencil_set["D2"]["closure_stencils"])
22 Δ = laplace(grid, inner_stencil,closure_stencils) 23 Δ = laplace(grid, inner_stencil,closure_stencils)
23 return Laplace(Δ,stencil_set) 24 return Laplace(Δ,stencil_set)
24 end 25 end