minlndos = 1e100
maxlndos = -1e100
numframes = 0
dataformat = 'data/lv/ww%.2f-ff%.2f-%gx%g-movie/%%06d' % (ww, ff, lenx, lenyz)
if 'tmi' in sys.argv:
dataformat = 'data/lv/ww%.2f-ff%.2f-%gx%g-tmi-movie/%%06d' % (ww, ff, lenx,
lenyz)
if 'toe' in sys.argv:
dataformat = 'data/lv/ww%.2f-ff%.2f-%gx%g-toe-movie/%%06d' % (ww, ff, lenx,
lenyz)
for frame in range(100000):
basename = dataformat % frame
try:
e, lndos = readandcompute.e_lndos(basename)
except:
break
numframes = frame + 1
minlndos = min(minlndos, lndos.min())
maxlndos = max(maxlndos, lndos.max())
try:
e, hist = readandcompute.e_and_total_init_histogram(basename)
except:
break
mine = min(mine, e.min() - 20)
maxe = max(maxe, e.max())
print('mine', mine)
print('maxe', maxe)
mine = 1e100
maxe = -1e100
numframes = 0
dataformat = 'data/%s/%s-%%s-movie/%%06d' % (subdirname, filename)
lastframe = -1
alldone = False
for frame in xrange(100000):
if alldone:
break
for suffix in suffixes:
basename = dataformat % (suffix, frame)
try:
e, lndos = readandcompute.e_lndos(basename)
except (KeyboardInterrupt, SystemExit):
raise
except:
alldone = True
break
numframes = frame + 1
mine = min(mine, e[lndos == lndos[-1]].max() - 20)
maxe = max(maxe, e[lndos == lndos[0]].min())
if numframes % 25 == 0 and frame != lastframe:
print 'counting %dth frame' % numframes
lastframe = frame
bestframe = sorted(
glob.glob('data/%s/%s-%s-movie/*-lndos.dat' %
(subdirname, filename, suffixes[0])))[-1]
fig, ax = plt.subplots()
mine = 1e100
maxe = -1e100
minlndos = 1e100
maxlndos = -1e100
numframes = 0
dataformat = 'data/%s/%s-%%s-movie/%%06d' % (subdirname, filename)
colors = ['k', 'b', 'r', 'g']
for frame in xrange(100000):
for suffix in suffixes:
basename = dataformat % (suffix, frame)
try:
e, lndos = readandcompute.e_lndos(basename)
except (KeyboardInterrupt, SystemExit):
raise
except:
# print "Failure in readandcompute.e_lndos:", sys.exc_info()[0]
break
numframes = frame+1
minlndos = min(minlndos, lndos.min())
maxlndos = max(maxlndos, lndos.max())
try:
e, hist = readandcompute.e_and_total_init_histogram(basename)
except (KeyboardInterrupt, SystemExit):
raise
except:
# print "Failure in readandcompute.e_and_total_init_histogram:", sys.exc_info()[0]
mine = 1e100
maxe = -1e100
numframes = 0
dataformat = 'data/%s/%s-%%s-movie/%%06d' % (subdirname, filename)
lastframe = -1
alldone = False
for frame in xrange(100000):
if alldone:
break
for suffix in suffixes:
basename = dataformat % (suffix, frame)
try:
e, lndos = readandcompute.e_lndos(basename)
except (KeyboardInterrupt, SystemExit):
raise
except:
alldone = True
break
numframes = frame+1
mine = min(mine, e[lndos == lndos[-1]].max() - 20)
maxe = max(maxe, e[lndos == lndos[0]].min())
if numframes % 25 == 0 and frame != lastframe:
print 'counting %dth frame' % numframes
lastframe = frame
bestframe = sorted(glob.glob('data/%s/%s-%s-movie/*-lndos.dat'
% (subdirname, filename, suffixes[0])))[-1]
best_e, best_lndos = readandcompute.e_lndos(bestframe)