30 |
from esys.finley.finleybench import * |
from esys.finley.finleybench import * |
31 |
from esys.escript.benchmark import BenchmarkSuite,Benchmark |
from esys.escript.benchmark import BenchmarkSuite,Benchmark |
32 |
|
|
33 |
thlist=[1,2,4,8,16] |
# thlist=[1,2,4,8,16] |
34 |
# thlist=[1,2,4,8,16,32] |
thlist=[1,2,4,8,16,32] |
35 |
# thlist=[1,2,4,8,16,32,64,128] |
# thlist=[1,2,4,8,16,32,64,128] |
36 |
# thlist=[1,2,4,8,16,32,64,128] |
# thlist=[1,2,4,8,16,32,64,128] |
37 |
show=False |
show=True |
38 |
ff=FinleyFilter() |
ff=FinleyFilter() |
39 |
|
|
40 |
opt1=FinleyOptions(solver_method=LinearPDE.PCG,preconditioner=LinearPDE.JACOBI,verbose=show) |
opt1=FinleyOptions(solver_method=LinearPDE.PCG,preconditioner=LinearPDE.JACOBI,verbose=show) |
41 |
opt2=FinleyOptions(solver_method=LinearPDE.PCG,preconditioner=LinearPDE.ILU0,verbose=show) |
opt2=FinleyOptions(solver_method=LinearPDE.PCG,preconditioner=LinearPDE.ILU0,verbose=show) |
42 |
opt3=FinleyOptions(solver_method=LinearPDE.DIRECT,verbose=show) |
opt3=FinleyOptions(solver_method=LinearPDE.PCG,preconditioner=LinearPDE.RILU,verbose=show) |
43 |
|
# opt4=FinleyOptions(solver_method=LinearPDE.DIRECT,verbose=show) |
44 |
|
|
45 |
bm_L2Do1=Benchmark(name="Laplace 2D (order 1)") |
bm_L2Do1=Benchmark(name="Laplace 2D (order 1)") |
46 |
bm_L2Do1.addProblem(Laplace2DOrder1_30k()) |
bm_L2Do1.addProblem(Laplace2DOrder1_30k()) |
113 |
bm_A2Do1g30.addProblem(Anisotropic2DOrder1Gamma30_3840k()) |
bm_A2Do1g30.addProblem(Anisotropic2DOrder1Gamma30_3840k()) |
114 |
bm_A2Do1g30.addProblem(Anisotropic2DOrder1Gamma30_7680k()) |
bm_A2Do1g30.addProblem(Anisotropic2DOrder1Gamma30_7680k()) |
115 |
bm_A2Do1g30.addProblem(Anisotropic2DOrder1Gamma30_15360k()) |
bm_A2Do1g30.addProblem(Anisotropic2DOrder1Gamma30_15360k()) |
116 |
|
bm_A2Do1g30.addOptions(opt1) |
117 |
|
bm_A2Do1g30.addOptions(opt2) |
118 |
|
bm_A2Do1g30.addOptions(opt3) |
119 |
|
|
120 |
bm_A2Do1g45=Benchmark("Anisotropic 2D (gamma=45, order=1)") |
bm_A2Do1g45=Benchmark("Anisotropic 2D (gamma=45, order=1)") |
121 |
bm_A2Do1g45.addProblem(Anisotropic2DOrder1Gamma45_30k()) |
bm_A2Do1g45.addProblem(Anisotropic2DOrder1Gamma45_30k()) |
128 |
bm_A2Do1g45.addProblem(Anisotropic2DOrder1Gamma45_3840k()) |
bm_A2Do1g45.addProblem(Anisotropic2DOrder1Gamma45_3840k()) |
129 |
bm_A2Do1g45.addProblem(Anisotropic2DOrder1Gamma45_7680k()) |
bm_A2Do1g45.addProblem(Anisotropic2DOrder1Gamma45_7680k()) |
130 |
bm_A2Do1g45.addProblem(Anisotropic2DOrder1Gamma45_15360k()) |
bm_A2Do1g45.addProblem(Anisotropic2DOrder1Gamma45_15360k()) |
131 |
|
bm_A2Do1g45.addOptions(opt1) |
132 |
|
bm_A2Do1g45.addOptions(opt2) |
133 |
|
bm_A2Do1g45.addOptions(opt3) |
134 |
|
|
135 |
bm_A2Do2g30=Benchmark("Anisotropic 2D (gamma=30, order=2)") |
bm_A2Do2g30=Benchmark("Anisotropic 2D (gamma=30, order=2)") |
136 |
bm_A2Do2g30.addProblem(Anisotropic2DOrder2Gamma30_30k()) |
bm_A2Do2g30.addProblem(Anisotropic2DOrder2Gamma30_30k()) |
143 |
bm_A2Do2g30.addProblem(Anisotropic2DOrder2Gamma30_3840k()) |
bm_A2Do2g30.addProblem(Anisotropic2DOrder2Gamma30_3840k()) |
144 |
bm_A2Do2g30.addProblem(Anisotropic2DOrder2Gamma30_7680k()) |
bm_A2Do2g30.addProblem(Anisotropic2DOrder2Gamma30_7680k()) |
145 |
bm_A2Do2g30.addProblem(Anisotropic2DOrder2Gamma30_15360k()) |
bm_A2Do2g30.addProblem(Anisotropic2DOrder2Gamma30_15360k()) |
146 |
|
bm_A2Do2g30.addOptions(opt1) |
147 |
|
bm_A2Do2g30.addOptions(opt2) |
148 |
|
bm_A2Do2g30.addOptions(opt3) |
149 |
|
|
150 |
|
|
151 |
bm_A2Do2g45=Benchmark("Anisotropic 2D (gamma=45, order=2)") |
bm_A2Do2g45=Benchmark("Anisotropic 2D (gamma=45, order=2)") |
159 |
bm_A2Do2g45.addProblem(Anisotropic2DOrder2Gamma45_3840k()) |
bm_A2Do2g45.addProblem(Anisotropic2DOrder2Gamma45_3840k()) |
160 |
bm_A2Do2g45.addProblem(Anisotropic2DOrder2Gamma45_7680k()) |
bm_A2Do2g45.addProblem(Anisotropic2DOrder2Gamma45_7680k()) |
161 |
bm_A2Do2g45.addProblem(Anisotropic2DOrder2Gamma45_15360k()) |
bm_A2Do2g45.addProblem(Anisotropic2DOrder2Gamma45_15360k()) |
162 |
|
bm_A2Do2g45.addOptions(opt1) |
163 |
|
bm_A2Do2g45.addOptions(opt2) |
164 |
|
bm_A2Do2g45.addOptions(opt3) |
165 |
|
|
166 |
|
|
167 |
bm_A3Do1g30=Benchmark("Anisotropic 3D (gamma=30, order=1)") |
bm_A3Do1g30=Benchmark("Anisotropic 3D (gamma=30, order=1)") |
175 |
bm_A3Do1g30.addProblem(Anisotropic3DOrder1Gamma30_3840k()) |
bm_A3Do1g30.addProblem(Anisotropic3DOrder1Gamma30_3840k()) |
176 |
bm_A3Do1g30.addProblem(Anisotropic3DOrder1Gamma30_7680k()) |
bm_A3Do1g30.addProblem(Anisotropic3DOrder1Gamma30_7680k()) |
177 |
bm_A3Do1g30.addProblem(Anisotropic3DOrder1Gamma30_15360k()) |
bm_A3Do1g30.addProblem(Anisotropic3DOrder1Gamma30_15360k()) |
178 |
|
bm_A3Do1g30.addOptions(opt1) |
179 |
|
bm_A3Do1g30.addOptions(opt2) |
180 |
|
bm_A3Do1g30.addOptions(opt3) |
181 |
|
|
182 |
bm_A3Do1g45=Benchmark("Anisotropic 3D (gamma=45, order=1)") |
bm_A3Do1g45=Benchmark("Anisotropic 3D (gamma=45, order=1)") |
183 |
bm_A3Do1g45.addProblem(Anisotropic3DOrder1Gamma45_30k()) |
bm_A3Do1g45.addProblem(Anisotropic3DOrder1Gamma45_30k()) |
190 |
bm_A3Do1g45.addProblem(Anisotropic3DOrder1Gamma45_3840k()) |
bm_A3Do1g45.addProblem(Anisotropic3DOrder1Gamma45_3840k()) |
191 |
bm_A3Do1g45.addProblem(Anisotropic3DOrder1Gamma45_7680k()) |
bm_A3Do1g45.addProblem(Anisotropic3DOrder1Gamma45_7680k()) |
192 |
bm_A3Do1g45.addProblem(Anisotropic3DOrder1Gamma45_15360k()) |
bm_A3Do1g45.addProblem(Anisotropic3DOrder1Gamma45_15360k()) |
193 |
|
bm_A3Do1g45.addOptions(opt1) |
194 |
|
bm_A3Do1g45.addOptions(opt2) |
195 |
|
bm_A3Do1g45.addOptions(opt3) |
196 |
|
|
197 |
bm_A3Do2g30=Benchmark("Anisotropic 3D (gamma=30, order=2)") |
bm_A3Do2g30=Benchmark("Anisotropic 3D (gamma=30, order=2)") |
198 |
bm_A3Do2g30.addProblem(Anisotropic3DOrder2Gamma30_30k()) |
bm_A3Do2g30.addProblem(Anisotropic3DOrder2Gamma30_30k()) |
205 |
bm_A3Do2g30.addProblem(Anisotropic3DOrder2Gamma30_3840k()) |
bm_A3Do2g30.addProblem(Anisotropic3DOrder2Gamma30_3840k()) |
206 |
bm_A3Do2g30.addProblem(Anisotropic3DOrder2Gamma30_7680k()) |
bm_A3Do2g30.addProblem(Anisotropic3DOrder2Gamma30_7680k()) |
207 |
bm_A3Do2g30.addProblem(Anisotropic3DOrder2Gamma30_15360k()) |
bm_A3Do2g30.addProblem(Anisotropic3DOrder2Gamma30_15360k()) |
208 |
|
bm_A3Do2g30.addOptions(opt1) |
209 |
|
bm_A3Do2g30.addOptions(opt2) |
210 |
|
bm_A3Do2g30.addOptions(opt3) |
211 |
|
|
212 |
bm_A3Do2g45=Benchmark("Anisotropic 3D (gamma=45, order=2)") |
bm_A3Do2g45=Benchmark("Anisotropic 3D (gamma=45, order=2)") |
213 |
bm_A3Do2g45.addProblem(Anisotropic3DOrder2Gamma45_30k()) |
bm_A3Do2g45.addProblem(Anisotropic3DOrder2Gamma45_30k()) |
220 |
bm_A3Do2g45.addProblem(Anisotropic3DOrder2Gamma45_3840k()) |
bm_A3Do2g45.addProblem(Anisotropic3DOrder2Gamma45_3840k()) |
221 |
bm_A3Do2g45.addProblem(Anisotropic3DOrder2Gamma45_7680k()) |
bm_A3Do2g45.addProblem(Anisotropic3DOrder2Gamma45_7680k()) |
222 |
bm_A3Do2g45.addProblem(Anisotropic3DOrder2Gamma45_15360k()) |
bm_A3Do2g45.addProblem(Anisotropic3DOrder2Gamma45_15360k()) |
223 |
|
bm_A3Do2g45.addOptions(opt1) |
224 |
|
bm_A3Do2g45.addOptions(opt2) |
225 |
|
bm_A3Do2g45.addOptions(opt3) |
226 |
|
|
227 |
|
|
228 |
bms=BenchmarkSuite("Paso/Finley") |
bms=BenchmarkSuite("Paso/Finley") |
230 |
bms.addBenchmark(bm_L2Do2) |
bms.addBenchmark(bm_L2Do2) |
231 |
bms.addBenchmark(bm_L3Do1) |
bms.addBenchmark(bm_L3Do1) |
232 |
bms.addBenchmark(bm_L3Do2) |
bms.addBenchmark(bm_L3Do2) |
233 |
# bms.addBenchmark(bm_A2Do1g30) |
bms.addBenchmark(bm_A2Do1g30) |
234 |
# bms.addBenchmark(bm_A2Do1g45) |
bms.addBenchmark(bm_A2Do1g45) |
235 |
# bms.addBenchmark(bm_A2Do2g30) |
bms.addBenchmark(bm_A2Do2g30) |
236 |
# bms.addBenchmark(bm_A2Do2g45) |
bms.addBenchmark(bm_A2Do2g45) |
237 |
# bms.addBenchmark(bm_A3Do1g30) |
bms.addBenchmark(bm_A3Do1g30) |
238 |
# bms.addBenchmark(bm_A3Do1g45) |
bms.addBenchmark(bm_A3Do1g45) |
239 |
# bms.addBenchmark(bm_A3Do2g30) |
bms.addBenchmark(bm_A3Do2g30) |
240 |
# bms.addBenchmark(bm_A3Do2g45) |
bms.addBenchmark(bm_A3Do2g45) |
241 |
bms.run(scale=thlist) |
bms.run(scale=thlist) |
242 |
out=bms.getHTML(filter=ff) |
out=bms.getHTML(filter=ff) |
243 |
print out |
print out |