Mercurial > repos > public > sbplib_julia

annotate test/SbpOperators/boundaryops/boundary_operator_test.jl @ 1153:f1bb1b6d85dd refactor/sbpoperators/inflation

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Review: Suggest changes to test and removal of conveninece constructor
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Tue, 25 Oct 2022 10:33:27 +0200
parents 6b24dc2d7b11
children ae006e844870
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
1153
f1bb1b6d85dd Review: Suggest changes to test and removal of conveninece constructor
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1108
diff changeset
17 # Review: Remove the first line below in case we remove convenience constructor in BoundaryOperator
1108
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
18 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
19 @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
20 @test op_l isa LazyTensor{T,0,1} where T
1153
f1bb1b6d85dd Review: Suggest changes to test and removal of conveninece constructor
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 1108
diff changeset
21 # Review: Remove the first line below in case we remove convenience constructor in BoundaryOperator
1108
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
22 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
23 @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
24 @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
25 end
1099
05a25a5063bb Try to remove volume_operator and boundary_operator methods
Jonatan Werpers <jonatan@werpers.com>
parents: 995
diff changeset
26
05a25a5063bb Try to remove volume_operator and boundary_operator methods
Jonatan Werpers <jonatan@werpers.com>
parents: 995
diff changeset
27 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
28 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
29
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
30 @testset "Sizes" begin
1108
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
31 @test domain_size(op_l) == (11,)
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
32 @test domain_size(op_r) == (11,)
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
33
1108
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
34 @test range_size(op_l) == ()
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
35 @test range_size(op_r) == ()
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
36 end
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
37
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
38 @testset "Application" begin
1108
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
39 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
40 u = fill(3.124)
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
41 @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
42 @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
43 @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
44 @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
45 @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
46
1108
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
47 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
48 @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
49
1108
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
50 u = fill(1. +im)
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
51 @test (op_l'*u)[1] isa ComplexF64
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
52 @test (op_l'*u)[5] isa ComplexF64
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
53 @test (op_l'*u)[11] isa ComplexF64
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
54
1108
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
55 u = fill(3.124)
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
56 @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
57 @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
58 @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
59 @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
60 @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
61
1108
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
62 @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
63 @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
64 @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
65 @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
66 @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
67 end
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
68
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
69 @testset "Inferred" begin
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
70 v = ones(Float64, 11)
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
71 u = fill(1.)
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
72
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
73 @inferred apply(op_l, v)
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
74 @inferred apply(op_r, v)
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
75
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
76 @inferred apply_transpose(op_l, u, 4)
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
77 @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
78 @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
79 @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
80 @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
81 @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
82
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
83 @inferred apply_transpose(op_r, u, 4)
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
84 @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
85 @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
86 @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
87 @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
88 @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
89 end
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
90 end