Ejemplo n.º 1
0
print "++Object Launched: %g" % (finish-start)


print "--Get Composition"
start = time.clock()
OT.properties()
n,X,Y,Z,dXc,dXo = OT.retrieve()
print u"Fixed composition at X=%.3g,Y=%.3g,Z=%.3g" % (X,Y,Z)
print u"Using table %d" % (n + 1)
finish = time.clock()
print "--Got Composition: %g" % (finish-start)


print "++Get Composition"
start = time.clock()
n,X,Y,Z,dXc,dXo = OO.properties()
print u"Fixed composition at X=%.3g,Y=%.3g,Z=%.3g" % (X,Y,Z)
print u"Using table %d" % (n + 1)
finish = time.clock()
print "++Got Composition: %g" % (finish-start)


ntest = 1e4
logTs   = np.linspace(4.0,6.3,ntest)
logrhos = np.linspace(-6.0,0.3,ntest)
logkappa = np.empty((ntest,2))
ntest = int(ntest)

print "--Doing Interpolation"
start = time.clock()
for i in xrange(ntest):