835 |
assert abs(float(M.y2) - 2.0437e-13) < 2e-14 |
assert abs(float(M.y2) - 2.0437e-13) < 2e-14 |
836 |
assert abs(float(M.y3) - 1.0) < 1e-5 |
assert abs(float(M.y3) - 1.0) < 1e-5 |
837 |
|
|
838 |
|
def testhires(self): |
839 |
|
self.L.load('test/hires.a4c') |
840 |
|
T = self.L.findType('hires') |
841 |
|
M = T.getSimulation('sim') |
842 |
|
M.setSolver(ascpy.Solver('QRSlv')) |
843 |
|
I = ascpy.Integrator(M) |
844 |
|
I.setEngine('IDA') |
845 |
|
I.setParameter('linsolver','DENSE') |
846 |
|
I.setParameter('rtol',1e-7) |
847 |
|
I.setParameter('atolvect',0) |
848 |
|
I.setParameter('atol',1e-7) |
849 |
|
I.setReporter(ascpy.IntegratorReporterConsole(I)) |
850 |
|
I.setLogTimesteps(ascpy.Units(""), 1, 321.8122, 5) |
851 |
|
I.setMaxSubStep(1); |
852 |
|
I.setInitialSubStep(1e-9) |
853 |
|
I.setMaxSubSteps(5000) |
854 |
|
I.analyse() |
855 |
|
I.solve() |
856 |
|
M.run(T.getMethod('self_test')) |
857 |
|
|
858 |
## @TODO fails during IDACalcIC (model too big?) |
## @TODO fails during IDACalcIC (model too big?) |
859 |
# def testkryx(self): |
# def testkryx(self): |
860 |
# self.L.load('johnpye/idakryx.a4c') |
# self.L.load('johnpye/idakryx.a4c') |