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

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

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

revision 2562 by jfenwick, Mon Jul 20 06:43:47 2009 UTC revision 2563 by gross, Tue Jul 28 03:50:45 2009 UTC
# Line 696  class Test_Mountains3D(unittest.TestCase Line 696  class Test_Mountains3D(unittest.TestCase
696         v[1]=a1*cos(pi*n0*x[0])* sin(pi*n1*x[1])* cos(pi*n2*x[2])         v[1]=a1*cos(pi*n0*x[0])* sin(pi*n1*x[1])* cos(pi*n2*x[2])
697         v[2]=a2*cos(pi*n0*x[0])* cos(pi*n1*x[1])* sin(pi*n2*x[2])         v[2]=a2*cos(pi*n0*x[0])* cos(pi*n1*x[1])* sin(pi*n2*x[2])
698    
699         H_t=Scalar(0.0, Solution(self.domain))         mts=Mountains(self.domain,eps=EPS)
700         mts=Mountains(self.domain,v,eps=EPS,z=1)         mts.setVelocity(v)
701         u,Z=mts.update(u=v,H_t=H_t)         Z=mts.update()
702                
703         error_int=integrate(Z)         error_int=integrate(Z)
704         self.failUnless(error_int<self.TOL, "Boundary intergral is too large.")         self.failUnless(error_int<self.TOL, "Boundary intergral is too large.")
# Line 723  class Test_Mountains2D(unittest.TestCase Line 723  class Test_Mountains2D(unittest.TestCase
723         v[1]=a1*cos(pi*n0*x[0])* sin(pi*n1*x[1])         v[1]=a1*cos(pi*n0*x[0])* sin(pi*n1*x[1])
724                
725         H_t=Scalar(0.0, Solution(self.domain))         H_t=Scalar(0.0, Solution(self.domain))
726         mts=Mountains(self.domain,v,eps=EPS,z=1)         mts=Mountains(self.domain,eps=EPS)
727         u,Z=mts.update(u=v,H_t=H_t)         mts.setVelocity(v)
728           Z=mts.update()
729                
730         error_int=integrate(Z)         error_int=integrate(Z)
731         self.failUnless(error_int<self.TOL, "Boundary intergral is too large.")         self.failUnless(error_int<self.TOL, "Boundary intergral is too large.")
Legend:
Removed from v.2562  
changed lines
  Added in v.2563
  ViewVC Help
Powered by ViewVC 1.1.26