/[escript]/trunk/finley/benchmarks/runbenchmark.py
ViewVC logotype

Diff of /trunk/finley/benchmarks/runbenchmark.py

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 415 by gross, Tue Dec 20 03:13:17 2005 UTC revision 416 by gross, Wed Jan 4 05:38:35 2006 UTC
# Line 12  Line 12 
12  #  #
13    
14  """  """
15  some benchmarks for tetsing the finley solver. The idea is to develop a set of standart benchmarks  some benchmarks for tetsing the finley solver.
   
   * Laplace2Dorder1_?k  
   * Laplace3Dorder2_?k  
   
 where ? is approximatively the number of unknowns in 1000.  
16    
17  @var __author__: name of author  @var __author__: name of author
18  @var __licence__: licence agreement  @var __licence__: licence agreement
# Line 35  __date__="$Date:$" Line 30  __date__="$Date:$"
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,32]  thlist=[1,2,4,8,16]
34    # 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=True  show=False
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)
43    
44  bm_L2Do1=Benchmark(name="Laplace 2D (order 1)")  bm_L2Do1=Benchmark(name="Laplace 2D (order 1)")
45  bm_L2Do1.addProblem(Laplace2DOrder1_30k())  bm_L2Do1.addProblem(Laplace2DOrder1_30k())
# Line 57  bm_L2Do1.addProblem(Laplace2DOrder1_7680 Line 54  bm_L2Do1.addProblem(Laplace2DOrder1_7680
54  bm_L2Do1.addProblem(Laplace2DOrder1_15360k())  bm_L2Do1.addProblem(Laplace2DOrder1_15360k())
55  bm_L2Do1.addOptions(opt1)  bm_L2Do1.addOptions(opt1)
56  bm_L2Do1.addOptions(opt2)  bm_L2Do1.addOptions(opt2)
57    bm_L2Do1.addOptions(opt3)
58    
59  bm_L2Do2=Benchmark("Laplace 2D (order 2)")  bm_L2Do2=Benchmark("Laplace 2D (order 2)")
60  bm_L2Do2.addProblem(Laplace2DOrder2_30k())  bm_L2Do2.addProblem(Laplace2DOrder2_30k())
# Line 71  bm_L2Do2.addProblem(Laplace2DOrder2_7680 Line 69  bm_L2Do2.addProblem(Laplace2DOrder2_7680
69  bm_L2Do2.addProblem(Laplace2DOrder2_15360k())  bm_L2Do2.addProblem(Laplace2DOrder2_15360k())
70  bm_L2Do2.addOptions(opt1)  bm_L2Do2.addOptions(opt1)
71  bm_L2Do2.addOptions(opt2)  bm_L2Do2.addOptions(opt2)
72    bm_L2Do2.addOptions(opt3)
73    
74  bm_L3Do1=Benchmark("Laplace 3D (order 1)")  bm_L3Do1=Benchmark("Laplace 3D (order 1)")
75  bm_L3Do1.addProblem(Laplace3DOrder1_30k())  bm_L3Do1.addProblem(Laplace3DOrder1_30k())
# Line 85  bm_L3Do1.addProblem(Laplace3DOrder1_7680 Line 84  bm_L3Do1.addProblem(Laplace3DOrder1_7680
84  bm_L3Do1.addProblem(Laplace3DOrder1_15360k())  bm_L3Do1.addProblem(Laplace3DOrder1_15360k())
85  bm_L3Do1.addOptions(opt1)  bm_L3Do1.addOptions(opt1)
86  bm_L3Do1.addOptions(opt2)  bm_L3Do1.addOptions(opt2)
87    bm_L3Do1.addOptions(opt3)
88    
89  bm_L3Do2=Benchmark("Laplace 3D (order 2)")  bm_L3Do2=Benchmark("Laplace 3D (order 2)")
90  bm_L3Do2.addProblem(Laplace3DOrder2_30k())  bm_L3Do2.addProblem(Laplace3DOrder2_30k())
# Line 99  bm_L3Do2.addProblem(Laplace3DOrder2_7680 Line 99  bm_L3Do2.addProblem(Laplace3DOrder2_7680
99  bm_L3Do2.addProblem(Laplace3DOrder2_15360k())  bm_L3Do2.addProblem(Laplace3DOrder2_15360k())
100  bm_L3Do2.addOptions(opt1)  bm_L3Do2.addOptions(opt1)
101  bm_L3Do2.addOptions(opt2)  bm_L3Do2.addOptions(opt2)
102    bm_L3Do2.addOptions(opt3)
103    
104    bm_A2Do1g30=Benchmark("Anisotropic 2D (gamma=30, order=1)")
105    bm_A2Do1g30.addProblem(Anisotropic2DOrder1Gamma30_30k())
106    bm_A2Do1g30.addProblem(Anisotropic2DOrder1Gamma30_60k())
107    bm_A2Do1g30.addProblem(Anisotropic2DOrder1Gamma30_120k())
108    bm_A2Do1g30.addProblem(Anisotropic2DOrder1Gamma30_240k())
109    bm_A2Do1g30.addProblem(Anisotropic2DOrder1Gamma30_480k())
110    bm_A2Do1g30.addProblem(Anisotropic2DOrder1Gamma30_960k())
111    bm_A2Do1g30.addProblem(Anisotropic2DOrder1Gamma30_1920k())
112    bm_A2Do1g30.addProblem(Anisotropic2DOrder1Gamma30_3840k())
113    bm_A2Do1g30.addProblem(Anisotropic2DOrder1Gamma30_7680k())
114    bm_A2Do1g30.addProblem(Anisotropic2DOrder1Gamma30_15360k())
115    
116    bm_A2Do1g45=Benchmark("Anisotropic 2D (gamma=45, order=1)")
117    bm_A2Do1g45.addProblem(Anisotropic2DOrder1Gamma45_30k())
118    bm_A2Do1g45.addProblem(Anisotropic2DOrder1Gamma45_60k())
119    bm_A2Do1g45.addProblem(Anisotropic2DOrder1Gamma45_120k())
120    bm_A2Do1g45.addProblem(Anisotropic2DOrder1Gamma45_240k())
121    bm_A2Do1g45.addProblem(Anisotropic2DOrder1Gamma45_480k())
122    bm_A2Do1g45.addProblem(Anisotropic2DOrder1Gamma45_960k())
123    bm_A2Do1g45.addProblem(Anisotropic2DOrder1Gamma45_1920k())
124    bm_A2Do1g45.addProblem(Anisotropic2DOrder1Gamma45_3840k())
125    bm_A2Do1g45.addProblem(Anisotropic2DOrder1Gamma45_7680k())
126    bm_A2Do1g45.addProblem(Anisotropic2DOrder1Gamma45_15360k())
127    
128    bm_A2Do2g30=Benchmark("Anisotropic 2D (gamma=30, order=2)")
129    bm_A2Do2g30.addProblem(Anisotropic2DOrder2Gamma30_30k())
130    bm_A2Do2g30.addProblem(Anisotropic2DOrder2Gamma30_60k())
131    bm_A2Do2g30.addProblem(Anisotropic2DOrder2Gamma30_120k())
132    bm_A2Do2g30.addProblem(Anisotropic2DOrder2Gamma30_240k())
133    bm_A2Do2g30.addProblem(Anisotropic2DOrder2Gamma30_480k())
134    bm_A2Do2g30.addProblem(Anisotropic2DOrder2Gamma30_960k())
135    bm_A2Do2g30.addProblem(Anisotropic2DOrder2Gamma30_1920k())
136    bm_A2Do2g30.addProblem(Anisotropic2DOrder2Gamma30_3840k())
137    bm_A2Do2g30.addProblem(Anisotropic2DOrder2Gamma30_7680k())
138    bm_A2Do2g30.addProblem(Anisotropic2DOrder2Gamma30_15360k())
139    
140    
141    bm_A2Do2g45=Benchmark("Anisotropic 2D (gamma=45, order=2)")
142    bm_A2Do2g45.addProblem(Anisotropic2DOrder2Gamma45_30k())
143    bm_A2Do2g45.addProblem(Anisotropic2DOrder2Gamma45_60k())
144    bm_A2Do2g45.addProblem(Anisotropic2DOrder2Gamma45_120k())
145    bm_A2Do2g45.addProblem(Anisotropic2DOrder2Gamma45_240k())
146    bm_A2Do2g45.addProblem(Anisotropic2DOrder2Gamma45_480k())
147    bm_A2Do2g45.addProblem(Anisotropic2DOrder2Gamma45_960k())
148    bm_A2Do2g45.addProblem(Anisotropic2DOrder2Gamma45_1920k())
149    bm_A2Do2g45.addProblem(Anisotropic2DOrder2Gamma45_3840k())
150    bm_A2Do2g45.addProblem(Anisotropic2DOrder2Gamma45_7680k())
151    bm_A2Do2g45.addProblem(Anisotropic2DOrder2Gamma45_15360k())
152    
153    
154    bm_A3Do1g30=Benchmark("Anisotropic 3D (gamma=30, order=1)")
155    bm_A3Do1g30.addProblem(Anisotropic3DOrder1Gamma30_30k())
156    bm_A3Do1g30.addProblem(Anisotropic3DOrder1Gamma30_60k())
157    bm_A3Do1g30.addProblem(Anisotropic3DOrder1Gamma30_120k())
158    bm_A3Do1g30.addProblem(Anisotropic3DOrder1Gamma30_240k())
159    bm_A3Do1g30.addProblem(Anisotropic3DOrder1Gamma30_480k())
160    bm_A3Do1g30.addProblem(Anisotropic3DOrder1Gamma30_960k())
161    bm_A3Do1g30.addProblem(Anisotropic3DOrder1Gamma30_1920k())
162    bm_A3Do1g30.addProblem(Anisotropic3DOrder1Gamma30_3840k())
163    bm_A3Do1g30.addProblem(Anisotropic3DOrder1Gamma30_7680k())
164    bm_A3Do1g30.addProblem(Anisotropic3DOrder1Gamma30_15360k())
165    
166    bm_A3Do1g45=Benchmark("Anisotropic 3D (gamma=45, order=1)")
167    bm_A3Do1g45.addProblem(Anisotropic3DOrder1Gamma45_30k())
168    bm_A3Do1g45.addProblem(Anisotropic3DOrder1Gamma45_60k())
169    bm_A3Do1g45.addProblem(Anisotropic3DOrder1Gamma45_120k())
170    bm_A3Do1g45.addProblem(Anisotropic3DOrder1Gamma45_240k())
171    bm_A3Do1g45.addProblem(Anisotropic3DOrder1Gamma45_480k())
172    bm_A3Do1g45.addProblem(Anisotropic3DOrder1Gamma45_960k())
173    bm_A3Do1g45.addProblem(Anisotropic3DOrder1Gamma45_1920k())
174    bm_A3Do1g45.addProblem(Anisotropic3DOrder1Gamma45_3840k())
175    bm_A3Do1g45.addProblem(Anisotropic3DOrder1Gamma45_7680k())
176    bm_A3Do1g45.addProblem(Anisotropic3DOrder1Gamma45_15360k())
177    
178    bm_A3Do2g30=Benchmark("Anisotropic 3D (gamma=30, order=2)")
179    bm_A3Do2g30.addProblem(Anisotropic3DOrder2Gamma30_30k())
180    bm_A3Do2g30.addProblem(Anisotropic3DOrder2Gamma30_60k())
181    bm_A3Do2g30.addProblem(Anisotropic3DOrder2Gamma30_120k())
182    bm_A3Do2g30.addProblem(Anisotropic3DOrder2Gamma30_240k())
183    bm_A3Do2g30.addProblem(Anisotropic3DOrder2Gamma30_480k())
184    bm_A3Do2g30.addProblem(Anisotropic3DOrder2Gamma30_960k())
185    bm_A3Do2g30.addProblem(Anisotropic3DOrder2Gamma30_1920k())
186    bm_A3Do2g30.addProblem(Anisotropic3DOrder2Gamma30_3840k())
187    bm_A3Do2g30.addProblem(Anisotropic3DOrder2Gamma30_7680k())
188    bm_A3Do2g30.addProblem(Anisotropic3DOrder2Gamma30_15360k())
189    
190    bm_A3Do2g45=Benchmark("Anisotropic 3D (gamma=45, order=2)")
191    bm_A3Do2g45.addProblem(Anisotropic3DOrder2Gamma45_30k())
192    bm_A3Do2g45.addProblem(Anisotropic3DOrder2Gamma45_60k())
193    bm_A3Do2g45.addProblem(Anisotropic3DOrder2Gamma45_120k())
194    bm_A3Do2g45.addProblem(Anisotropic3DOrder2Gamma45_240k())
195    bm_A3Do2g45.addProblem(Anisotropic3DOrder2Gamma45_480k())
196    bm_A3Do2g45.addProblem(Anisotropic3DOrder2Gamma45_960k())
197    bm_A3Do2g45.addProblem(Anisotropic3DOrder2Gamma45_1920k())
198    bm_A3Do2g45.addProblem(Anisotropic3DOrder2Gamma45_3840k())
199    bm_A3Do2g45.addProblem(Anisotropic3DOrder2Gamma45_7680k())
200    bm_A3Do2g45.addProblem(Anisotropic3DOrder2Gamma45_15360k())
201    
202    
203  bms=BenchmarkSuite("Paso/Finley")  bms=BenchmarkSuite("Paso/Finley")
# Line 107  bms.addBenchmark(bm_L2Do1) Line 205  bms.addBenchmark(bm_L2Do1)
205  bms.addBenchmark(bm_L2Do2)  bms.addBenchmark(bm_L2Do2)
206  bms.addBenchmark(bm_L3Do1)  bms.addBenchmark(bm_L3Do1)
207  bms.addBenchmark(bm_L3Do2)  bms.addBenchmark(bm_L3Do2)
208    # bms.addBenchmark(bm_A2Do1g30)
209    # bms.addBenchmark(bm_A2Do1g45)
210    # bms.addBenchmark(bm_A2Do2g30)
211    # bms.addBenchmark(bm_A2Do2g45)
212    # bms.addBenchmark(bm_A3Do1g30)
213    # bms.addBenchmark(bm_A3Do1g45)
214    # bms.addBenchmark(bm_A3Do2g30)
215    # bms.addBenchmark(bm_A3Do2g45)
216  bms.run(scale=thlist)  bms.run(scale=thlist)
217  print bms.getHTML(filter=ff)  out=bms.getHTML(filter=ff)
218    print out

Legend:
Removed from v.415  
changed lines
  Added in v.416

  ViewVC Help
Powered by ViewVC 1.1.26