Mercurial > repos > public > sbplib_julia

annotate test/Grids/mapped_grid_test.jl @ 1691:5bf4a35a78c5 feature/grids/curvilinear

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Factor out functions for the mappings used in tests of MappedGrid
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 23 Aug 2024 12:50:08 +0200
parents 5eabe1f560f0
children a4c52ae93b11 6eb5b48607e0
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
1426
9c689a627244 Add src and test file for curvilinear grid type
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
1 using Sbplib.Grids
1451
2e08f3444354 Implement boundary_identifiers, and a wishful implementation of boundary_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1450
diff changeset
2 using Sbplib.RegionIndices
1426
9c689a627244 Add src and test file for curvilinear grid type
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
3 using Test
1430
9fc3c1af33e5 Add testsets and a few tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1426
diff changeset
4 using StaticArrays
1688
72776d3d5fd6 Add min_spacing for 2D mapped grids
Jonatan Werpers <jonatan@werpers.com>
parents: 1684
diff changeset
5 using LinearAlgebra
1426
9c689a627244 Add src and test file for curvilinear grid type
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
6
1691
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
7
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
8 _skew_mapping(a,b) = (ξ̄->ξ̄[1]*a + ξ̄[2]*b, ξ̄->[a b])
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
9
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
10 function _partially_curved_mapping()
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
11 x̄((ξ, η)) = @SVector[ξ, η*(1+ξ*(ξ-1))]
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
12 J((ξ, η)) = @SMatrix[
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
13 1 0;
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
14 η*(2ξ-1) 1+ξ*(ξ-1);
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
15 ]
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
16
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
17 return x̄, J
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
18 end
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
19
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
20 function _fully_curved_mapping()
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
21 x̄((ξ, η)) = @SVector[2ξ + η*(1-η), 3η+(1+η/2)*ξ^2]
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
22 J((ξ, η)) = @SMatrix[
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
23 2 1-2η;
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
24 (2+η)*ξ 3+1/2*ξ^2;
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
25 ]
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
26
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
27 return x̄, J
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
28 end
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
29
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
30 @testset "MappedGrid" begin
1570
063a2bfb03da Fix bugs in calls to equidistant_grid after merge with default
Jonatan Werpers <jonatan@werpers.com>
parents: 1560
diff changeset
31 lg = equidistant_grid((0,0), (1,1), 11, 11) # TODO: Change dims of the grid to be different
1432
64b60b42d367 Implement indexing interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1431
diff changeset
32 x̄ = map(ξ̄ -> 2ξ̄, lg)
64b60b42d367 Implement indexing interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1431
diff changeset
33 J = map(ξ̄ -> @SArray(fill(2., 2, 2)), lg)
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
34 mg = MappedGrid(lg, x̄, J)
1430
9fc3c1af33e5 Add testsets and a few tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1426
diff changeset
35
1434
af73340a8f0e Add todos
Jonatan Werpers <jonatan@werpers.com>
parents: 1432
diff changeset
36 # TODO: Test constructor for different dims of range and domain for the coordinates
af73340a8f0e Add todos
Jonatan Werpers <jonatan@werpers.com>
parents: 1432
diff changeset
37 # TODO: Test constructor with different type than TensorGrid. a dummy type?
af73340a8f0e Add todos
Jonatan Werpers <jonatan@werpers.com>
parents: 1432
diff changeset
38
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
39 @test_broken false # @test_throws ArgumentError("Sizes must match") MappedGrid(lg, map(ξ̄ -> @SArray[ξ̄[1], ξ̄[2], -ξ̄[1]], lg), rand(SMatrix{2,3,Float64},15,11))
1434
af73340a8f0e Add todos
Jonatan Werpers <jonatan@werpers.com>
parents: 1432
diff changeset
40
af73340a8f0e Add todos
Jonatan Werpers <jonatan@werpers.com>
parents: 1432
diff changeset
41
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
42 @test mg isa Grid{SVector{2, Float64},2}
1430
9fc3c1af33e5 Add testsets and a few tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1426
diff changeset
43
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
44 @test jacobian(mg) isa Array{<:AbstractMatrix}
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
45 @test logicalgrid(mg) isa Grid
1430
9fc3c1af33e5 Add testsets and a few tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1426
diff changeset
46
1432
64b60b42d367 Implement indexing interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1431
diff changeset
47 @testset "Indexing Interface" begin
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
48 mg = MappedGrid(lg, x̄, J)
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
49 @test mg[1,1] == [0.0, 0.0]
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
50 @test mg[4,2] == [0.6, 0.2]
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
51 @test mg[6,10] == [1., 1.8]
1430
9fc3c1af33e5 Add testsets and a few tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1426
diff changeset
52
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
53 @test mg[begin, begin] == [0.0, 0.0]
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
54 @test mg[end,end] == [2.0, 2.0]
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
55 @test mg[begin,end] == [0., 2.]
1432
64b60b42d367 Implement indexing interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1431
diff changeset
56
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
57 @test eachindex(mg) == CartesianIndices((11,11))
1432
64b60b42d367 Implement indexing interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1431
diff changeset
58
64b60b42d367 Implement indexing interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1431
diff changeset
59 @testset "cartesian indexing" begin
64b60b42d367 Implement indexing interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1431
diff changeset
60 cases = [
64b60b42d367 Implement indexing interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1431
diff changeset
61 (1,1) ,
64b60b42d367 Implement indexing interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1431
diff changeset
62 (3,5) ,
64b60b42d367 Implement indexing interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1431
diff changeset
63 (10,6),
64b60b42d367 Implement indexing interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1431
diff changeset
64 (1,1) ,
64b60b42d367 Implement indexing interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1431
diff changeset
65 (3,2) ,
64b60b42d367 Implement indexing interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1431
diff changeset
66 ]
1430
9fc3c1af33e5 Add testsets and a few tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1426
diff changeset
67
1432
64b60b42d367 Implement indexing interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1431
diff changeset
68 @testset "i = $is" for (lg, is) ∈ cases
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
69 @test mg[CartesianIndex(is...)] == mg[is...]
1432
64b60b42d367 Implement indexing interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1431
diff changeset
70 end
64b60b42d367 Implement indexing interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1431
diff changeset
71 end
64b60b42d367 Implement indexing interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1431
diff changeset
72
64b60b42d367 Implement indexing interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1431
diff changeset
73 @testset "eachindex" begin
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
74 @test eachindex(mg) == CartesianIndices((11,11))
1432
64b60b42d367 Implement indexing interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1431
diff changeset
75 end
1430
9fc3c1af33e5 Add testsets and a few tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1426
diff changeset
76
1432
64b60b42d367 Implement indexing interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1431
diff changeset
77 @testset "firstindex" begin
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
78 @test firstindex(mg, 1) == 1
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
79 @test firstindex(mg, 2) == 1
1432
64b60b42d367 Implement indexing interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1431
diff changeset
80 end
64b60b42d367 Implement indexing interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1431
diff changeset
81
64b60b42d367 Implement indexing interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1431
diff changeset
82 @testset "lastindex" begin
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
83 @test lastindex(mg, 1) == 11
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
84 @test lastindex(mg, 2) == 11
1432
64b60b42d367 Implement indexing interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1431
diff changeset
85 end
1430
9fc3c1af33e5 Add testsets and a few tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1426
diff changeset
86 end
1432
64b60b42d367 Implement indexing interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1431
diff changeset
87 # TODO: Test with different types of logical grids
1430
9fc3c1af33e5 Add testsets and a few tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1426
diff changeset
88
9fc3c1af33e5 Add testsets and a few tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1426
diff changeset
89 @testset "Iterator interface" begin
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
90 sg = MappedGrid(
1570
063a2bfb03da Fix bugs in calls to equidistant_grid after merge with default
Jonatan Werpers <jonatan@werpers.com>
parents: 1560
diff changeset
91 equidistant_grid((0,0), (1,1), 15, 11),
1450
647c8b18b84f Implement iterator interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1434
diff changeset
92 map(ξ̄ -> @SArray[ξ̄[1], ξ̄[2], -ξ̄[1]], lg), rand(SMatrix{2,3,Float64},15,11)
647c8b18b84f Implement iterator interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1434
diff changeset
93 )
647c8b18b84f Implement iterator interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1434
diff changeset
94
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
95 @test eltype(mg) == SVector{2,Float64}
1450
647c8b18b84f Implement iterator interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1434
diff changeset
96 @test eltype(sg) == SVector{3,Float64}
647c8b18b84f Implement iterator interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1434
diff changeset
97
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
98 @test eltype(typeof(mg)) == SVector{2,Float64}
1450
647c8b18b84f Implement iterator interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1434
diff changeset
99 @test eltype(typeof(sg)) == SVector{3,Float64}
647c8b18b84f Implement iterator interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1434
diff changeset
100
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
101 @test size(mg) == (11,11)
1450
647c8b18b84f Implement iterator interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1434
diff changeset
102 @test size(sg) == (15,11)
647c8b18b84f Implement iterator interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1434
diff changeset
103
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
104 @test size(mg,2) == 11
1450
647c8b18b84f Implement iterator interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1434
diff changeset
105 @test size(sg,2) == 11
1430
9fc3c1af33e5 Add testsets and a few tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1426
diff changeset
106
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
107 @test length(mg) == 121
1450
647c8b18b84f Implement iterator interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1434
diff changeset
108 @test length(sg) == 165
647c8b18b84f Implement iterator interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1434
diff changeset
109
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
110 @test Base.IteratorSize(mg) == Base.HasShape{2}()
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
111 @test Base.IteratorSize(typeof(mg)) == Base.HasShape{2}()
1450
647c8b18b84f Implement iterator interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1434
diff changeset
112
647c8b18b84f Implement iterator interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1434
diff changeset
113 @test Base.IteratorSize(sg) == Base.HasShape{2}()
647c8b18b84f Implement iterator interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1434
diff changeset
114 @test Base.IteratorSize(typeof(sg)) == Base.HasShape{2}()
1430
9fc3c1af33e5 Add testsets and a few tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1426
diff changeset
115
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
116 element, state = iterate(mg)
1450
647c8b18b84f Implement iterator interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1434
diff changeset
117 @test element == lg[1,1].*2
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
118 element, _ = iterate(mg, state)
1450
647c8b18b84f Implement iterator interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1434
diff changeset
119 @test element == lg[2,1].*2
1430
9fc3c1af33e5 Add testsets and a few tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1426
diff changeset
120
1450
647c8b18b84f Implement iterator interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1434
diff changeset
121 element, state = iterate(sg)
647c8b18b84f Implement iterator interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1434
diff changeset
122 @test element == sg.physicalcoordinates[1,1]
647c8b18b84f Implement iterator interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1434
diff changeset
123 element, _ = iterate(sg, state)
647c8b18b84f Implement iterator interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1434
diff changeset
124 @test element == sg.physicalcoordinates[2,1]
647c8b18b84f Implement iterator interface
Jonatan Werpers <jonatan@werpers.com>
parents: 1434
diff changeset
125
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
126 @test collect(mg) == 2 .* lg
1430
9fc3c1af33e5 Add testsets and a few tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1426
diff changeset
127 end
9fc3c1af33e5 Add testsets and a few tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1426
diff changeset
128
9fc3c1af33e5 Add testsets and a few tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1426
diff changeset
129 @testset "Base" begin
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
130 @test ndims(mg) == 2
1430
9fc3c1af33e5 Add testsets and a few tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1426
diff changeset
131 end
9fc3c1af33e5 Add testsets and a few tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1426
diff changeset
132
9fc3c1af33e5 Add testsets and a few tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1426
diff changeset
133 @testset "boundary_identifiers" begin
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
134 @test boundary_identifiers(mg) == boundary_identifiers(lg)
1430
9fc3c1af33e5 Add testsets and a few tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1426
diff changeset
135 end
9fc3c1af33e5 Add testsets and a few tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1426
diff changeset
136
1493
58b8da9c7e56 Implement boundary_indices
Jonatan Werpers <jonatan@werpers.com>
parents: 1451
diff changeset
137 @testset "boundary_indices" begin
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
138 @test boundary_indices(mg, CartesianBoundary{1,Lower}()) == boundary_indices(lg,CartesianBoundary{1,Lower}())
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
139 @test boundary_indices(mg, CartesianBoundary{2,Lower}()) == boundary_indices(lg,CartesianBoundary{2,Lower}())
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
140 @test boundary_indices(mg, CartesianBoundary{1,Upper}()) == boundary_indices(lg,CartesianBoundary{1,Upper}())
1493
58b8da9c7e56 Implement boundary_indices
Jonatan Werpers <jonatan@werpers.com>
parents: 1451
diff changeset
141 end
58b8da9c7e56 Implement boundary_indices
Jonatan Werpers <jonatan@werpers.com>
parents: 1451
diff changeset
142
1430
9fc3c1af33e5 Add testsets and a few tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1426
diff changeset
143 @testset "boundary_grid" begin
1691
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
144 x̄, J = _partially_curved_mapping()
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
145 mg = mapped_grid(x̄, J, 10, 11)
1503
704a84eef8b6 Add tests for boundary_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1502
diff changeset
146 J1((ξ, η)) = @SMatrix[
704a84eef8b6 Add tests for boundary_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1502
diff changeset
147 1 ;
704a84eef8b6 Add tests for boundary_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1502
diff changeset
148 η*(2ξ-1);
704a84eef8b6 Add tests for boundary_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1502
diff changeset
149 ]
704a84eef8b6 Add tests for boundary_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1502
diff changeset
150 J2((ξ, η)) = @SMatrix[
704a84eef8b6 Add tests for boundary_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1502
diff changeset
151 0;
704a84eef8b6 Add tests for boundary_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1502
diff changeset
152 1+ξ*(ξ-1);
704a84eef8b6 Add tests for boundary_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1502
diff changeset
153 ]
704a84eef8b6 Add tests for boundary_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1502
diff changeset
154
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
155 function test_boundary_grid(mg, bId, Jb)
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
156 bg = boundary_grid(mg, bId)
1503
704a84eef8b6 Add tests for boundary_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1502
diff changeset
157
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
158 lg = logicalgrid(mg)
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
159 expected_bg = MappedGrid(
1503
704a84eef8b6 Add tests for boundary_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1502
diff changeset
160 boundary_grid(lg, bId),
704a84eef8b6 Add tests for boundary_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1502
diff changeset
161 map(x̄, boundary_grid(lg, bId)),
704a84eef8b6 Add tests for boundary_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1502
diff changeset
162 map(Jb, boundary_grid(lg, bId)),
704a84eef8b6 Add tests for boundary_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1502
diff changeset
163 )
704a84eef8b6 Add tests for boundary_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1502
diff changeset
164
704a84eef8b6 Add tests for boundary_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1502
diff changeset
165 @testset let bId=bId, bg=bg, expected_bg=expected_bg
704a84eef8b6 Add tests for boundary_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1502
diff changeset
166 @test collect(bg) == collect(expected_bg)
704a84eef8b6 Add tests for boundary_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1502
diff changeset
167 @test logicalgrid(bg) == logicalgrid(expected_bg)
704a84eef8b6 Add tests for boundary_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1502
diff changeset
168 @test jacobian(bg) == jacobian(expected_bg)
704a84eef8b6 Add tests for boundary_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1502
diff changeset
169 # TODO: Implement equality of a curvilinear grid and simlify the above
704a84eef8b6 Add tests for boundary_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1502
diff changeset
170 end
704a84eef8b6 Add tests for boundary_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1502
diff changeset
171 end
704a84eef8b6 Add tests for boundary_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1502
diff changeset
172
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
173 @testset test_boundary_grid(mg, TensorGridBoundary{1, Lower}(), J2)
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
174 @testset test_boundary_grid(mg, TensorGridBoundary{1, Upper}(), J2)
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
175 @testset test_boundary_grid(mg, TensorGridBoundary{2, Lower}(), J1)
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
176 @testset test_boundary_grid(mg, TensorGridBoundary{2, Upper}(), J1)
1430
9fc3c1af33e5 Add testsets and a few tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1426
diff changeset
177 end
1426
9c689a627244 Add src and test file for curvilinear grid type
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
178 end
1502
a2dc80396808 Add functions `curvilinear_grid`
Jonatan Werpers <jonatan@werpers.com>
parents: 1494
diff changeset
179
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
180 @testset "mapped_grid" begin
1691
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
181 x̄, J = _partially_curved_mapping()
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
182 mg = mapped_grid(x̄, J, 10, 11)
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
183 @test mg isa MappedGrid{SVector{2,Float64}, 2}
1502
a2dc80396808 Add functions `curvilinear_grid`
Jonatan Werpers <jonatan@werpers.com>
parents: 1494
diff changeset
184
1570
063a2bfb03da Fix bugs in calls to equidistant_grid after merge with default
Jonatan Werpers <jonatan@werpers.com>
parents: 1560
diff changeset
185 lg = equidistant_grid((0,0), (1,1), 10, 11)
1506
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
186 @test logicalgrid(mg) == lg
535f32316637 Rename from curvilinear to mapped
Jonatan Werpers <jonatan@werpers.com>
parents: 1503
diff changeset
187 @test collect(mg) == map(x̄, lg)
1690
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
188 end
1659
3bbcd496e021 Add function for computing the normal at the boundary of a mapped grid as a grid function
Jonatan Werpers <jonatan@werpers.com>
parents: 1570
diff changeset
189
1690
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
190 @testset "jacobian_determinant" begin
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
191 @test_broken false
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
192 end
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
193
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
194 @testset "metric_tensor" begin
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
195 @test_broken false
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
196 end
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
197
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
198 @testset "metric_tensor_inverse" begin
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
199 @test_broken false
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
200 end
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
201
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
202 @testset "min_spacing" begin
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
203 let g = mapped_grid(identity, x->@SMatrix[1], 11)
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
204 @test min_spacing(g) ≈ 0.1
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
205 end
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
206
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
207 let g = mapped_grid(x->x+x.^2/2, x->@SMatrix[1 .+ x], 11)
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
208 @test min_spacing(g) ≈ 0.105
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
209 end
1659
3bbcd496e021 Add function for computing the normal at the boundary of a mapped grid as a grid function
Jonatan Werpers <jonatan@werpers.com>
parents: 1570
diff changeset
210
1690
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
211 let g = mapped_grid(x->x + x.*(1 .- x)/2, x->@SMatrix[1.5 .- x], 11)
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
212 @test min_spacing(g) ≈ 0.055
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
213 end
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
214
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
215 let g = mapped_grid(identity, x->@SMatrix[1 0; 0 1], 11,11)
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
216 @test min_spacing(g) ≈ 0.1
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
217 end
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
218
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
219 let g = mapped_grid(identity, x->@SMatrix[1 0; 0 1], 11,21)
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
220 @test min_spacing(g) ≈ 0.05
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
221 end
1689
e11b5b6940a2 Add more tests for `normal`
Jonatan Werpers <jonatan@werpers.com>
parents: 1688
diff changeset
222
e11b5b6940a2 Add more tests for `normal`
Jonatan Werpers <jonatan@werpers.com>
parents: 1688
diff changeset
223
1690
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
224 @testset let a = @SVector[1,0], b = @SVector[1,1]/√2
1691
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
225 g = mapped_grid(_skew_mapping(a,b)...,11,11)
1690
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
226
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
227 @test min_spacing(g) ≈ 0.1*norm(b-a)
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
228 end
1689
e11b5b6940a2 Add more tests for `normal`
Jonatan Werpers <jonatan@werpers.com>
parents: 1688
diff changeset
229
1690
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
230 @testset let a = @SVector[1,0], b = @SVector[-1,1]/√2
1691
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
231 g = mapped_grid(_skew_mapping(a,b)...,11,11)
1690
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
232
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
233 @test min_spacing(g) ≈ 0.1*norm(a+b)
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
234 end
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
235 end
1689
e11b5b6940a2 Add more tests for `normal`
Jonatan Werpers <jonatan@werpers.com>
parents: 1688
diff changeset
236
1690
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
237 @testset "normal" begin
1691
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
238 g = mapped_grid(_partially_curved_mapping()...,10, 11)
1689
e11b5b6940a2 Add more tests for `normal`
Jonatan Werpers <jonatan@werpers.com>
parents: 1688
diff changeset
239
1690
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
240 @test normal(g, CartesianBoundary{1,Lower}()) == fill(@SVector[-1,0], 11)
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
241 @test normal(g, CartesianBoundary{1,Upper}()) == fill(@SVector[1,0], 11)
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
242 @test normal(g, CartesianBoundary{2,Lower}()) == fill(@SVector[0,-1], 10)
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
243 @test normal(g, CartesianBoundary{2,Upper}()) ≈ map(boundary_grid(g,CartesianBoundary{2,Upper}())|>logicalgrid) do ξ̄
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
244 α = 1-2ξ̄[1]
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
245 @SVector[α,1]/√(α^2 + 1)
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
246 end
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
247
1691
5bf4a35a78c5 Factor out functions for the mappings used in tests of MappedGrid
Jonatan Werpers <jonatan@werpers.com>
parents: 1690
diff changeset
248 g = mapped_grid(_fully_curved_mapping()...,5,4)
1690
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
249
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
250 unit(v) = v/norm(v)
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
251 @testset let bId = CartesianBoundary{1,Lower}()
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
252 lbg = boundary_grid(logicalgrid(g), bId)
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
253 @test normal(g, bId) ≈ map(lbg) do (ξ, η)
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
254 -unit(@SVector[1/2, η/3-1/6])
1689
e11b5b6940a2 Add more tests for `normal`
Jonatan Werpers <jonatan@werpers.com>
parents: 1688
diff changeset
255 end
1690
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
256 end
1689
e11b5b6940a2 Add more tests for `normal`
Jonatan Werpers <jonatan@werpers.com>
parents: 1688
diff changeset
257
1690
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
258 @testset let bId = CartesianBoundary{1,Upper}()
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
259 lbg = boundary_grid(logicalgrid(g), bId)
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
260 @test normal(g, bId) ≈ map(lbg) do (ξ, η)
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
261 unit(@SVector[7/2, 2η-1]/(5 + 3η + 2η^2))
1689
e11b5b6940a2 Add more tests for `normal`
Jonatan Werpers <jonatan@werpers.com>
parents: 1688
diff changeset
262 end
1690
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
263 end
1689
e11b5b6940a2 Add more tests for `normal`
Jonatan Werpers <jonatan@werpers.com>
parents: 1688
diff changeset
264
1690
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
265 @testset let bId = CartesianBoundary{2,Lower}()
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
266 lbg = boundary_grid(logicalgrid(g), bId)
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
267 @test normal(g, bId) ≈ map(lbg) do (ξ, η)
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
268 -unit(@SVector[-2ξ, 2]/(6 + ξ^2 - 2ξ))
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
269 end
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
270 end
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
271
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
272 @testset let bId = CartesianBoundary{2,Upper}()
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
273 lbg = boundary_grid(logicalgrid(g), bId)
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
274 @test normal(g, bId) ≈ map(lbg) do (ξ, η)
5eabe1f560f0 Reorganize nesting of tests for mapped_grid
Jonatan Werpers <jonatan@werpers.com>
parents: 1689
diff changeset
275 unit(@SVector[-3ξ, 2]/(6 + ξ^2 + 3ξ))
1689
e11b5b6940a2 Add more tests for `normal`
Jonatan Werpers <jonatan@werpers.com>
parents: 1688
diff changeset
276 end
1659
3bbcd496e021 Add function for computing the normal at the boundary of a mapped grid as a grid function
Jonatan Werpers <jonatan@werpers.com>
parents: 1570
diff changeset
277 end
1502
a2dc80396808 Add functions `curvilinear_grid`
Jonatan Werpers <jonatan@werpers.com>
parents: 1494
diff changeset
278 end