/[escript]/trunk/finley/test/python/run_darcy.py
ViewVC logotype

Diff of /trunk/finley/test/python/run_darcy.py

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

revision 3792 by jfenwick, Wed Jan 18 04:19:33 2012 UTC revision 3793 by gross, Wed Feb 1 07:39:43 2012 UTC
# Line 243  class Test_Darcy(unittest.TestCase): Line 243  class Test_Darcy(unittest.TestCase):
243          p=p_ref*mp          p=p_ref*mp
244          u=u_ref*mv          u=u_ref*mv
245          df=DarcyFlow(self.dom, verbose=VERBOSE, solver=self.SOLVER)          df=DarcyFlow(self.dom, verbose=VERBOSE, solver=self.SOLVER)
246            #df.getSolverOptionsPressure().setVerbosityOn()
247          df.setValue(g=f,          df.setValue(g=f,
248                        location_of_fixed_pressure=mp,                        location_of_fixed_pressure=mp,
249                        location_of_fixed_flux=mv,                        location_of_fixed_flux=mv,
250                        permeability=Scalar(k,Function(self.dom)))                        permeability=Scalar(k,Function(self.dom)))
251          v,p=df.solve(u,p)          v,p=df.solve(u,p)
   
   
252          self.assertTrue(Lsup(p-p_ref)<self.TEST_TOL*Lsup(p_ref), "pressure error too big.")          self.assertTrue(Lsup(p-p_ref)<self.TEST_TOL*Lsup(p_ref), "pressure error too big.")
253          self.assertTrue(Lsup(v-u_ref)<self.TEST_TOL*Lsup(u_ref), "flux error too big.")          self.assertTrue(Lsup(v-u_ref)<self.TEST_TOL*Lsup(u_ref), "flux error too big.")
254    
# Line 357  if __name__ == '__main__': Line 356  if __name__ == '__main__':
356     suite.addTest(unittest.makeSuite(Test_Darcy3D_SMOOTH))     suite.addTest(unittest.makeSuite(Test_Darcy3D_SMOOTH))
357     suite.addTest(unittest.makeSuite(Test_Darcy3D_POST))     suite.addTest(unittest.makeSuite(Test_Darcy3D_POST))
358     suite.addTest(unittest.makeSuite(Test_Darcy3D_EVAL))     suite.addTest(unittest.makeSuite(Test_Darcy3D_EVAL))
   
   
359     #suite.addTest(Test_Darcy2D_SMOOTH("testConstF_FreeBottom_largeK"))     #suite.addTest(Test_Darcy2D_SMOOTH("testConstF_FreeBottom_largeK"))
360     s=unittest.TextTestRunner(verbosity=2).run(suite)     s=unittest.TextTestRunner(verbosity=2).run(suite)
361     if not s.wasSuccessful(): sys.exit(1)     if not s.wasSuccessful(): sys.exit(1)

Legend:
Removed from v.3792  
changed lines
  Added in v.3793

  ViewVC Help
Powered by ViewVC 1.1.26