/[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 446 by gross, Mon Jan 23 00:36:07 2006 UTC revision 452 by gross, Mon Jan 23 23:48:50 2006 UTC
# Line 166  class FinleyProblem(BenchmarkProblem): Line 166  class FinleyProblem(BenchmarkProblem):
166         pde.setSolverMethod(options.solver_method,options.preconditioner)         pde.setSolverMethod(options.solver_method,options.preconditioner)
167         pde.setSolverPackage(options.package)         pde.setSolverPackage(options.package)
168         a=os.times()[4]         a=os.times()[4]
169         uh=pde.getSolution(verbose=options.verbose)         uh=pde.getSolution(verbose=options.verbose,iter_max=6000)
170         a=os.times()[4]-a         a=os.times()[4]-a
171         if u==None:         if u==None:
172            return {FinleyFilter.TIME : a , FinleyFilter.ERROR : None }            return {FinleyFilter.TIME : a , FinleyFilter.ERROR : None }
# Line 280  class AnisotropicProblem(RegularFinleyPr Line 280  class AnisotropicProblem(RegularFinleyPr
280              u*=(x[i]-i)              u*=(x[i]-i)
281    
282           gamma_rad=self.gamma/360.*8*math.atan(1.)           gamma_rad=self.gamma/360.*8*math.atan(1.)
283           cg=math.cos(self.gamma_rad)           cg=math.cos(gamma_rad)
284           sg=math.sin(self.gamma_rad)           sg=math.sin(gamma_rad)
285           C=kronecker(domain)           C=kronecker(domain)
286           C[0,0]=cg**2+self.c*sg**2           C[0,0]=cg**2+self.c*sg**2
287           C[1,0]=(self.c-1.)*cg*sg           C[1,0]=(self.c-1.)*cg*sg

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

  ViewVC Help
Powered by ViewVC 1.1.26