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

Diff of /trunk/test.py

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

revision 1394 by jpye, Fri Apr 6 10:50:41 2007 UTC revision 1395 by jpye, Sat Apr 21 15:11:45 2007 UTC
# Line 811  class TestSteam(AscendSelfTester): Line 811  class TestSteam(AscendSelfTester):
811          # here is the peturbation...          # here is the peturbation...
812          qdot_s = float(M.qdot_s)          qdot_s = float(M.qdot_s)
813          print "OLD QDOT_S = %f" % qdot_s          print "OLD QDOT_S = %f" % qdot_s
814          M.qdot_s.setRealValueWithUnits(250,"W/m")          M.qdot_s.setRealValueWithUnits(6000,"W/m")
815          # IDA has its own initial conditions solver, so no need to call QRSlv here          # IDA has its own initial conditions solver, so no need to call QRSlv here
816          I = ascpy.Integrator(M)          I = ascpy.Integrator(M)
817          I.setEngine('IDA')          I.setEngine('IDA')
818          I.setParameter('linsolver','DENSE')          I.setParameter('linsolver','DENSE')
         I.setParameter('safeeval',True)  
         I.setParameter('rtol',1e-5)  
         I.setParameter('atolvect',False)  
         I.setParameter('atol',1e-4)  
         I.setParameter('maxord',3)  
         I.setInitialSubStep(0.001)  
         #I.setMaxSubStep(1);  
819          I.setReporter(ascpy.IntegratorReporterConsole(I))          I.setReporter(ascpy.IntegratorReporterConsole(I))
820          #I.setLinearTimesteps(ascpy.Units("s"), 0,300,300)          #I.setLinearTimesteps(ascpy.Units("s"), 0,300,300)
821          I.setLogTimesteps(ascpy.Units("s"), 0.001, 3600, 40)          I.setLogTimesteps(ascpy.Units("s"), 0.009, 1200, 150)
822          I.analyse()          I.analyse()
823          F = file('ga.mm','w')          F = file('ga.mm','w')
824          I.writeMatrix(F,'dg/dz')          I.writeMatrix(F,'dg/dz')

Legend:
Removed from v.1394  
changed lines
  Added in v.1395

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