Mercurial > repos > public > sbplib_julia

diff test/Grids/mapped_grid_test.jl @ 1741:248304f88e49 feature/grids/curvilinear

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add TBD for deciding about what jacobian_determinant should mean for Rn->Rm grids
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 11 Sep 2024 12:09:41 +0200
parents 9f2bc8ea6a07
children c9c601678a14
line wrap: on
line diff
--- a/test/Grids/mapped_grid_test.jl	Wed Sep 11 09:03:55 2024 +0200
+++ b/test/Grids/mapped_grid_test.jl	Wed Sep 11 12:09:41 2024 +0200
@@ -258,6 +258,14 @@
         2η^2 - ξ
     end
     @test jacobian_determinant(g) ≈ J
+
+
+    lg = equidistant_grid((0,0), (1,1), 11, 21)
+    x̄ = map(ξ̄ -> @SVector[ξ̄[1],ξ̄[2], ξ̄[1] + ξ̄[2]], lg)
+    J = map(ξ̄ -> @SMatrix[1 0; 0 1; 1 1], lg)
+    mg = MappedGrid(lg, x̄, J)
+
+    @test_broken jacobian(mg) isa AbstractArray{2,Float64}
 end
 
 @testset "metric_tensor" begin