Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/volumeops/laplace/laplace.jl @ 1648:b7dcd3dd3181 feature/sbp_operators/laplace_curvilinear

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 26 Jun 2024 12:36:41 +0200
parents 4f6f5e5daa35
children 707fc9761c2b
comparison
equal deleted inserted replaced
1647:0685d97ebcb0 1648:b7dcd3dd3181
60 60
61 Jg = map(*, J, geometric_tensor_inverse(grid)) 61 Jg = map(*, J, geometric_tensor_inverse(grid))
62 lg = logicalgrid(grid) 62 lg = logicalgrid(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[1], I[2]) 66 Jgⁱʲ = componentview(Jg, i, j)
67 67
68 if i == j 68 if i == j
69 J⁻¹∘second_derivative_variable(lg, Jgⁱʲ, stencil_set, i) 69 J⁻¹∘second_derivative_variable(lg, Jgⁱʲ, stencil_set, i)
70 else 70 else
71 Dᵢ = first_derivative(lg, stencil_set, i) 71 Dᵢ = first_derivative(lg, stencil_set, i)