Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/volumeops/derivatives/second_derivative.jl @ 925:6b47a9ee1632 feature/laplace_opset

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add functionality for creating operators from a stencil set. Note: Tests are not updated yet!
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Mon, 21 Feb 2022 13:13:37 +0100
parents b4acd25943f4
children 22c80fb36400
comparison
equal deleted inserted replaced
924:12e8e431b43c 925:6b47a9ee1632
14 function second_derivative(grid::EquidistantGrid, inner_stencil, closure_stencils, direction) 14 function second_derivative(grid::EquidistantGrid, inner_stencil, closure_stencils, direction)
15 h_inv = inverse_spacing(grid)[direction] 15 h_inv = inverse_spacing(grid)[direction]
16 return SbpOperators.volume_operator(grid, scale(inner_stencil,h_inv^2), scale.(closure_stencils,h_inv^2), even, direction) 16 return SbpOperators.volume_operator(grid, scale(inner_stencil,h_inv^2), scale.(closure_stencils,h_inv^2), even, direction)
17 end 17 end
18 second_derivative(grid::EquidistantGrid{1}, inner_stencil, closure_stencils) = second_derivative(grid,inner_stencil,closure_stencils,1) 18 second_derivative(grid::EquidistantGrid{1}, inner_stencil, closure_stencils) = second_derivative(grid,inner_stencil,closure_stencils,1)
19 export second_derivative