225 |
P = I.getParameters() |
P = I.getParameters() |
226 |
for p in P: |
for p in P: |
227 |
print p.getName(),"=",p.getValue() |
print p.getName(),"=",p.getValue() |
228 |
assert len(P)==11 |
assert len(P)==12 |
229 |
assert P[0].isStr() |
assert P[0].isStr() |
230 |
assert P[0].getName()=="linsolver" |
assert P[0].getName()=="linsolver" |
231 |
assert P[0].getValue()=='SPGMR' |
assert P[0].getValue()=='SPGMR' |
232 |
assert P[2].getName()=="autodiff" |
assert P[2].getName()=="maxord" |
233 |
assert P[2].getValue()==True |
assert P[3].getName()=="autodiff" |
234 |
assert P[7].getName()=="atolvect" |
assert P[3].getValue()==True |
235 |
assert P[7].getBoolValue() == True |
assert P[8].getName()=="atolvect" |
236 |
P[2].setBoolValue(False) |
assert P[8].getBoolValue() == True |
237 |
assert P[2].getBoolValue()==False |
P[3].setBoolValue(False) |
238 |
|
assert P[3].getBoolValue()==False |
239 |
I.setParameters(P) |
I.setParameters(P) |
240 |
assert I.getParameterValue('autodiff')==False |
assert I.getParameterValue('autodiff')==False |
241 |
I.setParameter('autodiff',True) |
I.setParameter('autodiff',True) |
603 |
I = ascpy.Integrator(M) |
I = ascpy.Integrator(M) |
604 |
I.setEngine('LSODE') |
I.setEngine('LSODE') |
605 |
I.setReporter(ascpy.IntegratorReporterConsole(I)) |
I.setReporter(ascpy.IntegratorReporterConsole(I)) |
606 |
I.setLinearTimesteps(ascpy.Units("s"), 0, 5, 1) |
I.setLinearTimesteps(ascpy.Units("s"), 0, 3600, 10) |
607 |
I.analyse() |
I.analyse() |
608 |
I.solve() |
I.solve() |
609 |
|
|
631 |
I.setEngine('IDA') |
I.setEngine('IDA') |
632 |
I.setParameter('linsolver','DENSE') |
I.setParameter('linsolver','DENSE') |
633 |
I.setParameter('safeeval',True) |
I.setParameter('safeeval',True) |
634 |
I.setParameter('rtol',1e-5) |
I.setParameter('rtol',1e-4) |
635 |
I.setInitialSubStep(0.01) |
I.setParameter('atolvect',False) |
636 |
I.setMaxSubSteps(100) |
I.setParameter('atol',1e-4) |
637 |
|
I.setParameter('maxord',3) |
638 |
|
I.setInitialSubStep(0.001) |
639 |
I.setReporter(ascpy.IntegratorReporterConsole(I)) |
I.setReporter(ascpy.IntegratorReporterConsole(I)) |
640 |
I.setLinearTimesteps(ascpy.Units("s"), 0, 3600, 5) |
I.setLinearTimesteps(ascpy.Units("s"), 0, 3600, 10) |
641 |
I.analyse() |
I.analyse() |
642 |
I.solve() |
I.solve() |
643 |
self.assertAlmostEqual(float(M.T_w[2]),Tw1) |
self.assertAlmostEqual(float(M.T_w[2]),Tw1) |