Mercurial > repos > public > sbplib_julia

annotate test/SbpOperators/volumeops/laplace/laplace_test.jl @ 752:ff0ef711c388 feature/laplace_opset

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Remove obsolete file
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Fri, 19 Mar 2021 17:12:43 +0100
parents f94feb005e7d
children fc83d672be36
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
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
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
9 @testset "Laplace" begin
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
10 g_1D = EquidistantGrid(101, 0.0, 1.)
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
11 g_3D = EquidistantGrid((51,101,52), (0.0, -1.0, 0.0), (1., 1., 1.))
750
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
12 op = read_D2_operator(sbp_operators_path()*"standard_diagonal.toml"; order=4)
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
13 @testset "Constructors" begin
750
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
14
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
15 @testset "1D" begin
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
16 # Create all tensor mappings included in Laplace
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
17 Δ = laplace(g_1D, op.innerStencil, op.closureStencils)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
18 H = inner_product(g_1D, op.quadratureClosure)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
19 Hi = inverse_inner_product(g_1D, op.quadratureClosure)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
20
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
21 (id_l, id_r) = boundary_identifiers(g_1D)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
22
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
23 e_l = boundary_restriction(g_1D,op.eClosure,id_l)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
24 e_r = boundary_restriction(g_1D,op.eClosure,id_r)
751
f94feb005e7d Change from Dict to StaticDict in Laplace
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 750
diff changeset
25 e_dict = StaticDict(Pair(id_l,e_l),Pair(id_r,e_r))
750
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
26
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
27 d_l = normal_derivative(g_1D,op.dClosure,id_l)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
28 d_r = normal_derivative(g_1D,op.dClosure,id_r)
751
f94feb005e7d Change from Dict to StaticDict in Laplace
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 750
diff changeset
29 d_dict = StaticDict(Pair(id_l,d_l),Pair(id_r,d_r))
750
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
30
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
31 H_l = inner_product(boundary_grid(g_1D,id_l),op.quadratureClosure)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
32 H_r = inner_product(boundary_grid(g_1D,id_r),op.quadratureClosure)
751
f94feb005e7d Change from Dict to StaticDict in Laplace
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 750
diff changeset
33 Hb_dict = StaticDict(Pair(id_l,H_l),Pair(id_r,H_r))
750
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
34
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
35 L = Laplace(g_1D, sbp_operators_path()*"standard_diagonal.toml"; order=4)
751
f94feb005e7d Change from Dict to StaticDict in Laplace
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 750
diff changeset
36 @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
37 @test L isa TensorMapping{T,1,1} where T
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
38 @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
39 end
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
40 @testset "3D" begin
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
41 # Create all tensor mappings included in Laplace
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
42 Δ = laplace(g_3D, op.innerStencil, op.closureStencils)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
43 H = inner_product(g_3D, op.quadratureClosure)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
44 Hi = inverse_inner_product(g_3D, op.quadratureClosure)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
45
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
46 (id_l, id_r, id_s, id_n, id_b, id_t) = boundary_identifiers(g_3D)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
47
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
48 e_l = boundary_restriction(g_3D,op.eClosure,id_l)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
49 e_r = boundary_restriction(g_3D,op.eClosure,id_r)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
50 e_s = boundary_restriction(g_3D,op.eClosure,id_s)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
51 e_n = boundary_restriction(g_3D,op.eClosure,id_n)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
52 e_b = boundary_restriction(g_3D,op.eClosure,id_b)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
53 e_t = boundary_restriction(g_3D,op.eClosure,id_t)
751
f94feb005e7d Change from Dict to StaticDict in Laplace
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 750
diff changeset
54 e_dict = StaticDict(Pair(id_l,e_l),Pair(id_r,e_r),
750
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
55 Pair(id_s,e_s),Pair(id_n,e_n),
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
56 Pair(id_b,e_b),Pair(id_t,e_t))
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
57
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
58 d_l = normal_derivative(g_3D,op.dClosure,id_l)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
59 d_r = normal_derivative(g_3D,op.dClosure,id_r)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
60 d_s = normal_derivative(g_3D,op.dClosure,id_s)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
61 d_n = normal_derivative(g_3D,op.dClosure,id_n)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
62 d_b = normal_derivative(g_3D,op.dClosure,id_b)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
63 d_t = normal_derivative(g_3D,op.dClosure,id_t)
751
f94feb005e7d Change from Dict to StaticDict in Laplace
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 750
diff changeset
64 d_dict = StaticDict(Pair(id_l,d_l),Pair(id_r,d_r),
750
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
65 Pair(id_s,d_s),Pair(id_n,d_n),
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
66 Pair(id_b,d_b),Pair(id_t,d_t))
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
67
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
68 H_l = inner_product(boundary_grid(g_3D,id_l),op.quadratureClosure)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
69 H_r = inner_product(boundary_grid(g_3D,id_r),op.quadratureClosure)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
70 H_s = inner_product(boundary_grid(g_3D,id_s),op.quadratureClosure)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
71 H_n = inner_product(boundary_grid(g_3D,id_n),op.quadratureClosure)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
72 H_b = inner_product(boundary_grid(g_3D,id_b),op.quadratureClosure)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
73 H_t = inner_product(boundary_grid(g_3D,id_t),op.quadratureClosure)
751
f94feb005e7d Change from Dict to StaticDict in Laplace
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 750
diff changeset
74 Hb_dict = StaticDict(Pair(id_l,H_l),Pair(id_r,H_r),
750
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
75 Pair(id_s,H_s),Pair(id_n,H_n),
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
76 Pair(id_b,H_b),Pair(id_t,H_t))
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
77
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
78 L = Laplace(g_3D, sbp_operators_path()*"standard_diagonal.toml"; order=4)
751
f94feb005e7d Change from Dict to StaticDict in Laplace
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 750
diff changeset
79 @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
80 @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
81 @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
82 end
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
83 end
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
84
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
85 @testset "laplace" begin
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
86 @testset "1D" begin
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
87 L = laplace(g_1D, op.innerStencil, op.closureStencils)
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
88 @test L == second_derivative(g_1D, op.innerStencil, op.closureStencils)
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
89 @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
90 end
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
91 @testset "3D" begin
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
92 L = laplace(g_3D, op.innerStencil, op.closureStencils)
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
93 @test L isa TensorMapping{T,3,3} where T
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
94 Dxx = second_derivative(g_3D, op.innerStencil, op.closureStencils,1)
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
95 Dyy = second_derivative(g_3D, op.innerStencil, op.closureStencils,2)
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
96 Dzz = second_derivative(g_3D, op.innerStencil, op.closureStencils,3)
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
97 @test L == Dxx + Dyy + Dzz
750
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
98 @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
99 end
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
100 end
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
101
750
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
102 @testset "inner_product" begin
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
103 L = Laplace(g_3D, sbp_operators_path()*"standard_diagonal.toml"; order=4)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
104 @test inner_product(L) == inner_product(g_3D,op.quadratureClosure)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
105 end
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
106
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
107 @testset "inverse_inner_product" begin
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
108 L = Laplace(g_3D, sbp_operators_path()*"standard_diagonal.toml"; order=4)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
109 @test inverse_inner_product(L) == inverse_inner_product(g_3D,op.quadratureClosure)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
110 end
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
111
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
112 @testset "boundary_restriction" begin
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
113 L = Laplace(g_3D, sbp_operators_path()*"standard_diagonal.toml"; order=4)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
114 id_l = CartesianBoundary{1,Lower}()
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
115 id_r = CartesianBoundary{1,Upper}()
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
116 id_s = CartesianBoundary{2,Lower}()
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
117 id_n = CartesianBoundary{2,Upper}()
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
118 id_b = CartesianBoundary{3,Lower}()
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
119 id_t = CartesianBoundary{3,Upper}()
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
120 @test boundary_restriction(L,id_l) == boundary_restriction(g_3D,op.eClosure,id_l)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
121 @test boundary_restriction(L,id_r) == boundary_restriction(g_3D,op.eClosure,id_r)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
122 @test boundary_restriction(L,id_s) == boundary_restriction(g_3D,op.eClosure,id_s)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
123 @test boundary_restriction(L,id_n) == boundary_restriction(g_3D,op.eClosure,id_n)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
124 @test boundary_restriction(L,id_b) == boundary_restriction(g_3D,op.eClosure,id_b)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
125 @test boundary_restriction(L,id_t) == boundary_restriction(g_3D,op.eClosure,id_t)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
126 end
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
127
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
128 @testset "normal_derivative" begin
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
129 L = Laplace(g_3D, sbp_operators_path()*"standard_diagonal.toml"; order=4)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
130 id_l = CartesianBoundary{1,Lower}()
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
131 id_r = CartesianBoundary{1,Upper}()
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
132 id_s = CartesianBoundary{2,Lower}()
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
133 id_n = CartesianBoundary{2,Upper}()
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
134 id_b = CartesianBoundary{3,Lower}()
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
135 id_t = CartesianBoundary{3,Upper}()
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
136 @test normal_derivative(L,id_l) == normal_derivative(g_3D,op.dClosure,id_l)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
137 @test normal_derivative(L,id_r) == normal_derivative(g_3D,op.dClosure,id_r)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
138 @test normal_derivative(L,id_s) == normal_derivative(g_3D,op.dClosure,id_s)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
139 @test normal_derivative(L,id_n) == normal_derivative(g_3D,op.dClosure,id_n)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
140 @test normal_derivative(L,id_b) == normal_derivative(g_3D,op.dClosure,id_b)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
141 @test normal_derivative(L,id_t) == normal_derivative(g_3D,op.dClosure,id_t)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
142 end
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
143
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
144 @testset "boundary_quadrature" begin
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
145 L = Laplace(g_3D, sbp_operators_path()*"standard_diagonal.toml"; order=4)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
146 id_l = CartesianBoundary{1,Lower}()
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
147 id_r = CartesianBoundary{1,Upper}()
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
148 id_s = CartesianBoundary{2,Lower}()
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
149 id_n = CartesianBoundary{2,Upper}()
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
150 id_b = CartesianBoundary{3,Lower}()
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
151 id_t = CartesianBoundary{3,Upper}()
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
152 @test boundary_quadrature(L,id_l) == inner_product(boundary_grid(g_3D,id_l),op.quadratureClosure)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
153 @test boundary_quadrature(L,id_r) == inner_product(boundary_grid(g_3D,id_r),op.quadratureClosure)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
154 @test boundary_quadrature(L,id_s) == inner_product(boundary_grid(g_3D,id_s),op.quadratureClosure)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
155 @test boundary_quadrature(L,id_n) == inner_product(boundary_grid(g_3D,id_n),op.quadratureClosure)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
156 @test boundary_quadrature(L,id_b) == inner_product(boundary_grid(g_3D,id_b),op.quadratureClosure)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
157 @test boundary_quadrature(L,id_t) == inner_product(boundary_grid(g_3D,id_t),op.quadratureClosure)
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
158 end
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
159
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
160 # 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
161 # the error is measured in the l2-norm.
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
162 @testset "Accuracy" begin
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
163 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
164 polynomials = ()
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
165 maxOrder = 4;
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
166 for i = 0:maxOrder-1
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
167 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
168 polynomials = (polynomials...,evalOn(g_3D,f_i))
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
169 end
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
170 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
171 Δ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
172
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
173 # 2nd order interior stencil, 1st order boundary stencil,
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
174 # 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
175 @testset "2nd order" begin
750
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
176 L = Laplace(g_3D, sbp_operators_path()*"standard_diagonal.toml"; order=2)
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
177 @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
178 @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
179 @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
180 @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
181 end
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
182
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
183 # 4th order interior stencil, 2nd order boundary stencil,
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
184 # 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
185 @testset "4th order" begin
750
f88b2117dc69 Merge in default
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 732
diff changeset
186 L = Laplace(g_3D, sbp_operators_path()*"standard_diagonal.toml"; order=4)
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
187 # NOTE: high tolerances for checking the "exact" differentiation
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
188 # 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
189 @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
190 @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
191 @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
192 @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
193 @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
194 end
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
195 end
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
196 end