# Diff of /branches/diaplayground/poisson.py

revision 5095 by caltinay, Fri Jul 11 05:54:25 2014 UTC revision 5122 by caltinay, Thu Aug 21 04:32:39 2014 UTC
# Line 3  from esys.escript.linearPDEs import Line Line 3  from esys.escript.linearPDEs import Line
3  from esys.ripley import Rectangle, Brick  from esys.ripley import Rectangle, Brick
4  from time import time  from time import time
5
6  BLOCKSIZE=3  BLOCKSIZE=1
7
8  dom = Rectangle(l0=1.,l1=1.,n0=299, n1=299)  dom = Rectangle(l0=1.,l1=1.,n0=299, n1=299)
9  #dom = Brick(l0=1.,l1=1.,n0=51, n1=51, n2=17)  #dom = Brick(l0=1.,l1=1.,n0=51, n1=51, n2=17)
# Line 27  def solve(target): Line 27  def solve(target):
27              Y_reduced[i] = 1.              Y_reduced[i] = 1.
28
29      pde.setValue(A=A, Y_reduced=Y_reduced, q=q)      pde.setValue(A=A, Y_reduced=Y_reduced, q=q)
30        pde.getSolverOptions().setPackage(SolverOptions.PASO)
31        pde.getSolverOptions().setSolverTarget(target)
32      pde.getSolverOptions().setSolverMethod(SolverOptions.PRES20)      pde.getSolverOptions().setSolverMethod(SolverOptions.PRES20)
33      pde.getSolverOptions().setSolverMethod(SolverOptions.BICGSTAB)      pde.getSolverOptions().setSolverMethod(SolverOptions.BICGSTAB)
34      pde.getSolverOptions().setSolverMethod(SolverOptions.PCG)      pde.getSolverOptions().setSolverMethod(SolverOptions.PCG)
pde.getSolverOptions().setSolverTarget(target)
35      pde.getSolverOptions().setPreconditioner(SolverOptions.NO_PRECONDITIONER)      pde.getSolverOptions().setPreconditioner(SolverOptions.NO_PRECONDITIONER)
36      pde.getSolverOptions().setPreconditioner(SolverOptions.JACOBI)      pde.getSolverOptions().setPreconditioner(SolverOptions.JACOBI)
37      pde.getSolverOptions().setVerbosityOn()      pde.getSolverOptions().setVerbosityOn()

Legend:
 Removed from v.5095 changed lines Added in v.5122