changeset 1787:43c0bfc13de3 feature/grids/curvilinear

Delete metric_tensor_inverse
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 18 Sep 2024 22:39:38 +0200
parents ed1b43d07029
children 3e4ae9bfe86a
files src/Grids/Grids.jl src/Grids/mapped_grid.jl test/Grids/mapped_grid_test.jl
diffstat 3 files changed, 0 insertions(+), 31 deletions(-) [+]
line wrap: on
line diff
--- a/src/Grids/Grids.jl	Mon Sep 16 09:46:26 2024 +0200
+++ b/src/Grids/Grids.jl	Wed Sep 18 22:39:38 2024 +0200
@@ -44,7 +44,6 @@
 export logical_grid
 export mapped_grid
 export metric_tensor
-export metric_tensor_inverse
 
 include("grid.jl")
 include("tensor_grid.jl")
--- a/src/Grids/mapped_grid.jl	Mon Sep 16 09:46:26 2024 +0200
+++ b/src/Grids/mapped_grid.jl	Wed Sep 18 22:39:38 2024 +0200
@@ -133,17 +133,6 @@
     end
 end
 
-"""
-    metric_tensor_inverse(g::MappedGrid)
-
-The inverse of the metric tensor of `g` as a grid function.
-"""
-function metric_tensor_inverse(g::MappedGrid)
-    return map(jacobian(g)) do ∂x∂ξ
-        inv(∂x∂ξ'*∂x∂ξ)
-    end
-end
-
 function min_spacing(g::MappedGrid{T,1} where T)
     n, = size(g)
 
--- a/test/Grids/mapped_grid_test.jl	Mon Sep 16 09:46:26 2024 +0200
+++ b/test/Grids/mapped_grid_test.jl	Wed Sep 18 22:39:38 2024 +0200
@@ -302,25 +302,6 @@
     @test metric_tensor(g) ≈ G
 end
 
-@testset "metric_tensor_inverse" begin
-    x̄((ξ, η)) = @SVector[ξ + ξ^2/2, η + η^2 + ξ^2/2]
-    J((ξ, η)) = @SMatrix[
-        1+ξ   0;
-        ξ    1+η;
-    ]
-
-    g = mapped_grid(x̄, J, 10, 11)
-    G⁻¹ = map(logical_grid(g)) do (ξ,η)
-        @SMatrix[
-            (1+η)^2  -ξ*(1+η);
-            -ξ*(1+η) (1+ξ)^2+ξ^2;
-        ]/(((1+ξ)^2+ξ^2)*(1+η)^2 - ξ^2*(1+η)^2)
-
-    end
-
-    @test metric_tensor_inverse(g) ≈ G⁻¹
-end
-
 @testset "min_spacing" begin
     let g = mapped_grid(identity, x->@SMatrix[1], 11)
         @test min_spacing(g) ≈ 0.1