Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/volumeops/laplace/laplace.jl @ 1592:4f6f5e5daa35 feature/sbp_operators/laplace_curvilinear

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
very minor edit
author Vidar Stiernström <vidar.stiernstrom@gmail.com>
date Sat, 25 May 2024 15:47:01 -0700
parents b22b9097d907
children 707fc9761c2b
comparison
equal deleted inserted replaced
1590:62cb622cbe6b 1592:4f6f5e5daa35
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)