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

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

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

revision 1423 by gross, Wed Feb 27 04:28:48 2008 UTC revision 1513 by gross, Tue Apr 15 08:47:57 2008 UTC
# Line 2422  class TransportPDE(object): Line 2422  class TransportPDE(object):
2422       def solve(self,dt,**kwarg):       def solve(self,dt,**kwarg):
2423             if self.__useSUPG:             if self.__useSUPG:
2424                  if self.__reassemble:                  if self.__reassemble:
2425                      self.__pde.setValue(D=self.__M,d=self.__d,d_contact=self.__d_contact,q=self.__q,r=self.__r)                      self.__pde.setValue(D=self.__M,d=self.__d,d_contact=self.__d_contact,q=self.__q)
2426                      self.__reassemble=False                      self.__reassemble=False
2427                  dt2=self.getSafeTimeStepSize()                  dt2=self.getSafeTimeStepSize()
2428                  nn=max(math.ceil(dt/self.getSafeTimeStepSize()),1.)                  nn=max(math.ceil(dt/self.getSafeTimeStepSize()),1.)
# Line 2464  class TransportPDE(object): Line 2464  class TransportPDE(object):
2464                 self.__pde.setValue(y=dt*self.__y)                 self.__pde.setValue(y=dt*self.__y)
2465              if not self.__y_contact.isEmpty():              if not self.__y_contact.isEmpty():
2466                 self.__pde.setValue(y=dt*self.__y_contact)                 self.__pde.setValue(y=dt*self.__y_contact)
2467                  self.__pde.setValue(r=u0)

Legend:
Removed from v.1423  
changed lines
  Added in v.1513

  ViewVC Help
Powered by ViewVC 1.1.26