/[escript]/trunk/finley/test/python/SimpleSolve.py
ViewVC logotype

Diff of /trunk/finley/test/python/SimpleSolve.py

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

revision 102 by jgs, Wed Dec 15 07:08:39 2004 UTC revision 104 by jgs, Fri Dec 17 07:43:12 2004 UTC
# Line 48  print "norm_u_ex=u_ex.Lsup():" Line 48  print "norm_u_ex=u_ex.Lsup():"
48  norm_u_ex=u_ex.Lsup()  norm_u_ex=u_ex.Lsup()
49    
50  print "mypde=LinearPDE( A=[[1.,0.8],[0.4,1.]], D=alpha, Y=alpha, domain=mydomain)"  print "mypde=LinearPDE( A=[[1.,0.8],[0.4,1.]], D=alpha, Y=alpha, domain=mydomain)"
51  mypde=LinearPDE(A=[[1.,0.8],[0.4,1.]],D=alpha,Y=alpha,domain=mydomain)  mypde=LinearPDE(mydomain)
52  mypde.setDebugOn()  mypde.setDebugOn()
53  #mypde=LinearPDE(A=[[1.,0.],[0.,1.]],D=alpha,Y=alpha,domain=mydomain)  mypde.setValue(A=[[1.,0.8],[0.4,1.]],D=alpha,Y=alpha)
54  mypde.getOperator().saveMM("t.msh")  mypde.getOperator().saveMM("t.msh")
55    
56  # generate a test solution 1  # generate a test solution 1
# Line 87  print "msk=x[0].whereZero()+(x[0]-1.).wh Line 87  print "msk=x[0].whereZero()+(x[0]-1.).wh
87  msk=x[0].whereZero()+(x[0]-1.).whereZero()  msk=x[0].whereZero()+(x[0]-1.).whereZero()
88    
89  print "mypde=LinearPDE(A=[[1.,0.],[0.,1.]],q=msk,r=u_ex)"  print "mypde=LinearPDE(A=[[1.,0.],[0.,1.]],q=msk,r=u_ex)"
90  mypde=LinearPDE(A=[[1.,0.],[0.,1.]],q=msk,r=u_ex)  mypde=LinearPDE(mydomain)
91    mypde.setValue(A=[[1.,0.],[0.,1.]],q=msk,r=u_ex)
92  mypde.setDebugOn()  mypde.setDebugOn()
93    
94  # generate a test solution 2  # generate a test solution 2
# Line 98  print "----------------------------" Line 99  print "----------------------------"
99  print "\nDirect Solver (2)=>"  print "\nDirect Solver (2)=>"
100    
101  # mypde.setSymmetryOn() : is not woking yet!  # mypde.setSymmetryOn() : is not woking yet!
102  mypde.setSolverMethod(DIRECT)  mypde.setSolverMethod(mypde.DIRECT)
103  u_d=mypde.getSolution()  u_d=mypde.getSolution()
104    
105  print "\nIterative Solver (2)=>"  print "\nIterative Solver (2)=>"

Legend:
Removed from v.102  
changed lines
  Added in v.104

  ViewVC Help
Powered by ViewVC 1.1.26