Mercurial > repos > public > sbplib_julia
diff src/SbpOperators/volumeops/laplace/laplace.jl @ 926:47425442bbc5 feature/laplace_opset
Fix tests after refactoring
author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
---|---|
date | Mon, 21 Feb 2022 23:33:29 +0100 |
parents | 12e8e431b43c |
children | 22c80fb36400 |
line wrap: on
line diff
--- a/src/SbpOperators/volumeops/laplace/laplace.jl Mon Feb 21 13:13:37 2022 +0100 +++ b/src/SbpOperators/volumeops/laplace/laplace.jl Mon Feb 21 23:33:29 2022 +0100 @@ -1,12 +1,11 @@ """ - Laplace{T, DiffOp} <: TensorMapping{T,Dim,Dim} - Laplace(grid::Equidistant, stencil_set) + Laplace{T, Dim, DiffOp} <: TensorMapping{T, Dim, Dim} Implements the Laplace operator, approximating ∑d²/xᵢ² , i = 1,...,`Dim` as a `TensorMapping`. Additionally `Laplace` stores the stencil set (parsed from TOML) used to construct the `TensorMapping`. """ -struct Laplace{T, DiffOp<:TensorMapping{T,Dim,Dim}} <: TensorMapping{T,Dim,Dim} +struct Laplace{T, Dim, DiffOp<:TensorMapping{T, Dim, Dim}} <: TensorMapping{T, Dim, Dim} D::DiffOp# Differential operator stencil_set # Stencil set of the operator end @@ -17,7 +16,7 @@ Creates the `Laplace`` operator `Δ` on `grid` given a parsed TOML `stencil_set`. See also [`laplace`](@ref). """ -function Laplace(grid::Equidistant, stencil_set) +function Laplace(grid::EquidistantGrid, stencil_set) inner_stencil = parse_stencil(stencil_set["D2"]["inner_stencil"]) closure_stencils = parse_stencil.(stencil_set["D2"]["closure_stencils"]) Δ = laplace(grid, inner_stencil,closure_stencils) @@ -44,8 +43,8 @@ multi-dimensional `grid`, `Δ` is the sum of multi-dimensional `second_derivative`s where the sum is carried out lazily. See also [`second_derivative`](@ref). """ -function laplace(grid::Equidistant, inner_stencil, closure_stencils) - second_derivative(grid, inner_stencil, closure_stencils, 1) +function laplace(grid::EquidistantGrid, inner_stencil, closure_stencils) + Δ = second_derivative(grid, inner_stencil, closure_stencils, 1) for d = 2:dimension(grid) Δ += second_derivative(grid, inner_stencil, closure_stencils, d) end