/[ascend]/trunk/test.py
ViewVC logotype

Diff of /trunk/test.py

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

revision 956 by johnpye, Sat Dec 9 15:38:05 2006 UTC revision 961 by johnpye, Mon Dec 11 14:04:13 2006 UTC
# Line 87  class Ascend(unittest.TestCase): Line 87  class Ascend(unittest.TestCase):
87      def testLSODE(self):      def testLSODE(self):
88          self._testIntegrator('LSODE')          self._testIntegrator('LSODE')
89    
90        def testlotka(self):
91            self.L.load('johnpye/lotka.a4c')
92            M = self.L.findType('lotka').getSimulation('sim')
93            M.solve(ascpy.Solver("QRSlv"),ascpy.SolverReporter())  
94            I = ascpy.Integrator(M)
95            I.setEngine('LSODE')
96            I.setReporter(ascpy.IntegratorReporterConsole(I))
97            I.setLinearTimesteps(ascpy.Units("s"), 0, 200, 5);
98            I.analyse()
99            I.solve()
100            assert I.getNumObservedVars() == 3;
101            assert abs(float(M.sim.R) - 832) < 1.0
102            assert abs(float(M.sim.F) - 21.36) < 0.1
103            
104    
105      def testIDA(self):      def testIDA(self):
106          self._testIntegrator('IDA')          self._testIntegrator('IDA')
107    
108    
109      def testIDAparameters(self):      def testIDAparameters(self):
110          self.L.load('johnpye/shm.a4c')          self.L.load('johnpye/shm.a4c')
111          M = self.L.findType('shm').getSimulation('sim')          M = self.L.findType('shm').getSimulation('sim')
# Line 98  class Ascend(unittest.TestCase): Line 114  class Ascend(unittest.TestCase):
114          P = I.getParameters()          P = I.getParameters()
115          for p in P:          for p in P:
116              print p.getName(),"=",p.getValue()              print p.getName(),"=",p.getValue()
117          assert len(P)==5          assert len(P)==7
118          assert P[0].isStr()          assert P[0].isStr()
119          assert P[0].getName()=="linsolver"          assert P[0].getName()=="linsolver"
120          assert P[0].getValue()=='DENSE'          assert P[0].getValue()=='SPGMR'
121          assert P[1].getName()=="autodiff"          assert P[1].getName()=="autodiff"
122          assert P[1].getValue()==True          assert P[1].getValue()==True
123          assert P[4].getName()=="atolvect"          assert P[5].getName()=="atolvect"
124          assert P[4].getBoolValue() == True          assert P[5].getBoolValue() == True
125          P[1].setBoolValue(False)          P[1].setBoolValue(False)
126          assert P[1].getBoolValue()==False          assert P[1].getBoolValue()==False
127          I.setParameters(P)          I.setParameters(P)

Legend:
Removed from v.956  
changed lines
  Added in v.961

john.pye@anu.edu.au
ViewVC Help
Powered by ViewVC 1.1.22