/[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 1417 by gross, Mon Feb 25 04:45:48 2008 UTC revision 1423 by gross, Wed Feb 27 04:28:48 2008 UTC
# Line 2241  class TransportPDE(object): Line 2241  class TransportPDE(object):
2241          self.__trace=trace          self.__trace=trace
2242          self.__theta=theta          self.__theta=theta
2243          self.__matrix_type=0          self.__matrix_type=0
2244          self.__reduced=False          self.__reduced=True
2245          self.__reassemble=True          self.__reassemble=True
2246          if self.__useSUPG:          if self.__useSUPG:
2247             self.__pde=LinearPDE(domain,numEquations=num_equations,numSolutions=num_equations,debug=trace)             self.__pde=LinearPDE(domain,numEquations=num_equations,numSolutions=num_equations,debug=trace)
2248             self.__pde.setSymmetryOn()             self.__pde.setSymmetryOn()
2249               self.__pde.setReducedOrderOn()
2250          else:          else:
2251             self.__transport_problem=self.__getNewTransportProblem()             self.__transport_problem=self.__getNewTransportProblem()
2252          self.setTolerance()          self.setTolerance()
         self.setReducedOn()  
2253          self.__M=escript.Data()          self.__M=escript.Data()
2254          self.__A=escript.Data()          self.__A=escript.Data()
2255          self.__B=escript.Data()          self.__B=escript.Data()
# Line 2415  class TransportPDE(object): Line 2415  class TransportPDE(object):
2415    
2416       def setInitialSolution(self,u):       def setInitialSolution(self,u):
2417               if self.__useSUPG:               if self.__useSUPG:
2418                   self.__u=u                   self.__u=util.interpolate(u,self.getFunctionSpace())
2419               else:               else:
2420                   self.__transport_problem.setInitialValue(util.interpolate(u,self.getFunctionSpace()))                   self.__transport_problem.setInitialValue(util.interpolate(u,self.getFunctionSpace()))
2421    

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

  ViewVC Help
Powered by ViewVC 1.1.26