/[escript]/trunk/finley/py_src/finleybench.py
ViewVC logotype

Diff of /trunk/finley/py_src/finleybench.py

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

revision 389 by gross, Tue Dec 20 03:17:28 2005 UTC revision 446 by gross, Mon Jan 23 00:36:07 2006 UTC
# Line 18  some benchmarks for tetsing the finley s Line 18  some benchmarks for tetsing the finley s
18    * Laplace3Dorder1_?k    * Laplace3Dorder1_?k
19    * Laplace3Dorder2_?k    * Laplace3Dorder2_?k
20    
21    
22    
23  where ? is approximatively the number of unknowns in 1000.  where ? is approximatively the number of unknowns in 1000.
24    
25  @var __author__: name of author  @var __author__: name of author
# Line 38  from esys.escript.benchmark import Bench Line 40  from esys.escript.benchmark import Bench
40  import esys.finley  import esys.finley
41  from esys.escript.linearPDEs import LinearPDE  from esys.escript.linearPDEs import LinearPDE
42  import os  import os
43    import math
44    
45  class FinleyFilter(BenchmarkFilter):  class FinleyFilter(BenchmarkFilter):
46     """     """
# Line 110  class FinleyOptions(Options): Line 113  class FinleyOptions(Options):
113                        LinearPDE.SCSL:  "SCSL",                        LinearPDE.SCSL:  "SCSL",
114                        LinearPDE.MKL:  "MKL",                        LinearPDE.MKL:  "MKL",
115                        LinearPDE.UMFPACK: "UMFPACK",                        LinearPDE.UMFPACK: "UMFPACK",
116                        LinearPDE.PASO:  "PASO"                        LinearPDE.PASO:  "PASO",
117                          LinearPDE.RILU: "RILU",
118                          LinearPDE.AMG:  "AMG"
119                    }                    }
120         name=""         name=""
121         if solver_method==None:         if solver_method==None:
# Line 248  class LaplaceProblem(RegularFinleyProble Line 253  class LaplaceProblem(RegularFinleyProble
253           pde.setValue(A=kronecker(domain),q=msk,r=u)           pde.setValue(A=kronecker(domain),q=msk,r=u)
254           return pde,u           return pde,u
255    
256    class AnisotropicProblem(RegularFinleyProblem):
257        """
258        base class for the Anisotropic scalar problem on a rectangular mesh
259        """
260        def __init__(self,n,order,dim,gamma,c):
261            self.c=c
262            self.gamma=gamma
263            super(AnisotropicProblem,self).__init__(n,order,dim)
264    
265    
266        def getTestProblem(self,domain):
267             """
268             returns a PDE and a test solution on the given domain
269        
270             @param doamin: a domain
271             @type domain: L{escript.Domain}
272             @return: the Laplace equation and a test solution
273             @rtype: C{tuple} of C{LinearPDE} and C{escript.Data}
274             """
275             x=domain.getX()
276             msk=whereZero(x[0])+whereZero(x[0]-1.)
277             u=x[0]
278             for i in range(1,domain.getDim()):
279                msk+=whereZero(x[i])+whereZero(x[i]-1.)
280                u*=(x[i]-i)
281    
282             gamma_rad=self.gamma/360.*8*math.atan(1.)
283             cg=math.cos(self.gamma_rad)
284             sg=math.sin(self.gamma_rad)
285             C=kronecker(domain)
286             C[0,0]=cg**2+self.c*sg**2
287             C[1,0]=(self.c-1.)*cg*sg
288             C[0,1]=C[0,1]
289             C[1,1]=sg**2+self.c*cg**2
290             F=2*(1.-self.c)*cg*sg
291             if domain.getDim()==3: F*=x[2]-2.
292             pde=LinearPDE(domain)
293             pde.setSymmetryOn()
294             pde.setValue(A=C,Y=F,q=msk,r=u)
295             return pde,u
296    
297  class Laplace2DOrder1_30k(LaplaceProblem):  class Laplace2DOrder1_30k(LaplaceProblem):
298     def __init__(self):     def __init__(self):
299        super(Laplace2DOrder1_30k,self).__init__(n=172,order=1,dim=2)        super(Laplace2DOrder1_30k,self).__init__(n=172,order=1,dim=2)
# Line 369  class Laplace3DOrder2_15360k(LaplaceProb Line 415  class Laplace3DOrder2_15360k(LaplaceProb
415     def __init__(self):     def __init__(self):
416        super(Laplace3DOrder2_15360k,self).__init__(n=124,order=2,dim=3)        super(Laplace3DOrder2_15360k,self).__init__(n=124,order=2,dim=3)
417    
418    class Anisotropic2DOrder1Gamma30_30k(AnisotropicProblem):
419       def __init__(self):
420          super(Anisotropic2DOrder1Gamma30_30k,self).__init__(n=172,order=1,dim=2,gamma=30,c=0.001)
421    class Anisotropic2DOrder1Gamma30_60k(AnisotropicProblem):
422       def __init__(self):
423          super(Anisotropic2DOrder1Gamma30_60k,self).__init__(n=244,order=1,dim=2,gamma=30,c=0.001)
424    class Anisotropic2DOrder1Gamma30_120k(AnisotropicProblem):
425       def __init__(self):
426          super(Anisotropic2DOrder1Gamma30_120k,self).__init__(n=345,order=1,dim=2,gamma=30,c=0.001)
427    class Anisotropic2DOrder1Gamma30_240k(AnisotropicProblem):
428       def __init__(self):
429          super(Anisotropic2DOrder1Gamma30_240k,self).__init__(n=489,order=1,dim=2,gamma=30,c=0.001)
430    class Anisotropic2DOrder1Gamma30_480k(AnisotropicProblem):
431       def __init__(self):
432          super(Anisotropic2DOrder1Gamma30_480k,self).__init__(n=692,order=1,dim=2,gamma=30,c=0.001)
433    class Anisotropic2DOrder1Gamma30_960k(AnisotropicProblem):
434       def __init__(self):
435          super(Anisotropic2DOrder1Gamma30_960k,self).__init__(n=979,order=1,dim=2,gamma=30,c=0.001)
436    class Anisotropic2DOrder1Gamma30_1920k(AnisotropicProblem):
437       def __init__(self):
438          super(Anisotropic2DOrder1Gamma30_1920k,self).__init__(n=1385,order=1,dim=2,gamma=30,c=0.001)
439    class Anisotropic2DOrder1Gamma30_3840k(AnisotropicProblem):
440       def __init__(self):
441          super(Anisotropic2DOrder1Gamma30_3840k,self).__init__(n=1959,order=1,dim=2,gamma=30,c=0.001)
442    class Anisotropic2DOrder1Gamma30_7680k(AnisotropicProblem):
443       def __init__(self):
444          super(Anisotropic2DOrder1Gamma30_7680k,self).__init__(n=2770,order=1,dim=2,gamma=30,c=0.001)
445    class Anisotropic2DOrder1Gamma30_15360k(AnisotropicProblem):
446       def __init__(self):
447          super(Anisotropic2DOrder1Gamma30_15360k,self).__init__(n=3918,order=1,dim=2,gamma=30,c=0.001)
448    class Anisotropic2DOrder1Gamma45_30k(AnisotropicProblem):
449       def __init__(self):
450          super(Anisotropic2DOrder1Gamma45_30k,self).__init__(n=172,order=1,dim=2,gamma=45,c=0.001)
451    class Anisotropic2DOrder1Gamma45_60k(AnisotropicProblem):
452       def __init__(self):
453          super(Anisotropic2DOrder1Gamma45_60k,self).__init__(n=244,order=1,dim=2,gamma=45,c=0.001)
454    class Anisotropic2DOrder1Gamma45_120k(AnisotropicProblem):
455       def __init__(self):
456          super(Anisotropic2DOrder1Gamma45_120k,self).__init__(n=345,order=1,dim=2,gamma=45,c=0.001)
457    class Anisotropic2DOrder1Gamma45_240k(AnisotropicProblem):
458       def __init__(self):
459          super(Anisotropic2DOrder1Gamma45_240k,self).__init__(n=489,order=1,dim=2,gamma=45,c=0.001)
460    class Anisotropic2DOrder1Gamma45_480k(AnisotropicProblem):
461       def __init__(self):
462          super(Anisotropic2DOrder1Gamma45_480k,self).__init__(n=692,order=1,dim=2,gamma=45,c=0.001)
463    class Anisotropic2DOrder1Gamma45_960k(AnisotropicProblem):
464       def __init__(self):
465          super(Anisotropic2DOrder1Gamma45_960k,self).__init__(n=979,order=1,dim=2,gamma=45,c=0.001)
466    class Anisotropic2DOrder1Gamma45_1920k(AnisotropicProblem):
467       def __init__(self):
468          super(Anisotropic2DOrder1Gamma45_1920k,self).__init__(n=1385,order=1,dim=2,gamma=45,c=0.001)
469    class Anisotropic2DOrder1Gamma45_3840k(AnisotropicProblem):
470       def __init__(self):
471          super(Anisotropic2DOrder1Gamma45_3840k,self).__init__(n=1959,order=1,dim=2,gamma=45,c=0.001)
472    class Anisotropic2DOrder1Gamma45_7680k(AnisotropicProblem):
473       def __init__(self):
474          super(Anisotropic2DOrder1Gamma45_7680k,self).__init__(n=2770,order=1,dim=2,gamma=45,c=0.001)
475    class Anisotropic2DOrder1Gamma45_15360k(AnisotropicProblem):
476       def __init__(self):
477          super(Anisotropic2DOrder1Gamma45_15360k,self).__init__(n=3918,order=1,dim=2,gamma=45,c=0.001)
478    class Anisotropic2DOrder2Gamma30_30k(AnisotropicProblem):
479       def __init__(self):
480          super(Anisotropic2DOrder2Gamma30_30k,self).__init__(n=86,order=2,dim=2,gamma=30,c=0.001)
481    class Anisotropic2DOrder2Gamma30_60k(AnisotropicProblem):
482       def __init__(self):
483          super(Anisotropic2DOrder2Gamma30_60k,self).__init__(n=122,order=2,dim=2,gamma=30,c=0.001)
484    class Anisotropic2DOrder2Gamma30_120k(AnisotropicProblem):
485       def __init__(self):
486          super(Anisotropic2DOrder2Gamma30_120k,self).__init__(n=173,order=2,dim=2,gamma=30,c=0.001)
487    class Anisotropic2DOrder2Gamma30_240k(AnisotropicProblem):
488       def __init__(self):
489          super(Anisotropic2DOrder2Gamma30_240k,self).__init__(n=244,order=2,dim=2,gamma=30,c=0.001)
490    class Anisotropic2DOrder2Gamma30_480k(AnisotropicProblem):
491       def __init__(self):
492          super(Anisotropic2DOrder2Gamma30_480k,self).__init__(n=346,order=2,dim=2,gamma=30,c=0.001)
493    class Anisotropic2DOrder2Gamma30_960k(AnisotropicProblem):
494       def __init__(self):
495          super(Anisotropic2DOrder2Gamma30_960k,self).__init__(n=489,order=2,dim=2,gamma=30,c=0.001)
496    class Anisotropic2DOrder2Gamma30_1920k(AnisotropicProblem):
497       def __init__(self):
498          super(Anisotropic2DOrder2Gamma30_1920k,self).__init__(n=692,order=2,dim=2,gamma=30,c=0.001)
499    class Anisotropic2DOrder2Gamma30_3840k(AnisotropicProblem):
500       def __init__(self):
501          super(Anisotropic2DOrder2Gamma30_3840k,self).__init__(n=979,order=2,dim=2,gamma=30,c=0.001)
502    class Anisotropic2DOrder2Gamma30_7680k(AnisotropicProblem):
503       def __init__(self):
504          super(Anisotropic2DOrder2Gamma30_7680k,self).__init__(n=1385,order=2,dim=2,gamma=30,c=0.001)
505    class Anisotropic2DOrder2Gamma30_15360k(AnisotropicProblem):
506       def __init__(self):
507          super(Anisotropic2DOrder2Gamma30_15360k,self).__init__(n=1959,order=2,dim=2,gamma=30,c=0.001)
508    class Anisotropic2DOrder2Gamma45_30k(AnisotropicProblem):
509       def __init__(self):
510          super(Anisotropic2DOrder2Gamma45_30k,self).__init__(n=86,order=2,dim=2,gamma=45,c=0.001)
511    class Anisotropic2DOrder2Gamma45_60k(AnisotropicProblem):
512       def __init__(self):
513          super(Anisotropic2DOrder2Gamma45_60k,self).__init__(n=122,order=2,dim=2,gamma=45,c=0.001)
514    class Anisotropic2DOrder2Gamma45_120k(AnisotropicProblem):
515       def __init__(self):
516          super(Anisotropic2DOrder2Gamma45_120k,self).__init__(n=173,order=2,dim=2,gamma=45,c=0.001)
517    class Anisotropic2DOrder2Gamma45_240k(AnisotropicProblem):
518       def __init__(self):
519          super(Anisotropic2DOrder2Gamma45_240k,self).__init__(n=244,order=2,dim=2,gamma=45,c=0.001)
520    class Anisotropic2DOrder2Gamma45_480k(AnisotropicProblem):
521       def __init__(self):
522          super(Anisotropic2DOrder2Gamma45_480k,self).__init__(n=346,order=2,dim=2,gamma=45,c=0.001)
523    class Anisotropic2DOrder2Gamma45_960k(AnisotropicProblem):
524       def __init__(self):
525          super(Anisotropic2DOrder2Gamma45_960k,self).__init__(n=489,order=2,dim=2,gamma=45,c=0.001)
526    class Anisotropic2DOrder2Gamma45_1920k(AnisotropicProblem):
527       def __init__(self):
528          super(Anisotropic2DOrder2Gamma45_1920k,self).__init__(n=692,order=2,dim=2,gamma=45,c=0.001)
529    class Anisotropic2DOrder2Gamma45_3840k(AnisotropicProblem):
530       def __init__(self):
531          super(Anisotropic2DOrder2Gamma45_3840k,self).__init__(n=979,order=2,dim=2,gamma=45,c=0.001)
532    class Anisotropic2DOrder2Gamma45_7680k(AnisotropicProblem):
533       def __init__(self):
534          super(Anisotropic2DOrder2Gamma45_7680k,self).__init__(n=1385,order=2,dim=2,gamma=45,c=0.001)
535    class Anisotropic2DOrder2Gamma45_15360k(AnisotropicProblem):
536       def __init__(self):
537          super(Anisotropic2DOrder2Gamma45_15360k,self).__init__(n=1959,order=2,dim=2,gamma=45,c=0.001)
538    class Anisotropic3DOrder1Gamma30_30k(AnisotropicProblem):
539       def __init__(self):
540          super(Anisotropic3DOrder1Gamma30_30k,self).__init__(n=30,order=1,dim=3,gamma=30,c=0.001)
541    class Anisotropic3DOrder1Gamma30_60k(AnisotropicProblem):
542       def __init__(self):
543          super(Anisotropic3DOrder1Gamma30_60k,self).__init__(n=38,order=1,dim=3,gamma=30,c=0.001)
544    class Anisotropic3DOrder1Gamma30_120k(AnisotropicProblem):
545       def __init__(self):
546          super(Anisotropic3DOrder1Gamma30_120k,self).__init__(n=48,order=1,dim=3,gamma=30,c=0.001)
547    class Anisotropic3DOrder1Gamma30_240k(AnisotropicProblem):
548       def __init__(self):
549          super(Anisotropic3DOrder1Gamma30_240k,self).__init__(n=61,order=1,dim=3,gamma=30,c=0.001)
550    class Anisotropic3DOrder1Gamma30_480k(AnisotropicProblem):
551       def __init__(self):
552          super(Anisotropic3DOrder1Gamma30_480k,self).__init__(n=77,order=1,dim=3,gamma=30,c=0.001)
553    class Anisotropic3DOrder1Gamma30_960k(AnisotropicProblem):
554       def __init__(self):
555          super(Anisotropic3DOrder1Gamma30_960k,self).__init__(n=98,order=1,dim=3,gamma=30,c=0.001)
556    class Anisotropic3DOrder1Gamma30_1920k(AnisotropicProblem):
557       def __init__(self):
558          super(Anisotropic3DOrder1Gamma30_1920k,self).__init__(n=123,order=1,dim=3,gamma=30,c=0.001)
559    class Anisotropic3DOrder1Gamma30_3840k(AnisotropicProblem):
560       def __init__(self):
561          super(Anisotropic3DOrder1Gamma30_3840k,self).__init__(n=156,order=1,dim=3,gamma=30,c=0.001)
562    class Anisotropic3DOrder1Gamma30_7680k(AnisotropicProblem):
563       def __init__(self):
564          super(Anisotropic3DOrder1Gamma30_7680k,self).__init__(n=196,order=1,dim=3,gamma=30,c=0.001)
565    class Anisotropic3DOrder1Gamma30_15360k(AnisotropicProblem):
566       def __init__(self):
567          super(Anisotropic3DOrder1Gamma30_15360k,self).__init__(n=248,order=1,dim=3,gamma=30,c=0.001)
568    class Anisotropic3DOrder1Gamma45_30k(AnisotropicProblem):
569       def __init__(self):
570          super(Anisotropic3DOrder1Gamma45_30k,self).__init__(n=30,order=1,dim=3,gamma=45,c=0.001)
571    class Anisotropic3DOrder1Gamma45_60k(AnisotropicProblem):
572       def __init__(self):
573          super(Anisotropic3DOrder1Gamma45_60k,self).__init__(n=38,order=1,dim=3,gamma=45,c=0.001)
574    class Anisotropic3DOrder1Gamma45_120k(AnisotropicProblem):
575       def __init__(self):
576          super(Anisotropic3DOrder1Gamma45_120k,self).__init__(n=48,order=1,dim=3,gamma=45,c=0.001)
577    class Anisotropic3DOrder1Gamma45_240k(AnisotropicProblem):
578       def __init__(self):
579          super(Anisotropic3DOrder1Gamma45_240k,self).__init__(n=61,order=1,dim=3,gamma=45,c=0.001)
580    class Anisotropic3DOrder1Gamma45_480k(AnisotropicProblem):
581       def __init__(self):
582          super(Anisotropic3DOrder1Gamma45_480k,self).__init__(n=77,order=1,dim=3,gamma=45,c=0.001)
583    class Anisotropic3DOrder1Gamma45_960k(AnisotropicProblem):
584       def __init__(self):
585          super(Anisotropic3DOrder1Gamma45_960k,self).__init__(n=98,order=1,dim=3,gamma=45,c=0.001)
586    class Anisotropic3DOrder1Gamma45_1920k(AnisotropicProblem):
587       def __init__(self):
588          super(Anisotropic3DOrder1Gamma45_1920k,self).__init__(n=123,order=1,dim=3,gamma=45,c=0.001)
589    class Anisotropic3DOrder1Gamma45_3840k(AnisotropicProblem):
590       def __init__(self):
591          super(Anisotropic3DOrder1Gamma45_3840k,self).__init__(n=156,order=1,dim=3,gamma=45,c=0.001)
592    class Anisotropic3DOrder1Gamma45_7680k(AnisotropicProblem):
593       def __init__(self):
594          super(Anisotropic3DOrder1Gamma45_7680k,self).__init__(n=196,order=1,dim=3,gamma=45,c=0.001)
595    class Anisotropic3DOrder1Gamma45_15360k(AnisotropicProblem):
596       def __init__(self):
597          super(Anisotropic3DOrder1Gamma45_15360k,self).__init__(n=248,order=1,dim=3,gamma=45,c=0.001)
598    class Anisotropic3DOrder2Gamma30_30k(AnisotropicProblem):
599       def __init__(self):
600          super(Anisotropic3DOrder2Gamma30_30k,self).__init__(n=15,order=2,dim=3,gamma=30,c=0.001)
601    class Anisotropic3DOrder2Gamma30_60k(AnisotropicProblem):
602       def __init__(self):
603          super(Anisotropic3DOrder2Gamma30_60k,self).__init__(n=19,order=2,dim=3,gamma=30,c=0.001)
604    class Anisotropic3DOrder2Gamma30_120k(AnisotropicProblem):
605       def __init__(self):
606          super(Anisotropic3DOrder2Gamma30_120k,self).__init__(n=24,order=2,dim=3,gamma=30,c=0.001)
607    class Anisotropic3DOrder2Gamma30_240k(AnisotropicProblem):
608       def __init__(self):
609          super(Anisotropic3DOrder2Gamma30_240k,self).__init__(n=31,order=2,dim=3,gamma=30,c=0.001)
610    class Anisotropic3DOrder2Gamma30_480k(AnisotropicProblem):
611       def __init__(self):
612          super(Anisotropic3DOrder2Gamma30_480k,self).__init__(n=39,order=2,dim=3,gamma=30,c=0.001)
613    class Anisotropic3DOrder2Gamma30_960k(AnisotropicProblem):
614       def __init__(self):
615          super(Anisotropic3DOrder2Gamma30_960k,self).__init__(n=49,order=2,dim=3,gamma=30,c=0.001)
616    class Anisotropic3DOrder2Gamma30_1920k(AnisotropicProblem):
617       def __init__(self):
618          super(Anisotropic3DOrder2Gamma30_1920k,self).__init__(n=62,order=2,dim=3,gamma=30,c=0.001)
619    class Anisotropic3DOrder2Gamma30_3840k(AnisotropicProblem):
620       def __init__(self):
621          super(Anisotropic3DOrder2Gamma30_3840k,self).__init__(n=78,order=2,dim=3,gamma=30,c=0.001)
622    class Anisotropic3DOrder2Gamma30_7680k(AnisotropicProblem):
623       def __init__(self):
624          super(Anisotropic3DOrder2Gamma30_7680k,self).__init__(n=98,order=2,dim=3,gamma=30,c=0.001)
625    class Anisotropic3DOrder2Gamma30_15360k(AnisotropicProblem):
626       def __init__(self):
627          super(Anisotropic3DOrder2Gamma30_15360k,self).__init__(n=124,order=2,dim=3,gamma=30,c=0.001)
628    class Anisotropic3DOrder2Gamma45_30k(AnisotropicProblem):
629       def __init__(self):
630          super(Anisotropic3DOrder2Gamma45_30k,self).__init__(n=15,order=2,dim=3,gamma=45,c=0.001)
631    class Anisotropic3DOrder2Gamma45_60k(AnisotropicProblem):
632       def __init__(self):
633          super(Anisotropic3DOrder2Gamma45_60k,self).__init__(n=19,order=2,dim=3,gamma=45,c=0.001)
634    class Anisotropic3DOrder2Gamma45_120k(AnisotropicProblem):
635       def __init__(self):
636          super(Anisotropic3DOrder2Gamma45_120k,self).__init__(n=24,order=2,dim=3,gamma=45,c=0.001)
637    class Anisotropic3DOrder2Gamma45_240k(AnisotropicProblem):
638       def __init__(self):
639          super(Anisotropic3DOrder2Gamma45_240k,self).__init__(n=31,order=2,dim=3,gamma=45,c=0.001)
640    class Anisotropic3DOrder2Gamma45_480k(AnisotropicProblem):
641       def __init__(self):
642          super(Anisotropic3DOrder2Gamma45_480k,self).__init__(n=39,order=2,dim=3,gamma=45,c=0.001)
643    class Anisotropic3DOrder2Gamma45_960k(AnisotropicProblem):
644       def __init__(self):
645          super(Anisotropic3DOrder2Gamma45_960k,self).__init__(n=49,order=2,dim=3,gamma=45,c=0.001)
646    class Anisotropic3DOrder2Gamma45_1920k(AnisotropicProblem):
647       def __init__(self):
648          super(Anisotropic3DOrder2Gamma45_1920k,self).__init__(n=62,order=2,dim=3,gamma=45,c=0.001)
649    class Anisotropic3DOrder2Gamma45_3840k(AnisotropicProblem):
650       def __init__(self):
651          super(Anisotropic3DOrder2Gamma45_3840k,self).__init__(n=78,order=2,dim=3,gamma=45,c=0.001)
652    class Anisotropic3DOrder2Gamma45_7680k(AnisotropicProblem):
653       def __init__(self):
654          super(Anisotropic3DOrder2Gamma45_7680k,self).__init__(n=98,order=2,dim=3,gamma=45,c=0.001)
655    class Anisotropic3DOrder2Gamma45_15360k(AnisotropicProblem):
656       def __init__(self):
657          super(Anisotropic3DOrder2Gamma45_15360k,self).__init__(n=124,order=2,dim=3,gamma=45,c=0.001)
658    
659  if __name__=="__main__":  if __name__=="__main__":
660     test=""     test=""
661     n0=30000     n0=30000
662     for d in [2,3]:     for d in [2,3]:
663       for o in [1,2]:      for o in [1,2]:
664          for g in [30,45]:
665          for i in range(10):          for i in range(10):
666               dofs=n0*2**i               dofs=n0*2**i
667               n=int((float(dofs)**(1./float(d))-1)/o+0.5)               n=int((float(dofs)**(1./float(d))-1)/o+0.5)
668               name="Laplace%sDOrder%s_%sk"%(d,o,dofs/1000)               name="Anisotropic%sDOrder%sGamma%s_%sk"%(d,o,g,dofs/1000)
669               print "class %s(LaplaceProblem):"%name               print "class %s(AnisotropicProblem):"%name
670               print "   def __init__(self):"               print "   def __init__(self):"
671               print "      super(%s,self).__init__(n=%s,order=%s,dim=%s)"%(name,n,o,d)               print "      super(%s,self).__init__(n=%s,order=%s,dim=%s,gamma=%s,c=0.001)"%(name,n,o,d,g)
672               test+="addProblem(%s())\n"%name               test+="addProblem(%s())\n"%name
673     print test     print test
674    

Legend:
Removed from v.389  
changed lines
  Added in v.446

  ViewVC Help
Powered by ViewVC 1.1.26