Mercurial > repos > public > sbplib_julia

annotate test/SbpOperators/boundaryops/boundary_operator_test.jl @ 1108:6b24dc2d7b11 refactor/sbpoperators/inflation

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Clean up testsets in boundary_operator_test
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 09 Jun 2022 07:30:20 +0200
parents f80e69b0566b
children f1bb1b6d85dd
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.LazyTensors
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
4 using Sbplib.SbpOperators
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
5 using Sbplib.Grids
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
6 using Sbplib.RegionIndices
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
7 import Sbplib.SbpOperators.Stencil
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
8 import Sbplib.SbpOperators.BoundaryOperator
1099
05a25a5063bb Try to remove volume_operator and boundary_operator methods
Jonatan Werpers <jonatan@werpers.com>
parents: 995
diff changeset
9
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
10
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
11 @testset "BoundaryOperator" begin
1062
a76830879c63 Fix bugs from stencil refactor
Jonatan Werpers <jonatan@werpers.com>
parents: 1017
diff changeset
12 closure_stencil = Stencil(2.,1.,3.; center = 1)
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
13 g_1D = EquidistantGrid(11, 0.0, 1.0)
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
14 g_2D = EquidistantGrid((11,15), (0.0, 0.0), (1.0,1.0))
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
15
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
16 @testset "Constructors" begin
1108
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
17 op_l = BoundaryOperator{Lower}(closure_stencil,size(g_1D)[1])
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
18 @test op_l == BoundaryOperator(g_1D,closure_stencil,Lower())
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
19 @test op_l isa LazyTensor{T,0,1} where T
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
20
1108
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
21 op_r = BoundaryOperator{Upper}(closure_stencil,size(g_1D)[1])
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
22 @test op_r == BoundaryOperator(g_1D,closure_stencil,Upper())
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
23 @test op_r isa LazyTensor{T,0,1} where T
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
24 end
1099
05a25a5063bb Try to remove volume_operator and boundary_operator methods
Jonatan Werpers <jonatan@werpers.com>
parents: 995
diff changeset
25
05a25a5063bb Try to remove volume_operator and boundary_operator methods
Jonatan Werpers <jonatan@werpers.com>
parents: 995
diff changeset
26 op_l = BoundaryOperator(g_1D, closure_stencil, Lower())
05a25a5063bb Try to remove volume_operator and boundary_operator methods
Jonatan Werpers <jonatan@werpers.com>
parents: 995
diff changeset
27 op_r = BoundaryOperator(g_1D, closure_stencil, Upper())
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
28
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
29 @testset "Sizes" begin
1108
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
30 @test domain_size(op_l) == (11,)
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
31 @test domain_size(op_r) == (11,)
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
32
1108
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
33 @test range_size(op_l) == ()
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
34 @test range_size(op_r) == ()
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
35 end
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
36
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
37 @testset "Application" begin
1108
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
38 v = evalOn(g_1D,x->1+x^2)
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
39 u = fill(3.124)
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
40 @test (op_l*v)[] == 2*v[1] + v[2] + 3*v[3]
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
41 @test (op_r*v)[] == 2*v[end] + v[end-1] + 3*v[end-2]
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
42 @test (op_r*v)[1] == 2*v[end] + v[end-1] + 3*v[end-2]
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
43 @test op_l'*u == [2*u[]; u[]; 3*u[]; zeros(8)]
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
44 @test op_r'*u == [zeros(8); 3*u[]; u[]; 2*u[]]
946
469ed954b493 Allow volume_operator, boundary_operator, and constant_interior_scaling_operator to act on arbitrary arrays
Jonatan Werpers <jonatan@werpers.com>
parents: 728
diff changeset
45
1108
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
46 v = evalOn(g_1D, x->1. +x*im)
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
47 @test (op_l*v)[] isa ComplexF64
946
469ed954b493 Allow volume_operator, boundary_operator, and constant_interior_scaling_operator to act on arbitrary arrays
Jonatan Werpers <jonatan@werpers.com>
parents: 728
diff changeset
48
1108
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
49 u = fill(1. +im)
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
50 @test (op_l'*u)[1] isa ComplexF64
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
51 @test (op_l'*u)[5] isa ComplexF64
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
52 @test (op_l'*u)[11] isa ComplexF64
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
53
1108
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
54 u = fill(3.124)
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
55 @test (op_l'*u)[Index(1,Lower)] == 2*u[]
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
56 @test (op_l'*u)[Index(2,Lower)] == u[]
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
57 @test (op_l'*u)[Index(6,Interior)] == 0
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
58 @test (op_l'*u)[Index(10,Upper)] == 0
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
59 @test (op_l'*u)[Index(11,Upper)] == 0
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
60
1108
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
61 @test (op_r'*u)[Index(1,Lower)] == 0
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
62 @test (op_r'*u)[Index(2,Lower)] == 0
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
63 @test (op_r'*u)[Index(6,Interior)] == 0
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
64 @test (op_r'*u)[Index(10,Upper)] == u[]
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
65 @test (op_r'*u)[Index(11,Upper)] == 2*u[]
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
66 end
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
67
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
68 @testset "Inferred" begin
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
69 v = ones(Float64, 11)
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
70 u = fill(1.)
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
71
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
72 @inferred apply(op_l, v)
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
73 @inferred apply(op_r, v)
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
74
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
75 @inferred apply_transpose(op_l, u, 4)
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
76 @inferred apply_transpose(op_l, u, Index(1,Lower))
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
77 @inferred apply_transpose(op_l, u, Index(2,Lower))
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
78 @inferred apply_transpose(op_l, u, Index(6,Interior))
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
79 @inferred apply_transpose(op_l, u, Index(10,Upper))
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
80 @inferred apply_transpose(op_l, u, Index(11,Upper))
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
81
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
82 @inferred apply_transpose(op_r, u, 4)
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
83 @inferred apply_transpose(op_r, u, Index(1,Lower))
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
84 @inferred apply_transpose(op_r, u, Index(2,Lower))
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
85 @inferred apply_transpose(op_r, u, Index(6,Interior))
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
86 @inferred apply_transpose(op_r, u, Index(10,Upper))
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
87 @inferred apply_transpose(op_r, u, Index(11,Upper))
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
88 end
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
89 end