Mercurial > repos > public > sbplib_julia
diff src/SbpOperators/volumeops/laplace/laplace.jl @ 1786:1f42944d4a72 feature/sbp_operators/laplace_curvilinear
Fix errors introduced by merge
author | Jonatan Werpers <jonatan@werpers.com> |
---|---|
date | Mon, 16 Sep 2024 11:03:37 +0200 |
parents | f3d7e2d7a43f |
children | 835b1dcee38e |
line wrap: on
line diff
--- a/src/SbpOperators/volumeops/laplace/laplace.jl Mon Sep 16 10:33:47 2024 +0200 +++ b/src/SbpOperators/volumeops/laplace/laplace.jl Mon Sep 16 11:03:37 2024 +0200 @@ -55,11 +55,11 @@ laplace(g::EquidistantGrid, stencil_set) = second_derivative(g, stencil_set) function laplace(grid::MappedGrid, stencil_set) - J = jacobian_determinant(grid) + J = map(det,jacobian(grid)) J⁻¹ = DiagonalTensor(map(inv, J)) Jg = map(*, J, metric_tensor_inverse(grid)) - lg = logicalgrid(grid) + lg = logical_grid(grid) return mapreduce(+, CartesianIndices(first(Jg))) do I i, j = I[1], I[2]