/[escript]/trunk/escript/py_src/flows.py
ViewVC logotype

Diff of /trunk/escript/py_src/flows.py

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

revision 1562 by gross, Wed May 21 13:04:40 2008 UTC revision 1639 by gross, Mon Jul 14 08:55:25 2008 UTC
# Line 63  class StokesProblemCartesian(Homogeneous Line 63  class StokesProblemCartesian(Homogeneous
63           self.vol=util.integrate(1.,Function(self.domain))           self.vol=util.integrate(1.,Function(self.domain))
64           self.__pde_u=LinearPDE(domain,numEquations=self.domain.getDim(),numSolutions=self.domain.getDim())           self.__pde_u=LinearPDE(domain,numEquations=self.domain.getDim(),numSolutions=self.domain.getDim())
65           self.__pde_u.setSymmetryOn()           self.__pde_u.setSymmetryOn()
66           self.__pde_u.setSolverMethod(preconditioner=LinearPDE.ILU0)           # self.__pde_u.setSolverMethod(preconditioner=LinearPDE.ILU0)
67                            
68           self.__pde_prec=LinearPDE(domain)           self.__pde_prec=LinearPDE(domain)
69           self.__pde_prec.setReducedOrderOn()           self.__pde_prec.setReducedOrderOn()
# Line 110  class StokesProblemCartesian(Homogeneous Line 110  class StokesProblemCartesian(Homogeneous
110        def getStress(self,u):        def getStress(self,u):
111           mg=util.grad(u)           mg=util.grad(u)
112           return 2.*self.eta*util.symmetric(mg)           return 2.*self.eta*util.symmetric(mg)
113          def getEtaEffective(self):
114             return self.eta
115    
116        def solve_A(self,u,p):        def solve_A(self,u,p):
117           """           """

Legend:
Removed from v.1562  
changed lines
  Added in v.1639

  ViewVC Help
Powered by ViewVC 1.1.26