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

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

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

revision 2549 by jfenwick, Mon Jul 20 06:43:47 2009 UTC revision 2620 by gross, Thu Aug 20 06:24:00 2009 UTC
# Line 465  class StokesProblemCartesian(Homogeneous Line 465  class StokesProblemCartesian(Homogeneous
465           self.getSolverOptionsVelocity().setAbsoluteTolerance(0.)           self.getSolverOptionsVelocity().setAbsoluteTolerance(0.)
466                    
467    
468       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(), restoration_factor=0):
469          """          """
470          assigns values to the model parameters          assigns values to the model parameters
471    
# Line 488  class StokesProblemCartesian(Homogeneous Line 488  class StokesProblemCartesian(Homogeneous
488          for j in range(self.domain.getDim()):          for j in range(self.domain.getDim()):
489              A[i,j,j,i] += 1.              A[i,j,j,i] += 1.
490              A[i,j,i,j] += 1.              A[i,j,i,j] += 1.
491            n=self.domain.getNormal()
492      self.__pde_prec.setValue(D=1/self.eta)      self.__pde_prec.setValue(D=1/self.eta)
493          self.__pde_u.setValue(A=A*self.eta,q=fixed_u_mask)          self.__pde_u.setValue(A=A*self.eta,q=fixed_u_mask, d=restoration_factor*util.outer(n,n))
494          self.__f=f          self.__f=f
495          self.__surface_stress=surface_stress          self.__surface_stress=surface_stress
496          self.__stress=stress          self.__stress=stress

Legend:
Removed from v.2549  
changed lines
  Added in v.2620

  ViewVC Help
Powered by ViewVC 1.1.26