/[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 1639 by gross, Mon Jul 14 08:55:25 2008 UTC
# Line 490  class LinearPDE(object): Line 490  class LinearPDE(object):
490     AMG= 22     AMG= 22
491     RILU = 23     RILU = 23
492     TRILINOS = 24     TRILINOS = 24
493       NONLINEAR_GMRES = 25
494    
495     SMALL_TOLERANCE=1.e-13     SMALL_TOLERANCE=1.e-13
496     __PACKAGE_KEY="package"     __PACKAGE_KEY="package"
# Line 2422  class TransportPDE(object): Line 2423  class TransportPDE(object):
2423       def solve(self,dt,**kwarg):       def solve(self,dt,**kwarg):
2424             if self.__useSUPG:             if self.__useSUPG:
2425                  if self.__reassemble:                  if self.__reassemble:
2426                      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)
2427                      self.__reassemble=False                      self.__reassemble=False
2428                  dt2=self.getSafeTimeStepSize()                  dt2=self.getSafeTimeStepSize()
2429                  nn=max(math.ceil(dt/self.getSafeTimeStepSize()),1.)                  nn=max(math.ceil(dt/self.getSafeTimeStepSize()),1.)
# Line 2447  class TransportPDE(object): Line 2448  class TransportPDE(object):
2448              X=0              X=0
2449              Y=self.__M*u0              Y=self.__M*u0
2450              X=0              X=0
2451                self.__pde.setValue(r=u0)
2452              if not self.__A.isEmpty():              if not self.__A.isEmpty():
2453                 X=X+dt*util.matrixmult(self.__A,g)                 X=X+dt*util.matrixmult(self.__A,g)
2454              if not self.__B.isEmpty():              if not self.__B.isEmpty():

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

  ViewVC Help
Powered by ViewVC 1.1.26