Mercurial > repos > public > sbplib_julia

changeset 1561:6fdc81860b0c feature/sbp_operators/laplace_curvilinear

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 25 Apr 2024 09:03:35 +0200
parents 43e6acbefdd1 (current diff) 5d32ecb98db8 (diff)
children efa994405c38
files
diffstat 3 files changed, 34 insertions(+), 0 deletions(-) [+]
line wrap: on
line diff
--- a/src/Grids/Grids.jl	Thu Apr 25 08:59:45 2024 +0200
+++ b/src/Grids/Grids.jl	Thu Apr 25 09:03:35 2024 +0200
@@ -38,6 +38,9 @@
 export jacobian
 export logicalgrid
 export mapped_grid
+export jacobian_determinant
+export geometric_tensor
+export geometric_tensor_inverse
 
 abstract type BoundaryIdentifier end
 
--- a/src/Grids/mapped_grid.jl	Thu Apr 25 08:59:45 2024 +0200
+++ b/src/Grids/mapped_grid.jl	Thu Apr 25 09:03:35 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
+
--- a/test/Grids/mapped_grid_test.jl	Thu Apr 25 08:59:45 2024 +0200
+++ b/test/Grids/mapped_grid_test.jl	Thu Apr 25 09:03:35 2024 +0200
@@ -157,6 +157,18 @@
         @testset test_boundary_grid(mg, TensorGridBoundary{2, Upper}(), J1)
     end
 
+    @testset "jacobian_determinant" begin
+        @test_broken false
+    end
+
+    @testset "geometric_tensor" begin
+        @test_broken false
+    end
+
+    @testset "geometric_tensor_inverse" begin
+        @test_broken false
+    end
+
 end
 
 @testset "mapped_grid" begin