sbplib_julia: test/Grids/mapped_grid

Simplify tests of boundary_grid

author	Jonatan Werpers <jonatan@werpers.com>
date	Tue, 10 Sep 2024 22:17:01 +0200
parents	863385aae454
children	7573c9244e47

comparison

equal deleted inserted replaced

-:863385aae454
+:a4e1721a7109
 J2((ξ, η)) = @SMatrix[
 0;
 1+ξ*(ξ-1);
 ]
-function test_boundary_grid(mg, bId, Jb)
+function expected_bg(mg, bId, Jb)
-bg = boundary_grid(mg, bId)
 lg = logicalgrid(mg)
-expected_bg = MappedGrid(
+return MappedGrid(
 boundary_grid(lg, bId),
 map(x̄, boundary_grid(lg, bId)),
 map(Jb, boundary_grid(lg, bId)),
 )
+end
-@testset let bId=bId, bg=bg, expected_bg=expected_bg
-@test collect(bg) == collect(expected_bg)
+let bid = TensorGridBoundary{1, LowerBoundary}()
-@test logicalgrid(bg) == logicalgrid(expected_bg)
+@test boundary_grid(mg, bid) == expected_bg(mg, bid, J2)
-@test jacobian(bg) == jacobian(expected_bg)
+end
-# TODO: Implement equality of a curvilinear grid and simlify the above
-end
+let bid = TensorGridBoundary{1, UpperBoundary}()
-end
+@test boundary_grid(mg, bid) == expected_bg(mg, bid, J2)
+end
-@testset test_boundary_grid(mg, TensorGridBoundary{1, LowerBoundary}(), J2)
-@testset test_boundary_grid(mg, TensorGridBoundary{1, UpperBoundary}(), J2)
+let bid = TensorGridBoundary{2, LowerBoundary}()
-@testset test_boundary_grid(mg, TensorGridBoundary{2, LowerBoundary}(), J1)
+@test boundary_grid(mg, bid) == expected_bg(mg, bid, J1)
-@testset test_boundary_grid(mg, TensorGridBoundary{2, UpperBoundary}(), J1)
+end
+let bid = TensorGridBoundary{2, UpperBoundary}()
+@test boundary_grid(mg, bid) == expected_bg(mg, bid, J1)
+end
 end
 end
 @testset "mapped_grid" begin
 x̄, J = _partially_curved_mapping()

Mercurial > repos > public > sbplib_julia