/[escript]/trunk/escript/test/python/test_pdetools.py
ViewVC logotype

Diff of /trunk/escript/test/python/test_pdetools.py

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

revision 798 by gross, Fri Aug 4 01:05:36 2006 UTC revision 854 by gross, Thu Sep 21 05:29:42 2006 UTC
# Line 270  class Test_pdetools(unittest.TestCase): Line 270  class Test_pdetools(unittest.TestCase):
270        p=NoPDE(self.domain)        p=NoPDE(self.domain)
271        x=self.domain.getX()        x=self.domain.getX()
272        msk=whereZero(x[0])*[1.,0.]        msk=whereZero(x[0])*[1.,0.]
273        p.setValue(D=1.,q=msk,r=2.)        p.setValue(D=numarray.ones([2]),q=msk,r=2.)
274        u=p.getSolution()        u=p.getSolution()
275        u_ex=msk*2.        u_ex=msk*2.
276        self.failUnless(Lsup(u_ex-u)<Lsup(u_ex)*self.RES_TOL,"value wrong")        self.failUnless(Lsup(u_ex-u)<Lsup(u_ex)*self.RES_TOL,"value wrong")
# Line 279  class Test_pdetools(unittest.TestCase): Line 279  class Test_pdetools(unittest.TestCase):
279        p=NoPDE(self.domain)        p=NoPDE(self.domain)
280        x=self.domain.getX()        x=self.domain.getX()
281        msk=whereZero(x[0])*[1.,0.]        msk=whereZero(x[0])*[1.,0.]
282        p.setValue(D=1.,Y=1.,q=msk)        p.setValue(D=numarray.ones([2]),Y=numarray.ones([2]),q=msk)
283        u=p.getSolution()        u=p.getSolution()
284        u_ex=(1.-msk)        u_ex=(1.-msk)
285        self.failUnless(Lsup(u_ex-u)<Lsup(u_ex)*self.RES_TOL,"value wrong")        self.failUnless(Lsup(u_ex-u)<Lsup(u_ex)*self.RES_TOL,"value wrong")
# Line 288  class Test_pdetools(unittest.TestCase): Line 288  class Test_pdetools(unittest.TestCase):
288        p=NoPDE(self.domain)        p=NoPDE(self.domain)
289        x=self.domain.getX()        x=self.domain.getX()
290        msk=whereZero(x[0])*[1.,0.]        msk=whereZero(x[0])*[1.,0.]
291        p.setValue(D=1.,Y=1.,q=msk,r=2.)        p.setValue(D=numarray.ones([2]),Y=numarray.ones([2]),q=msk,r=2.)
292        u=p.getSolution()        u=p.getSolution()
293        u_ex=(1.-msk)+msk*2.        u_ex=(1.-msk)+msk*2.
294        self.failUnless(Lsup(u_ex-u)<Lsup(u_ex)*self.RES_TOL,"value wrong")        self.failUnless(Lsup(u_ex-u)<Lsup(u_ex)*self.RES_TOL,"value wrong")

Legend:
Removed from v.798  
changed lines
  Added in v.854

  ViewVC Help
Powered by ViewVC 1.1.26