Mercurial > repos > public > sbplib_julia

annotate test/testDiffOps.jl @ 338:2b0c9b30ea3b refactor/combine_to_one_package

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add test sets for each submodule to make the test output nicer
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 25 Sep 2020 13:48:23 +0200
parents f4e3e71a4ff4
children ffddaf053085
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
235
a5fdc00d5070 Fix a bunch of compilation errors
Jonatan Werpers <jonatan@werpers.com>
parents: 216
diff changeset
1 using Test
335
f4e3e71a4ff4 Fix `using` commands to refer to local modules within the Sbplib package/module
Jonatan Werpers <jonatan@werpers.com>
parents: 333
diff changeset
2 using Sbplib.DiffOps
f4e3e71a4ff4 Fix `using` commands to refer to local modules within the Sbplib package/module
Jonatan Werpers <jonatan@werpers.com>
parents: 333
diff changeset
3 using Sbplib.Grids
f4e3e71a4ff4 Fix `using` commands to refer to local modules within the Sbplib package/module
Jonatan Werpers <jonatan@werpers.com>
parents: 333
diff changeset
4 using Sbplib.SbpOperators
f4e3e71a4ff4 Fix `using` commands to refer to local modules within the Sbplib package/module
Jonatan Werpers <jonatan@werpers.com>
parents: 333
diff changeset
5 using Sbplib.RegionIndices
f4e3e71a4ff4 Fix `using` commands to refer to local modules within the Sbplib package/module
Jonatan Werpers <jonatan@werpers.com>
parents: 333
diff changeset
6 using Sbplib.LazyTensors
216
30112f73555c Add Manifest.toml, test stub and dependencies to DiffOps
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
7
338
2b0c9b30ea3b Add test sets for each submodule to make the test output nicer
Jonatan Werpers <jonatan@werpers.com>
parents: 335
diff changeset
8 @testset "DiffOps" begin
2b0c9b30ea3b Add test sets for each submodule to make the test output nicer
Jonatan Werpers <jonatan@werpers.com>
parents: 335
diff changeset
9
282
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
10 @testset "Laplace2D" begin
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
11 op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
12 Lx = 3.5
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
13 Ly = 7.2
283
12a12a5cd973 Fix tests for Laplace2D.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 282
diff changeset
14 g = EquidistantGrid((42,41), (0.0, 0.0), (Lx,Ly))
282
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
15 L = Laplace(g, 1., op)
283
12a12a5cd973 Fix tests for Laplace2D.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 282
diff changeset
16 H = quadrature(L)
282
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
17
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
18 f0(x::Float64,y::Float64) = 2.
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
19 f1(x::Float64,y::Float64) = x+y
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
20 f2(x::Float64,y::Float64) = 1/2*x^2 + 1/2*y^2
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
21 f3(x::Float64,y::Float64) = 1/6*x^3 + 1/6*y^3
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
22 f4(x::Float64,y::Float64) = 1/24*x^4 + 1/24*y^4
283
12a12a5cd973 Fix tests for Laplace2D.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 282
diff changeset
23 f5(x::Float64,y::Float64) = sin(x) + cos(y)
12a12a5cd973 Fix tests for Laplace2D.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 282
diff changeset
24 f5ₓₓ(x::Float64,y::Float64) = -f5(x,y)
282
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
25
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
26 v0 = evalOn(g,f0)
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
27 v1 = evalOn(g,f1)
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
28 v2 = evalOn(g,f2)
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
29 v3 = evalOn(g,f3)
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
30 v4 = evalOn(g,f4)
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
31 v5 = evalOn(g,f5)
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
32 v5ₓₓ = evalOn(g,f5ₓₓ)
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
33
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
34 @test L isa TensorOperator{T,2} where T
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
35 @test L' isa TensorMapping{T,2,2} where T
283
12a12a5cd973 Fix tests for Laplace2D.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 282
diff changeset
36
282
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
37 # TODO: Should perhaps set tolerance level for isapporx instead?
283
12a12a5cd973 Fix tests for Laplace2D.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 282
diff changeset
38 # Are these tolerance levels resonable or should tests be constructed
12a12a5cd973 Fix tests for Laplace2D.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 282
diff changeset
39 # differently?
12a12a5cd973 Fix tests for Laplace2D.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 282
diff changeset
40 equalitytol = 0.5*1e-10
12a12a5cd973 Fix tests for Laplace2D.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 282
diff changeset
41 accuracytol = 0.5*1e-3
12a12a5cd973 Fix tests for Laplace2D.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 282
diff changeset
42 # 4th order interior stencil, 2nd order boundary stencil,
12a12a5cd973 Fix tests for Laplace2D.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 282
diff changeset
43 # implies that L*v should be exact for v - monomial up to order 3.
12a12a5cd973 Fix tests for Laplace2D.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 282
diff changeset
44 # Exact differentiation is measured point-wise. For other grid functions
12a12a5cd973 Fix tests for Laplace2D.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 282
diff changeset
45 # the error is measured in the H-norm.
12a12a5cd973 Fix tests for Laplace2D.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 282
diff changeset
46 @test all(abs.(collect(L*v0)) .<= equalitytol)
12a12a5cd973 Fix tests for Laplace2D.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 282
diff changeset
47 @test all(abs.(collect(L*v1)) .<= equalitytol)
282
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
48 @test all(collect(L*v2) .≈ v0) # Seems to be more accurate
283
12a12a5cd973 Fix tests for Laplace2D.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 282
diff changeset
49 @test all(abs.((collect(L*v3) - v1)) .<= equalitytol)
12a12a5cd973 Fix tests for Laplace2D.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 282
diff changeset
50 e4 = collect(L*v4) - v2
12a12a5cd973 Fix tests for Laplace2D.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 282
diff changeset
51 e5 = collect(L*v5) - v5ₓₓ
12a12a5cd973 Fix tests for Laplace2D.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 282
diff changeset
52 @test sum(collect(H*e4.^2)) <= accuracytol
12a12a5cd973 Fix tests for Laplace2D.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 282
diff changeset
53 @test sum(collect(H*e5.^2)) <= accuracytol
282
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
54 end
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
55
262
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
56 @testset "Quadrature" begin
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
57 op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
58 Lx = 2.3
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
59 Ly = 5.2
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
60 g = EquidistantGrid((77,66), (0.0, 0.0), (Lx,Ly))
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
61 H = Quadrature(op,g)
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
62 v = ones(Float64, size(g))
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
63
282
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
64 @test H isa TensorOperator{T,2} where T
262
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
65 @test H' isa TensorMapping{T,2,2} where T
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
66 @test sum(collect(H*v)) ≈ (Lx*Ly)
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
67 @test collect(H*v) == collect(H'*v)
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
68 end
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
69
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
70 @testset "InverseQuadrature" begin
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
71 op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
72 Lx = 7.3
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
73 Ly = 8.2
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
74 g = EquidistantGrid((77,66), (0.0, 0.0), (Lx,Ly))
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
75 H = Quadrature(op,g)
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
76 Hinv = InverseQuadrature(op,g)
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
77 v = evalOn(g, (x,y)-> x^2 + (y-1)^2 + x*y)
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
78
282
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
79 @test Hinv isa TensorOperator{T,2} where T
262
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
80 @test Hinv' isa TensorMapping{T,2,2} where T
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
81 @test collect(Hinv*H*v) ≈ v
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
82 @test collect(Hinv*v) == collect(Hinv'*v)
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
83 end
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
84
235
a5fdc00d5070 Fix a bunch of compilation errors
Jonatan Werpers <jonatan@werpers.com>
parents: 216
diff changeset
85 @testset "BoundaryValue" begin
a5fdc00d5070 Fix a bunch of compilation errors
Jonatan Werpers <jonatan@werpers.com>
parents: 216
diff changeset
86 op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
239
60011a10e17d Add tests for BoundaryValue and fix index types
Jonatan Werpers <jonatan@werpers.com>
parents: 235
diff changeset
87 g = EquidistantGrid((4,5), (0.0, 0.0), (1.0,1.0))
235
a5fdc00d5070 Fix a bunch of compilation errors
Jonatan Werpers <jonatan@werpers.com>
parents: 216
diff changeset
88
a5fdc00d5070 Fix a bunch of compilation errors
Jonatan Werpers <jonatan@werpers.com>
parents: 216
diff changeset
89 e_w = BoundaryValue(op, g, CartesianBoundary{1,Lower}())
a5fdc00d5070 Fix a bunch of compilation errors
Jonatan Werpers <jonatan@werpers.com>
parents: 216
diff changeset
90 e_e = BoundaryValue(op, g, CartesianBoundary{1,Upper}())
a5fdc00d5070 Fix a bunch of compilation errors
Jonatan Werpers <jonatan@werpers.com>
parents: 216
diff changeset
91 e_s = BoundaryValue(op, g, CartesianBoundary{2,Lower}())
a5fdc00d5070 Fix a bunch of compilation errors
Jonatan Werpers <jonatan@werpers.com>
parents: 216
diff changeset
92 e_n = BoundaryValue(op, g, CartesianBoundary{2,Upper}())
a5fdc00d5070 Fix a bunch of compilation errors
Jonatan Werpers <jonatan@werpers.com>
parents: 216
diff changeset
93
239
60011a10e17d Add tests for BoundaryValue and fix index types
Jonatan Werpers <jonatan@werpers.com>
parents: 235
diff changeset
94 v = zeros(Float64, 4, 5)
60011a10e17d Add tests for BoundaryValue and fix index types
Jonatan Werpers <jonatan@werpers.com>
parents: 235
diff changeset
95 v[:,5] = [1, 2, 3,4]
60011a10e17d Add tests for BoundaryValue and fix index types
Jonatan Werpers <jonatan@werpers.com>
parents: 235
diff changeset
96 v[:,4] = [1, 2, 3,4]
60011a10e17d Add tests for BoundaryValue and fix index types
Jonatan Werpers <jonatan@werpers.com>
parents: 235
diff changeset
97 v[:,3] = [4, 5, 6, 7]
60011a10e17d Add tests for BoundaryValue and fix index types
Jonatan Werpers <jonatan@werpers.com>
parents: 235
diff changeset
98 v[:,2] = [7, 8, 9, 10]
60011a10e17d Add tests for BoundaryValue and fix index types
Jonatan Werpers <jonatan@werpers.com>
parents: 235
diff changeset
99 v[:,1] = [10, 11, 12, 13]
235
a5fdc00d5070 Fix a bunch of compilation errors
Jonatan Werpers <jonatan@werpers.com>
parents: 216
diff changeset
100
a5fdc00d5070 Fix a bunch of compilation errors
Jonatan Werpers <jonatan@werpers.com>
parents: 216
diff changeset
101 @test e_w isa TensorMapping{T,2,1} where T
a5fdc00d5070 Fix a bunch of compilation errors
Jonatan Werpers <jonatan@werpers.com>
parents: 216
diff changeset
102 @test e_w' isa TensorMapping{T,1,2} where T
a5fdc00d5070 Fix a bunch of compilation errors
Jonatan Werpers <jonatan@werpers.com>
parents: 216
diff changeset
103
a5fdc00d5070 Fix a bunch of compilation errors
Jonatan Werpers <jonatan@werpers.com>
parents: 216
diff changeset
104 @test domain_size(e_w, (3,2)) == (2,)
a5fdc00d5070 Fix a bunch of compilation errors
Jonatan Werpers <jonatan@werpers.com>
parents: 216
diff changeset
105 @test domain_size(e_e, (3,2)) == (2,)
a5fdc00d5070 Fix a bunch of compilation errors
Jonatan Werpers <jonatan@werpers.com>
parents: 216
diff changeset
106 @test domain_size(e_s, (3,2)) == (3,)
a5fdc00d5070 Fix a bunch of compilation errors
Jonatan Werpers <jonatan@werpers.com>
parents: 216
diff changeset
107 @test domain_size(e_n, (3,2)) == (3,)
a5fdc00d5070 Fix a bunch of compilation errors
Jonatan Werpers <jonatan@werpers.com>
parents: 216
diff changeset
108
239
60011a10e17d Add tests for BoundaryValue and fix index types
Jonatan Werpers <jonatan@werpers.com>
parents: 235
diff changeset
109 @test size(e_w'*v) == (5,)
60011a10e17d Add tests for BoundaryValue and fix index types
Jonatan Werpers <jonatan@werpers.com>
parents: 235
diff changeset
110 @test size(e_e'*v) == (5,)
60011a10e17d Add tests for BoundaryValue and fix index types
Jonatan Werpers <jonatan@werpers.com>
parents: 235
diff changeset
111 @test size(e_s'*v) == (4,)
60011a10e17d Add tests for BoundaryValue and fix index types
Jonatan Werpers <jonatan@werpers.com>
parents: 235
diff changeset
112 @test size(e_n'*v) == (4,)
235
a5fdc00d5070 Fix a bunch of compilation errors
Jonatan Werpers <jonatan@werpers.com>
parents: 216
diff changeset
113
239
60011a10e17d Add tests for BoundaryValue and fix index types
Jonatan Werpers <jonatan@werpers.com>
parents: 235
diff changeset
114 @test collect(e_w'*v) == [10,7,4,1.0,1]
60011a10e17d Add tests for BoundaryValue and fix index types
Jonatan Werpers <jonatan@werpers.com>
parents: 235
diff changeset
115 @test collect(e_e'*v) == [13,10,7,4,4.0]
60011a10e17d Add tests for BoundaryValue and fix index types
Jonatan Werpers <jonatan@werpers.com>
parents: 235
diff changeset
116 @test collect(e_s'*v) == [10,11,12,13.0]
60011a10e17d Add tests for BoundaryValue and fix index types
Jonatan Werpers <jonatan@werpers.com>
parents: 235
diff changeset
117 @test collect(e_n'*v) == [1,2,3,4.0]
235
a5fdc00d5070 Fix a bunch of compilation errors
Jonatan Werpers <jonatan@werpers.com>
parents: 216
diff changeset
118
242
9819243102dd Add test for and fix apply(::BoundaryValue)
Jonatan Werpers <jonatan@werpers.com>
parents: 239
diff changeset
119 g_x = [1,2,3,4.0]
9819243102dd Add test for and fix apply(::BoundaryValue)
Jonatan Werpers <jonatan@werpers.com>
parents: 239
diff changeset
120 g_y = [5,4,3,2,1.0]
9819243102dd Add test for and fix apply(::BoundaryValue)
Jonatan Werpers <jonatan@werpers.com>
parents: 239
diff changeset
121
9819243102dd Add test for and fix apply(::BoundaryValue)
Jonatan Werpers <jonatan@werpers.com>
parents: 239
diff changeset
122 G_w = zeros(Float64, (4,5))
9819243102dd Add test for and fix apply(::BoundaryValue)
Jonatan Werpers <jonatan@werpers.com>
parents: 239
diff changeset
123 G_w[1,:] = g_y
9819243102dd Add test for and fix apply(::BoundaryValue)
Jonatan Werpers <jonatan@werpers.com>
parents: 239
diff changeset
124
9819243102dd Add test for and fix apply(::BoundaryValue)
Jonatan Werpers <jonatan@werpers.com>
parents: 239
diff changeset
125 G_e = zeros(Float64, (4,5))
9819243102dd Add test for and fix apply(::BoundaryValue)
Jonatan Werpers <jonatan@werpers.com>
parents: 239
diff changeset
126 G_e[4,:] = g_y
9819243102dd Add test for and fix apply(::BoundaryValue)
Jonatan Werpers <jonatan@werpers.com>
parents: 239
diff changeset
127
9819243102dd Add test for and fix apply(::BoundaryValue)
Jonatan Werpers <jonatan@werpers.com>
parents: 239
diff changeset
128 G_s = zeros(Float64, (4,5))
9819243102dd Add test for and fix apply(::BoundaryValue)
Jonatan Werpers <jonatan@werpers.com>
parents: 239
diff changeset
129 G_s[:,1] = g_x
9819243102dd Add test for and fix apply(::BoundaryValue)
Jonatan Werpers <jonatan@werpers.com>
parents: 239
diff changeset
130
9819243102dd Add test for and fix apply(::BoundaryValue)
Jonatan Werpers <jonatan@werpers.com>
parents: 239
diff changeset
131 G_n = zeros(Float64, (4,5))
9819243102dd Add test for and fix apply(::BoundaryValue)
Jonatan Werpers <jonatan@werpers.com>
parents: 239
diff changeset
132 G_n[:,5] = g_x
9819243102dd Add test for and fix apply(::BoundaryValue)
Jonatan Werpers <jonatan@werpers.com>
parents: 239
diff changeset
133
285
e21dcda55163 Add type representing an unknown dimension in the domain of tensormappings
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 284
diff changeset
134 @test size(e_w*g_y) == (UnknownDim,5)
e21dcda55163 Add type representing an unknown dimension in the domain of tensormappings
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 284
diff changeset
135 @test size(e_e*g_y) == (UnknownDim,5)
e21dcda55163 Add type representing an unknown dimension in the domain of tensormappings
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 284
diff changeset
136 @test size(e_s*g_x) == (4,UnknownDim)
e21dcda55163 Add type representing an unknown dimension in the domain of tensormappings
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 284
diff changeset
137 @test size(e_n*g_x) == (4,UnknownDim)
244
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
138
284
0b8e041a1873 Change how range_size and domain_size work with BoundaryValues and NormalDerivative
Jonatan Werpers <jonatan@werpers.com>
parents: 283
diff changeset
139 # These tests should be moved to where they are possible (i.e we know what the grid should be)
0b8e041a1873 Change how range_size and domain_size work with BoundaryValues and NormalDerivative
Jonatan Werpers <jonatan@werpers.com>
parents: 283
diff changeset
140 @test_broken collect(e_w*g_y) == G_w
0b8e041a1873 Change how range_size and domain_size work with BoundaryValues and NormalDerivative
Jonatan Werpers <jonatan@werpers.com>
parents: 283
diff changeset
141 @test_broken collect(e_e*g_y) == G_e
0b8e041a1873 Change how range_size and domain_size work with BoundaryValues and NormalDerivative
Jonatan Werpers <jonatan@werpers.com>
parents: 283
diff changeset
142 @test_broken collect(e_s*g_x) == G_s
0b8e041a1873 Change how range_size and domain_size work with BoundaryValues and NormalDerivative
Jonatan Werpers <jonatan@werpers.com>
parents: 283
diff changeset
143 @test_broken collect(e_n*g_x) == G_n
244
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
144 end
242
9819243102dd Add test for and fix apply(::BoundaryValue)
Jonatan Werpers <jonatan@werpers.com>
parents: 239
diff changeset
145
244
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
146 @testset "NormalDerivative" begin
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
147 op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
148 g = EquidistantGrid((5,6), (0.0, 0.0), (4.0,5.0))
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
149
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
150 d_w = NormalDerivative(op, g, CartesianBoundary{1,Lower}())
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
151 d_e = NormalDerivative(op, g, CartesianBoundary{1,Upper}())
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
152 d_s = NormalDerivative(op, g, CartesianBoundary{2,Lower}())
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
153 d_n = NormalDerivative(op, g, CartesianBoundary{2,Upper}())
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
154
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
155
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
156 v = evalOn(g, (x,y)-> x^2 + (y-1)^2 + x*y)
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
157 v∂x = evalOn(g, (x,y)-> 2*x + y)
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
158 v∂y = evalOn(g, (x,y)-> 2*(y-1) + x)
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
159
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
160 @test d_w isa TensorMapping{T,2,1} where T
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
161 @test d_w' isa TensorMapping{T,1,2} where T
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
162
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
163 @test domain_size(d_w, (3,2)) == (2,)
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
164 @test domain_size(d_e, (3,2)) == (2,)
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
165 @test domain_size(d_s, (3,2)) == (3,)
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
166 @test domain_size(d_n, (3,2)) == (3,)
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
167
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
168 @test size(d_w'*v) == (6,)
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
169 @test size(d_e'*v) == (6,)
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
170 @test size(d_s'*v) == (5,)
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
171 @test size(d_n'*v) == (5,)
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
172
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
173 @test collect(d_w'*v) ≈ v∂x[1,:]
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
174 @test collect(d_e'*v) ≈ v∂x[5,:]
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
175 @test collect(d_s'*v) ≈ v∂y[:,1]
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
176 @test collect(d_n'*v) ≈ v∂y[:,6]
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
177
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
178
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
179 d_x_l = zeros(Float64, 5)
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
180 d_x_u = zeros(Float64, 5)
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
181 for i ∈ eachindex(d_x_l)
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
182 d_x_l[i] = op.dClosure[i-1]
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
183 d_x_u[i] = -op.dClosure[length(d_x_u)-i]
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
184 end
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
185
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
186 d_y_l = zeros(Float64, 6)
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
187 d_y_u = zeros(Float64, 6)
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
188 for i ∈ eachindex(d_y_l)
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
189 d_y_l[i] = op.dClosure[i-1]
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
190 d_y_u[i] = -op.dClosure[length(d_y_u)-i]
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
191 end
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
192
282
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
193 function prod_matrix(x,y)
244
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
194 G = zeros(Float64, length(x), length(y))
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
195 for I ∈ CartesianIndices(G)
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
196 G[I] = x[I[1]]*y[I[2]]
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
197 end
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
198
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
199 return G
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
200 end
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
201
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
202 g_x = [1,2,3,4.0,5]
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
203 g_y = [5,4,3,2,1.0,11]
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
204
282
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
205 G_w = prod_matrix(d_x_l, g_y)
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
206 G_e = prod_matrix(d_x_u, g_y)
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
207 G_s = prod_matrix(g_x, d_y_l)
ce6a2f3f732a Make Laplace a TensorOperator and add tests. NOTE: Two of the tests for Laplace2D are currently failing.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 262
diff changeset
208 G_n = prod_matrix(g_x, d_y_u)
244
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
209
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
210
285
e21dcda55163 Add type representing an unknown dimension in the domain of tensormappings
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 284
diff changeset
211 @test size(d_w*g_y) == (UnknownDim,6)
e21dcda55163 Add type representing an unknown dimension in the domain of tensormappings
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 284
diff changeset
212 @test size(d_e*g_y) == (UnknownDim,6)
e21dcda55163 Add type representing an unknown dimension in the domain of tensormappings
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 284
diff changeset
213 @test size(d_s*g_x) == (5,UnknownDim)
e21dcda55163 Add type representing an unknown dimension in the domain of tensormappings
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 284
diff changeset
214 @test size(d_n*g_x) == (5,UnknownDim)
244
a827568fc251 Fix NormalDerivative and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 242
diff changeset
215
284
0b8e041a1873 Change how range_size and domain_size work with BoundaryValues and NormalDerivative
Jonatan Werpers <jonatan@werpers.com>
parents: 283
diff changeset
216 # These tests should be moved to where they are possible (i.e we know what the grid should be)
0b8e041a1873 Change how range_size and domain_size work with BoundaryValues and NormalDerivative
Jonatan Werpers <jonatan@werpers.com>
parents: 283
diff changeset
217 @test_broken collect(d_w*g_y) ≈ G_w
0b8e041a1873 Change how range_size and domain_size work with BoundaryValues and NormalDerivative
Jonatan Werpers <jonatan@werpers.com>
parents: 283
diff changeset
218 @test_broken collect(d_e*g_y) ≈ G_e
0b8e041a1873 Change how range_size and domain_size work with BoundaryValues and NormalDerivative
Jonatan Werpers <jonatan@werpers.com>
parents: 283
diff changeset
219 @test_broken collect(d_s*g_x) ≈ G_s
0b8e041a1873 Change how range_size and domain_size work with BoundaryValues and NormalDerivative
Jonatan Werpers <jonatan@werpers.com>
parents: 283
diff changeset
220 @test_broken collect(d_n*g_x) ≈ G_n
235
a5fdc00d5070 Fix a bunch of compilation errors
Jonatan Werpers <jonatan@werpers.com>
parents: 216
diff changeset
221 end
259
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
222
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
223 @testset "BoundaryQuadrature" begin
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
224 op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
225 g = EquidistantGrid((10,11), (0.0, 0.0), (1.0,1.0))
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
226
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
227 H_w = BoundaryQuadrature(op, g, CartesianBoundary{1,Lower}())
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
228 H_e = BoundaryQuadrature(op, g, CartesianBoundary{1,Upper}())
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
229 H_s = BoundaryQuadrature(op, g, CartesianBoundary{2,Lower}())
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
230 H_n = BoundaryQuadrature(op, g, CartesianBoundary{2,Upper}())
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
231
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
232 v = evalOn(g, (x,y)-> x^2 + (y-1)^2 + x*y)
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
233
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
234 function get_quadrature(N)
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
235 qc = op.quadratureClosure
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
236 q = (qc..., ones(N-2*closuresize(op))..., reverse(qc)...)
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
237 @assert length(q) == N
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
238 return q
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
239 end
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
240
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
241 v_w = v[1,:]
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
242 v_e = v[10,:]
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
243 v_s = v[:,1]
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
244 v_n = v[:,11]
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
245
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
246 q_x = spacing(g)[1].*get_quadrature(10)
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
247 q_y = spacing(g)[2].*get_quadrature(11)
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
248
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
249 @test H_w isa TensorOperator{T,1} where T
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
250
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
251 @test domain_size(H_w, (3,)) == (3,)
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
252 @test domain_size(H_n, (3,)) == (3,)
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
253
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
254 @test range_size(H_w, (3,)) == (3,)
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
255 @test range_size(H_n, (3,)) == (3,)
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
256
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
257 @test size(H_w*v_w) == (11,)
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
258 @test size(H_e*v_e) == (11,)
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
259 @test size(H_s*v_s) == (10,)
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
260 @test size(H_n*v_n) == (10,)
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
261
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
262 @test collect(H_w*v_w) ≈ q_y.*v_w
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
263 @test collect(H_e*v_e) ≈ q_y.*v_e
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
264 @test collect(H_s*v_s) ≈ q_x.*v_s
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
265 @test collect(H_n*v_n) ≈ q_x.*v_n
262
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
266
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
267 @test collect(H_w'*v_w) == collect(H_w'*v_w)
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
268 @test collect(H_e'*v_e) == collect(H_e'*v_e)
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
269 @test collect(H_s'*v_s) == collect(H_s'*v_s)
f1e90a92ad74 Add Quadrature and InverseQuadrature for Laplace as TensorMappings. Implement and test the 2D case. Fix implementation of apply_transpose for BoundaryQuadrature and add tests.
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 259
diff changeset
270 @test collect(H_n'*v_n) == collect(H_n'*v_n)
259
5571d2c5bf0f Implement BaoundaryQuadrature for Laplace
Jonatan Werpers <jonatan@werpers.com>
parents: 244
diff changeset
271 end
338
2b0c9b30ea3b Add test sets for each submodule to make the test output nicer
Jonatan Werpers <jonatan@werpers.com>
parents: 335
diff changeset
272
2b0c9b30ea3b Add test sets for each submodule to make the test output nicer
Jonatan Werpers <jonatan@werpers.com>
parents: 335
diff changeset
273 end