Mercurial > repos > public > sbplib_julia
annotate test/SbpOperators/volumeops/laplace/laplace_test.jl @ 872:6a4d36eccf39 feature/laplace_opset
REVIEW: Add some comments about the tests
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Tue, 25 Jan 2022 10:36:13 +0100 |
| parents | 1784b1c0af3e |
| children | 12e8e431b43c |
| rev | line source |
|---|---|
|
728
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
1 using Test |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
2 |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
3 using Sbplib.SbpOperators |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
4 using Sbplib.Grids |
|
732
6114274447f5
Add missing using and imports in test files
Jonatan Werpers <jonatan@werpers.com>
parents:
728
diff
changeset
|
5 using Sbplib.LazyTensors |
|
750
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
6 using Sbplib.RegionIndices |
|
751
f94feb005e7d
Change from Dict to StaticDict in Laplace
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
750
diff
changeset
|
7 using Sbplib.StaticDicts |
|
728
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
8 |
| 866 | 9 operator_path = sbp_operators_path()*"standard_diagonal.toml" |
| 10 # Default stencils (4th order) | |
| 11 stencil_set = read_stencil_set(operator_path; order=4) | |
| 12 inner_stencil = parse_stencil(stencil_set["D2"]["inner_stencil"]) | |
| 13 closure_stencils = parse_stencil.(stencil_set["D2"]["closure_stencils"]) | |
| 14 e_closure = parse_stencil(stencil_set["e"]["closure"]) | |
| 15 d_closure = parse_stencil(stencil_set["d1"]["closure"]) | |
| 16 quadrature_interior = parse_scalar(stencil_set["H"]["inner"]) | |
| 17 quadrature_closure = parse_tuple(stencil_set["H"]["closure"]) | |
|
728
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
18 |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
19 @testset "Laplace" begin |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
20 g_1D = EquidistantGrid(101, 0.0, 1.) |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
21 g_3D = EquidistantGrid((51,101,52), (0.0, -1.0, 0.0), (1., 1., 1.)) |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
22 @testset "Constructors" begin |
|
750
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
23 |
|
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
24 @testset "1D" begin |
| 866 | 25 Δ = laplace(g_1D, inner_stencil, closure_stencils) |
| 26 H = inner_product(g_1D, quadrature_interior, quadrature_closure) | |
| 27 Hi = inverse_inner_product(g_1D, quadrature_interior, quadrature_closure) | |
|
750
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
28 |
|
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
29 (id_l, id_r) = boundary_identifiers(g_1D) |
|
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
30 |
| 866 | 31 e_l = boundary_restriction(g_1D, e_closure,id_l) |
| 32 e_r = boundary_restriction(g_1D, e_closure,id_r) | |
|
753
fc83d672be36
Minor cleanup of code
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
751
diff
changeset
|
33 e_dict = StaticDict(id_l => e_l, id_r => e_r) |
|
750
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
34 |
| 866 | 35 d_l = normal_derivative(g_1D, d_closure,id_l) |
| 36 d_r = normal_derivative(g_1D, d_closure,id_r) | |
|
753
fc83d672be36
Minor cleanup of code
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
751
diff
changeset
|
37 d_dict = StaticDict(id_l => d_l, id_r => d_r) |
|
750
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
38 |
| 866 | 39 H_l = inner_product(boundary_grid(g_1D,id_l), quadrature_interior, quadrature_closure) |
| 40 H_r = inner_product(boundary_grid(g_1D,id_r), quadrature_interior, quadrature_closure) | |
|
753
fc83d672be36
Minor cleanup of code
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
751
diff
changeset
|
41 Hb_dict = StaticDict(id_l => H_l, id_r => H_r) |
|
750
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
42 |
| 866 | 43 L = Laplace(g_1D, operator_path; order=4) |
| 44 @test L == Laplace(Δ, H, Hi, e_dict, d_dict, Hb_dict) | |
|
750
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
45 @test L isa TensorMapping{T,1,1} where T |
| 866 | 46 @inferred Laplace(Δ, H, Hi, e_dict, d_dict, Hb_dict) |
|
872
6a4d36eccf39
REVIEW: Add some comments about the tests
Jonatan Werpers <jonatan@werpers.com>
parents:
866
diff
changeset
|
47 # REVIEW: The tests above seem very tied to the implementation. Is |
|
6a4d36eccf39
REVIEW: Add some comments about the tests
Jonatan Werpers <jonatan@werpers.com>
parents:
866
diff
changeset
|
48 # it important that the components of the operator set are stored |
|
6a4d36eccf39
REVIEW: Add some comments about the tests
Jonatan Werpers <jonatan@werpers.com>
parents:
866
diff
changeset
|
49 # in static dicts? Is something like below better? |
|
6a4d36eccf39
REVIEW: Add some comments about the tests
Jonatan Werpers <jonatan@werpers.com>
parents:
866
diff
changeset
|
50 # |
|
6a4d36eccf39
REVIEW: Add some comments about the tests
Jonatan Werpers <jonatan@werpers.com>
parents:
866
diff
changeset
|
51 # ``` |
|
6a4d36eccf39
REVIEW: Add some comments about the tests
Jonatan Werpers <jonatan@werpers.com>
parents:
866
diff
changeset
|
52 # L = Laplace(g_1D, operator_path; order=4) |
|
6a4d36eccf39
REVIEW: Add some comments about the tests
Jonatan Werpers <jonatan@werpers.com>
parents:
866
diff
changeset
|
53 # @test L isa TensorMapping{T,1,1} where T |
|
6a4d36eccf39
REVIEW: Add some comments about the tests
Jonatan Werpers <jonatan@werpers.com>
parents:
866
diff
changeset
|
54 # @test boundary_restriction(L,id_l) == boundary_restriction(g_1D, e_closure,id_l) |
|
6a4d36eccf39
REVIEW: Add some comments about the tests
Jonatan Werpers <jonatan@werpers.com>
parents:
866
diff
changeset
|
55 # ... |
|
6a4d36eccf39
REVIEW: Add some comments about the tests
Jonatan Werpers <jonatan@werpers.com>
parents:
866
diff
changeset
|
56 # ``` |
|
6a4d36eccf39
REVIEW: Add some comments about the tests
Jonatan Werpers <jonatan@werpers.com>
parents:
866
diff
changeset
|
57 # I guess this is more or less simply a reorganization of the test and skipping testing for the struct layout |
|
750
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
58 end |
|
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
59 @testset "3D" begin |
| 866 | 60 Δ = laplace(g_3D, inner_stencil, closure_stencils) |
| 61 H = inner_product(g_3D, quadrature_interior, quadrature_closure) | |
| 62 Hi = inverse_inner_product(g_3D, quadrature_interior, quadrature_closure) | |
|
750
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
63 |
|
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
64 (id_l, id_r, id_s, id_n, id_b, id_t) = boundary_identifiers(g_3D) |
| 866 | 65 e_l = boundary_restriction(g_3D, e_closure,id_l) |
| 66 e_r = boundary_restriction(g_3D, e_closure,id_r) | |
| 67 e_s = boundary_restriction(g_3D, e_closure,id_s) | |
| 68 e_n = boundary_restriction(g_3D, e_closure,id_n) | |
| 69 e_b = boundary_restriction(g_3D, e_closure,id_b) | |
| 70 e_t = boundary_restriction(g_3D, e_closure,id_t) | |
|
753
fc83d672be36
Minor cleanup of code
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
751
diff
changeset
|
71 e_dict = StaticDict(id_l => e_l, id_r => e_r, |
|
fc83d672be36
Minor cleanup of code
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
751
diff
changeset
|
72 id_s => e_s, id_n => e_n, |
|
fc83d672be36
Minor cleanup of code
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
751
diff
changeset
|
73 id_b => e_b, id_t => e_t) |
|
750
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
74 |
| 866 | 75 d_l = normal_derivative(g_3D, d_closure,id_l) |
| 76 d_r = normal_derivative(g_3D, d_closure,id_r) | |
| 77 d_s = normal_derivative(g_3D, d_closure,id_s) | |
| 78 d_n = normal_derivative(g_3D, d_closure,id_n) | |
| 79 d_b = normal_derivative(g_3D, d_closure,id_b) | |
| 80 d_t = normal_derivative(g_3D, d_closure,id_t) | |
|
753
fc83d672be36
Minor cleanup of code
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
751
diff
changeset
|
81 d_dict = StaticDict(id_l => d_l, id_r => d_r, |
|
fc83d672be36
Minor cleanup of code
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
751
diff
changeset
|
82 id_s => d_s, id_n => d_n, |
|
fc83d672be36
Minor cleanup of code
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
751
diff
changeset
|
83 id_b => d_b, id_t => d_t) |
|
750
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
84 |
| 866 | 85 H_l = inner_product(boundary_grid(g_3D,id_l), quadrature_interior, quadrature_closure) |
| 86 H_r = inner_product(boundary_grid(g_3D,id_r), quadrature_interior, quadrature_closure) | |
| 87 H_s = inner_product(boundary_grid(g_3D,id_s), quadrature_interior, quadrature_closure) | |
| 88 H_n = inner_product(boundary_grid(g_3D,id_n), quadrature_interior, quadrature_closure) | |
| 89 H_b = inner_product(boundary_grid(g_3D,id_b), quadrature_interior, quadrature_closure) | |
| 90 H_t = inner_product(boundary_grid(g_3D,id_t), quadrature_interior, quadrature_closure) | |
|
753
fc83d672be36
Minor cleanup of code
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
751
diff
changeset
|
91 Hb_dict = StaticDict(id_l => H_l, id_r => H_r, |
|
fc83d672be36
Minor cleanup of code
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
751
diff
changeset
|
92 id_s => H_s, id_n => H_n, |
|
fc83d672be36
Minor cleanup of code
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
751
diff
changeset
|
93 id_b => H_b, id_t => H_t) |
|
750
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
94 |
| 866 | 95 L = Laplace(g_3D, operator_path; order=4) |
|
751
f94feb005e7d
Change from Dict to StaticDict in Laplace
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
750
diff
changeset
|
96 @test L == Laplace(Δ,H,Hi,e_dict,d_dict,Hb_dict) |
|
750
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
97 @test L isa TensorMapping{T,3,3} where T |
|
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
98 @inferred Laplace(Δ,H,Hi,e_dict,d_dict,Hb_dict) |
|
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
99 end |
|
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
100 end |
|
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
101 |
|
872
6a4d36eccf39
REVIEW: Add some comments about the tests
Jonatan Werpers <jonatan@werpers.com>
parents:
866
diff
changeset
|
102 # REVIEW: Is this testset misplaced? Should it really be inside the "Laplace" testset? |
|
750
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
103 @testset "laplace" begin |
|
728
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
104 @testset "1D" begin |
| 799 | 105 L = laplace(g_1D, inner_stencil, closure_stencils) |
| 106 @test L == second_derivative(g_1D, inner_stencil, closure_stencils) | |
|
728
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
107 @test L isa TensorMapping{T,1,1} where T |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
108 end |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
109 @testset "3D" begin |
| 799 | 110 L = laplace(g_3D, inner_stencil, closure_stencils) |
|
728
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
111 @test L isa TensorMapping{T,3,3} where T |
| 799 | 112 Dxx = second_derivative(g_3D, inner_stencil, closure_stencils, 1) |
| 113 Dyy = second_derivative(g_3D, inner_stencil, closure_stencils, 2) | |
| 114 Dzz = second_derivative(g_3D, inner_stencil, closure_stencils, 3) | |
|
728
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
115 @test L == Dxx + Dyy + Dzz |
|
750
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
116 @test L isa TensorMapping{T,3,3} where T |
|
728
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
117 end |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
118 end |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
119 |
|
750
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
120 @testset "inner_product" begin |
| 866 | 121 L = Laplace(g_3D, operator_path; order=4) |
| 122 @test inner_product(L) == inner_product(g_3D, quadrature_interior, quadrature_closure) | |
|
750
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
123 end |
|
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
124 |
|
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
125 @testset "inverse_inner_product" begin |
| 866 | 126 L = Laplace(g_3D, operator_path; order=4) |
| 127 @test inverse_inner_product(L) == inverse_inner_product(g_3D, quadrature_interior, quadrature_closure) | |
|
750
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
128 end |
|
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
129 |
|
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
130 @testset "boundary_restriction" begin |
| 866 | 131 L = Laplace(g_3D, operator_path; order=4) |
|
753
fc83d672be36
Minor cleanup of code
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
751
diff
changeset
|
132 (id_l, id_r, id_s, id_n, id_b, id_t) = boundary_identifiers(g_3D) |
| 866 | 133 @test boundary_restriction(L, id_l) == boundary_restriction(g_3D, e_closure,id_l) |
| 134 @test boundary_restriction(L, id_r) == boundary_restriction(g_3D, e_closure,id_r) | |
| 135 @test boundary_restriction(L, id_s) == boundary_restriction(g_3D, e_closure,id_s) | |
| 136 @test boundary_restriction(L, id_n) == boundary_restriction(g_3D, e_closure,id_n) | |
| 137 @test boundary_restriction(L, id_b) == boundary_restriction(g_3D, e_closure,id_b) | |
| 138 @test boundary_restriction(L, id_t) == boundary_restriction(g_3D, e_closure,id_t) | |
|
754
dc38e57ebd1b
Add convenience functions for returning multiple boundary operators from Laplace
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
753
diff
changeset
|
139 |
|
dc38e57ebd1b
Add convenience functions for returning multiple boundary operators from Laplace
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
753
diff
changeset
|
140 ids = boundary_identifiers(g_3D) |
| 866 | 141 es = boundary_restriction(L, ids) |
| 142 @test es == (boundary_restriction(L, id_l), | |
| 143 boundary_restriction(L, id_r), | |
| 144 boundary_restriction(L, id_s), | |
| 145 boundary_restriction(L, id_n), | |
| 146 boundary_restriction(L, id_b), | |
| 147 boundary_restriction(L, id_t)); | |
| 148 @test es == boundary_restriction(L, ids...) | |
|
750
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
149 end |
|
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
150 |
|
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
151 @testset "normal_derivative" begin |
| 866 | 152 L = Laplace(g_3D, operator_path; order=4) |
|
753
fc83d672be36
Minor cleanup of code
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
751
diff
changeset
|
153 (id_l, id_r, id_s, id_n, id_b, id_t) = boundary_identifiers(g_3D) |
| 866 | 154 @test normal_derivative(L, id_l) == normal_derivative(g_3D, d_closure,id_l) |
| 155 @test normal_derivative(L, id_r) == normal_derivative(g_3D, d_closure,id_r) | |
| 156 @test normal_derivative(L, id_s) == normal_derivative(g_3D, d_closure,id_s) | |
| 157 @test normal_derivative(L, id_n) == normal_derivative(g_3D, d_closure,id_n) | |
| 158 @test normal_derivative(L, id_b) == normal_derivative(g_3D, d_closure,id_b) | |
| 159 @test normal_derivative(L, id_t) == normal_derivative(g_3D, d_closure,id_t) | |
|
754
dc38e57ebd1b
Add convenience functions for returning multiple boundary operators from Laplace
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
753
diff
changeset
|
160 |
|
dc38e57ebd1b
Add convenience functions for returning multiple boundary operators from Laplace
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
753
diff
changeset
|
161 ids = boundary_identifiers(g_3D) |
| 866 | 162 ds = normal_derivative(L, ids) |
| 163 @test ds == (normal_derivative(L, id_l), | |
| 164 normal_derivative(L, id_r), | |
| 165 normal_derivative(L, id_s), | |
| 166 normal_derivative(L, id_n), | |
| 167 normal_derivative(L, id_b), | |
| 168 normal_derivative(L, id_t)); | |
| 169 @test ds == normal_derivative(L, ids...) | |
|
750
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
170 end |
|
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
171 |
|
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
172 @testset "boundary_quadrature" begin |
| 866 | 173 L = Laplace(g_3D, operator_path; order=4) |
|
753
fc83d672be36
Minor cleanup of code
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
751
diff
changeset
|
174 (id_l, id_r, id_s, id_n, id_b, id_t) = boundary_identifiers(g_3D) |
| 866 | 175 @test boundary_quadrature(L, id_l) == inner_product(boundary_grid(g_3D, id_l), quadrature_interior, quadrature_closure) |
| 176 @test boundary_quadrature(L, id_r) == inner_product(boundary_grid(g_3D, id_r), quadrature_interior, quadrature_closure) | |
| 177 @test boundary_quadrature(L, id_s) == inner_product(boundary_grid(g_3D, id_s), quadrature_interior, quadrature_closure) | |
| 178 @test boundary_quadrature(L, id_n) == inner_product(boundary_grid(g_3D, id_n), quadrature_interior, quadrature_closure) | |
| 179 @test boundary_quadrature(L, id_b) == inner_product(boundary_grid(g_3D, id_b), quadrature_interior, quadrature_closure) | |
| 180 @test boundary_quadrature(L, id_t) == inner_product(boundary_grid(g_3D, id_t), quadrature_interior, quadrature_closure) | |
|
754
dc38e57ebd1b
Add convenience functions for returning multiple boundary operators from Laplace
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
753
diff
changeset
|
181 |
|
dc38e57ebd1b
Add convenience functions for returning multiple boundary operators from Laplace
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
753
diff
changeset
|
182 ids = boundary_identifiers(g_3D) |
| 866 | 183 H_gammas = boundary_quadrature(L, ids) |
| 184 @test H_gammas == (boundary_quadrature(L, id_l), | |
| 185 boundary_quadrature(L, id_r), | |
| 186 boundary_quadrature(L, id_s), | |
| 187 boundary_quadrature(L, id_n), | |
| 188 boundary_quadrature(L, id_b), | |
| 189 boundary_quadrature(L, id_t)); | |
| 190 @test H_gammas == boundary_quadrature(L, ids...) | |
|
750
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
191 end |
|
f88b2117dc69
Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
732
diff
changeset
|
192 |
|
728
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
193 # Exact differentiation is measured point-wise. In other cases |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
194 # the error is measured in the l2-norm. |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
195 @testset "Accuracy" begin |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
196 l2(v) = sqrt(prod(spacing(g_3D))*sum(v.^2)); |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
197 polynomials = () |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
198 maxOrder = 4; |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
199 for i = 0:maxOrder-1 |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
200 f_i(x,y,z) = 1/factorial(i)*(y^i + x^i + z^i) |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
201 polynomials = (polynomials...,evalOn(g_3D,f_i)) |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
202 end |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
203 v = evalOn(g_3D, (x,y,z) -> sin(x) + cos(y) + exp(z)) |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
204 Δv = evalOn(g_3D,(x,y,z) -> -sin(x) - cos(y) + exp(z)) |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
205 |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
206 # 2nd order interior stencil, 1st order boundary stencil, |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
207 # implies that L*v should be exact for binomials up to order 2. |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
208 @testset "2nd order" begin |
| 866 | 209 stencil_set = read_stencil_set(operator_path; order=2) |
| 799 | 210 inner_stencil = parse_stencil(stencil_set["D2"]["inner_stencil"]) |
| 211 closure_stencils = parse_stencil.(stencil_set["D2"]["closure_stencils"]) | |
| 212 L = laplace(g_3D, inner_stencil, closure_stencils) | |
|
728
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
213 @test L*polynomials[1] ≈ zeros(Float64, size(g_3D)...) atol = 5e-9 |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
214 @test L*polynomials[2] ≈ zeros(Float64, size(g_3D)...) atol = 5e-9 |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
215 @test L*polynomials[3] ≈ polynomials[1] atol = 5e-9 |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
216 @test L*v ≈ Δv rtol = 5e-2 norm = l2 |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
217 end |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
218 |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
219 # 4th order interior stencil, 2nd order boundary stencil, |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
220 # implies that L*v should be exact for binomials up to order 3. |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
221 @testset "4th order" begin |
| 866 | 222 stencil_set = read_stencil_set(operator_path; order=4) |
| 799 | 223 inner_stencil = parse_stencil(stencil_set["D2"]["inner_stencil"]) |
| 224 closure_stencils = parse_stencil.(stencil_set["D2"]["closure_stencils"]) | |
| 225 L = laplace(g_3D, inner_stencil, closure_stencils) | |
|
728
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
226 # NOTE: high tolerances for checking the "exact" differentiation |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
227 # due to accumulation of round-off errors/cancellation errors? |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
228 @test L*polynomials[1] ≈ zeros(Float64, size(g_3D)...) atol = 5e-9 |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
229 @test L*polynomials[2] ≈ zeros(Float64, size(g_3D)...) atol = 5e-9 |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
230 @test L*polynomials[3] ≈ polynomials[1] atol = 5e-9 |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
231 @test L*polynomials[4] ≈ polynomials[2] atol = 5e-9 |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
232 @test L*v ≈ Δv rtol = 5e-4 norm = l2 |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
233 end |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
234 end |
|
45966c77cb20
Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
235 end |