1722 |
|
|
1723 |
class Test_TransportPDE(Test_linearPDEs): |
class Test_TransportPDE(Test_linearPDEs): |
1724 |
N=4 |
N=4 |
1725 |
def test_init_theta(self): |
def test_init_useBackwardEuler(self): |
1726 |
mypde=TransportPDE(self.domain,debug=self.DEBUG, theta=1.) |
mypde=TransportPDE(self.domain,debug=self.DEBUG, useBackwardEuler=True) |
1727 |
self.failUnless(mypde.getTheta()==1.,'wrong theta') |
self.failUnless(mypde.useBackwardEuler()==True,'backward Euler should be used') |
1728 |
def test_init_invalid_theta(self): |
def test_init_donntUseBackwardEuler(self): |
1729 |
self.failUnlessRaises(ValueError,TransportPDE,self.domain,debug=self.DEBUG, theta=-1.) |
mypde=TransportPDE(self.domain,debug=self.DEBUG, useBackwardEuler=False) |
1730 |
|
self.failUnless(mypde.useBackwardEuler()==False,'backward Euler should not be used') |
1731 |
def test_setCoefficient_WithWrongName(self): |
def test_setCoefficient_WithWrongName(self): |
1732 |
mypde=TransportPDE(self.domain,debug=self.DEBUG) |
mypde=TransportPDE(self.domain,debug=self.DEBUG) |
1733 |
self.failUnlessRaises(IllegalCoefficient,mypde.setValue, ROMA=Vector(0.,FunctionOnBoundary(self.domain))) |
self.failUnlessRaises(IllegalCoefficient,mypde.setValue, ROMA=Vector(0.,FunctionOnBoundary(self.domain))) |