Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/volumeops/laplace/laplace.jl @ 1786:1f42944d4a72 feature/sbp_operators/laplace_curvilinear

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix errors introduced by merge
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 16 Sep 2024 11:03:37 +0200
parents f3d7e2d7a43f
children 835b1dcee38e
comparison
equal deleted inserted replaced
1785:96f8cad255b4 1786:1f42944d4a72
53 end 53 end
54 54
55 laplace(g::EquidistantGrid, stencil_set) = second_derivative(g, stencil_set) 55 laplace(g::EquidistantGrid, stencil_set) = second_derivative(g, stencil_set)
56 56
57 function laplace(grid::MappedGrid, stencil_set) 57 function laplace(grid::MappedGrid, stencil_set)
58 J = jacobian_determinant(grid) 58 J = map(det,jacobian(grid))
59 J⁻¹ = DiagonalTensor(map(inv, J)) 59 J⁻¹ = DiagonalTensor(map(inv, J))
60 60
61 Jg = map(*, J, metric_tensor_inverse(grid)) 61 Jg = map(*, J, metric_tensor_inverse(grid))
62 lg = logicalgrid(grid) 62 lg = logical_grid(grid)
63 63
64 return mapreduce(+, CartesianIndices(first(Jg))) do I 64 return mapreduce(+, CartesianIndices(first(Jg))) do I
65 i, j = I[1], I[2] 65 i, j = I[1], I[2]
66 Jgⁱʲ = componentview(Jg, i, j) 66 Jgⁱʲ = componentview(Jg, i, j)
67 67