diff test/Grids/mapped_grid_test.jl @ 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 ecec2b0eea0f
children 2b5f81e288f1 d91a9f47380f
line wrap: on
line diff
--- 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