Mercurial > repos > public > sbplib_julia

annotate test/SbpOperators/stencil_test.jl @ 1064:a601427023e3 feature/nested_stencils

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add type stability tests for nested stencil application (grafted from 47f5451cdfc4ad3d070c6f017c88644746e17669)
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 15 Feb 2022 07:54:52 +0100
parents f98893154e22
children a3bc90c59e8e
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
1 using Test
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
2 using Sbplib.SbpOperators
732
6114274447f5 Add missing using and imports in test files
Jonatan Werpers <jonatan@werpers.com>
parents: 728
diff changeset
3 import Sbplib.SbpOperators.Stencil
1052
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
4 import Sbplib.SbpOperators.NestedStencil
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
5 import Sbplib.SbpOperators.scale
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
6
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
7 @testset "Stencil" begin
1060
ff0e819f2075 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 1059
diff changeset
8 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
9 @test s isa Stencil{Float64, 5}
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
10
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
11 @test eltype(s) == Float64
1057
0728e84af0dc Add length method for stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1055
diff changeset
12
0728e84af0dc Add length method for stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1055
diff changeset
13 @test length(s) == 5
1060
ff0e819f2075 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 1059
diff changeset
14 @test length(Stencil(-1:2, (1,2,3,4))) == 4
1057
0728e84af0dc Add length method for stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1055
diff changeset
15
1060
ff0e819f2075 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 1059
diff changeset
16 @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
17
1060
ff0e819f2075 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 1059
diff changeset
18 @test Stencil(1,2,3,4; center=1) == Stencil(0:3,(1,2,3,4))
ff0e819f2075 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 1059
diff changeset
19 @test Stencil(1,2,3,4; center=2) == Stencil(-1:2,(1,2,3,4))
ff0e819f2075 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 1059
diff changeset
20 @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
21
1060
ff0e819f2075 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 1059
diff changeset
22 @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
23 @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
24
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
25 # 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
26 @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
27 @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
28 @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
29 @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
30
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 @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
32 @test convert(Stencil{Float64}, Stencil(1,2,3,4,5; center=2)) == Stencil(1.,2.,3.,4.,5.; center=2)
1060
ff0e819f2075 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 1059
diff changeset
33 @test convert(Stencil{Float64,5}, CenteredStencil(1,2,3,4,5)) == CenteredStencil(1.,2.,3.,4.,5.)
ff0e819f2075 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 1059
diff changeset
34 @test convert(Stencil{Int,5}, Stencil(1.,2.,3.,4.,5.; center=2)) == Stencil(1,2,3,4,5; center=2)
ff0e819f2075 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 1059
diff changeset
35 @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)
ff0e819f2075 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 1059
diff changeset
36 end
ff0e819f2075 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 1059
diff changeset
37
ff0e819f2075 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 1059
diff changeset
38 @testset "promotion of weights" begin
ff0e819f2075 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 1059
diff changeset
39 @test Stencil(1.,2; center = 1) isa Stencil{Float64, 2}
ff0e819f2075 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 1059
diff changeset
40 @test Stencil(1,2//2; center = 1) isa Stencil{Rational{Int64}, 2}
ff0e819f2075 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 1059
diff changeset
41 end
ff0e819f2075 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 1059
diff changeset
42
ff0e819f2075 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 1059
diff changeset
43 @testset "promotion" begin
ff0e819f2075 Refactor code for regular stencils to use fewer type parameters and allow promotion
Jonatan Werpers <jonatan@werpers.com>
parents: 1059
diff changeset
44 @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
45 end
1063
f98893154e22 Fix type instability for stencil application when weights and vector didn't match
Jonatan Werpers <jonatan@werpers.com>
parents: 1061
diff changeset
46
1064
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
47 @testset "type stability" begin
1063
f98893154e22 Fix type instability for stencil application when weights and vector didn't match
Jonatan Werpers <jonatan@werpers.com>
parents: 1061
diff changeset
48 s_int = CenteredStencil(1,2,3)
f98893154e22 Fix type instability for stencil application when weights and vector didn't match
Jonatan Werpers <jonatan@werpers.com>
parents: 1061
diff changeset
49 s_float = CenteredStencil(1.,2.,3.)
f98893154e22 Fix type instability for stencil application when weights and vector didn't match
Jonatan Werpers <jonatan@werpers.com>
parents: 1061
diff changeset
50 v_int = rand(1:10,10);
f98893154e22 Fix type instability for stencil application when weights and vector didn't match
Jonatan Werpers <jonatan@werpers.com>
parents: 1061
diff changeset
51 v_float = rand(10);
f98893154e22 Fix type instability for stencil application when weights and vector didn't match
Jonatan Werpers <jonatan@werpers.com>
parents: 1061
diff changeset
52
f98893154e22 Fix type instability for stencil application when weights and vector didn't match
Jonatan Werpers <jonatan@werpers.com>
parents: 1061
diff changeset
53 @inferred SbpOperators.apply_stencil(s_int, v_int, 2)
f98893154e22 Fix type instability for stencil application when weights and vector didn't match
Jonatan Werpers <jonatan@werpers.com>
parents: 1061
diff changeset
54 @inferred SbpOperators.apply_stencil(s_float, v_float, 2)
f98893154e22 Fix type instability for stencil application when weights and vector didn't match
Jonatan Werpers <jonatan@werpers.com>
parents: 1061
diff changeset
55 @inferred SbpOperators.apply_stencil(s_int, v_float, 2)
f98893154e22 Fix type instability for stencil application when weights and vector didn't match
Jonatan Werpers <jonatan@werpers.com>
parents: 1061
diff changeset
56 @inferred SbpOperators.apply_stencil(s_float, v_int, 2)
1064
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
57
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
58 # TODO: apply backwards
1063
f98893154e22 Fix type instability for stencil application when weights and vector didn't match
Jonatan Werpers <jonatan@werpers.com>
parents: 1061
diff changeset
59 end
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
60 end
1052
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
61
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
62 @testset "NestedStencil" begin
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
63
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
64 @testset "Constructors" begin
1053
cb0a5f41be02 Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 1052
diff changeset
65 s1 = CenteredStencil(-1, 1, 0)
cb0a5f41be02 Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 1052
diff changeset
66 s2 = CenteredStencil(-1, 0, 1)
cb0a5f41be02 Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 1052
diff changeset
67 s3 = CenteredStencil( 0,-1, 1)
1052
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
68
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
69 ns = NestedStencil(CenteredStencil(s1,s2,s3))
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
70 @test ns isa NestedStencil{Int,3}
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
71
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
72 @test CenteredNestedStencil(s1,s2,s3) == ns
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
73
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
74 @test NestedStencil(s1,s2,s3, center = 2) == ns
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
75 @test NestedStencil(s1,s2,s3, center = 1) == NestedStencil(Stencil(s1,s2,s3, center=1))
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
76
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
77 @test NestedStencil((-1,1,0),(-1,0,1),(0,-1,1), center=2) == ns
1055
df498ce0cf52 Add methods for CenteredNestedStencil with tuple input
Jonatan Werpers <jonatan@werpers.com>
parents: 1054
diff changeset
78 @test CenteredNestedStencil((-1,1,0),(-1,0,1),(0,-1,1)) == ns
1054
ce16519aa7ff Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1053
diff changeset
79 @test NestedStencil((-1,1,0),(-1,0,1),(0,-1,1), center=1) == NestedStencil(Stencil(
ce16519aa7ff Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1053
diff changeset
80 Stencil(-1, 1, 0; center=1),
ce16519aa7ff Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1053
diff changeset
81 Stencil(-1, 0, 1; center=1),
ce16519aa7ff Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1053
diff changeset
82 Stencil( 0,-1, 1; center=1);
ce16519aa7ff Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1053
diff changeset
83 center=1
ce16519aa7ff Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1053
diff changeset
84 ))
1052
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
85
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
86 @testset "Error handling" begin
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
87 end
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
88 end
1053
cb0a5f41be02 Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 1052
diff changeset
89
1059
4d06642174ec Add scale method for nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1058
diff changeset
90 @testset "scale" begin
4d06642174ec Add scale method for nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1058
diff changeset
91 ns = NestedStencil((-1,1,0),(-1,0,1),(0,-1,1), center=2)
4d06642174ec Add scale method for nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1058
diff changeset
92 @test SbpOperators.scale(ns, 2) == NestedStencil((-2,2,0),(-2,0,2),(0,-2,2), center=2)
4d06642174ec Add scale method for nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1058
diff changeset
93 end
4d06642174ec Add scale method for nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1058
diff changeset
94
1058
c4f71d6f2d63 Add functions for converting the element type of nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1057
diff changeset
95 @testset "conversion" begin
c4f71d6f2d63 Add functions for converting the element type of nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1057
diff changeset
96 ns = NestedStencil((-1,1,0),(-1,0,1),(0,-1,1), center=2)
c4f71d6f2d63 Add functions for converting the element type of nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1057
diff changeset
97 @test NestedStencil{Float64}(ns) == NestedStencil((-1.,1.,0.),(-1.,0.,1.),(0.,-1.,1.), center=2)
c4f71d6f2d63 Add functions for converting the element type of nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1057
diff changeset
98 @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)
c4f71d6f2d63 Add functions for converting the element type of nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1057
diff changeset
99
c4f71d6f2d63 Add functions for converting the element type of nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1057
diff changeset
100 @test convert(NestedStencil{Float64}, ns) == NestedStencil((-1.,1.,0.),(-1.,0.,1.),(0.,-1.,1.), center=2)
c4f71d6f2d63 Add functions for converting the element type of nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1057
diff changeset
101 @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)
c4f71d6f2d63 Add functions for converting the element type of nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1057
diff changeset
102 end
c4f71d6f2d63 Add functions for converting the element type of nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1057
diff changeset
103
1061
d1ddcdc098bf Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 1060
diff changeset
104 @testset "promotion of weights" begin
d1ddcdc098bf Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 1060
diff changeset
105 @test NestedStencil((-1,1,0),(-1.,0.,1.),(0,-1,1), center=2) isa NestedStencil{Float64,3,3}
d1ddcdc098bf Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 1060
diff changeset
106 @test NestedStencil((-1,1,0),(-1,0,1),(0//1,-1,1), center=2) isa NestedStencil{Rational{Int64},3,3}
d1ddcdc098bf Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 1060
diff changeset
107 end
d1ddcdc098bf Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 1060
diff changeset
108
d1ddcdc098bf Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 1060
diff changeset
109 @testset "promotion" begin
d1ddcdc098bf Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 1060
diff changeset
110 promote(
d1ddcdc098bf Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 1060
diff changeset
111 CenteredNestedStencil((-1,1,0),(-1,0,1),(0,-1,1)),
d1ddcdc098bf Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 1060
diff changeset
112 CenteredNestedStencil((-1.,1.,0.),(-1.,0.,1.),(0.,-1.,1.))
d1ddcdc098bf Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 1060
diff changeset
113 ) == (
d1ddcdc098bf Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 1060
diff changeset
114 CenteredNestedStencil((-1.,1.,0.),(-1.,0.,1.),(0.,-1.,1.)),
d1ddcdc098bf Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 1060
diff changeset
115 CenteredNestedStencil((-1.,1.,0.),(-1.,0.,1.),(0.,-1.,1.))
d1ddcdc098bf Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 1060
diff changeset
116 )
d1ddcdc098bf Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 1060
diff changeset
117 end
d1ddcdc098bf Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 1060
diff changeset
118
1053
cb0a5f41be02 Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 1052
diff changeset
119 @testset "apply" begin
cb0a5f41be02 Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 1052
diff changeset
120 c = [ 1, 3, 6, 10, 15, 21, 28, 36, 45, 55]
cb0a5f41be02 Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 1052
diff changeset
121 v = [ 2, 3, 5, 7, 11, 13, 17, 19, 23, 29]
cb0a5f41be02 Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 1052
diff changeset
122
1054
ce16519aa7ff Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1053
diff changeset
123 # Centered
1053
cb0a5f41be02 Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 1052
diff changeset
124 ns = NestedStencil((-1,1,0),(-1,0,1),(0,-2,2), center=2)
cb0a5f41be02 Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 1052
diff changeset
125 @test SbpOperators.apply_inner_stencils(ns, c, 4) == Stencil(4,9,10; center=2)
cb0a5f41be02 Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 1052
diff changeset
126 @test SbpOperators.apply_inner_stencils_backwards(ns, c, 4) == Stencil(-5,-9,-8; center=2)
cb0a5f41be02 Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 1052
diff changeset
127
cb0a5f41be02 Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 1052
diff changeset
128 @test SbpOperators.apply_stencil(ns, c, v, 4) == 4*5 + 9*7 + 10*11
cb0a5f41be02 Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 1052
diff changeset
129 @test SbpOperators.apply_stencil_backwards(ns, c, v, 4) == -8*5 - 9*7 - 5*11
1054
ce16519aa7ff Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1053
diff changeset
130
ce16519aa7ff Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1053
diff changeset
131 # Non-centered
ce16519aa7ff Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1053
diff changeset
132 ns = NestedStencil((-1,1,0),(-1,0,1),(0,-1,1), center=1)
ce16519aa7ff Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1053
diff changeset
133 @test SbpOperators.apply_inner_stencils(ns, c, 4) == Stencil(5,11,6; center=1)
ce16519aa7ff Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1053
diff changeset
134 @test SbpOperators.apply_inner_stencils_backwards(ns, c, 4) == Stencil(-4,-7,-3; center=1)
ce16519aa7ff Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1053
diff changeset
135
ce16519aa7ff Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1053
diff changeset
136 @test SbpOperators.apply_stencil(ns, c, v, 4) == 5*7 + 11*11 + 6*13
ce16519aa7ff Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1053
diff changeset
137 @test SbpOperators.apply_stencil_backwards(ns, c, v, 4) == -3*3 - 7*5 - 4*7
1053
cb0a5f41be02 Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 1052
diff changeset
138 end
1064
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
139
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
140 @testset "type stability" begin
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
141 s_int = CenteredNestedStencil((1,2,3),(1,2,3),(1,2,3))
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
142 s_float = CenteredNestedStencil((1.,2.,3.),(1.,2.,3.),(1.,2.,3.))
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
143
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
144 v_int = rand(1:10,10);
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
145 v_float = rand(10);
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
146
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
147 c_int = rand(1:10,10);
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
148 c_float = rand(10);
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
149
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
150 @inferred SbpOperators.apply_stencil(s_int, c_int, v_int, 2)
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
151 @inferred SbpOperators.apply_stencil(s_float, c_int, v_float, 2)
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
152 @inferred SbpOperators.apply_stencil(s_int, c_int, v_float, 2)
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
153 @inferred SbpOperators.apply_stencil(s_float, c_int, v_int, 2)
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
154
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
155 @inferred SbpOperators.apply_stencil(s_int, c_float, v_int, 2)
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
156 @inferred SbpOperators.apply_stencil(s_float, c_float, v_float, 2)
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
157 @inferred SbpOperators.apply_stencil(s_int, c_float, v_float, 2)
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
158 @inferred SbpOperators.apply_stencil(s_float, c_float, v_int, 2)
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
159
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
160 # TODO: apply backwards
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
161 end
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
162
1052
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
163 end