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

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

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

revision 877 by gross, Wed Oct 25 03:06:58 2006 UTC revision 878 by gross, Wed Oct 25 03:42:56 2006 UTC
# Line 102  class StokesProblem(SaddlePointProblem): Line 102  class StokesProblem(SaddlePointProblem):
102           dp=self.__pde_p.getSolution()           dp=self.__pde_p.getSolution()
103           return  dp           return  dp
104    
105  NE=1  NE=50
106  dom=Rectangle(NE,NE,order=2)  dom=Rectangle(NE,NE,order=2)
107  # prop=SimpleStokesProblem(dom)  # prop=SimpleStokesProblem(dom)
108  prop=StokesProblem(dom)  prop=StokesProblem(dom)
# Line 118  u,p=prop.solve(u0,p0,tolerance=0.01) Line 118  u,p=prop.solve(u0,p0,tolerance=0.01)
118    
119  eta=whereNegative(x[1]-0.5)*1.e6+whereNonNegative(x[1]-0.5)  eta=whereNegative(x[1]-0.5)*1.e6+whereNonNegative(x[1]-0.5)
120  prop.initialize(fixed_u_mask=mask,eta=eta)  prop.initialize(fixed_u_mask=mask,eta=eta)
121  u,p=prop.solve(u0,p0,tolerance=0.01,tolerance_u=0.1,relaxation=1.)  u,p=prop.solve(u0,p0,tolerance=0.01,tolerance_u=0.1,accepted_reduction=0.8)
122  saveVTK("stokes.xml",u=u,p=p,m=mask,u0=u0)  saveVTK("stokes.xml",u=u,p=p,m=mask,u0=u0)
123                        
124  # vim: expandtab shiftwidth=4:  # vim: expandtab shiftwidth=4:

Legend:
Removed from v.877  
changed lines
  Added in v.878

  ViewVC Help
Powered by ViewVC 1.1.26