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

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

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

revision 2485 by gross, Tue Jun 16 06:32:15 2009 UTC revision 2486 by gross, Tue Jun 23 03:38:54 2009 UTC
# Line 62  class DarcyFlow(object): Line 62  class DarcyFlow(object):
62          if weight == None:          if weight == None:
63             s=self.domain.getSize()             s=self.domain.getSize()
64             self.__l=(3.*util.longestEdge(self.domain)*s/util.sup(s))**2             self.__l=(3.*util.longestEdge(self.domain)*s/util.sup(s))**2
65               # self.__l=(3.*util.longestEdge(self.domain))**2
66               # self.__l=(0.1*util.longestEdge(self.domain)*s/util.sup(s))**2
67          else:          else:
68             self.__l=weight             self.__l=weight
69          self.__pde_v=LinearPDESystem(domain)          self.__pde_v=LinearPDESystem(domain)
# Line 296  class DarcyFlow(object): Line 298  class DarcyFlow(object):
298                 if self.verbose:                 if self.verbose:
299                      print "DarcyFlux: L2 norm of v = %e."%norm_v                      print "DarcyFlux: L2 norm of v = %e."%norm_v
300                      print "DarcyFlux: L2 norm of k.grad(p) = %e."%norm_Qp                      print "DarcyFlux: L2 norm of k.grad(p) = %e."%norm_Qp
301                        print "DarcyFlux: L2 defect u = %e."%(util.integrate(util.length(self.__g-util.interpolate(v,Function(self.domain))-Qp)**2)**(0.5),)
302                        print "DarcyFlux: L2 defect div(v) = %e."%(util.integrate((self.__f-util.div(v))**2)**(0.5),)
303                      print "DarcyFlux: absolute tolerance ATOL = %e."%ATOL                      print "DarcyFlux: absolute tolerance ATOL = %e."%ATOL
304                 if norm_r == None or norm_r>ATOL:                 if norm_r == None or norm_r>ATOL:
305                     if num_corrections>max_num_corrections:                     if num_corrections>max_num_corrections:

Legend:
Removed from v.2485  
changed lines
  Added in v.2486

  ViewVC Help
Powered by ViewVC 1.1.26