/[escript]/trunk/escript/py_src/flows.py
ViewVC logotype

Diff of /trunk/escript/py_src/flows.py

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 2107 by gross, Wed Nov 26 08:13:00 2008 UTC revision 2108 by gross, Fri Nov 28 05:09:23 2008 UTC
# Line 201  class DarcyFlow(object): Line 201  class DarcyFlow(object):
201           if not ATOL>0:           if not ATOL>0:
202               raise ValueError,"Negative absolute tolerance (rtol = %e, norm right hand side =%, atol =%e)."%(rtol, util.sqrt(norm_rhs), atol)               raise ValueError,"Negative absolute tolerance (rtol = %e, norm right hand side =%, atol =%e)."%(rtol, util.sqrt(norm_rhs), atol)
203           rhs=ArithmeticTuple(g,dv)           rhs=ArithmeticTuple(g,dv)
204           dp,r=PCG(rhs,self.__Aprod_PCG,self.__Msolve_PCG,self.__inner_PCG,atol=ATOL, rtol=0.,iter_max=max_iter, x=p0, verbose=self.verbose, initial_guess=False)           dp,r=PCG(rhs,self.__Aprod_PCG,self.__Msolve_PCG,self.__inner_PCG,atol=ATOL, rtol=0.,iter_max=max_iter, x=p0-p2, verbose=self.verbose, initial_guess=True)
205           return u2+r[1],p2+dp           return u2+r[1],p2+dp
206                    
207      def __Aprod_PCG(self,p):      def __Aprod_PCG(self,p):
Legend:
Removed from v.2107  
changed lines
  Added in v.2108
  ViewVC Help
Powered by ViewVC 1.1.26