21 |
|
|
22 |
from esys.finley import Rectangle |
from esys.finley import Rectangle |
23 |
class Test_LinearPDEOnFinley(Test_LinearPDE): |
class Test_LinearPDEOnFinley(Test_LinearPDE): |
24 |
RES_TOL=1-8 |
RES_TOL=1.e-8 |
25 |
def setUp(self): |
def setUp(self): |
26 |
self.domain = Rectangle(10,10,2) |
self.domain = Rectangle(10,10,2) |
27 |
def tearDown(self): |
def tearDown(self): |
329 |
p=NoPDE(self.domain) |
p=NoPDE(self.domain) |
330 |
x=self.domain.getX() |
x=self.domain.getX() |
331 |
msk=whereZero(x[0])*[1.,0.] |
msk=whereZero(x[0])*[1.,0.] |
332 |
p.setValue(D=x[:2],Y=2*x[:2],q=msk,r=2.) |
p.setValue(D=x[:2]+1,Y=2*(x[:2]+1),q=msk,r=2.) |
333 |
u=p.getSolution() |
u=p.getSolution() |
334 |
u_ex=2. |
u_ex=2. |
335 |
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") |