Mercurial > repos > public > sbplib_julia

diff src/Grids/mapped_grid.jl @ 1560:5d32ecb98db8 feature/grids/curvilinear

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add functions for jacobian_determiant, geometric_tensor, geometric_tensor_inverse
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 25 Apr 2024 09:03:04 +0200
parents 69790e9d1652
children 64baaf29ae4e 063a2bfb03da
line wrap: on
line diff
--- a/src/Grids/mapped_grid.jl	Sat Apr 13 23:52:40 2024 +0200
+++ b/src/Grids/mapped_grid.jl	Thu Apr 25 09:03:04 2024 +0200
@@ -60,3 +60,22 @@
         map(J,lg),
     )
 end
+
+function jacobian_determinant(g::MappedGrid)
+    return map(jacobian(g)) do ∂x∂ξ
+        det(∂x∂ξ)
+    end
+end
+
+function geometric_tensor(g::MappedGrid)
+    return map(jacobian(g)) do ∂x∂ξ
+        ∂x∂ξ'*∂x∂ξ
+    end
+end
+
+function geometric_tensor_inverse(g::MappedGrid)
+    return map(jacobian(g)) do ∂x∂ξ
+        inv(∂x∂ξ'*∂x∂ξ)
+    end
+end
+