/[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 720 by gross, Thu Apr 27 10:16:05 2006 UTC revision 772 by gross, Fri Jul 7 05:47:28 2006 UTC
# Line 664  class LinearPDE(object): Line 664  class LinearPDE(object):
664       @rtype: L{Data<escript.Data>}       @rtype: L{Data<escript.Data>}
665       """       """
666       if u==None:       if u==None:
667            return self.getOperator()*self.getSolution()          return self.getOperator()*self.getSolution()
668       else:       else:
669          self.getOperator()*escript.Data(u,self.getFunctionSpaceForSolution())          return self.getOperator()*escript.Data(u,self.getFunctionSpaceForSolution())
670    
671     def getResidual(self,u=None):     def getResidual(self,u=None):
672       """       """
# Line 890  class LinearPDE(object): Line 890  class LinearPDE(object):
890         """         """
891         if package==None: package=self.DEFAULT         if package==None: package=self.DEFAULT
892         if not package==self.getSolverPackage():         if not package==self.getSolverPackage():
893             self.__solver_method=solver             self.__solver_package=package
894             self.__checkMatrixType()             self.__checkMatrixType()
895             self.trace("New solver is %s"%self.getSolverMethodName())             self.trace("New solver is %s"%self.getSolverMethodName())
896    
# Line 1513  class LinearPDE(object): Line 1513  class LinearPDE(object):
1513           r=self.getCoefficientOfGeneralPDE("r")           r=self.getCoefficientOfGeneralPDE("r")
1514           homogeneous_constraint=True           homogeneous_constraint=True
1515           if not q.isEmpty() and not r.isEmpty():           if not q.isEmpty() and not r.isEmpty():
1516               if util.Lsup(q*r)>=1.e-13*util.Lsup(r):               if util.Lsup(q*r)>0.:
1517                 self.trace("Inhomogeneous constraint detected.")                 self.trace("Inhomogeneous constraint detected.")
1518                 self.__invalidateSystem()                 self.__invalidateSystem()
1519    

Legend:
Removed from v.720  
changed lines
  Added in v.772

  ViewVC Help
Powered by ViewVC 1.1.26