/[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 791 by bcumming, Thu Jul 27 00:37:57 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    
# Line 1538  class LinearPDE(object): Line 1538  class LinearPDE(object):
1538                   D=self.getCoefficientOfGeneralPDE("D")                   D=self.getCoefficientOfGeneralPDE("D")
1539                   if not D.isEmpty():                   if not D.isEmpty():
1540                       if self.getNumSolutions()>1:                       if self.getNumSolutions()>1:
1541                          D_times_e=util.matrixmult(D,numarray.ones((self.getNumSolutions(),)))                          #D_times_e=util.matrixmult(D,numarray.ones((self.getNumSolutions(),)))
1542                            D_times_e=util.matrix_mult(D,numarray.ones((self.getNumSolutions(),)))
1543                       else:                       else:
1544                          D_times_e=D                          D_times_e=D
1545                   else:                   else:
# Line 1546  class LinearPDE(object): Line 1547  class LinearPDE(object):
1547                   d=self.getCoefficientOfGeneralPDE("d")                   d=self.getCoefficientOfGeneralPDE("d")
1548                   if not d.isEmpty():                   if not d.isEmpty():
1549                       if self.getNumSolutions()>1:                       if self.getNumSolutions()>1:
1550                          d_times_e=util.matrixmult(d,numarray.ones((self.getNumSolutions(),)))                          #d_times_e=util.matrixmult(d,numarray.ones((self.getNumSolutions(),)))
1551                            d_times_e=util.matrix_mult(d,numarray.ones((self.getNumSolutions(),)))
1552                       else:                       else:
1553                          d_times_e=d                          d_times_e=d
1554                   else:                   else:

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

  ViewVC Help
Powered by ViewVC 1.1.26