Mercurial > repos > public > sbplib_julia

annotate test/SbpOperators/stencil_test.jl @ 1956:b0fcb29e3620 feature/grids/multiblock_boundaries

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge default
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 11 Feb 2025 08:54:18 +0100
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
1726
471a948cd2b2 Rename project from Sbplib to Diffinitive
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents: 1452
diff changeset
2 using Diffinitive.SbpOperators
1446
d744b01be520 Add broken tests for applying stencils to vectors of vectors
Jonatan Werpers <jonatan@werpers.com>
parents: 1445
diff changeset
3 using StaticArrays
1726
471a948cd2b2 Rename project from Sbplib to Diffinitive
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents: 1452
diff changeset
4 import Diffinitive.SbpOperators.Stencil
471a948cd2b2 Rename project from Sbplib to Diffinitive
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents: 1452
diff changeset
5 import Diffinitive.SbpOperators.NestedStencil
471a948cd2b2 Rename project from Sbplib to Diffinitive
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents: 1452
diff changeset
6 import Diffinitive.SbpOperators.scale
471a948cd2b2 Rename project from Sbplib to Diffinitive
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents: 1452
diff changeset
7 import Diffinitive.SbpOperators: apply_stencil, apply_stencil_backwards
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
8
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
9 @testset "Stencil" begin
897
737cd68318c7 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 896
diff changeset
10 s = Stencil(-2:2, (1.,2.,2.,3.,4.))
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
11 @test s isa Stencil{Float64, 5}
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
12
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
13 @test eltype(s) == Float64
887
0be29e65521e Add length method for stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 879
diff changeset
14
0be29e65521e Add length method for stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 879
diff changeset
15 @test length(s) == 5
897
737cd68318c7 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 896
diff changeset
16 @test length(Stencil(-1:2, (1,2,3,4))) == 4
887
0be29e65521e Add length method for stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 879
diff changeset
17
897
737cd68318c7 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 896
diff changeset
18 @test SbpOperators.scale(s, 2) == Stencil(-2:2, (2.,4.,4.,6.,8.))
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
19
897
737cd68318c7 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 896
diff changeset
20 @test Stencil(1,2,3,4; center=1) == Stencil(0:3,(1,2,3,4))
737cd68318c7 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 896
diff changeset
21 @test Stencil(1,2,3,4; center=2) == Stencil(-1:2,(1,2,3,4))
737cd68318c7 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 896
diff changeset
22 @test Stencil(1,2,3,4; center=4) == Stencil(-3:0,(1,2,3,4))
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
23
897
737cd68318c7 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 896
diff changeset
24 @test CenteredStencil(1,2,3,4,5) == Stencil(-2:2, (1,2,3,4,5))
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
25 @test_throws ArgumentError CenteredStencil(1,2,3,4)
826
4433be383840 Add stencil constructor to change the type of the weights along with a convert method
Jonatan Werpers <jonatan@werpers.com>
parents: 732
diff changeset
26
4433be383840 Add stencil constructor to change the type of the weights along with a convert method
Jonatan Werpers <jonatan@werpers.com>
parents: 732
diff changeset
27 # Changing the type of the weights
4433be383840 Add stencil constructor to change the type of the weights along with a convert method
Jonatan Werpers <jonatan@werpers.com>
parents: 732
diff changeset
28 @test Stencil{Float64}(Stencil(1,2,3,4,5; center=2)) == Stencil(1.,2.,3.,4.,5.; center=2)
4433be383840 Add stencil constructor to change the type of the weights along with a convert method
Jonatan Werpers <jonatan@werpers.com>
parents: 732
diff changeset
29 @test Stencil{Float64}(CenteredStencil(1,2,3,4,5)) == CenteredStencil(1.,2.,3.,4.,5.)
4433be383840 Add stencil constructor to change the type of the weights along with a convert method
Jonatan Werpers <jonatan@werpers.com>
parents: 732
diff changeset
30 @test Stencil{Int}(Stencil(1.,2.,3.,4.,5.; center=2)) == Stencil(1,2,3,4,5; center=2)
4433be383840 Add stencil constructor to change the type of the weights along with a convert method
Jonatan Werpers <jonatan@werpers.com>
parents: 732
diff changeset
31 @test Stencil{Rational}(Stencil(1.,2.,3.,4.,5.; center=2)) == Stencil(1//1,2//1,3//1,4//1,5//1; center=2)
4433be383840 Add stencil constructor to change the type of the weights along with a convert method
Jonatan Werpers <jonatan@werpers.com>
parents: 732
diff changeset
32
4433be383840 Add stencil constructor to change the type of the weights along with a convert method
Jonatan Werpers <jonatan@werpers.com>
parents: 732
diff changeset
33 @testset "convert" begin
4433be383840 Add stencil constructor to change the type of the weights along with a convert method
Jonatan Werpers <jonatan@werpers.com>
parents: 732
diff changeset
34 @test convert(Stencil{Float64}, Stencil(1,2,3,4,5; center=2)) == Stencil(1.,2.,3.,4.,5.; center=2)
897
737cd68318c7 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 896
diff changeset
35 @test convert(Stencil{Float64,5}, CenteredStencil(1,2,3,4,5)) == CenteredStencil(1.,2.,3.,4.,5.)
737cd68318c7 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 896
diff changeset
36 @test convert(Stencil{Int,5}, Stencil(1.,2.,3.,4.,5.; center=2)) == Stencil(1,2,3,4,5; center=2)
737cd68318c7 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 896
diff changeset
37 @test convert(Stencil{Rational,5}, Stencil(1.,2.,3.,4.,5.; center=2)) == Stencil(1//1,2//1,3//1,4//1,5//1; center=2)
737cd68318c7 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 896
diff changeset
38 end
737cd68318c7 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 896
diff changeset
39
737cd68318c7 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 896
diff changeset
40 @testset "promotion of weights" begin
737cd68318c7 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 896
diff changeset
41 @test Stencil(1.,2; center = 1) isa Stencil{Float64, 2}
737cd68318c7 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 896
diff changeset
42 @test Stencil(1,2//2; center = 1) isa Stencil{Rational{Int64}, 2}
737cd68318c7 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 896
diff changeset
43 end
737cd68318c7 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 896
diff changeset
44
737cd68318c7 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 896
diff changeset
45 @testset "promotion" begin
737cd68318c7 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 896
diff changeset
46 @test promote(Stencil(1,1;center=1), Stencil(2.,2.;center=2)) == (Stencil(1.,1.;center=1), Stencil(2.,2.;center=2))
826
4433be383840 Add stencil constructor to change the type of the weights along with a convert method
Jonatan Werpers <jonatan@werpers.com>
parents: 732
diff changeset
47 end
903
a7f898b1ce1e Fix type instability for stencil application when weights and vector didn't match
Jonatan Werpers <jonatan@werpers.com>
parents: 898
diff changeset
48
1444
c779af4d02e5 Add tests for apply_stencil(::Stencil) and apply_stencil_backwards(::Stencil)
Jonatan Werpers <jonatan@werpers.com>
parents: 1072
diff changeset
49 @testset "apply_stencil" begin
c779af4d02e5 Add tests for apply_stencil(::Stencil) and apply_stencil_backwards(::Stencil)
Jonatan Werpers <jonatan@werpers.com>
parents: 1072
diff changeset
50 v = [1, 2, 4, 8, 16, 32, 64, 128]
c779af4d02e5 Add tests for apply_stencil(::Stencil) and apply_stencil_backwards(::Stencil)
Jonatan Werpers <jonatan@werpers.com>
parents: 1072
diff changeset
51 s = Stencil(1,2,3,4, center = 2)
c779af4d02e5 Add tests for apply_stencil(::Stencil) and apply_stencil_backwards(::Stencil)
Jonatan Werpers <jonatan@werpers.com>
parents: 1072
diff changeset
52 @test apply_stencil(s,v, 2) == v[1] + 2*v[2] + 3*v[3] + 4*v[4]
c779af4d02e5 Add tests for apply_stencil(::Stencil) and apply_stencil_backwards(::Stencil)
Jonatan Werpers <jonatan@werpers.com>
parents: 1072
diff changeset
53 @test apply_stencil(s,v, 4) == v[3] + 2*v[4] + 3*v[5] + 4*v[6]
c779af4d02e5 Add tests for apply_stencil(::Stencil) and apply_stencil_backwards(::Stencil)
Jonatan Werpers <jonatan@werpers.com>
parents: 1072
diff changeset
54 @test apply_stencil_backwards(s,v, 3) == 4*v[1] + 3*v[2] + 2*v[3] + 1*v[4]
c779af4d02e5 Add tests for apply_stencil(::Stencil) and apply_stencil_backwards(::Stencil)
Jonatan Werpers <jonatan@werpers.com>
parents: 1072
diff changeset
55 @test apply_stencil_backwards(s,v, 7) == 4*v[5] + 3*v[6] + 2*v[7] + 1*v[8]
c779af4d02e5 Add tests for apply_stencil(::Stencil) and apply_stencil_backwards(::Stencil)
Jonatan Werpers <jonatan@werpers.com>
parents: 1072
diff changeset
56 @test apply_stencil(s,v, 2) isa Int
c779af4d02e5 Add tests for apply_stencil(::Stencil) and apply_stencil_backwards(::Stencil)
Jonatan Werpers <jonatan@werpers.com>
parents: 1072
diff changeset
57 @test apply_stencil_backwards(s,v, 7) isa Int
c779af4d02e5 Add tests for apply_stencil(::Stencil) and apply_stencil_backwards(::Stencil)
Jonatan Werpers <jonatan@werpers.com>
parents: 1072
diff changeset
58
c779af4d02e5 Add tests for apply_stencil(::Stencil) and apply_stencil_backwards(::Stencil)
Jonatan Werpers <jonatan@werpers.com>
parents: 1072
diff changeset
59 v = [1, 2, 4, 8, 16, 32, 64, 128]
c779af4d02e5 Add tests for apply_stencil(::Stencil) and apply_stencil_backwards(::Stencil)
Jonatan Werpers <jonatan@werpers.com>
parents: 1072
diff changeset
60 s = Stencil(1.,2.,3.,4., center = 2)
c779af4d02e5 Add tests for apply_stencil(::Stencil) and apply_stencil_backwards(::Stencil)
Jonatan Werpers <jonatan@werpers.com>
parents: 1072
diff changeset
61 @test apply_stencil(s,v, 4) == v[3] + 2. *v[4] + 3. *v[5] + 4. *v[6]
c779af4d02e5 Add tests for apply_stencil(::Stencil) and apply_stencil_backwards(::Stencil)
Jonatan Werpers <jonatan@werpers.com>
parents: 1072
diff changeset
62 @test apply_stencil_backwards(s,v, 7) == 4. *v[5] + 3. *v[6] + 2. *v[7] + v[8]
c779af4d02e5 Add tests for apply_stencil(::Stencil) and apply_stencil_backwards(::Stencil)
Jonatan Werpers <jonatan@werpers.com>
parents: 1072
diff changeset
63 @test apply_stencil(s,v, 2) isa Float64
c779af4d02e5 Add tests for apply_stencil(::Stencil) and apply_stencil_backwards(::Stencil)
Jonatan Werpers <jonatan@werpers.com>
parents: 1072
diff changeset
64 @test apply_stencil_backwards(s,v, 7) isa Float64
c779af4d02e5 Add tests for apply_stencil(::Stencil) and apply_stencil_backwards(::Stencil)
Jonatan Werpers <jonatan@werpers.com>
parents: 1072
diff changeset
65
c779af4d02e5 Add tests for apply_stencil(::Stencil) and apply_stencil_backwards(::Stencil)
Jonatan Werpers <jonatan@werpers.com>
parents: 1072
diff changeset
66 v = [1., 2., 4., 8., 16., 32., 64., 128.]
c779af4d02e5 Add tests for apply_stencil(::Stencil) and apply_stencil_backwards(::Stencil)
Jonatan Werpers <jonatan@werpers.com>
parents: 1072
diff changeset
67 s = Stencil(1,2,3,4, center = 2)
c779af4d02e5 Add tests for apply_stencil(::Stencil) and apply_stencil_backwards(::Stencil)
Jonatan Werpers <jonatan@werpers.com>
parents: 1072
diff changeset
68 @test apply_stencil(s,v, 2) == v[1] + 2*v[2] + 3*v[3] + 4*v[4]
c779af4d02e5 Add tests for apply_stencil(::Stencil) and apply_stencil_backwards(::Stencil)
Jonatan Werpers <jonatan@werpers.com>
parents: 1072
diff changeset
69 @test apply_stencil_backwards(s,v, 3) == 4*v[1] + 3*v[2] + 2*v[3] + 1*v[4]
c779af4d02e5 Add tests for apply_stencil(::Stencil) and apply_stencil_backwards(::Stencil)
Jonatan Werpers <jonatan@werpers.com>
parents: 1072
diff changeset
70 @test apply_stencil(s,v, 2) isa Float64
c779af4d02e5 Add tests for apply_stencil(::Stencil) and apply_stencil_backwards(::Stencil)
Jonatan Werpers <jonatan@werpers.com>
parents: 1072
diff changeset
71 @test apply_stencil_backwards(s,v, 3) isa Float64
1446
d744b01be520 Add broken tests for applying stencils to vectors of vectors
Jonatan Werpers <jonatan@werpers.com>
parents: 1445
diff changeset
72
d744b01be520 Add broken tests for applying stencils to vectors of vectors
Jonatan Werpers <jonatan@werpers.com>
parents: 1445
diff changeset
73 v = [@SVector[1, 2], @SVector[3, 4], @SVector[5, 6], @SVector[7, 8]]
d744b01be520 Add broken tests for applying stencils to vectors of vectors
Jonatan Werpers <jonatan@werpers.com>
parents: 1445
diff changeset
74 s = Stencil(1,2, center = 1)
1452
4f79ab676ebc Fix broken tests using sum(f,...). Needs to be evaluated for performance
Jonatan Werpers <jonatan@werpers.com>
parents: 1446
diff changeset
75 @test apply_stencil(s,v,1) == @SVector[7, 10]
4f79ab676ebc Fix broken tests using sum(f,...). Needs to be evaluated for performance
Jonatan Werpers <jonatan@werpers.com>
parents: 1446
diff changeset
76 @test apply_stencil_backwards(s,v,3) == @SVector[11, 14]
4f79ab676ebc Fix broken tests using sum(f,...). Needs to be evaluated for performance
Jonatan Werpers <jonatan@werpers.com>
parents: 1446
diff changeset
77 @test apply_stencil(s,v,1) isa SVector{2, Int}
4f79ab676ebc Fix broken tests using sum(f,...). Needs to be evaluated for performance
Jonatan Werpers <jonatan@werpers.com>
parents: 1446
diff changeset
78 @test apply_stencil_backwards(s,v,3) isa SVector{2, Int}
1446
d744b01be520 Add broken tests for applying stencils to vectors of vectors
Jonatan Werpers <jonatan@werpers.com>
parents: 1445
diff changeset
79
d744b01be520 Add broken tests for applying stencils to vectors of vectors
Jonatan Werpers <jonatan@werpers.com>
parents: 1445
diff changeset
80 v = [@SVector[1., 2.], @SVector[3., 4.], @SVector[5., 6.], @SVector[7., 8.]]
d744b01be520 Add broken tests for applying stencils to vectors of vectors
Jonatan Werpers <jonatan@werpers.com>
parents: 1445
diff changeset
81 s = Stencil(1,2, center = 1)
1452
4f79ab676ebc Fix broken tests using sum(f,...). Needs to be evaluated for performance
Jonatan Werpers <jonatan@werpers.com>
parents: 1446
diff changeset
82 @test apply_stencil(s,v,1) == @SVector[7., 10.]
4f79ab676ebc Fix broken tests using sum(f,...). Needs to be evaluated for performance
Jonatan Werpers <jonatan@werpers.com>
parents: 1446
diff changeset
83 @test apply_stencil_backwards(s,v,3) == @SVector[11., 14.]
4f79ab676ebc Fix broken tests using sum(f,...). Needs to be evaluated for performance
Jonatan Werpers <jonatan@werpers.com>
parents: 1446
diff changeset
84 @test apply_stencil(s,v,1) isa SVector{2, Float64}
4f79ab676ebc Fix broken tests using sum(f,...). Needs to be evaluated for performance
Jonatan Werpers <jonatan@werpers.com>
parents: 1446
diff changeset
85 @test apply_stencil_backwards(s,v,3) isa SVector{2, Float64}
1446
d744b01be520 Add broken tests for applying stencils to vectors of vectors
Jonatan Werpers <jonatan@werpers.com>
parents: 1445
diff changeset
86
d744b01be520 Add broken tests for applying stencils to vectors of vectors
Jonatan Werpers <jonatan@werpers.com>
parents: 1445
diff changeset
87 v = [@SVector[1, 2], @SVector[3, 4], @SVector[5, 6], @SVector[7, 8]]
d744b01be520 Add broken tests for applying stencils to vectors of vectors
Jonatan Werpers <jonatan@werpers.com>
parents: 1445
diff changeset
88 s = Stencil(1.,2., center = 1)
1452
4f79ab676ebc Fix broken tests using sum(f,...). Needs to be evaluated for performance
Jonatan Werpers <jonatan@werpers.com>
parents: 1446
diff changeset
89 @test apply_stencil(s,v,1) == @SVector[7., 10.]
4f79ab676ebc Fix broken tests using sum(f,...). Needs to be evaluated for performance
Jonatan Werpers <jonatan@werpers.com>
parents: 1446
diff changeset
90 @test apply_stencil_backwards(s,v,3) == @SVector[11., 14.]
4f79ab676ebc Fix broken tests using sum(f,...). Needs to be evaluated for performance
Jonatan Werpers <jonatan@werpers.com>
parents: 1446
diff changeset
91 @test apply_stencil(s,v,1) isa SVector{2, Float64}
4f79ab676ebc Fix broken tests using sum(f,...). Needs to be evaluated for performance
Jonatan Werpers <jonatan@werpers.com>
parents: 1446
diff changeset
92 @test apply_stencil_backwards(s,v,3) isa SVector{2, Float64}
1444
c779af4d02e5 Add tests for apply_stencil(::Stencil) and apply_stencil_backwards(::Stencil)
Jonatan Werpers <jonatan@werpers.com>
parents: 1072
diff changeset
93 end
c779af4d02e5 Add tests for apply_stencil(::Stencil) and apply_stencil_backwards(::Stencil)
Jonatan Werpers <jonatan@werpers.com>
parents: 1072
diff changeset
94
904
47f5451cdfc4 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 903
diff changeset
95 @testset "type stability" begin
903
a7f898b1ce1e Fix type instability for stencil application when weights and vector didn't match
Jonatan Werpers <jonatan@werpers.com>
parents: 898
diff changeset
96 s_int = CenteredStencil(1,2,3)
a7f898b1ce1e Fix type instability for stencil application when weights and vector didn't match
Jonatan Werpers <jonatan@werpers.com>
parents: 898
diff changeset
97 s_float = CenteredStencil(1.,2.,3.)
a7f898b1ce1e Fix type instability for stencil application when weights and vector didn't match
Jonatan Werpers <jonatan@werpers.com>
parents: 898
diff changeset
98 v_int = rand(1:10,10);
a7f898b1ce1e Fix type instability for stencil application when weights and vector didn't match
Jonatan Werpers <jonatan@werpers.com>
parents: 898
diff changeset
99 v_float = rand(10);
a7f898b1ce1e Fix type instability for stencil application when weights and vector didn't match
Jonatan Werpers <jonatan@werpers.com>
parents: 898
diff changeset
100
a7f898b1ce1e Fix type instability for stencil application when weights and vector didn't match
Jonatan Werpers <jonatan@werpers.com>
parents: 898
diff changeset
101 @inferred SbpOperators.apply_stencil(s_int, v_int, 2)
a7f898b1ce1e Fix type instability for stencil application when weights and vector didn't match
Jonatan Werpers <jonatan@werpers.com>
parents: 898
diff changeset
102 @inferred SbpOperators.apply_stencil(s_float, v_float, 2)
a7f898b1ce1e Fix type instability for stencil application when weights and vector didn't match
Jonatan Werpers <jonatan@werpers.com>
parents: 898
diff changeset
103 @inferred SbpOperators.apply_stencil(s_int, v_float, 2)
a7f898b1ce1e Fix type instability for stencil application when weights and vector didn't match
Jonatan Werpers <jonatan@werpers.com>
parents: 898
diff changeset
104 @inferred SbpOperators.apply_stencil(s_float, v_int, 2)
904
47f5451cdfc4 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 903
diff changeset
105
905
a3bc90c59e8e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 904
diff changeset
106 @inferred SbpOperators.apply_stencil_backwards(s_int, v_int, 5)
a3bc90c59e8e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 904
diff changeset
107 @inferred SbpOperators.apply_stencil_backwards(s_float, v_float, 5)
a3bc90c59e8e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 904
diff changeset
108 @inferred SbpOperators.apply_stencil_backwards(s_int, v_float, 5)
a3bc90c59e8e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 904
diff changeset
109 @inferred SbpOperators.apply_stencil_backwards(s_float, v_int, 5)
903
a7f898b1ce1e Fix type instability for stencil application when weights and vector didn't match
Jonatan Werpers <jonatan@werpers.com>
parents: 898
diff changeset
110 end
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
111 end
867
313648b01504 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
112
1032
11767fbb29f4 Add padding functions for stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
113 @testset "left_pad" begin
11767fbb29f4 Add padding functions for stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
114 @test SbpOperators.left_pad(Stencil(1,1, center = 1), 2) == Stencil(1,1, center=1)
11767fbb29f4 Add padding functions for stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
115 @test SbpOperators.left_pad(Stencil(1,1, center = 1), 3) == Stencil(0,1,1, center=2)
11767fbb29f4 Add padding functions for stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
116 @test SbpOperators.left_pad(Stencil(2,3, center = 2), 4) == Stencil(0,0,2,3, center=4)
11767fbb29f4 Add padding functions for stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
117
11767fbb29f4 Add padding functions for stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
118 @test SbpOperators.left_pad(Stencil(2.,3., center = 2), 4) == Stencil(0.,0.,2.,3., center=4)
11767fbb29f4 Add padding functions for stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
119 end
11767fbb29f4 Add padding functions for stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
120
11767fbb29f4 Add padding functions for stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
121 @testset "right_pad" begin
11767fbb29f4 Add padding functions for stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
122 @test SbpOperators.right_pad(Stencil(1,1, center = 1), 2) == Stencil(1,1, center=1)
11767fbb29f4 Add padding functions for stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
123 @test SbpOperators.right_pad(Stencil(1,1, center = 1), 3) == Stencil(1,1,0, center=1)
11767fbb29f4 Add padding functions for stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
124 @test SbpOperators.right_pad(Stencil(2,3, center = 2), 4) == Stencil(2,3,0,0, center=2)
11767fbb29f4 Add padding functions for stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
125
11767fbb29f4 Add padding functions for stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
126 @test SbpOperators.right_pad(Stencil(2.,3., center = 2), 4) == Stencil(2.,3.,0.,0., center=2)
11767fbb29f4 Add padding functions for stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
127 end
1072
14cb97284373 Merge default
Jonatan Werpers <jonatan@werpers.com>
parents: 1032 905
diff changeset
128
14cb97284373 Merge default
Jonatan Werpers <jonatan@werpers.com>
parents: 1032 905
diff changeset
129
867
313648b01504 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
130 @testset "NestedStencil" begin
313648b01504 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
131
313648b01504 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
132 @testset "Constructors" begin
868
e37ee63bf9ac Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 867
diff changeset
133 s1 = CenteredStencil(-1, 1, 0)
e37ee63bf9ac Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 867
diff changeset
134 s2 = CenteredStencil(-1, 0, 1)
e37ee63bf9ac Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 867
diff changeset
135 s3 = CenteredStencil( 0,-1, 1)
867
313648b01504 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
136
313648b01504 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
137 ns = NestedStencil(CenteredStencil(s1,s2,s3))
313648b01504 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
138 @test ns isa NestedStencil{Int,3}
313648b01504 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
139
313648b01504 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
140 @test CenteredNestedStencil(s1,s2,s3) == ns
313648b01504 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
141
313648b01504 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
142 @test NestedStencil(s1,s2,s3, center = 2) == ns
313648b01504 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
143 @test NestedStencil(s1,s2,s3, center = 1) == NestedStencil(Stencil(s1,s2,s3, center=1))
313648b01504 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
144
313648b01504 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
145 @test NestedStencil((-1,1,0),(-1,0,1),(0,-1,1), center=2) == ns
879
011d3ff65e9f Add methods for CenteredNestedStencil with tuple input
Jonatan Werpers <jonatan@werpers.com>
parents: 878
diff changeset
146 @test CenteredNestedStencil((-1,1,0),(-1,0,1),(0,-1,1)) == ns
878
0fe46190a81a Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 868
diff changeset
147 @test NestedStencil((-1,1,0),(-1,0,1),(0,-1,1), center=1) == NestedStencil(Stencil(
0fe46190a81a Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 868
diff changeset
148 Stencil(-1, 1, 0; center=1),
0fe46190a81a Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 868
diff changeset
149 Stencil(-1, 0, 1; center=1),
0fe46190a81a Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 868
diff changeset
150 Stencil( 0,-1, 1; center=1);
0fe46190a81a Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 868
diff changeset
151 center=1
0fe46190a81a Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 868
diff changeset
152 ))
867
313648b01504 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
153
313648b01504 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
154 @testset "Error handling" begin
313648b01504 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
155 end
313648b01504 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
156 end
868
e37ee63bf9ac Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 867
diff changeset
157
896
004324d7ed35 Add scale method for nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 894
diff changeset
158 @testset "scale" begin
004324d7ed35 Add scale method for nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 894
diff changeset
159 ns = NestedStencil((-1,1,0),(-1,0,1),(0,-1,1), center=2)
004324d7ed35 Add scale method for nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 894
diff changeset
160 @test SbpOperators.scale(ns, 2) == NestedStencil((-2,2,0),(-2,0,2),(0,-2,2), center=2)
004324d7ed35 Add scale method for nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 894
diff changeset
161 end
004324d7ed35 Add scale method for nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 894
diff changeset
162
894
54e36688dab8 Add functions for converting the element type of nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 887
diff changeset
163 @testset "conversion" begin
54e36688dab8 Add functions for converting the element type of nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 887
diff changeset
164 ns = NestedStencil((-1,1,0),(-1,0,1),(0,-1,1), center=2)
54e36688dab8 Add functions for converting the element type of nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 887
diff changeset
165 @test NestedStencil{Float64}(ns) == NestedStencil((-1.,1.,0.),(-1.,0.,1.),(0.,-1.,1.), center=2)
54e36688dab8 Add functions for converting the element type of nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 887
diff changeset
166 @test NestedStencil{Rational}(ns) == NestedStencil((-1//1,1//1,0//1),(-1//1,0//1,1//1),(0//1,-1//1,1//1), center=2)
54e36688dab8 Add functions for converting the element type of nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 887
diff changeset
167
54e36688dab8 Add functions for converting the element type of nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 887
diff changeset
168 @test convert(NestedStencil{Float64}, ns) == NestedStencil((-1.,1.,0.),(-1.,0.,1.),(0.,-1.,1.), center=2)
54e36688dab8 Add functions for converting the element type of nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 887
diff changeset
169 @test convert(NestedStencil{Rational}, ns) == NestedStencil((-1//1,1//1,0//1),(-1//1,0//1,1//1),(0//1,-1//1,1//1), center=2)
54e36688dab8 Add functions for converting the element type of nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 887
diff changeset
170 end
54e36688dab8 Add functions for converting the element type of nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 887
diff changeset
171
898
cd6d71781137 Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 897
diff changeset
172 @testset "promotion of weights" begin
cd6d71781137 Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 897
diff changeset
173 @test NestedStencil((-1,1,0),(-1.,0.,1.),(0,-1,1), center=2) isa NestedStencil{Float64,3,3}
cd6d71781137 Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 897
diff changeset
174 @test NestedStencil((-1,1,0),(-1,0,1),(0//1,-1,1), center=2) isa NestedStencil{Rational{Int64},3,3}
cd6d71781137 Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 897
diff changeset
175 end
cd6d71781137 Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 897
diff changeset
176
cd6d71781137 Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 897
diff changeset
177 @testset "promotion" begin
cd6d71781137 Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 897
diff changeset
178 promote(
cd6d71781137 Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 897
diff changeset
179 CenteredNestedStencil((-1,1,0),(-1,0,1),(0,-1,1)),
cd6d71781137 Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 897
diff changeset
180 CenteredNestedStencil((-1.,1.,0.),(-1.,0.,1.),(0.,-1.,1.))
cd6d71781137 Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 897
diff changeset
181 ) == (
cd6d71781137 Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 897
diff changeset
182 CenteredNestedStencil((-1.,1.,0.),(-1.,0.,1.),(0.,-1.,1.)),
cd6d71781137 Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 897
diff changeset
183 CenteredNestedStencil((-1.,1.,0.),(-1.,0.,1.),(0.,-1.,1.))
cd6d71781137 Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 897
diff changeset
184 )
cd6d71781137 Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 897
diff changeset
185 end
cd6d71781137 Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 897
diff changeset
186
868
e37ee63bf9ac Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 867
diff changeset
187 @testset "apply" begin
e37ee63bf9ac Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 867
diff changeset
188 c = [ 1, 3, 6, 10, 15, 21, 28, 36, 45, 55]
e37ee63bf9ac Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 867
diff changeset
189 v = [ 2, 3, 5, 7, 11, 13, 17, 19, 23, 29]
e37ee63bf9ac Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 867
diff changeset
190
878
0fe46190a81a Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 868
diff changeset
191 # Centered
868
e37ee63bf9ac Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 867
diff changeset
192 ns = NestedStencil((-1,1,0),(-1,0,1),(0,-2,2), center=2)
e37ee63bf9ac Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 867
diff changeset
193 @test SbpOperators.apply_inner_stencils(ns, c, 4) == Stencil(4,9,10; center=2)
e37ee63bf9ac Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 867
diff changeset
194 @test SbpOperators.apply_inner_stencils_backwards(ns, c, 4) == Stencil(-5,-9,-8; center=2)
e37ee63bf9ac Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 867
diff changeset
195
e37ee63bf9ac Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 867
diff changeset
196 @test SbpOperators.apply_stencil(ns, c, v, 4) == 4*5 + 9*7 + 10*11
e37ee63bf9ac Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 867
diff changeset
197 @test SbpOperators.apply_stencil_backwards(ns, c, v, 4) == -8*5 - 9*7 - 5*11
878
0fe46190a81a Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 868
diff changeset
198
0fe46190a81a Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 868
diff changeset
199 # Non-centered
0fe46190a81a Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 868
diff changeset
200 ns = NestedStencil((-1,1,0),(-1,0,1),(0,-1,1), center=1)
0fe46190a81a Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 868
diff changeset
201 @test SbpOperators.apply_inner_stencils(ns, c, 4) == Stencil(5,11,6; center=1)
0fe46190a81a Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 868
diff changeset
202 @test SbpOperators.apply_inner_stencils_backwards(ns, c, 4) == Stencil(-4,-7,-3; center=1)
0fe46190a81a Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 868
diff changeset
203
0fe46190a81a Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 868
diff changeset
204 @test SbpOperators.apply_stencil(ns, c, v, 4) == 5*7 + 11*11 + 6*13
0fe46190a81a Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 868
diff changeset
205 @test SbpOperators.apply_stencil_backwards(ns, c, v, 4) == -3*3 - 7*5 - 4*7
1445
76a121ee7b80 Add tests for apply with NestedStencils and differing types in the stencils and arrays
Jonatan Werpers <jonatan@werpers.com>
parents: 1444
diff changeset
206
76a121ee7b80 Add tests for apply with NestedStencils and differing types in the stencils and arrays
Jonatan Werpers <jonatan@werpers.com>
parents: 1444
diff changeset
207 # Different types in vector and stencil
76a121ee7b80 Add tests for apply with NestedStencils and differing types in the stencils and arrays
Jonatan Werpers <jonatan@werpers.com>
parents: 1444
diff changeset
208 ns = NestedStencil((-1.,1.,0.),(-1.,0.,1.),(0.,-2.,2.), center=2)
76a121ee7b80 Add tests for apply with NestedStencils and differing types in the stencils and arrays
Jonatan Werpers <jonatan@werpers.com>
parents: 1444
diff changeset
209 @test SbpOperators.apply_inner_stencils(ns, c, 4) isa Stencil{Float64, 3}
76a121ee7b80 Add tests for apply with NestedStencils and differing types in the stencils and arrays
Jonatan Werpers <jonatan@werpers.com>
parents: 1444
diff changeset
210 @test SbpOperators.apply_inner_stencils(ns, c, 4) == Stencil(4.,9.,10.; center=2)
76a121ee7b80 Add tests for apply with NestedStencils and differing types in the stencils and arrays
Jonatan Werpers <jonatan@werpers.com>
parents: 1444
diff changeset
211 @test SbpOperators.apply_inner_stencils_backwards(ns, c, 4) isa Stencil{Float64, 3}
76a121ee7b80 Add tests for apply with NestedStencils and differing types in the stencils and arrays
Jonatan Werpers <jonatan@werpers.com>
parents: 1444
diff changeset
212 @test SbpOperators.apply_inner_stencils_backwards(ns, c, 4) == Stencil(-5.,-9.,-8.; center=2)
76a121ee7b80 Add tests for apply with NestedStencils and differing types in the stencils and arrays
Jonatan Werpers <jonatan@werpers.com>
parents: 1444
diff changeset
213
76a121ee7b80 Add tests for apply with NestedStencils and differing types in the stencils and arrays
Jonatan Werpers <jonatan@werpers.com>
parents: 1444
diff changeset
214 @test SbpOperators.apply_stencil(ns, c, v, 4) isa Float64
76a121ee7b80 Add tests for apply with NestedStencils and differing types in the stencils and arrays
Jonatan Werpers <jonatan@werpers.com>
parents: 1444
diff changeset
215 @test SbpOperators.apply_stencil(ns, c, v, 4) == 193.
76a121ee7b80 Add tests for apply with NestedStencils and differing types in the stencils and arrays
Jonatan Werpers <jonatan@werpers.com>
parents: 1444
diff changeset
216 @test SbpOperators.apply_stencil_backwards(ns, c, v, 4) isa Float64
76a121ee7b80 Add tests for apply with NestedStencils and differing types in the stencils and arrays
Jonatan Werpers <jonatan@werpers.com>
parents: 1444
diff changeset
217 @test SbpOperators.apply_stencil_backwards(ns, c, v, 4) == -158.
1446
d744b01be520 Add broken tests for applying stencils to vectors of vectors
Jonatan Werpers <jonatan@werpers.com>
parents: 1445
diff changeset
218
d744b01be520 Add broken tests for applying stencils to vectors of vectors
Jonatan Werpers <jonatan@werpers.com>
parents: 1445
diff changeset
219 # Arrays of vectors
d744b01be520 Add broken tests for applying stencils to vectors of vectors
Jonatan Werpers <jonatan@werpers.com>
parents: 1445
diff changeset
220 ns = NestedStencil((-1.,1.,0.),(-1.,0.,1.),(0.,-2.,2.), center=2)
d744b01be520 Add broken tests for applying stencils to vectors of vectors
Jonatan Werpers <jonatan@werpers.com>
parents: 1445
diff changeset
221 c = [ 1, 3, 6, 10]
d744b01be520 Add broken tests for applying stencils to vectors of vectors
Jonatan Werpers <jonatan@werpers.com>
parents: 1445
diff changeset
222 v = [@SVector[1, 2], @SVector[3, 4], @SVector[5, 6], @SVector[7, 8]]
1452
4f79ab676ebc Fix broken tests using sum(f,...). Needs to be evaluated for performance
Jonatan Werpers <jonatan@werpers.com>
parents: 1446
diff changeset
223 @test SbpOperators.apply_stencil(ns, c, v, 2) isa SVector{2,Float64}
4f79ab676ebc Fix broken tests using sum(f,...). Needs to be evaluated for performance
Jonatan Werpers <jonatan@werpers.com>
parents: 1446
diff changeset
224 @test SbpOperators.apply_stencil(ns, c, v, 2) == 2v[1] + 5v[2] + 6v[3]
4f79ab676ebc Fix broken tests using sum(f,...). Needs to be evaluated for performance
Jonatan Werpers <jonatan@werpers.com>
parents: 1446
diff changeset
225 @test SbpOperators.apply_stencil_backwards(ns, c, v, 2) isa SVector{2,Float64}
4f79ab676ebc Fix broken tests using sum(f,...). Needs to be evaluated for performance
Jonatan Werpers <jonatan@werpers.com>
parents: 1446
diff changeset
226 @test SbpOperators.apply_stencil_backwards(ns, c, v, 2) == -4v[1] - 5v[2] - 3v[3]
868
e37ee63bf9ac Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 867
diff changeset
227 end
904
47f5451cdfc4 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 903
diff changeset
228
47f5451cdfc4 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 903
diff changeset
229 @testset "type stability" begin
47f5451cdfc4 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 903
diff changeset
230 s_int = CenteredNestedStencil((1,2,3),(1,2,3),(1,2,3))
47f5451cdfc4 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 903
diff changeset
231 s_float = CenteredNestedStencil((1.,2.,3.),(1.,2.,3.),(1.,2.,3.))
47f5451cdfc4 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 903
diff changeset
232
47f5451cdfc4 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 903
diff changeset
233 v_int = rand(1:10,10);
47f5451cdfc4 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 903
diff changeset
234 v_float = rand(10);
47f5451cdfc4 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 903
diff changeset
235
47f5451cdfc4 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 903
diff changeset
236 c_int = rand(1:10,10);
47f5451cdfc4 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 903
diff changeset
237 c_float = rand(10);
47f5451cdfc4 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 903
diff changeset
238
47f5451cdfc4 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 903
diff changeset
239 @inferred SbpOperators.apply_stencil(s_int, c_int, v_int, 2)
47f5451cdfc4 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 903
diff changeset
240 @inferred SbpOperators.apply_stencil(s_float, c_int, v_float, 2)
47f5451cdfc4 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 903
diff changeset
241 @inferred SbpOperators.apply_stencil(s_int, c_int, v_float, 2)
47f5451cdfc4 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 903
diff changeset
242 @inferred SbpOperators.apply_stencil(s_float, c_int, v_int, 2)
47f5451cdfc4 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 903
diff changeset
243
47f5451cdfc4 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 903
diff changeset
244 @inferred SbpOperators.apply_stencil(s_int, c_float, v_int, 2)
47f5451cdfc4 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 903
diff changeset
245 @inferred SbpOperators.apply_stencil(s_float, c_float, v_float, 2)
47f5451cdfc4 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 903
diff changeset
246 @inferred SbpOperators.apply_stencil(s_int, c_float, v_float, 2)
47f5451cdfc4 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 903
diff changeset
247 @inferred SbpOperators.apply_stencil(s_float, c_float, v_int, 2)
47f5451cdfc4 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 903
diff changeset
248
905
a3bc90c59e8e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 904
diff changeset
249 @inferred SbpOperators.apply_stencil_backwards(s_int, c_int, v_int, 2)
a3bc90c59e8e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 904
diff changeset
250 @inferred SbpOperators.apply_stencil_backwards(s_float, c_int, v_float, 2)
a3bc90c59e8e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 904
diff changeset
251 @inferred SbpOperators.apply_stencil_backwards(s_int, c_int, v_float, 2)
a3bc90c59e8e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 904
diff changeset
252 @inferred SbpOperators.apply_stencil_backwards(s_float, c_int, v_int, 2)
a3bc90c59e8e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 904
diff changeset
253
a3bc90c59e8e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 904
diff changeset
254 @inferred SbpOperators.apply_stencil_backwards(s_int, c_float, v_int, 2)
a3bc90c59e8e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 904
diff changeset
255 @inferred SbpOperators.apply_stencil_backwards(s_float, c_float, v_float, 2)
a3bc90c59e8e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 904
diff changeset
256 @inferred SbpOperators.apply_stencil_backwards(s_int, c_float, v_float, 2)
a3bc90c59e8e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 904
diff changeset
257 @inferred SbpOperators.apply_stencil_backwards(s_float, c_float, v_int, 2)
904
47f5451cdfc4 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 903
diff changeset
258 end
867
313648b01504 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
259 end