/[escript]/trunk/esys2/modellib/py_src/flow.py
ViewVC logotype

Diff of /trunk/esys2/modellib/py_src/flow.py

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

revision 147 by jgs, Fri Aug 12 01:45:47 2005 UTC revision 148 by jgs, Tue Aug 23 01:24:31 2005 UTC
# Line 84  class SteadyIncompressibleFlow(Model): Line 84  class SteadyIncompressibleFlow(Model):
84            self.__pde.setTolerance(self.rel_tol/10.)            self.__pde.setTolerance(self.rel_tol/10.)
85            self.velocity=self.__pde.getSolution()            self.velocity=self.__pde.getSolution()
86            update=penalty*div(self.velocity)            update=penalty*div(self.velocity)
           print "f=",inf(update),sup(update)  
87            self.pressure=self.pressure-update            self.pressure=self.pressure-update
88            self.__diff,diff_old=Lsup(update),self.__diff            self.__diff,diff_old=Lsup(update),self.__diff
89            self.__iter+=1            self.__iter+=1

Legend:
Removed from v.147  
changed lines
  Added in v.148

  ViewVC Help
Powered by ViewVC 1.1.26