Mercurial > repos > public > sbplib_julia

diff src/SbpOperators/volumeops/laplace/laplace.jl @ 1953:835b1dcee38e feature/sbp_operators/laplace_curvilinear

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Replace metric_tensor_inverse calls
author Jonatan Werpers <jonatan@werpers.com>
date Sat, 08 Feb 2025 09:35:13 +0100
parents 1f42944d4a72
children
line wrap: on
line diff
--- 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