Mercurial > repos > public > sbplib_julia

annotate test/SbpOperators/stencil_test.jl @ 1889:c1519a9b4903 allocation_testing

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Close branch.
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 28 Jan 2025 22:07:06 +0100
parents fdd594b2a15e
children 14cb97284373
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
1065
fdd594b2a15e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1064
diff changeset
58 @inferred SbpOperators.apply_stencil_backwards(s_int, v_int, 5)
fdd594b2a15e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1064
diff changeset
59 @inferred SbpOperators.apply_stencil_backwards(s_float, v_float, 5)
fdd594b2a15e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1064
diff changeset
60 @inferred SbpOperators.apply_stencil_backwards(s_int, v_float, 5)
fdd594b2a15e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1064
diff changeset
61 @inferred SbpOperators.apply_stencil_backwards(s_float, v_int, 5)
1063
f98893154e22 Fix type instability for stencil application when weights and vector didn't match
Jonatan Werpers <jonatan@werpers.com>
parents: 1061
diff changeset
62 end
728
45966c77cb20 Split tests for SbpOperators over several files
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
63 end
1052
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
64
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
65 @testset "NestedStencil" begin
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
66
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
67 @testset "Constructors" begin
1053
cb0a5f41be02 Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 1052
diff changeset
68 s1 = CenteredStencil(-1, 1, 0)
cb0a5f41be02 Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 1052
diff changeset
69 s2 = CenteredStencil(-1, 0, 1)
cb0a5f41be02 Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 1052
diff changeset
70 s3 = CenteredStencil( 0,-1, 1)
1052
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 ns = NestedStencil(CenteredStencil(s1,s2,s3))
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
73 @test ns isa NestedStencil{Int,3}
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
74
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
75 @test CenteredNestedStencil(s1,s2,s3) == ns
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(s1,s2,s3, center = 2) == ns
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
78 @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
79
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
80 @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
81 @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
82 @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
83 Stencil(-1, 1, 0; center=1),
ce16519aa7ff Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1053
diff changeset
84 Stencil(-1, 0, 1; center=1),
ce16519aa7ff Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1053
diff changeset
85 Stencil( 0,-1, 1; center=1);
ce16519aa7ff Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1053
diff changeset
86 center=1
ce16519aa7ff Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1053
diff changeset
87 ))
1052
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
88
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
89 @testset "Error handling" begin
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
90 end
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
91 end
1053
cb0a5f41be02 Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 1052
diff changeset
92
1059
4d06642174ec Add scale method for nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1058
diff changeset
93 @testset "scale" begin
4d06642174ec Add scale method for nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1058
diff changeset
94 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
95 @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
96 end
4d06642174ec Add scale method for nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1058
diff changeset
97
1058
c4f71d6f2d63 Add functions for converting the element type of nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1057
diff changeset
98 @testset "conversion" begin
c4f71d6f2d63 Add functions for converting the element type of nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1057
diff changeset
99 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
100 @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
101 @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
102
c4f71d6f2d63 Add functions for converting the element type of nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1057
diff changeset
103 @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
104 @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
105 end
c4f71d6f2d63 Add functions for converting the element type of nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1057
diff changeset
106
1061
d1ddcdc098bf Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 1060
diff changeset
107 @testset "promotion of weights" begin
d1ddcdc098bf Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 1060
diff changeset
108 @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
109 @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
110 end
d1ddcdc098bf Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 1060
diff changeset
111
d1ddcdc098bf Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 1060
diff changeset
112 @testset "promotion" begin
d1ddcdc098bf Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 1060
diff changeset
113 promote(
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 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
118 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
119 )
d1ddcdc098bf Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 1060
diff changeset
120 end
d1ddcdc098bf Add promotion functionality for NestedTuples
Jonatan Werpers <jonatan@werpers.com>
parents: 1060
diff changeset
121
1053
cb0a5f41be02 Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 1052
diff changeset
122 @testset "apply" begin
cb0a5f41be02 Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 1052
diff changeset
123 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
124 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
125
1054
ce16519aa7ff Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1053
diff changeset
126 # Centered
1053
cb0a5f41be02 Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 1052
diff changeset
127 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
128 @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
129 @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
130
cb0a5f41be02 Add some tests and fix some bugs
Jonatan Werpers <jonatan@werpers.com>
parents: 1052
diff changeset
131 @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
132 @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
133
ce16519aa7ff Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1053
diff changeset
134 # Non-centered
ce16519aa7ff Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1053
diff changeset
135 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
136 @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
137 @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
138
ce16519aa7ff Tests for non centered nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 1053
diff changeset
139 @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
140 @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
141 end
1064
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
142
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
143 @testset "type stability" begin
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
144 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
145 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
146
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
147 v_int = rand(1:10,10);
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
148 v_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 c_int = rand(1:10,10);
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
151 c_float = rand(10);
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
152
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
153 @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
154 @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
155 @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
156 @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
157
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
158 @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
159 @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
160 @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
161 @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
162
1065
fdd594b2a15e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1064
diff changeset
163 @inferred SbpOperators.apply_stencil_backwards(s_int, c_int, v_int, 2)
fdd594b2a15e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1064
diff changeset
164 @inferred SbpOperators.apply_stencil_backwards(s_float, c_int, v_float, 2)
fdd594b2a15e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1064
diff changeset
165 @inferred SbpOperators.apply_stencil_backwards(s_int, c_int, v_float, 2)
fdd594b2a15e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1064
diff changeset
166 @inferred SbpOperators.apply_stencil_backwards(s_float, c_int, v_int, 2)
fdd594b2a15e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1064
diff changeset
167
fdd594b2a15e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1064
diff changeset
168 @inferred SbpOperators.apply_stencil_backwards(s_int, c_float, v_int, 2)
fdd594b2a15e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1064
diff changeset
169 @inferred SbpOperators.apply_stencil_backwards(s_float, c_float, v_float, 2)
fdd594b2a15e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1064
diff changeset
170 @inferred SbpOperators.apply_stencil_backwards(s_int, c_float, v_float, 2)
fdd594b2a15e Fix type instability in apply_stencil_backwards and add tests
Jonatan Werpers <jonatan@werpers.com>
parents: 1064
diff changeset
171 @inferred SbpOperators.apply_stencil_backwards(s_float, c_float, v_int, 2)
1064
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
172 end
a601427023e3 Add type stability tests for nested stencil application
Jonatan Werpers <jonatan@werpers.com>
parents: 1063
diff changeset
173
1052
ca718fd4e816 Start implementing nested stencils
Jonatan Werpers <jonatan@werpers.com>
parents: 826
diff changeset
174 end