Mercurial > repos > public > sbplib_julia

annotate test/SbpOperators/boundaryops/boundary_operator_test.jl @ 2005:52e5ab4a96d5

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge feature/grids/parameter_spaces/in
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 01 May 2025 14:57:09 +0200
parents 471a948cd2b2
children
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
1726
471a948cd2b2 Rename project from Sbplib to Diffinitive
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents: 1672
diff changeset
3 using Diffinitive.LazyTensors
471a948cd2b2 Rename project from Sbplib to Diffinitive
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents: 1672
diff changeset
4 using Diffinitive.SbpOperators
471a948cd2b2 Rename project from Sbplib to Diffinitive
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents: 1672
diff changeset
5 using Diffinitive.Grids
471a948cd2b2 Rename project from Sbplib to Diffinitive
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents: 1672
diff changeset
6 using Diffinitive.RegionIndices
471a948cd2b2 Rename project from Sbplib to Diffinitive
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents: 1672
diff changeset
7 import Diffinitive.SbpOperators.Stencil
471a948cd2b2 Rename project from Sbplib to Diffinitive
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents: 1672
diff changeset
8 import Diffinitive.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)
1281
1cc45207817e Fix tests for boundary operator and boundary restriction
Jonatan Werpers <jonatan@werpers.com>
parents: 1154
diff changeset
13 g_1D = EquidistantGrid(range(0,1,length=11))
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
14
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
15 @testset "Constructors" begin
1672
3714a391545a Make the boundary identifiers for EquidistantGrid subtype BoundaryIdentifer
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents: 1365
diff changeset
16 @test BoundaryOperator(g_1D, closure_stencil, LowerBoundary()) isa LazyTensor{T,0,1} where T
3714a391545a Make the boundary identifiers for EquidistantGrid subtype BoundaryIdentifer
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents: 1365
diff changeset
17 @test BoundaryOperator(g_1D, closure_stencil, UpperBoundary()) isa LazyTensor{T,0,1} where T
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
18 end
1099
05a25a5063bb Try to remove volume_operator and boundary_operator methods
Jonatan Werpers <jonatan@werpers.com>
parents: 995
diff changeset
19
1672
3714a391545a Make the boundary identifiers for EquidistantGrid subtype BoundaryIdentifer
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents: 1365
diff changeset
20 op_l = BoundaryOperator(g_1D, closure_stencil, LowerBoundary())
3714a391545a Make the boundary identifiers for EquidistantGrid subtype BoundaryIdentifer
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents: 1365
diff changeset
21 op_r = BoundaryOperator(g_1D, closure_stencil, UpperBoundary())
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
22
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
23 @testset "Sizes" begin
1108
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
24 @test domain_size(op_l) == (11,)
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
25 @test domain_size(op_r) == (11,)
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
26
1108
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
27 @test range_size(op_l) == ()
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
28 @test range_size(op_r) == ()
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
29 end
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
30
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
31 @testset "Application" begin
1281
1cc45207817e Fix tests for boundary operator and boundary restriction
Jonatan Werpers <jonatan@werpers.com>
parents: 1154
diff changeset
32 v = eval_on(g_1D,x->1+x^2)
1108
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
33 u = fill(3.124)
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
34 @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
35 @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
36 @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
37 @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
38 @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
39
1281
1cc45207817e Fix tests for boundary operator and boundary restriction
Jonatan Werpers <jonatan@werpers.com>
parents: 1154
diff changeset
40 v = eval_on(g_1D, x->1. +x*im)
1108
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
41 @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
42
1108
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
43 u = fill(1. +im)
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
44 @test (op_l'*u)[1] isa ComplexF64
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
45 @test (op_l'*u)[5] isa ComplexF64
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
46 @test (op_l'*u)[11] isa ComplexF64
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
47
1108
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
48 u = fill(3.124)
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
49 @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
50 @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
51 @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
52 @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
53 @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
54
1108
6b24dc2d7b11 Clean up testsets in boundary_operator_test
Jonatan Werpers <jonatan@werpers.com>
parents: 1107
diff changeset
55 @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
56 @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
57 @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
58 @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
59 @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
60 end
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
61
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
62 @testset "Inferred" begin
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
63 v = ones(Float64, 11)
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
64 u = fill(1.)
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
65
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
66 @inferred apply(op_l, v)
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
67 @inferred apply(op_r, v)
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 @inferred apply_transpose(op_l, u, 4)
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
70 @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
71 @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
72 @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
73 @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
74 @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
75
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
76 @inferred apply_transpose(op_r, u, 4)
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
77 @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
78 @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
79 @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
80 @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
81 @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
82 end
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
83 end