/[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 2354 by gross, Wed Apr 1 04:01:58 2009 UTC revision 2370 by gross, Mon Apr 6 06:41:49 2009 UTC
# Line 333  class DarcyFlow(object): Line 333  class DarcyFlow(object):
333            Qdp=self.__Q(dp)            Qdp=self.__Q(dp)
334            self.__pde_v.setValue(Y=-Qdp,X=Data(), r=Data())            self.__pde_v.setValue(Y=-Qdp,X=Data(), r=Data())
335            du=self.__pde_v.getSolution(verbose=self.show_details, iter_max = 100000)            du=self.__pde_v.getSolution(verbose=self.show_details, iter_max = 100000)
336              # self.__pde_v.getOperator().saveMM("proj.mm")
337            return Qdp+du            return Qdp+du
338      def __inner_GMRES(self,r,s):      def __inner_GMRES(self,r,s):
339           return util.integrate(util.inner(r,s))           return util.integrate(util.inner(r,s))
# Line 344  class DarcyFlow(object): Line 345  class DarcyFlow(object):
345            self.__pde_p.setTolerance(self.getSubProblemTolerance())            self.__pde_p.setTolerance(self.getSubProblemTolerance())
346            if self.show_details: print "DarcyFlux: Applying preconditioner"            if self.show_details: print "DarcyFlux: Applying preconditioner"
347            self.__pde_p.setValue(X=util.transposed_tensor_mult(self.__permeability,r), Y=Data(), r=Data())            self.__pde_p.setValue(X=util.transposed_tensor_mult(self.__permeability,r), Y=Data(), r=Data())
348              # self.__pde_p.getOperator().saveMM("prec.mm")
349            return self.__pde_p.getSolution(verbose=self.show_details, iter_max = 100000)            return self.__pde_p.getSolution(verbose=self.show_details, iter_max = 100000)
350    
351      def getFlux(self,p=None, fixed_flux=Data(), show_details=False):      def getFlux(self,p=None, fixed_flux=Data(), show_details=False):

Legend:
Removed from v.2354  
changed lines
  Added in v.2370

  ViewVC Help
Powered by ViewVC 1.1.26