Nseg = 2 * int(duration / Tseg)
else:
Nseg = 2 * (int(duration / Tseg) + 1)
Ajoins, Ejoins, Tjoins = [], [], []
leftover = None
for i, seg in enumerate(range(0, Nseg, 2)):
t0l = T0 + Tseg / 2 * seg
t0r = T0 + Tseg / 2 * (seg + 1)
print ""
print "~~~~ Segment pair number: %d ~~~~" % (i + 1)
print "Initial times:", t0l, ',', t0r
print ""
if options.seed == 'None':
shortftl = AS.CSpectra_to_ShortTermFT(cspec_dict, seed=None)
shortftr = AS.CSpectra_to_ShortTermFT(cspec_dict, seed=None)
else:
seedl = int(options.seed) + (day - 1) * Nseg + seg
seedr = int(options.seed) + (day - 1) * Nseg + (seg + 1)
print 'seeds: left, right', seedl, seedr
shortftl = AS.CSpectra_to_ShortTermFT(cspec_dict, seed=seedl)
shortftr = AS.CSpectra_to_ShortTermFT(cspec_dict, seed=seedr)
stsl = AS.InverseFT(shortftl, Neven=True, tOffset=t0l)
stsr = AS.InverseFT(shortftr, Neven=True, tOffset=t0r)
if i == 0:
Offset = stsl.s1.Offset1
Cadence = stsl.s1.Cadence1
joinrems = [
AS.WindowAndJoin(getattr(stsl, s),
getattr(stsr, s),
