87 |
def testLSODE(self): |
def testLSODE(self): |
88 |
self._testIntegrator('LSODE') |
self._testIntegrator('LSODE') |
89 |
|
|
90 |
|
def testzill(self): |
91 |
|
self.L.load('johnpye/zill.a4c') |
92 |
|
T = self.L.findType('zill') |
93 |
|
M = T.getSimulation('sim') |
94 |
|
M.solve(ascpy.Solver("QRSlv"),ascpy.SolverReporter()) |
95 |
|
I = ascpy.Integrator(M) |
96 |
|
I.setEngine('LSODE') |
97 |
|
I.setReporter(ascpy.IntegratorReporterConsole(I)) |
98 |
|
I.setLinearTimesteps(ascpy.Units(), 0, 1.5, 2); |
99 |
|
I.analyse() |
100 |
|
I.solve() |
101 |
|
M.run(T.getMethod('self_test')) |
102 |
|
|
103 |
|
|
104 |
def testnewton(self): |
def testnewton(self): |
105 |
self.L.load('johnpye/newton.a4c') |
self.L.load('johnpye/newton.a4c') |
106 |
T = self.L.findType('newton') |
T = self.L.findType('newton') |
117 |
I.setMaxSubSteps(10000) |
I.setMaxSubSteps(10000) |
118 |
|
|
119 |
I.setReporter(ascpy.IntegratorReporterConsole(I)) |
I.setReporter(ascpy.IntegratorReporterConsole(I)) |
120 |
I.setLinearTimesteps(ascpy.Units("s"), 0, 2*float(M.sim.v)/float(M.sim.g), 3); |
I.setLinearTimesteps(ascpy.Units("s"), 0, 2*float(M.sim.v)/float(M.sim.g), 2); |
121 |
I.analyse() |
I.analyse() |
122 |
I.solve() |
I.solve() |
123 |
print "At end of simulation," |
print "At end of simulation," |