1 |
johnpye |
1008 |
#!/usr/bin/env python |
2 |
johnpye |
1102 |
# ASCEND modelling environment |
3 |
|
|
# Copyright (C) 2006 Carnegie Mellon University |
4 |
|
|
# |
5 |
|
|
# This program is free software; you can redistribute it and/or modify |
6 |
|
|
# it under the terms of the GNU General Public License as published by |
7 |
|
|
# the Free Software Foundation; either version 2, or (at your option) |
8 |
|
|
# any later version. |
9 |
|
|
# |
10 |
|
|
# This program is distributed in the hope that it will be useful, |
11 |
|
|
# but WITHOUT ANY WARRANTY; without even the implied warranty of |
12 |
|
|
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
13 |
|
|
# GNU General Public License for more details. |
14 |
|
|
# |
15 |
|
|
# You should have received a copy of the GNU General Public License |
16 |
|
|
# along with this program; if not, write to the Free Software |
17 |
|
|
# Foundation, Inc., 59 Temple Place - Suite 330, |
18 |
|
|
# Boston, MA 02111-1307, USA. |
19 |
|
|
|
20 |
|
|
# This script gives a test suite for the high-level interface of ASCEND via |
21 |
|
|
# Python. It is also planned to be a wrapper for the CUnit test suite, although |
22 |
|
|
# this is still experimental. |
23 |
|
|
|
24 |
johnpye |
669 |
import unittest |
25 |
johnpye |
1104 |
import os, sys |
26 |
johnpye |
1098 |
import math |
27 |
|
|
import atexit |
28 |
johnpye |
1091 |
|
29 |
|
|
import platform |
30 |
johnpye |
1028 |
if platform.system() != "Windows": |
31 |
|
|
import dl |
32 |
|
|
sys.setdlopenflags(dl.RTLD_GLOBAL|dl.RTLD_NOW) |
33 |
|
|
|
34 |
johnpye |
956 |
class Ascend(unittest.TestCase): |
35 |
johnpye |
669 |
|
36 |
johnpye |
933 |
def setUp(self): |
37 |
|
|
import ascpy |
38 |
johnpye |
1118 |
self.L = ascpy.Library() |
39 |
johnpye |
933 |
|
40 |
|
|
def tearDown(self): |
41 |
|
|
self.L.clear() |
42 |
|
|
del self.L |
43 |
|
|
|
44 |
johnpye |
1024 |
class AscendSelfTester(Ascend): |
45 |
|
|
|
46 |
|
|
def _run(self,modelname,solvername="QRSlv",filename=None): |
47 |
|
|
if filename==None: |
48 |
|
|
filename = 'johnpye/%s.a4c' % modelname |
49 |
|
|
self.L.load(filename) |
50 |
|
|
T = self.L.findType(modelname) |
51 |
|
|
M = T.getSimulation('sim') |
52 |
|
|
M.build() |
53 |
|
|
M.solve(ascpy.Solver(solvername),ascpy.SolverReporter()) |
54 |
|
|
M.run(T.getMethod('self_test')) |
55 |
|
|
return M |
56 |
|
|
|
57 |
johnpye |
966 |
class TestCompiler(Ascend): |
58 |
|
|
|
59 |
johnpye |
941 |
def testloading(self): |
60 |
|
|
pass |
61 |
|
|
|
62 |
|
|
def testsystema4l(self): |
63 |
|
|
self.L.load('system.a4l') |
64 |
|
|
|
65 |
|
|
def testatomsa4l(self): |
66 |
|
|
self.L.load('atoms.a4l') |
67 |
|
|
|
68 |
johnpye |
1024 |
class TestSolver(AscendSelfTester): |
69 |
johnpye |
966 |
|
70 |
|
|
def testlog10(self): |
71 |
|
|
self._run('testlog10') |
72 |
|
|
|
73 |
johnpye |
1156 |
def testrootsofpoly(self): |
74 |
|
|
self._run('roots_of_poly',filename="roots_of_poly.a4c") |
75 |
|
|
|
76 |
johnpye |
1152 |
def testcollapsingcan(self): |
77 |
|
|
self._run('collapsingcan',filename="collapsingcan.a4c") |
78 |
|
|
|
79 |
johnpye |
1165 |
def testdistancecalc(self): |
80 |
|
|
self._run('distance_calc',filename="distance_calc.a4c") |
81 |
|
|
|
82 |
johnpye |
966 |
def testconopt(self): |
83 |
|
|
self._run('testconopt',"CONOPT") |
84 |
|
|
|
85 |
|
|
def testcmslv2(self): |
86 |
johnpye |
974 |
self._run('testcmslv2',"CMSlv") |
87 |
johnpye |
966 |
|
88 |
|
|
def testsunpos1(self): |
89 |
|
|
self._run('example_1_6_1',"QRSlv","johnpye/sunpos.a4c") |
90 |
|
|
|
91 |
|
|
def testsunpos2(self): |
92 |
|
|
self._run('example_1_6_2',"QRSlv","johnpye/sunpos.a4c") |
93 |
|
|
|
94 |
|
|
def testsunpos3(self): |
95 |
|
|
self._run('example_1_7_1',"QRSlv","johnpye/sunpos.a4c") |
96 |
|
|
|
97 |
|
|
def testsunpos4(self): |
98 |
|
|
self._run('example_1_7_2',"QRSlv","johnpye/sunpos.a4c") |
99 |
|
|
|
100 |
|
|
def testsunpos5(self): |
101 |
|
|
self._run('example_1_7_3',"QRSlv","johnpye/sunpos.a4c") |
102 |
|
|
|
103 |
|
|
def testsunpos6(self): |
104 |
|
|
self._run('example_1_8_1',"QRSlv","johnpye/sunpos.a4c") |
105 |
|
|
|
106 |
johnpye |
1073 |
def testinstanceas(self): |
107 |
|
|
M = self._run('example_1_6_1',"QRSlv","johnpye/sunpos.a4c") |
108 |
|
|
self.assertAlmostEqual( float(M.t_solar), M.t_solar.as("s")) |
109 |
|
|
self.assertAlmostEqual( float(M.t_solar)/3600, M.t_solar.as("h")) |
110 |
|
|
|
111 |
johnpye |
1126 |
class TestMatrix(AscendSelfTester): |
112 |
|
|
def testlog10(self): |
113 |
|
|
M = self._run('testlog10') |
114 |
|
|
print M.getMatrix().write(sys.stderr,"mmio") |
115 |
|
|
|
116 |
|
|
|
117 |
johnpye |
966 |
class TestIntegrator(Ascend): |
118 |
|
|
|
119 |
johnpye |
941 |
def testListIntegrators(self): |
120 |
|
|
I = ascpy.Integrator.getEngines() |
121 |
|
|
s1 = sorted([str(i) for i in I.values()]) |
122 |
johnpye |
972 |
s2 = sorted(['IDA','LSODE','AWW']) |
123 |
johnpye |
941 |
assert s1==s2 |
124 |
|
|
|
125 |
johnpye |
942 |
# this routine is reused by both testIDA and testLSODE |
126 |
johnpye |
941 |
def _testIntegrator(self,integratorname): |
127 |
johnpye |
940 |
self.L.load('johnpye/shm.a4c') |
128 |
|
|
M = self.L.findType('shm').getSimulation('sim') |
129 |
johnpye |
972 |
M.setSolver(ascpy.Solver('QRSlv')) |
130 |
johnpye |
1133 |
P = M.getParameters() |
131 |
|
|
M.setParameter('feastol',1e-12) |
132 |
johnpye |
979 |
print M.getChildren() |
133 |
|
|
assert float(M.x) == 10.0 |
134 |
|
|
assert float(M.v) == 0.0 |
135 |
johnpye |
941 |
t_end = math.pi |
136 |
johnpye |
940 |
|
137 |
|
|
I = ascpy.Integrator(M) |
138 |
|
|
I.setReporter(ascpy.IntegratorReporterNull(I)) |
139 |
johnpye |
941 |
I.setEngine(integratorname); |
140 |
johnpye |
940 |
I.setLinearTimesteps(ascpy.Units("s"), 0.0, t_end, 100); |
141 |
johnpye |
1133 |
I.setMinSubStep(0.0001); # these limits are required by IDA at present (numeric diff) |
142 |
|
|
I.setMaxSubStep(0.1); |
143 |
johnpye |
941 |
I.setInitialSubStep(0.001); |
144 |
|
|
I.setMaxSubSteps(200); |
145 |
johnpye |
944 |
if(integratorname=='IDA'): |
146 |
|
|
I.setParameter('autodiff',False) |
147 |
johnpye |
1133 |
for p in M.getParameters(): |
148 |
|
|
print p.getName(),"=",p.getValue() |
149 |
johnpye |
940 |
I.analyse(); |
150 |
|
|
I.solve(); |
151 |
johnpye |
941 |
print "At end of simulation," |
152 |
johnpye |
979 |
print "x = %f" % M.x |
153 |
|
|
print "v = %f" % M.v |
154 |
|
|
assert abs(float(M.x) + 10) < 1e-2 |
155 |
|
|
assert abs(float(M.v)) < 1e-2 |
156 |
johnpye |
940 |
assert I.getNumObservedVars() == 3 |
157 |
|
|
|
158 |
johnpye |
941 |
def testInvalidIntegrator(self): |
159 |
johnpye |
966 |
self.L.load('johnpye/shm.a4c') |
160 |
johnpye |
941 |
M = self.L.findType('shm').getSimulation('sim') |
161 |
johnpye |
972 |
M.setSolver(ascpy.Solver('QRSlv')) |
162 |
johnpye |
941 |
I = ascpy.Integrator(M) |
163 |
|
|
try: |
164 |
|
|
I.setEngine('___NONEXISTENT____') |
165 |
johnpye |
972 |
except RuntimeError: |
166 |
johnpye |
941 |
return |
167 |
|
|
self.fail("setEngine did not raise error!") |
168 |
|
|
|
169 |
|
|
def testLSODE(self): |
170 |
|
|
self._testIntegrator('LSODE') |
171 |
|
|
|
172 |
johnpye |
972 |
def testIDA(self): |
173 |
|
|
self._testIntegrator('IDA') |
174 |
|
|
|
175 |
johnpye |
1016 |
def testparameters(self): |
176 |
|
|
self.L.load('johnpye/shm.a4c') |
177 |
|
|
M = self.L.findType('shm').getSimulation('sim') |
178 |
|
|
M.build() |
179 |
|
|
I = ascpy.Integrator(M) |
180 |
|
|
I.setEngine('IDA') |
181 |
|
|
P = I.getParameters() |
182 |
|
|
for p in P: |
183 |
|
|
print p.getName(),"=",p.getValue() |
184 |
|
|
assert len(P)==11 |
185 |
|
|
assert P[0].isStr() |
186 |
|
|
assert P[0].getName()=="linsolver" |
187 |
|
|
assert P[0].getValue()=='SPGMR' |
188 |
|
|
assert P[2].getName()=="autodiff" |
189 |
|
|
assert P[2].getValue()==True |
190 |
|
|
assert P[7].getName()=="atolvect" |
191 |
|
|
assert P[7].getBoolValue() == True |
192 |
|
|
P[2].setBoolValue(False) |
193 |
|
|
assert P[2].getBoolValue()==False |
194 |
|
|
I.setParameters(P) |
195 |
|
|
assert I.getParameterValue('autodiff')==False |
196 |
|
|
I.setParameter('autodiff',True) |
197 |
|
|
try: |
198 |
|
|
v = I.getParameterValue('nonexist') |
199 |
|
|
except KeyError: |
200 |
|
|
pass |
201 |
|
|
else: |
202 |
|
|
self.fail('Failed to trip invalid Integrator parameter') |
203 |
|
|
|
204 |
johnpye |
972 |
class TestLSODE(Ascend): |
205 |
|
|
|
206 |
johnpye |
964 |
def testzill(self): |
207 |
|
|
self.L.load('johnpye/zill.a4c') |
208 |
|
|
T = self.L.findType('zill') |
209 |
|
|
M = T.getSimulation('sim') |
210 |
johnpye |
972 |
M.setSolver(ascpy.Solver('QRSlv')) |
211 |
johnpye |
964 |
I = ascpy.Integrator(M) |
212 |
johnpye |
966 |
I.setEngine('LSODE') |
213 |
|
|
I.setMinSubStep(1e-7) |
214 |
|
|
I.setMaxSubStep(0.001) |
215 |
|
|
I.setMaxSubSteps(10000) |
216 |
johnpye |
964 |
I.setReporter(ascpy.IntegratorReporterConsole(I)) |
217 |
johnpye |
1017 |
I.setLinearTimesteps(ascpy.Units(), 1.0, 1.5, 5) |
218 |
johnpye |
964 |
I.analyse() |
219 |
|
|
I.solve() |
220 |
|
|
M.run(T.getMethod('self_test')) |
221 |
|
|
|
222 |
johnpye |
962 |
def testnewton(self): |
223 |
johnpye |
973 |
sys.stderr.write("STARTING TESTNEWTON\n") |
224 |
johnpye |
962 |
self.L.load('johnpye/newton.a4c') |
225 |
|
|
T = self.L.findType('newton') |
226 |
|
|
M = T.getSimulation('sim') |
227 |
|
|
M.solve(ascpy.Solver("QRSlv"),ascpy.SolverReporter()) |
228 |
|
|
I = ascpy.Integrator(M) |
229 |
|
|
I.setEngine('LSODE') |
230 |
johnpye |
963 |
I.setParameter('rtolvect',False) |
231 |
|
|
I.setParameter('rtol',1e-7) |
232 |
|
|
I.setParameter('atolvect',False) |
233 |
|
|
I.setParameter('atol',1e-7) |
234 |
|
|
I.setMinSubStep(1e-7) |
235 |
|
|
I.setMaxSubStep(0.001) |
236 |
|
|
I.setMaxSubSteps(10000) |
237 |
|
|
|
238 |
johnpye |
962 |
I.setReporter(ascpy.IntegratorReporterConsole(I)) |
239 |
johnpye |
1017 |
I.setLinearTimesteps(ascpy.Units("s"), 0, 2*float(M.v)/float(M.g), 2) |
240 |
johnpye |
962 |
I.analyse() |
241 |
|
|
I.solve() |
242 |
|
|
print "At end of simulation," |
243 |
johnpye |
979 |
print "x = %f" % M.x |
244 |
|
|
print "v = %f" % M.v |
245 |
johnpye |
962 |
M.run(T.getMethod('self_test')) |
246 |
|
|
|
247 |
johnpye |
961 |
def testlotka(self): |
248 |
|
|
self.L.load('johnpye/lotka.a4c') |
249 |
|
|
M = self.L.findType('lotka').getSimulation('sim') |
250 |
johnpye |
980 |
M.setSolver(ascpy.Solver("QRSlv")) |
251 |
johnpye |
961 |
I = ascpy.Integrator(M) |
252 |
|
|
I.setEngine('LSODE') |
253 |
|
|
I.setReporter(ascpy.IntegratorReporterConsole(I)) |
254 |
johnpye |
1017 |
I.setLinearTimesteps(ascpy.Units("s"), 0, 200, 5) |
255 |
johnpye |
961 |
I.analyse() |
256 |
johnpye |
979 |
print "Number of vars = %d" % I.getNumVars() |
257 |
|
|
assert I.getNumVars()==2 |
258 |
johnpye |
961 |
I.solve() |
259 |
|
|
assert I.getNumObservedVars() == 3; |
260 |
johnpye |
979 |
assert abs(M.R - 832) < 1.0 |
261 |
|
|
assert abs(M.F - 21.36) < 0.1 |
262 |
johnpye |
1017 |
|
263 |
|
|
#------------------------------------------------------------------------------- |
264 |
johnpye |
1032 |
# Testing of a external blackbox functions |
265 |
johnpye |
1021 |
|
266 |
johnpye |
1032 |
class TestBlackBox(AscendSelfTester): |
267 |
|
|
def testparsefail0(self): |
268 |
|
|
try: |
269 |
|
|
self.L.load('test/blackbox/parsefail0.a4c') |
270 |
|
|
self.fail("parsefail0 should not have loaded without errors") |
271 |
|
|
except: |
272 |
|
|
pass |
273 |
|
|
|
274 |
|
|
def testparsefail1(self): |
275 |
|
|
try: |
276 |
|
|
self.L.load('test/blackbox/parsefail1.a4c') |
277 |
|
|
self.fail("parsefail1 should not have loaded without errors") |
278 |
|
|
except: |
279 |
|
|
pass |
280 |
|
|
|
281 |
|
|
def testparsefail2(self): |
282 |
|
|
try: |
283 |
|
|
self.L.load('test/blackbox/parsefail2.a4c') |
284 |
|
|
self.fail("parsefail2 should not have loaded without errors") |
285 |
|
|
except: |
286 |
|
|
pass |
287 |
|
|
|
288 |
|
|
def testparsefail3(self): |
289 |
|
|
try: |
290 |
|
|
self.L.load('test/blackbox/parsefail3.a4c') |
291 |
|
|
self.fail("parsefail3 should not have loaded without errors") |
292 |
|
|
except: |
293 |
|
|
pass |
294 |
|
|
|
295 |
|
|
def testparsefail4(self): |
296 |
|
|
try: |
297 |
|
|
self.L.load('test/blackbox/parsefail4.a4c') |
298 |
|
|
self.fail("parsefail4 should not have loaded") |
299 |
|
|
except: |
300 |
|
|
pass |
301 |
|
|
|
302 |
|
|
def testfail1(self): |
303 |
|
|
"""Mismatched arg counts check-- tests bbox, not ascend.""" |
304 |
|
|
self.L.load('test/blackbox/fail1.a4c') |
305 |
johnpye |
1034 |
try: |
306 |
johnpye |
1035 |
M = self.L.findType('fail1').getSimulation('sim') |
307 |
johnpye |
1034 |
self.fail("expected exception was not raised") |
308 |
|
|
except RuntimeError,e: |
309 |
|
|
print "Caught exception '%s', assumed ok" % e |
310 |
johnpye |
1032 |
|
311 |
|
|
def testfail2(self): |
312 |
|
|
"""Incorrect data arg check -- tests bbox, not ascend""" |
313 |
|
|
self.L.load('test/blackbox/fail2.a4c') |
314 |
johnpye |
1035 |
try: |
315 |
|
|
M = self.L.findType('fail2').getSimulation('sim') |
316 |
|
|
self.fail("expected exception was not raised") |
317 |
|
|
except RuntimeError,e: |
318 |
|
|
print "Caught exception '%s', assumed ok (should mention errors during instantiation)" % e |
319 |
johnpye |
1032 |
|
320 |
|
|
def testpass1(self): |
321 |
|
|
"""simple single bbox forward solve""" |
322 |
|
|
M = self._run('pass1',filename='test/blackbox/pass.a4c') |
323 |
|
|
|
324 |
|
|
def testpass2(self): |
325 |
|
|
"""simple single bbox reverse solve""" |
326 |
|
|
M = self._run('pass2',filename='test/blackbox/pass.a4c') |
327 |
|
|
|
328 |
|
|
def testpass3(self): |
329 |
|
|
"""simple double bbox solve""" |
330 |
|
|
M = self._run('pass3',filename='test/blackbox/pass3.a4c') |
331 |
|
|
|
332 |
|
|
def testpass4(self): |
333 |
|
|
"""simple double bbox reverse solve""" |
334 |
|
|
M = self._run('pass4',filename='test/blackbox/pass3.a4c') |
335 |
|
|
|
336 |
|
|
def testpass5(self): |
337 |
|
|
M = self._run('pass5',filename='test/blackbox/pass5.a4c') |
338 |
|
|
|
339 |
|
|
def testpass6(self): |
340 |
|
|
M = self._run('pass6',filename='test/blackbox/pass5.a4c') |
341 |
|
|
|
342 |
|
|
def testpass7(self): |
343 |
|
|
M = self._run('pass7',filename='test/blackbox/passmerge.a4c') |
344 |
|
|
|
345 |
|
|
def testpass8(self): |
346 |
|
|
M = self._run('pass8',filename='test/blackbox/passmerge.a4c') |
347 |
|
|
|
348 |
|
|
def testpass9(self): |
349 |
|
|
M = self._run('pass9',filename='test/blackbox/passmerge.a4c') |
350 |
|
|
|
351 |
|
|
def testpass10(self): |
352 |
|
|
M = self._run('pass10',filename='test/blackbox/passmerge.a4c') |
353 |
|
|
|
354 |
|
|
def testpass11(self): |
355 |
|
|
M = self._run('pass11',filename='test/blackbox/passmerge.a4c') |
356 |
|
|
|
357 |
|
|
def testpass12(self): |
358 |
|
|
M = self._run('pass12',filename='test/blackbox/passmerge.a4c') |
359 |
|
|
|
360 |
johnpye |
1037 |
# this test doesn't work: 'system is inconsistent' -- and structurally singular |
361 |
|
|
# def testpass13(self): |
362 |
|
|
# """cross-merged input/output solve""" |
363 |
|
|
# M = self._run('pass13',filename='test/blackbox/passmerge.a4c') |
364 |
johnpye |
1032 |
|
365 |
|
|
def testpass14(self): |
366 |
|
|
"""cross-merged input/output reverse solve""" |
367 |
|
|
M = self._run('pass14',filename='test/blackbox/passmerge.a4c') |
368 |
|
|
|
369 |
|
|
def testpass20(self): |
370 |
|
|
M = self._run('pass20',filename='test/blackbox/passarray.a4c') |
371 |
|
|
|
372 |
|
|
def testparsefail21(self): |
373 |
|
|
"""dense array of black boxes wrong syntax""" |
374 |
|
|
try: |
375 |
|
|
self.L.load('test/blackbox/parsefail21.a4c') |
376 |
|
|
self.fail("parsefail21 should not have loaded without errors") |
377 |
|
|
except: |
378 |
|
|
pass |
379 |
|
|
|
380 |
|
|
def testpass22(self): |
381 |
|
|
M = self._run('pass22',filename='test/blackbox/passarray.a4c') |
382 |
|
|
|
383 |
|
|
def testpass23(self): |
384 |
|
|
M = self._run('pass23',filename='test/blackbox/passarray.a4c') |
385 |
|
|
|
386 |
|
|
def testpass61(self): |
387 |
|
|
M = self._run('pass61',filename='test/blackbox/reinstantiate.a4c') |
388 |
|
|
|
389 |
|
|
def testpass62(self): |
390 |
|
|
M = self._run('pass62',filename='test/blackbox/reinstantiate.a4c') |
391 |
|
|
|
392 |
|
|
def testpass64(self): |
393 |
|
|
M = self._run('pass64',filename='test/blackbox/reinstantiate.a4c') |
394 |
|
|
|
395 |
|
|
def testpass65(self): |
396 |
|
|
M = self._run('pass65',filename='test/blackbox/reinstantiate.a4c') |
397 |
|
|
|
398 |
|
|
def testpass66(self): |
399 |
|
|
M = self._run('pass66',filename='test/blackbox/reinstantiate.a4c') |
400 |
|
|
|
401 |
|
|
def testpass67(self): |
402 |
|
|
M = self._run('pass67',filename='test/blackbox/reinstantiate.a4c') |
403 |
|
|
|
404 |
johnpye |
1024 |
class TestExtFn(AscendSelfTester): |
405 |
johnpye |
1021 |
def testextfntest(self): |
406 |
johnpye |
1024 |
M = self._run('extfntest',filename='johnpye/extfn/extfntest.a4c') |
407 |
|
|
self.assertAlmostEqual(M.y, 2); |
408 |
|
|
self.assertAlmostEqual(M.x, 1); |
409 |
johnpye |
1021 |
self.assertAlmostEqual(M.y, M.x + 1); |
410 |
|
|
|
411 |
johnpye |
1039 |
def testextrelfor(self): |
412 |
|
|
M = self._run('extrelfor',filename='johnpye/extfn/extrelfor.a4c') |
413 |
johnpye |
1024 |
|
414 |
johnpye |
1056 |
## @TODO fix bug with badly-named bbox rel in a loop (Ben, maybe) |
415 |
|
|
# def testextrelforbadnaming(self): |
416 |
|
|
# self.L.load('johnpye/extfn/extrelforbadnaming.a4c') |
417 |
|
|
# T = self.L.findType('extrelfor') |
418 |
|
|
# M = T.getSimulation('sim') |
419 |
|
|
# M.solve(ascpy.Solver('QRSlv'),ascpy.SolverReporter()) |
420 |
|
|
# print "x[1] = %f" % M.x[1] |
421 |
|
|
# print "x[2] = %f" % M.x[2] |
422 |
|
|
# print "x[3] = %f" % M.x[3] |
423 |
|
|
# print "x[4] = %f" % M.x[4] |
424 |
|
|
# print "x[5] = %f" % M.x[5] |
425 |
|
|
# M.run(T.getMethod('self_test')) |
426 |
johnpye |
1039 |
|
427 |
johnpye |
1024 |
def testextrelrepeat(self): |
428 |
|
|
M = self._run('extrelrepeat',filename='johnpye/extfn/extrelrepeat.a4c') |
429 |
|
|
|
430 |
johnpye |
1021 |
#------------------------------------------------------------------------------- |
431 |
johnpye |
1162 |
# Testing of Sensitivity module |
432 |
|
|
|
433 |
|
|
class TestSensitivity(AscendSelfTester): |
434 |
|
|
def test1(self): |
435 |
|
|
self.L.load('sensitivity_test.a4c') |
436 |
|
|
T = self.L.findType('sensitivity_test') |
437 |
|
|
M = T.getSimulation('sim',False) |
438 |
|
|
M.run(T.getMethod('on_load')) |
439 |
|
|
M.solve(ascpy.Solver('QRSlv'),ascpy.SolverReporter()) |
440 |
|
|
M.run(T.getMethod('analyse')) |
441 |
johnpye |
1163 |
M.run(T.getMethod('self_test')) |
442 |
johnpye |
1167 |
|
443 |
johnpye |
1168 |
# def testall(self): |
444 |
|
|
# self.L.load('sensitivity_test.a4c') |
445 |
|
|
# T = self.L.findType('sensitivity_test_all') |
446 |
|
|
# M = T.getSimulation('sim',False) |
447 |
|
|
# M.run(T.getMethod('on_load')) |
448 |
|
|
# M.solve(ascpy.Solver('QRSlv'),ascpy.SolverReporter()) |
449 |
|
|
# M.run(T.getMethod('analyse')) |
450 |
|
|
# M.run(T.getMethod('self_test')) |
451 |
|
|
# CAUSES CRASH |
452 |
johnpye |
1162 |
|
453 |
|
|
#------------------------------------------------------------------------------- |
454 |
johnpye |
1024 |
# Testing of a ExtPy - external python methods |
455 |
|
|
|
456 |
|
|
class TestExtPy(AscendSelfTester): |
457 |
johnpye |
1055 |
def test1(self): |
458 |
|
|
self.L.load('johnpye/extpy/extpytest.a4c') |
459 |
|
|
T = self.L.findType('extpytest') |
460 |
|
|
M = T.getSimulation('sim') |
461 |
|
|
M.run(T.getMethod('self_test')) |
462 |
|
|
|
463 |
|
|
def test2(self): |
464 |
|
|
self.L.load('johnpye/extpy/extpytest.a4c') |
465 |
|
|
T = self.L.findType('extpytest') |
466 |
|
|
M = T.getSimulation('sim') |
467 |
|
|
M.run(T.getMethod('pythonthing')) |
468 |
|
|
M.run(T.getMethod('pythonthing')) |
469 |
|
|
M.run(T.getMethod('pythonthing')) |
470 |
|
|
M.run(T.getMethod('pythonthing')) |
471 |
|
|
# causes crash! |
472 |
johnpye |
1024 |
|
473 |
|
|
#------------------------------------------------------------------------------- |
474 |
johnpye |
1042 |
# Testing of saturated steam properties library (iapwssatprops.a4c) |
475 |
|
|
|
476 |
|
|
class TestSteam(AscendSelfTester): |
477 |
johnpye |
1161 |
|
478 |
johnpye |
1042 |
def testiapwssatprops1(self): |
479 |
|
|
M = self._run('testiapwssatprops1',filename='steam/iapwssatprops.a4c') |
480 |
|
|
def testiapwssatprops2(self): |
481 |
|
|
M = self._run('testiapwssatprops2',filename='steam/iapwssatprops.a4c') |
482 |
|
|
def testiapwssatprops3(self): |
483 |
|
|
M = self._run('testiapwssatprops3',filename='steam/iapwssatprops.a4c') |
484 |
|
|
|
485 |
johnpye |
1161 |
# test the stream model basically works |
486 |
johnpye |
1043 |
def testsatsteamstream(self): |
487 |
|
|
M = self._run('satsteamstream',filename='steam/satsteamstream.a4c') |
488 |
johnpye |
1042 |
|
489 |
johnpye |
1161 |
# test that we can solve in terms of various (rho,u) |
490 |
|
|
def testsatuv(self): |
491 |
|
|
self.L.load('steam/iapwssat.a4c') |
492 |
|
|
T = self.L.findType('testiapwssatuv') |
493 |
|
|
M = T.getSimulation('sim',False) |
494 |
|
|
M.run(T.getMethod('on_load')) |
495 |
|
|
M.solve(ascpy.Solver('QRSlv'),ascpy.SolverReporter()) |
496 |
|
|
print "p = %f bar" % M.p.as('bar'); |
497 |
|
|
print "T = %f C" % (M.T.as('K') - 273.15); |
498 |
|
|
print "x = %f" % M.x; |
499 |
|
|
M.run(T.getMethod('self_test')) |
500 |
|
|
M.run(T.getMethod('values2')) |
501 |
|
|
# M.v.setRealValueWithUnits(1.0/450,"m^3/kg"); |
502 |
|
|
# M.solve(ascpy.Solver('QRSlv'),ascpy.SolverReporter()) |
503 |
|
|
M.solve(ascpy.Solver('QRSlv'),ascpy.SolverReporter()) |
504 |
|
|
print "p = %f bar" % M.p.as('bar'); |
505 |
|
|
print "T = %f C" % (M.T.as('K') - 273.15); |
506 |
|
|
print "x = %f" % M.x; |
507 |
|
|
M.run(T.getMethod('self_test2')) |
508 |
|
|
|
509 |
|
|
|
510 |
johnpye |
1056 |
## @TODO fix error capture from bounds checking during initialisation |
511 |
|
|
# def testiapwssat1(self): |
512 |
|
|
# M = self._run('testiapwssat1',filename='steam/iapwssat.a4c') |
513 |
johnpye |
1042 |
|
514 |
johnpye |
1100 |
def testdsgsat(self): |
515 |
|
|
self.L.load('steam/dsgsat2.a4c') |
516 |
|
|
T = self.L.findType('dsgsat2') |
517 |
|
|
M = T.getSimulation('sim',False) |
518 |
johnpye |
1127 |
M.run(T.getMethod('on_load')) |
519 |
johnpye |
1100 |
M.solve(ascpy.Solver('QRSlv'),ascpy.SolverReporter()) |
520 |
johnpye |
1152 |
self.assertAlmostEqual(M.dTw_dt[2],0.0); |
521 |
johnpye |
1127 |
M.run(T.getMethod('configure_dynamic')) |
522 |
|
|
M.solve(ascpy.Solver('QRSlv'),ascpy.SolverReporter()) |
523 |
johnpye |
1132 |
return M |
524 |
|
|
|
525 |
|
|
def testintegLSODE(self): |
526 |
|
|
M = self.testdsgsat() |
527 |
|
|
M.qdot_s.setRealValueWithUnits(1000,"W/m") |
528 |
|
|
M.solve(ascpy.Solver('QRSlv'),ascpy.SolverReporter()) |
529 |
johnpye |
1144 |
#M.setParameter(' |
530 |
johnpye |
1132 |
I = ascpy.Integrator(M) |
531 |
|
|
I.setEngine('LSODE') |
532 |
johnpye |
1142 |
I.setParameter('meth','AM') |
533 |
|
|
I.setParameter('maxord',12) |
534 |
johnpye |
1132 |
I.setReporter(ascpy.IntegratorReporterConsole(I)) |
535 |
|
|
I.setLinearTimesteps(ascpy.Units("s"), 0, 5, 1) |
536 |
|
|
I.analyse() |
537 |
|
|
I.solve() |
538 |
|
|
|
539 |
|
|
def testintegIDA(self): |
540 |
johnpye |
1137 |
M = self.testdsgsat() |
541 |
johnpye |
1152 |
self.assertAlmostEqual(M.dTw_dt[2],0.0) |
542 |
|
|
Tw1 = float(M.T_w[2]) |
543 |
johnpye |
1137 |
T = self.L.findType('dsgsat2') |
544 |
|
|
M.run(T.getMethod('free_states')) |
545 |
johnpye |
1100 |
I = ascpy.Integrator(M) |
546 |
johnpye |
1132 |
I.setEngine('IDA') |
547 |
|
|
I.setParameter('linsolver','DENSE') |
548 |
|
|
I.setParameter('safeeval',True) |
549 |
|
|
I.setParameter('rtol',1e-8) |
550 |
|
|
I.setInitialSubStep(0.01) |
551 |
johnpye |
1137 |
I.setMinSubStep(0.001) |
552 |
johnpye |
1132 |
I.setMaxSubSteps(100) |
553 |
johnpye |
1100 |
I.setReporter(ascpy.IntegratorReporterConsole(I)) |
554 |
johnpye |
1132 |
I.setLinearTimesteps(ascpy.Units("s"), 0, 3600, 100) |
555 |
johnpye |
1100 |
I.analyse() |
556 |
|
|
I.solve() |
557 |
johnpye |
1132 |
I.analyse() |
558 |
|
|
I.solve() |
559 |
johnpye |
1152 |
self.assertAlmostEqual(float(M.T_w[2]),Tw1) |
560 |
johnpye |
1127 |
M.qdot_s.setRealValueWithUnits(1000,"W/m") |
561 |
|
|
self.assertAlmostEqual(M.qdot_s.as("W/m"),1000) |
562 |
|
|
M.solve(ascpy.Solver('QRSlv'),ascpy.SolverReporter()) |
563 |
johnpye |
1152 |
self.assertNotAlmostEqual(M.dTw_dt[2],0.0) |
564 |
johnpye |
1127 |
# I = ascpy.Integrator(M) |
565 |
|
|
# I.setEngine('LSODE') |
566 |
|
|
# I.setReporter(ascpy.IntegratorReporterConsole(I)) |
567 |
|
|
# I.setReporter(ascpy.IntegratorReporterConsole(I)) |
568 |
|
|
# I.setLinearTimesteps(ascpy.Units("s"), 0, 5, 100) |
569 |
|
|
# I.setMinSubStep(0.0001) |
570 |
|
|
# I.setMaxSubStep(100) |
571 |
|
|
# I.setInitialSubStep(0.1) |
572 |
|
|
# I.analyse() |
573 |
|
|
# I.solve() |
574 |
johnpye |
1046 |
|
575 |
johnpye |
1042 |
#------------------------------------------------------------------------------- |
576 |
johnpye |
1017 |
# Testing of freesteam external steam properties functions |
577 |
|
|
|
578 |
johnpye |
1032 |
with_freesteam = True |
579 |
johnpye |
1017 |
try: |
580 |
johnpye |
1039 |
# we assume that if the freesteam python module is installed, the ASCEND |
581 |
|
|
# external library will also be. |
582 |
johnpye |
1017 |
import freesteam |
583 |
|
|
have_freesteam = True |
584 |
johnpye |
1018 |
except ImportError,e: |
585 |
johnpye |
1017 |
have_freesteam = False |
586 |
|
|
|
587 |
johnpye |
1024 |
if with_freesteam and have_freesteam: |
588 |
johnpye |
1039 |
class TestFreesteam(AscendSelfTester): |
589 |
johnpye |
1119 |
# def testfreesteamtest(self): |
590 |
|
|
# """run the self-test cases bundled with freesteam""" |
591 |
|
|
# self._run('testfreesteam',filename='testfreesteam.a4c') |
592 |
johnpye |
1039 |
|
593 |
johnpye |
1021 |
def testload(self): |
594 |
johnpye |
1039 |
"""check that we can load 'thermalequilibrium2' (IMPORT "freesteam", etc)""" |
595 |
johnpye |
1017 |
self.L.load('johnpye/thermalequilibrium2.a4c') |
596 |
johnpye |
1018 |
|
597 |
johnpye |
1021 |
def testinstantiate(self): |
598 |
johnpye |
1039 |
"""load an instantiate 'thermalequilibrium2'""" |
599 |
johnpye |
1021 |
self.testload() |
600 |
johnpye |
1018 |
M = self.L.findType('thermalequilibrium2').getSimulation('sim') |
601 |
johnpye |
1039 |
return M |
602 |
johnpye |
1021 |
|
603 |
johnpye |
1039 |
def testintegrate(self): |
604 |
|
|
"""integrate transfer of heat from one mass of water/steam to another |
605 |
|
|
according to Newton's law of cooling""" |
606 |
|
|
M = self.testinstantiate() |
607 |
johnpye |
1018 |
M.setSolver(ascpy.Solver("QRSlv")) |
608 |
johnpye |
1039 |
I = ascpy.Integrator(M) |
609 |
|
|
I.setEngine('LSODE') |
610 |
|
|
I.setReporter(ascpy.IntegratorReporterConsole(I)) |
611 |
|
|
I.setLinearTimesteps(ascpy.Units("s"), 0, 3000, 30) |
612 |
johnpye |
1055 |
I.setMinSubStep(0.01) |
613 |
|
|
I.setInitialSubStep(1) |
614 |
johnpye |
1039 |
I.analyse() |
615 |
|
|
print "Number of vars = %d" % I.getNumVars() |
616 |
|
|
assert I.getNumVars()==2 |
617 |
|
|
I.solve() |
618 |
|
|
assert I.getNumObservedVars() == 3; |
619 |
|
|
print "S[1].T = %f K" % M.S[1].T |
620 |
|
|
print "S[2].T = %f K" % M.S[2].T |
621 |
|
|
print "Q = %f W" % M.Q |
622 |
johnpye |
1119 |
self.assertAlmostEqual(float(M.S[1].T),506.77225109,4); |
623 |
johnpye |
1056 |
self.assertAlmostEqual(float(M.S[2].T),511.605173967,5); |
624 |
|
|
self.assertAlmostEqual(float(M.Q),-48.32922877329,3); |
625 |
johnpye |
1039 |
self.assertAlmostEqual(float(M.t),3000); |
626 |
|
|
print "Note that the above values have not been verified analytically" |
627 |
johnpye |
1017 |
|
628 |
johnpye |
1170 |
def testcollapsingcan2(self): |
629 |
|
|
""" solve the collapsing can model using IAPWS-IF97 steam props """ |
630 |
|
|
M = self._run("collapsingcan2",filename="collapsingcan2.a4c"); |
631 |
johnpye |
1121 |
|
632 |
johnpye |
1017 |
#------------------------------------------------------------------------------- |
633 |
|
|
# Testing of IDA models using DENSE linear solver |
634 |
|
|
|
635 |
johnpye |
1016 |
class TestIDADENSE(Ascend): |
636 |
johnpye |
1017 |
"""IDA DAE integrator, DENSE linear solver""" |
637 |
johnpye |
961 |
|
638 |
johnpye |
979 |
def testnewton(self): |
639 |
|
|
sys.stderr.write("STARTING TESTNEWTON\n") |
640 |
|
|
self.L.load('johnpye/newton.a4c') |
641 |
|
|
T = self.L.findType('newton') |
642 |
|
|
M = T.getSimulation('sim') |
643 |
|
|
M.solve(ascpy.Solver("QRSlv"),ascpy.SolverReporter()) |
644 |
|
|
I = ascpy.Integrator(M) |
645 |
|
|
I.setEngine('IDA') |
646 |
johnpye |
1016 |
I.setParameter('linsolver','DENSE') |
647 |
johnpye |
979 |
I.setParameter('safeeval',True) |
648 |
|
|
I.setParameter('rtol',1e-8) |
649 |
|
|
I.setMaxSubStep(0.001) |
650 |
|
|
I.setMaxSubSteps(10000) |
651 |
|
|
|
652 |
|
|
I.setReporter(ascpy.IntegratorReporterConsole(I)) |
653 |
johnpye |
1017 |
I.setLinearTimesteps(ascpy.Units("s"), 0, 2*float(M.v)/float(M.g), 2) |
654 |
johnpye |
979 |
I.analyse() |
655 |
|
|
I.solve() |
656 |
|
|
print "At end of simulation," |
657 |
|
|
print "x = %f" % M.x |
658 |
|
|
print "v = %f" % M.v |
659 |
|
|
M.run(T.getMethod('self_test')) |
660 |
|
|
|
661 |
johnpye |
1042 |
def testlotka(self): |
662 |
johnpye |
991 |
self.L.load('johnpye/lotka.a4c') |
663 |
|
|
M = self.L.findType('lotka').getSimulation('sim') |
664 |
|
|
M.setSolver(ascpy.Solver("QRSlv")) |
665 |
|
|
I = ascpy.Integrator(M) |
666 |
|
|
I.setEngine('IDA') |
667 |
|
|
I.setReporter(ascpy.IntegratorReporterConsole(I)) |
668 |
|
|
I.setLinearTimesteps(ascpy.Units("s"), 0, 200, 5); |
669 |
|
|
I.setParameter('linsolver','DENSE') |
670 |
|
|
I.setParameter('rtol',1e-8); |
671 |
|
|
I.analyse() |
672 |
|
|
assert I.getNumVars()==2 |
673 |
|
|
assert abs(M.R - 1000) < 1e-300 |
674 |
|
|
I.solve() |
675 |
johnpye |
1017 |
assert I.getNumObservedVars() == 3 |
676 |
johnpye |
991 |
assert abs(M.R - 832) < 1.0 |
677 |
|
|
assert abs(M.F - 21.36) < 0.1 |
678 |
johnpye |
972 |
|
679 |
johnpye |
975 |
def testdenx(self): |
680 |
johnpye |
1026 |
print "-----------------------------=====" |
681 |
johnpye |
942 |
self.L.load('johnpye/idadenx.a4c') |
682 |
|
|
M = self.L.findType('idadenx').getSimulation('sim') |
683 |
johnpye |
1017 |
M.setSolver(ascpy.Solver("QRSlv")) |
684 |
johnpye |
942 |
I = ascpy.Integrator(M) |
685 |
|
|
I.setEngine('IDA') |
686 |
johnpye |
1024 |
I.setParameter('calcic','YA_YPD') |
687 |
johnpye |
972 |
I.setParameter('linsolver','DENSE') |
688 |
johnpye |
944 |
I.setReporter(ascpy.IntegratorReporterConsole(I)) |
689 |
johnpye |
1017 |
I.setLogTimesteps(ascpy.Units("s"), 0.4, 4e10, 11) |
690 |
johnpye |
950 |
I.setMaxSubStep(0); |
691 |
johnpye |
1017 |
I.setInitialSubStep(0) |
692 |
johnpye |
950 |
I.setMaxSubSteps(0); |
693 |
johnpye |
944 |
I.setParameter('autodiff',True) |
694 |
|
|
I.analyse() |
695 |
|
|
I.solve() |
696 |
johnpye |
1022 |
assert abs(float(M.y1) - 5.1091e-08) < 2e-9 |
697 |
|
|
assert abs(float(M.y2) - 2.0437e-13) < 2e-14 |
698 |
johnpye |
1017 |
assert abs(float(M.y3) - 1.0) < 1e-5 |
699 |
johnpye |
942 |
|
700 |
johnpye |
1058 |
## @TODO fails during IDACalcIC (model too big?) |
701 |
|
|
# def testkryx(self): |
702 |
|
|
# self.L.load('johnpye/idakryx.a4c') |
703 |
|
|
# ascpy.getCompiler().setUseRelationSharing(False) |
704 |
|
|
# M = self.L.findType('idakryx').getSimulation('sim') |
705 |
|
|
# M.setSolver(ascpy.Solver('QRSlv')) |
706 |
|
|
# M.build() |
707 |
|
|
# I = ascpy.Integrator(M) |
708 |
|
|
# I.setEngine('IDA') |
709 |
|
|
# I.setReporter(ascpy.IntegratorReporterConsole(I)) |
710 |
|
|
# I.setParameter('linsolver','DENSE') |
711 |
|
|
# I.setParameter('maxl',8) |
712 |
|
|
# I.setParameter('gsmodified',False) |
713 |
|
|
# I.setParameter('autodiff',True) |
714 |
|
|
# I.setParameter('rtol',0) |
715 |
|
|
# I.setParameter('atol',1e-3); |
716 |
|
|
# I.setParameter('atolvect',False) |
717 |
|
|
# I.setParameter('calcic','YA_YDP') |
718 |
|
|
# I.analyse() |
719 |
|
|
# I.setLogTimesteps(ascpy.Units("s"), 0.01, 10.24, 11) |
720 |
|
|
# I.solve() |
721 |
|
|
# assert abs(M.u[2][2].getValue()) < 1e-5 |
722 |
johnpye |
1017 |
|
723 |
|
|
#------------------------------------------------------------------------------- |
724 |
|
|
# Testing of IDA models using SPGMR linear solver (Krylov) |
725 |
|
|
|
726 |
johnpye |
1016 |
# these tests are disabled until SPGMR preconditioning has been implemented |
727 |
|
|
class TestIDASPGMR:#(Ascend): |
728 |
|
|
def testlotka(self): |
729 |
|
|
self.L.load('johnpye/lotka.a4c') |
730 |
|
|
M = self.L.findType('lotka').getSimulation('sim') |
731 |
|
|
M.setSolver(ascpy.Solver("QRSlv")) |
732 |
johnpye |
951 |
I = ascpy.Integrator(M) |
733 |
|
|
I.setEngine('IDA') |
734 |
|
|
I.setReporter(ascpy.IntegratorReporterConsole(I)) |
735 |
johnpye |
1017 |
I.setLinearTimesteps(ascpy.Units("s"), 0, 200, 5) |
736 |
|
|
I.setParameter('rtol',1e-8) |
737 |
johnpye |
951 |
I.analyse() |
738 |
johnpye |
1016 |
assert I.getNumVars()==2 |
739 |
|
|
assert abs(M.R - 1000) < 1e-300 |
740 |
johnpye |
951 |
I.solve() |
741 |
johnpye |
1017 |
assert I.getNumObservedVars() == 3 |
742 |
johnpye |
1016 |
assert abs(M.R - 832) < 1.0 |
743 |
|
|
assert abs(M.F - 21.36) < 0.1 |
744 |
johnpye |
951 |
|
745 |
johnpye |
1016 |
|
746 |
johnpye |
991 |
def testkryx(self): |
747 |
johnpye |
951 |
self.L.load('johnpye/idakryx.a4c') |
748 |
|
|
M = self.L.findType('idakryx').getSimulation('sim') |
749 |
johnpye |
952 |
M.build() |
750 |
johnpye |
951 |
I = ascpy.Integrator(M) |
751 |
|
|
I.setEngine('IDA') |
752 |
|
|
I.setReporter(ascpy.IntegratorReporterConsole(I)) |
753 |
johnpye |
992 |
I.setParameter('linsolver','SPGMR') |
754 |
johnpye |
993 |
I.setParameter('prec','JACOBI') |
755 |
johnpye |
970 |
I.setParameter('maxl',8) |
756 |
johnpye |
952 |
I.setParameter('gsmodified',False) |
757 |
|
|
I.setParameter('autodiff',True) |
758 |
johnpye |
993 |
I.setParameter('gsmodified',True) |
759 |
johnpye |
952 |
I.setParameter('rtol',0) |
760 |
|
|
I.setParameter('atol',1e-3); |
761 |
|
|
I.setParameter('atolvect',False) |
762 |
johnpye |
993 |
I.setParameter('calcic','Y') |
763 |
johnpye |
952 |
I.analyse() |
764 |
|
|
I.setLogTimesteps(ascpy.Units("s"), 0.01, 10.24, 10); |
765 |
johnpye |
1017 |
print M.udot[1][3] |
766 |
johnpye |
952 |
I.solve() |
767 |
|
|
assert 0 |
768 |
johnpye |
967 |
|
769 |
johnpye |
1016 |
def testzill(self): |
770 |
|
|
self.L.load('johnpye/zill.a4c') |
771 |
|
|
T = self.L.findType('zill') |
772 |
|
|
M = T.getSimulation('sim') |
773 |
|
|
M.setSolver(ascpy.Solver('QRSlv')) |
774 |
|
|
I = ascpy.Integrator(M) |
775 |
|
|
I.setEngine('IDA') |
776 |
|
|
I.setParameter('safeeval',False) |
777 |
|
|
I.setMinSubStep(1e-7) |
778 |
|
|
I.setMaxSubStep(0.001) |
779 |
|
|
I.setMaxSubSteps(10000) |
780 |
|
|
I.setReporter(ascpy.IntegratorReporterConsole(I)) |
781 |
johnpye |
1017 |
I.setLinearTimesteps(ascpy.Units(), 1.0, 1.5, 5) |
782 |
johnpye |
1016 |
I.analyse() |
783 |
|
|
I.solve() |
784 |
|
|
M.run(T.getMethod('self_test')) |
785 |
|
|
|
786 |
|
|
def testdenxSPGMR(self): |
787 |
|
|
self.L.load('johnpye/idadenx.a4c') |
788 |
|
|
M = self.L.findType('idadenx').getSimulation('sim') |
789 |
|
|
M.setSolver(ascpy.Solver('QRSlv')) |
790 |
|
|
I = ascpy.Integrator(M) |
791 |
|
|
I.setEngine('IDA') |
792 |
|
|
I.setReporter(ascpy.IntegratorReporterConsole(I)) |
793 |
johnpye |
1017 |
I.setLogTimesteps(ascpy.Units("s"), 0.4, 4e10, 11) |
794 |
johnpye |
1016 |
I.setMaxSubStep(0); |
795 |
|
|
I.setInitialSubStep(0); |
796 |
|
|
I.setMaxSubSteps(0); |
797 |
|
|
I.setParameter('autodiff',True) |
798 |
|
|
I.setParameter('linsolver','SPGMR') |
799 |
|
|
I.setParameter('gsmodified',False) |
800 |
|
|
I.setParameter('maxncf',10) |
801 |
|
|
I.analyse() |
802 |
|
|
I.solve() |
803 |
johnpye |
1017 |
assert abs(float(M.y1) - 5.1091e-08) < 1e-10 |
804 |
|
|
assert abs(float(M.y2) - 2.0437e-13) < 1e-15 |
805 |
|
|
assert abs(float(M.y3) - 1.0) < 1e-5 |
806 |
johnpye |
1016 |
|
807 |
johnpye |
943 |
# move code above down here if you want to temporarily avoid testing it |
808 |
johnpye |
932 |
class NotToBeTested: |
809 |
|
|
def nothing(self): |
810 |
|
|
pass |
811 |
johnpye |
1016 |
|
812 |
johnpye |
669 |
if __name__=='__main__': |
813 |
johnpye |
1118 |
# a whole bag of tricks to make sure we get the necessary dirs in our ascend, python and ld path vars |
814 |
johnpye |
1098 |
restart = 0 |
815 |
|
|
|
816 |
|
|
if platform.system()=="Windows": |
817 |
johnpye |
1120 |
LD_LIBRARY_PATH="PATH" |
818 |
johnpye |
1098 |
SEP = ";" |
819 |
|
|
else: |
820 |
|
|
LD_LIBRARY_PATH="LD_LIBRARY_PATH" |
821 |
|
|
SEP = ":" |
822 |
|
|
|
823 |
johnpye |
1119 |
freesteamdir = os.path.expanduser("~/freesteam/ascend") |
824 |
|
|
modeldirs = [os.path.abspath(os.path.join(sys.path[0],"models")),os.path.abspath(freesteamdir)] |
825 |
johnpye |
1118 |
if not os.environ.get('ASCENDLIBRARY'): |
826 |
johnpye |
1119 |
os.environ['ASCENDLIBRARY'] = SEP.join(modeldirs) |
827 |
johnpye |
1118 |
restart = 1 |
828 |
|
|
else: |
829 |
|
|
envmodelsdir = [os.path.abspath(i) for i in os.environ['ASCENDLIBRARY'].split(SEP)] |
830 |
johnpye |
1119 |
for l in modeldirs: |
831 |
|
|
if l in envmodelsdir[len(modeldirs):]: |
832 |
|
|
envmodelsdir.remove(l) |
833 |
|
|
restart = 1 |
834 |
|
|
for l in modeldirs: |
835 |
|
|
if l not in envmodelsdir: |
836 |
|
|
envmodelsdir.insert(0,l) |
837 |
|
|
restart = 1 |
838 |
|
|
os.environ['ASCENDLIBRARY'] = SEP.join(envmodelsdir) |
839 |
johnpye |
1118 |
|
840 |
johnpye |
1102 |
libdirs = ["pygtk","."] |
841 |
johnpye |
1098 |
libdirs = [os.path.normpath(os.path.join(sys.path[0],l)) for l in libdirs] |
842 |
|
|
if not os.environ.get(LD_LIBRARY_PATH): |
843 |
johnpye |
1105 |
os.environ[LD_LIBRARY_PATH]=SEP.join(libdirs) |
844 |
johnpye |
1106 |
restart = 1 |
845 |
johnpye |
1098 |
else: |
846 |
|
|
envlibdirs = [os.path.normpath(i) for i in os.environ[LD_LIBRARY_PATH].split(SEP)] |
847 |
|
|
for l in libdirs: |
848 |
johnpye |
1106 |
if l in envlibdirs[len(libdirs):]: |
849 |
|
|
envlibdirs.remove(l) |
850 |
|
|
restart = 1 |
851 |
|
|
for l in libdirs: |
852 |
johnpye |
1098 |
if l not in envlibdirs: |
853 |
|
|
envlibdirs.insert(0,l) |
854 |
johnpye |
1106 |
restart = 1 |
855 |
johnpye |
1098 |
os.environ[LD_LIBRARY_PATH] = SEP.join(envlibdirs) |
856 |
|
|
|
857 |
johnpye |
1102 |
pypath = os.path.normpath(os.path.join(sys.path[0],"pygtk")) |
858 |
|
|
if not os.environ.get('PYTHONPATH'): |
859 |
|
|
os.environ['PYTHONPATH']=pypath |
860 |
|
|
else: |
861 |
|
|
envpypath = os.environ['PYTHONPATH'].split(SEP) |
862 |
|
|
if pypath not in envpypath: |
863 |
|
|
envpypath.insert(0,pypath) |
864 |
johnpye |
1137 |
os.environ['PYTHONPATH']=SEP.join(envpypath) |
865 |
johnpye |
1102 |
restart = 1 |
866 |
|
|
|
867 |
johnpye |
1098 |
if restart: |
868 |
johnpye |
1167 |
script = os.path.join(sys.path[0],"test.py") |
869 |
johnpye |
1119 |
print "Restarting with..." |
870 |
johnpye |
1194 |
print " export LD_LIBRARY_PATH=%s" % os.environ.get(LD_LIBRARY_PATH) |
871 |
|
|
print " export PYTHONPATH=%s" % os.environ.get('PYTHONPATH') |
872 |
|
|
print " export ASCENDLIBRARY=%s" % os.environ.get('ASCENDLIBRARY') |
873 |
johnpye |
1098 |
|
874 |
johnpye |
1167 |
if sys.argv[0] !="/usr/bin/gdb": |
875 |
|
|
os.execvp("python",[script] + sys.argv) |
876 |
|
|
else: |
877 |
|
|
print "Can't restart from inside GDB" |
878 |
|
|
|
879 |
johnpye |
1098 |
import ascpy |
880 |
|
|
|
881 |
|
|
try: |
882 |
|
|
import cunit |
883 |
|
|
except: |
884 |
|
|
pass |
885 |
|
|
|
886 |
johnpye |
966 |
atexit.register(ascpy.shutdown) |
887 |
johnpye |
1008 |
#suite = unittest.TestSuite() |
888 |
johnpye |
1003 |
#suite = unittest.defaultTestLoader.loadTestsFromName('__main__') |
889 |
johnpye |
1008 |
#unittest.TextTestRunner(verbosity=2).run(suite) |
890 |
|
|
unittest.main() |