/[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 2386 by gross, Wed Apr 15 03:54:25 2009 UTC revision 2415 by gross, Wed May 13 02:48:39 2009 UTC
# Line 367  class StokesProblemCartesian(Homogeneous Line 367  class StokesProblemCartesian(Homogeneous
367           # self.__pde_prec.setSolverMethod(self.__pde_prec.LUMPING)           # self.__pde_prec.setSolverMethod(self.__pde_prec.LUMPING)
368           self.__pde_prec.setSymmetryOn()           self.__pde_prec.setSymmetryOn()
369    
370             self.__pde_proj=LinearPDE(domain)
371             self.__pde_proj.setReducedOrderOn()
372         self.__pde_proj.setValue(D=1)
373             self.__pde_proj.setSymmetryOn()
374    
375    
376       def initialize(self,f=Data(),fixed_u_mask=Data(),eta=1,surface_stress=Data(),stress=Data()):       def initialize(self,f=Data(),fixed_u_mask=Data(),eta=1,surface_stress=Data(),stress=Data()):
377          """          """
378          assigns values to the model parameters          assigns values to the model parameters
# Line 459  class StokesProblemCartesian(Homogeneous Line 465  class StokesProblemCartesian(Homogeneous
465          @rtype: equal to the type of p          @rtype: equal to the type of p
466          @note: boundary conditions on p should be zero!          @note: boundary conditions on p should be zero!
467          """          """
468          return util.sqrt(util.integrate(util.div(v)**2))          self.__pde_proj.setValue(Y=util.div(v))
469            self.__pde_prec.setTolerance(self.getSubProblemTolerance())
470            return util.sqrt(util.integrate(util.interpolate(self.__pde_proj.getSolution(),Function(self.domain))**2))
471    
472       def solve_AinvBt(self,p):       def solve_AinvBt(self,p):
473           """           """

Legend:
Removed from v.2386  
changed lines
  Added in v.2415

  ViewVC Help
Powered by ViewVC 1.1.26