Mercurial > repos > public > sbplib_julia
changeset 1585:d359d0d7fb12 feature/sbp_operators/laplace_curvilinear
Fix some typos and rename g to grid
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Fri, 26 Apr 2024 08:53:32 +0200 |
| parents | d7483e8af705 |
| children | aef3827ef522 |
| files | src/SbpOperators/volumeops/laplace/laplace.jl |
| diffstat | 1 files changed, 7 insertions(+), 7 deletions(-) [+] |
line wrap: on
line diff
--- a/src/SbpOperators/volumeops/laplace/laplace.jl Fri Apr 26 08:45:54 2024 +0200 +++ b/src/SbpOperators/volumeops/laplace/laplace.jl Fri Apr 26 08:53:32 2024 +0200 @@ -54,16 +54,16 @@ laplace(g::EquidistantGrid, stencil_set) = second_derivative(g, stencil_set) -function laplace(g::MappedGrid, stencil_set) - J = jacobian_determinant(g) - J⁻¹ = map(inv, J) +function laplace(grid::MappedGrid, stencil_set) + J = jacobian_determinant(grid) + J⁻¹ = DiagonalTensor(map(inv, J)) - Jḡ = map(*, J, ggeometric_tensor_inverse(g)) - lg = logicalgrid(g) + Jg = map(*, J, ggeometric_tensor_inverse(grid)) + lg = logicalgrid(grid) - return mapreduce(+, CartesianIndices(first(ḡ))) do I + return mapreduce(+, CartesianIndices(first(Jg))) do I i,j = I[1], I[2] - Jgⁱʲ = componentview(Jḡ, I[1], I[2]) + Jgⁱʲ = componentview(Jg, I[1], I[2]) if i == j J⁻¹∘second_derivative_variable(lg, Jgⁱʲ, stencil_set, i)