Mercurial > repos > public > sbplib_julia
annotate benchmark/benchmarks.jl @ 1713:b5690ab5f0b8 refactor/grids/boundary_identifiers_1d
Merge default
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Thu, 05 Sep 2024 08:41:38 +0200 |
| parents | 3714a391545a |
| children | 471a948cd2b2 |
| rev | line source |
|---|---|
| 1171 | 1 using BenchmarkTools |
| 1375 | 2 |
|
1198
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
3 using Sbplib |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
4 using Sbplib.Grids |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
5 using Sbplib.SbpOperators |
|
1314
6ae9b98b15b9
Add missing import of LazyTensors in benchmarks.jl
Jonatan Werpers <jonatan@werpers.com>
parents:
1304
diff
changeset
|
6 using Sbplib.LazyTensors |
| 1171 | 7 |
| 1375 | 8 using LinearAlgebra |
| 9 | |
| 1171 | 10 const SUITE = BenchmarkGroup() |
| 11 | |
|
1198
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
12 |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
13 sz(d) = ntuple(i->100, d) |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
14 ll(d) = ntuple(i->0., d) |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
15 lu(d) = ntuple(i->1., d) |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
16 |
|
1634
29c7a787e67d
Fix bug in benchmarks.jl
Jonatan Werpers <jonatan@werpers.com>
parents:
1529
diff
changeset
|
17 g1 = equidistant_grid(ll(1)[1], lu(1)[1], sz(1)...) |
|
29c7a787e67d
Fix bug in benchmarks.jl
Jonatan Werpers <jonatan@werpers.com>
parents:
1529
diff
changeset
|
18 g2 = equidistant_grid(ll(2), lu(2), sz(2)...) |
|
29c7a787e67d
Fix bug in benchmarks.jl
Jonatan Werpers <jonatan@werpers.com>
parents:
1529
diff
changeset
|
19 g3 = equidistant_grid(ll(3), lu(3), sz(3)...) |
|
1198
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
20 |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
21 v1 = rand(sz(1)...) |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
22 v2 = rand(sz(2)...) |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
23 v3 = rand(sz(3)...) |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
24 |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
25 u1 = rand(sz(1)...) |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
26 u2 = rand(sz(2)...) |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
27 u3 = rand(sz(3)...) |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
28 |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
29 stencil_set = read_stencil_set(joinpath(sbp_operators_path(),"standard_diagonal.toml"); order=4) |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
30 |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
31 SUITE["derivatives"] = BenchmarkGroup() |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
32 |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
33 |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
34 SUITE["derivatives"]["first_derivative"] = BenchmarkGroup() |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
35 |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
36 D₁ = first_derivative(g1,stencil_set) |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
37 SUITE["derivatives"]["first_derivative"]["1D"] = @benchmarkable $u1 .= $D₁*$v1 |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
38 |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
39 Dx = first_derivative(g2,stencil_set,1) |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
40 Dy = first_derivative(g2,stencil_set,2) |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
41 SUITE["derivatives"]["first_derivative"]["2D"] = BenchmarkGroup() |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
42 SUITE["derivatives"]["first_derivative"]["2D"]["x"] = @benchmarkable $u2 .= $Dx*$v2 |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
43 SUITE["derivatives"]["first_derivative"]["2D"]["y"] = @benchmarkable $u2 .= $Dy*$v2 |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
44 |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
45 Dx = first_derivative(g3,stencil_set,1) |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
46 Dy = first_derivative(g3,stencil_set,2) |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
47 Dz = first_derivative(g3,stencil_set,3) |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
48 SUITE["derivatives"]["first_derivative"]["3D"] = BenchmarkGroup() |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
49 SUITE["derivatives"]["first_derivative"]["3D"]["x"] = @benchmarkable $u3 .= $Dx*$v3 |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
50 SUITE["derivatives"]["first_derivative"]["3D"]["y"] = @benchmarkable $u3 .= $Dy*$v3 |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
51 SUITE["derivatives"]["first_derivative"]["3D"]["z"] = @benchmarkable $u3 .= $Dz*$v3 |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
52 |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
53 |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
54 SUITE["derivatives"]["second_derivative"] = BenchmarkGroup() |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
55 |
| 1199 | 56 D₂ = second_derivative(g1,stencil_set) |
| 57 SUITE["derivatives"]["second_derivative"]["1D"] = @benchmarkable $u1 .= $D₂*$v1 | |
|
1198
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
58 |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
59 Dx = second_derivative(g2,stencil_set,1) |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
60 Dy = second_derivative(g2,stencil_set,2) |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
61 SUITE["derivatives"]["second_derivative"]["2D"] = BenchmarkGroup() |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
62 SUITE["derivatives"]["second_derivative"]["2D"]["x"] = @benchmarkable $u2 .= $Dx*$v2 |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
63 SUITE["derivatives"]["second_derivative"]["2D"]["y"] = @benchmarkable $u2 .= $Dy*$v2 |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
64 |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
65 Dx = second_derivative(g3,stencil_set,1) |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
66 Dy = second_derivative(g3,stencil_set,2) |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
67 Dz = second_derivative(g3,stencil_set,3) |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
68 SUITE["derivatives"]["second_derivative"]["3D"] = BenchmarkGroup() |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
69 SUITE["derivatives"]["second_derivative"]["3D"]["x"] = @benchmarkable $u3 .= $Dx*$v3 |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
70 SUITE["derivatives"]["second_derivative"]["3D"]["y"] = @benchmarkable $u3 .= $Dy*$v3 |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
71 SUITE["derivatives"]["second_derivative"]["3D"]["z"] = @benchmarkable $u3 .= $Dz*$v3 |
|
5d4a15d6118f
Add some actual benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1190
diff
changeset
|
72 |
| 1171 | 73 |
| 1375 | 74 SUITE["derivatives"]["second_derivative_variable"] = BenchmarkGroup() |
| 75 | |
| 76 c1 = map(x->sin(x)+2, g1) | |
| 77 D₂ = second_derivative_variable(g1, c1, stencil_set) | |
| 78 SUITE["derivatives"]["second_derivative_variable"]["1D"] = @benchmarkable $u1 .= $D₂*$v1 | |
| 79 | |
| 80 c2 = map(x->sin(x[1] + x[2])+2, g2) | |
| 81 Dx = second_derivative_variable(g2, c2, stencil_set, 1) | |
| 82 Dy = second_derivative_variable(g2, c2, stencil_set, 2) | |
| 83 SUITE["derivatives"]["second_derivative_variable"]["2D"] = BenchmarkGroup() | |
| 84 SUITE["derivatives"]["second_derivative_variable"]["2D"]["x"] = @benchmarkable $u2 .= $Dx*$v2 | |
| 85 SUITE["derivatives"]["second_derivative_variable"]["2D"]["y"] = @benchmarkable $u2 .= $Dy*$v2 | |
| 86 | |
| 87 c3 = map(x->sin(norm(x))+2, g3) | |
| 88 Dx = second_derivative_variable(g3, c3, stencil_set, 1) | |
| 89 Dy = second_derivative_variable(g3, c3, stencil_set, 2) | |
| 90 Dz = second_derivative_variable(g3, c3, stencil_set, 3) | |
| 91 SUITE["derivatives"]["second_derivative_variable"]["3D"] = BenchmarkGroup() | |
| 92 SUITE["derivatives"]["second_derivative_variable"]["3D"]["x"] = @benchmarkable $u3 .= $Dx*$v3 | |
| 93 SUITE["derivatives"]["second_derivative_variable"]["3D"]["y"] = @benchmarkable $u3 .= $Dy*$v3 | |
| 94 SUITE["derivatives"]["second_derivative_variable"]["3D"]["z"] = @benchmarkable $u3 .= $Dz*$v3 | |
| 95 | |
| 96 | |
| 97 | |
| 1199 | 98 |
| 99 SUITE["derivatives"]["addition"] = BenchmarkGroup() | |
| 100 | |
| 101 D₁ = first_derivative(g1,stencil_set) | |
| 102 D₂ = second_derivative(g1,stencil_set) | |
| 103 SUITE["derivatives"]["addition"]["1D"] = BenchmarkGroup() | |
| 104 SUITE["derivatives"]["addition"]["1D"]["apply,add"] = @benchmarkable $u1 .= $D₁*$v1 + $D₂*$v1 | |
| 105 SUITE["derivatives"]["addition"]["1D"]["add,apply"] = @benchmarkable $u1 .= ($D₁ + $D₂)*$v1 | |
| 106 | |
| 107 Dxx = second_derivative(g2,stencil_set,1) | |
| 108 Dyy = second_derivative(g2,stencil_set,2) | |
| 109 SUITE["derivatives"]["addition"]["2D"] = BenchmarkGroup() | |
| 110 SUITE["derivatives"]["addition"]["2D"]["apply,add"] = @benchmarkable $u2 .= $Dxx*$v2 + $Dyy*$v2 | |
| 111 SUITE["derivatives"]["addition"]["2D"]["add,apply"] = @benchmarkable $u2 .= ($Dxx + $Dyy)*$v2 | |
| 112 | |
| 113 Dxx = second_derivative(g3,stencil_set,1) | |
| 114 Dyy = second_derivative(g3,stencil_set,2) | |
| 115 Dzz = second_derivative(g3,stencil_set,3) | |
| 116 SUITE["derivatives"]["addition"]["3D"] = BenchmarkGroup() | |
| 117 SUITE["derivatives"]["addition"]["3D"]["apply,add"] = @benchmarkable $u3 .= $Dxx*$v3 + $Dyy*$v3 + $Dzz*$v3 | |
| 118 SUITE["derivatives"]["addition"]["3D"]["add,apply"] = @benchmarkable $u3 .= ($Dxx + $Dyy + $Dzz)*$v3 | |
| 119 | |
| 120 | |
| 121 SUITE["derivatives"]["composition"] = BenchmarkGroup() | |
| 122 | |
| 123 Dx = first_derivative(g1,stencil_set) | |
| 124 SUITE["derivatives"]["composition"]["1D"] = BenchmarkGroup() | |
| 125 SUITE["derivatives"]["composition"]["1D"]["apply,apply"] = @benchmarkable $u1 .= $Dx*($Dx*$v1) | |
| 126 SUITE["derivatives"]["composition"]["1D"]["compose,apply"] = @benchmarkable $u1 .= ($Dx∘$Dx)*$v1 | |
| 127 | |
| 128 Dx = first_derivative(g2,stencil_set,1) | |
| 129 Dy = first_derivative(g2,stencil_set,2) | |
| 130 SUITE["derivatives"]["composition"]["2D"] = BenchmarkGroup() | |
| 131 SUITE["derivatives"]["composition"]["2D"]["apply,apply"] = @benchmarkable $u2 .= $Dy*($Dx*$v2) | |
| 132 SUITE["derivatives"]["composition"]["2D"]["compose,apply"] = @benchmarkable $u2 .= ($Dy∘$Dx)*$v2 | |
| 133 | |
| 134 Dx = first_derivative(g3,stencil_set,1) | |
| 135 Dy = first_derivative(g3,stencil_set,2) | |
| 136 Dz = first_derivative(g3,stencil_set,3) | |
| 137 SUITE["derivatives"]["composition"]["3D"] = BenchmarkGroup() | |
| 138 SUITE["derivatives"]["composition"]["3D"]["xy"] = BenchmarkGroup() | |
| 139 SUITE["derivatives"]["composition"]["3D"]["xy"]["apply,apply"] = @benchmarkable $u3 .= $Dx*($Dy*$v3) | |
| 140 SUITE["derivatives"]["composition"]["3D"]["xy"]["compose,apply"] = @benchmarkable $u3 .= ($Dx∘$Dy)*$v3 | |
| 141 | |
| 142 SUITE["derivatives"]["composition"]["3D"]["yz"] = BenchmarkGroup() | |
| 143 SUITE["derivatives"]["composition"]["3D"]["yz"]["apply,apply"] = @benchmarkable $u3 .= $Dy*($Dz*$v3) | |
| 144 SUITE["derivatives"]["composition"]["3D"]["yz"]["compose,apply"] = @benchmarkable $u3 .= ($Dy∘$Dz)*$v3 | |
| 145 | |
| 146 SUITE["derivatives"]["composition"]["3D"]["xz"] = BenchmarkGroup() | |
| 147 SUITE["derivatives"]["composition"]["3D"]["xz"]["apply,apply"] = @benchmarkable $u3 .= $Dx*($Dz*$v3) | |
| 148 SUITE["derivatives"]["composition"]["3D"]["xz"]["compose,apply"] = @benchmarkable $u3 .= ($Dx∘$Dz)*$v3 | |
| 149 | |
| 150 SUITE["derivatives"]["composition"]["3D"]["xx"] = BenchmarkGroup() | |
| 151 SUITE["derivatives"]["composition"]["3D"]["xx"]["apply,apply"] = @benchmarkable $u3 .= $Dx*($Dx*$v3) | |
| 152 SUITE["derivatives"]["composition"]["3D"]["xx"]["compose,apply"] = @benchmarkable $u3 .= ($Dx∘$Dx)*$v3 | |
| 153 | |
| 154 SUITE["derivatives"]["composition"]["3D"]["yy"] = BenchmarkGroup() | |
| 155 SUITE["derivatives"]["composition"]["3D"]["yy"]["apply,apply"] = @benchmarkable $u3 .= $Dy*($Dy*$v3) | |
| 156 SUITE["derivatives"]["composition"]["3D"]["yy"]["compose,apply"] = @benchmarkable $u3 .= ($Dy∘$Dy)*$v3 | |
| 157 | |
| 158 SUITE["derivatives"]["composition"]["3D"]["zz"] = BenchmarkGroup() | |
| 159 SUITE["derivatives"]["composition"]["3D"]["zz"]["apply,apply"] = @benchmarkable $u3 .= $Dz*($Dz*$v3) | |
| 160 SUITE["derivatives"]["composition"]["3D"]["zz"]["compose,apply"] = @benchmarkable $u3 .= ($Dz∘$Dz)*$v3 | |
| 161 | |
|
1211
58a7aa15ce68
Add boundary term benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1199
diff
changeset
|
162 |
|
58a7aa15ce68
Add boundary term benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1199
diff
changeset
|
163 SUITE["boundary_terms"] = BenchmarkGroup() |
|
58a7aa15ce68
Add boundary term benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1199
diff
changeset
|
164 |
|
58a7aa15ce68
Add boundary term benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1199
diff
changeset
|
165 H = inner_product(g2, stencil_set) |
|
58a7aa15ce68
Add boundary term benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1199
diff
changeset
|
166 H⁻¹ = inverse_inner_product(g2, stencil_set) |
|
58a7aa15ce68
Add boundary term benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1199
diff
changeset
|
167 Dxx = second_derivative(g2, stencil_set, 1) |
|
58a7aa15ce68
Add boundary term benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1199
diff
changeset
|
168 Dyy = second_derivative(g2, stencil_set, 2) |
|
58a7aa15ce68
Add boundary term benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1199
diff
changeset
|
169 |
|
1672
3714a391545a
Make the boundary identifiers for EquidistantGrid subtype BoundaryIdentifer
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents:
1635
diff
changeset
|
170 e₁ₗ = boundary_restriction(g2, stencil_set, CartesianBoundary{1,LowerBoundary}()) |
|
3714a391545a
Make the boundary identifiers for EquidistantGrid subtype BoundaryIdentifer
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents:
1635
diff
changeset
|
171 e₁ᵤ = boundary_restriction(g2, stencil_set, CartesianBoundary{1,UpperBoundary}()) |
|
3714a391545a
Make the boundary identifiers for EquidistantGrid subtype BoundaryIdentifer
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents:
1635
diff
changeset
|
172 e₂ₗ = boundary_restriction(g2, stencil_set, CartesianBoundary{2,LowerBoundary}()) |
|
3714a391545a
Make the boundary identifiers for EquidistantGrid subtype BoundaryIdentifer
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents:
1635
diff
changeset
|
173 e₂ᵤ = boundary_restriction(g2, stencil_set, CartesianBoundary{2,UpperBoundary}()) |
|
1211
58a7aa15ce68
Add boundary term benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1199
diff
changeset
|
174 |
|
1672
3714a391545a
Make the boundary identifiers for EquidistantGrid subtype BoundaryIdentifer
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents:
1635
diff
changeset
|
175 d₁ₗ = normal_derivative(g2, stencil_set, CartesianBoundary{1,LowerBoundary}()) |
|
3714a391545a
Make the boundary identifiers for EquidistantGrid subtype BoundaryIdentifer
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents:
1635
diff
changeset
|
176 d₁ᵤ = normal_derivative(g2, stencil_set, CartesianBoundary{1,UpperBoundary}()) |
|
3714a391545a
Make the boundary identifiers for EquidistantGrid subtype BoundaryIdentifer
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents:
1635
diff
changeset
|
177 d₂ₗ = normal_derivative(g2, stencil_set, CartesianBoundary{2,LowerBoundary}()) |
|
3714a391545a
Make the boundary identifiers for EquidistantGrid subtype BoundaryIdentifer
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents:
1635
diff
changeset
|
178 d₂ᵤ = normal_derivative(g2, stencil_set, CartesianBoundary{2,UpperBoundary}()) |
|
1211
58a7aa15ce68
Add boundary term benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1199
diff
changeset
|
179 |
|
1672
3714a391545a
Make the boundary identifiers for EquidistantGrid subtype BoundaryIdentifer
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents:
1635
diff
changeset
|
180 H₁ₗ = inner_product(boundary_grid(g2, CartesianBoundary{1,LowerBoundary}()), stencil_set) |
|
3714a391545a
Make the boundary identifiers for EquidistantGrid subtype BoundaryIdentifer
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents:
1635
diff
changeset
|
181 H₁ᵤ = inner_product(boundary_grid(g2, CartesianBoundary{1,UpperBoundary}()), stencil_set) |
|
3714a391545a
Make the boundary identifiers for EquidistantGrid subtype BoundaryIdentifer
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents:
1635
diff
changeset
|
182 H₂ₗ = inner_product(boundary_grid(g2, CartesianBoundary{2,LowerBoundary}()), stencil_set) |
|
3714a391545a
Make the boundary identifiers for EquidistantGrid subtype BoundaryIdentifer
Vidar Stiernström <vidar.stiernstrom@gmail.com>
parents:
1635
diff
changeset
|
183 H₂ᵤ = inner_product(boundary_grid(g2, CartesianBoundary{2,UpperBoundary}()), stencil_set) |
|
1211
58a7aa15ce68
Add boundary term benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1199
diff
changeset
|
184 |
|
58a7aa15ce68
Add boundary term benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1199
diff
changeset
|
185 SUITE["boundary_terms"]["pre_composition"] = @benchmarkable $u2 .= $(H⁻¹∘e₁ₗ'∘H₁ₗ∘d₁ₗ)*$v2 |
|
58a7aa15ce68
Add boundary term benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1199
diff
changeset
|
186 SUITE["boundary_terms"]["composition"] = @benchmarkable $u2 .= ($H⁻¹∘$e₁ₗ'∘$H₁ₗ∘$d₁ₗ)*$v2 |
|
58a7aa15ce68
Add boundary term benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1199
diff
changeset
|
187 SUITE["boundary_terms"]["application"] = @benchmarkable $u2 .= $H⁻¹*$e₁ₗ'*$H₁ₗ* $d₁ₗ*$v2 |
|
58a7aa15ce68
Add boundary term benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1199
diff
changeset
|
188 # An investigation of these allocations can be found in the branch `allocation_test` |
|
58a7aa15ce68
Add boundary term benchmarks
Jonatan Werpers <jonatan@werpers.com>
parents:
1199
diff
changeset
|
189 |
| 1199 | 190 #TODO: Reorg with dimension as first level? To reduce operator creation? |
| 191 | |
|
1304
37d615522c2f
Add a simple benchmark demonstrating allocations with compositions
Jonatan Werpers <jonatan@werpers.com>
parents:
1211
diff
changeset
|
192 |
|
37d615522c2f
Add a simple benchmark demonstrating allocations with compositions
Jonatan Werpers <jonatan@werpers.com>
parents:
1211
diff
changeset
|
193 |
|
37d615522c2f
Add a simple benchmark demonstrating allocations with compositions
Jonatan Werpers <jonatan@werpers.com>
parents:
1211
diff
changeset
|
194 SUITE["lazy_tensors"] = BenchmarkGroup() |
|
37d615522c2f
Add a simple benchmark demonstrating allocations with compositions
Jonatan Werpers <jonatan@werpers.com>
parents:
1211
diff
changeset
|
195 |
|
37d615522c2f
Add a simple benchmark demonstrating allocations with compositions
Jonatan Werpers <jonatan@werpers.com>
parents:
1211
diff
changeset
|
196 SUITE["lazy_tensors"]["compositions"] = BenchmarkGroup() |
|
37d615522c2f
Add a simple benchmark demonstrating allocations with compositions
Jonatan Werpers <jonatan@werpers.com>
parents:
1211
diff
changeset
|
197 s = ScalingTensor(1.,(10,)) |
|
37d615522c2f
Add a simple benchmark demonstrating allocations with compositions
Jonatan Werpers <jonatan@werpers.com>
parents:
1211
diff
changeset
|
198 u = rand(10) |
|
37d615522c2f
Add a simple benchmark demonstrating allocations with compositions
Jonatan Werpers <jonatan@werpers.com>
parents:
1211
diff
changeset
|
199 v = similar(u) |
|
37d615522c2f
Add a simple benchmark demonstrating allocations with compositions
Jonatan Werpers <jonatan@werpers.com>
parents:
1211
diff
changeset
|
200 s3 = s∘s∘s |
|
37d615522c2f
Add a simple benchmark demonstrating allocations with compositions
Jonatan Werpers <jonatan@werpers.com>
parents:
1211
diff
changeset
|
201 s4 = s∘s∘s∘s |
|
37d615522c2f
Add a simple benchmark demonstrating allocations with compositions
Jonatan Werpers <jonatan@werpers.com>
parents:
1211
diff
changeset
|
202 SUITE["lazy_tensors"]["compositions"]["s∘s∘s"] = @benchmarkable $v .= $s3*$u |
|
37d615522c2f
Add a simple benchmark demonstrating allocations with compositions
Jonatan Werpers <jonatan@werpers.com>
parents:
1211
diff
changeset
|
203 SUITE["lazy_tensors"]["compositions"]["s∘s∘s∘s"] = @benchmarkable $v .= $s4*$u |
|
37d615522c2f
Add a simple benchmark demonstrating allocations with compositions
Jonatan Werpers <jonatan@werpers.com>
parents:
1211
diff
changeset
|
204 |
|
37d615522c2f
Add a simple benchmark demonstrating allocations with compositions
Jonatan Werpers <jonatan@werpers.com>
parents:
1211
diff
changeset
|
205 |
| 1171 | 206 SUITE |