Mercurial > repos > public > sbplib_julia

annotate test/SbpOperators/volumeops/laplace/laplace_test.jl @ 750:f88b2117dc69 feature/laplace_opset

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