/[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 1513 by gross, Tue Apr 15 08:47:57 2008 UTC revision 1552 by gross, Thu May 8 08:52:41 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)                      self.__pde.setValue(D=self.__M,d=self.__d,d_contact=self.__d_contact,q=self.__q) # ,r=self.__r)
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 2447  class TransportPDE(object): Line 2447  class TransportPDE(object):
2447              X=0              X=0
2448              Y=self.__M*u0              Y=self.__M*u0
2449              X=0              X=0
2450                self.__pde.setValue(r=u0)
2451              if not self.__A.isEmpty():              if not self.__A.isEmpty():
2452                 X=X+dt*util.matrixmult(self.__A,g)                 X=X+dt*util.matrixmult(self.__A,g)
2453              if not self.__B.isEmpty():              if not self.__B.isEmpty():

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

  ViewVC Help
Powered by ViewVC 1.1.26