Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/volumeops/laplace/laplace.jl @ 947:38d1752a9aff feature/laplace_opset

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Reformat "See also:"
author Jonatan Werpers <jonatan@werpers.com>
date Sun, 13 Mar 2022 21:01:09 +0100
parents 22c80fb36400
children 1484073dfe27
comparison
equal deleted inserted replaced
936:22c80fb36400 947:38d1752a9aff
39 the stencil `inner_stencil` in the interior and a set of stencils `closure_stencils` 39 the stencil `inner_stencil` in the interior and a set of stencils `closure_stencils`
40 for the points in the closure regions. 40 for the points in the closure regions.
41 41
42 On a one-dimensional `grid`, `Δ` is equivalent to `second_derivative`. On a 42 On a one-dimensional `grid`, `Δ` is equivalent to `second_derivative`. On a
43 multi-dimensional `grid`, `Δ` is the sum of multi-dimensional `second_derivative`s 43 multi-dimensional `grid`, `Δ` is the sum of multi-dimensional `second_derivative`s
44 where the sum is carried out lazily. See also [`second_derivative`](@ref). 44 where the sum is carried out lazily.
45
46 See also: [`second_derivative`](@ref).
45 """ 47 """
46 function laplace(grid::EquidistantGrid, inner_stencil, closure_stencils) 48 function laplace(grid::EquidistantGrid, inner_stencil, closure_stencils)
47 Δ = second_derivative(grid, inner_stencil, closure_stencils, 1) 49 Δ = second_derivative(grid, inner_stencil, closure_stencils, 1)
48 for d = 2:dimension(grid) 50 for d = 2:dimension(grid)
49 Δ += second_derivative(grid, inner_stencil, closure_stencils, d) 51 Δ += second_derivative(grid, inner_stencil, closure_stencils, d)