Mercurial > repos > public > sbplib_julia
changeset 1953:835b1dcee38e feature/sbp_operators/laplace_curvilinear
Replace metric_tensor_inverse calls
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Sat, 08 Feb 2025 09:35:13 +0100 |
| parents | 847286c70d7c |
| children | b0915f43b122 |
| files | src/SbpOperators/boundaryops/normal_derivative.jl src/SbpOperators/volumeops/laplace/laplace.jl |
| diffstat | 2 files changed, 5 insertions(+), 2 deletions(-) [+] |
line wrap: on
line diff
--- a/src/SbpOperators/boundaryops/normal_derivative.jl Fri Feb 07 23:56:48 2025 +0100 +++ b/src/SbpOperators/boundaryops/normal_derivative.jl Sat Feb 08 09:35:13 2025 +0100 @@ -34,7 +34,7 @@ b_indices = boundary_indices(g, boundary) # Compute the weights for the logical derivatives - g⁻¹ = metric_tensor_inverse(g) + g⁻¹ = map(inv, metric_tensor(g)) α = map(CartesianIndices(g⁻¹)[b_indices...]) do I # TODO: Fix iterator here gᵏⁱ = g⁻¹[I][k,:] gᵏᵏ = g⁻¹[I][k,k]
--- a/src/SbpOperators/volumeops/laplace/laplace.jl Fri Feb 07 23:56:48 2025 +0100 +++ b/src/SbpOperators/volumeops/laplace/laplace.jl Sat Feb 08 09:35:13 2025 +0100 @@ -58,7 +58,10 @@ J = map(det,jacobian(grid)) J⁻¹ = DiagonalTensor(map(inv, J)) - Jg = map(*, J, metric_tensor_inverse(grid)) + Jg = map(J, metric_tensor(grid)) do Jₓ, gₓ + Jₓ*inv(gₓ) + end + lg = logical_grid(grid) return mapreduce(+, CartesianIndices(first(Jg))) do I