# if not opts.Ethr:
# opts.Ethr = np.logspace(-21, -15, 100)
# if not opts.collE:
# opts.collE = np.logspace(-21, -15, 100)
# =================================================
### load data
# =================================================
if opts.verbose:
print "loading data from", opts.outfilename
t_P, q, N_m = nmu.load_out(opts.outfilename, tmin=opts.time_xmin, tmax=opts.time_xmax, downsample=opts.downsample)
if opts.verbose:
print "loading system from", opts.logfilename
system = nmu.load_log(opts.logfilename)
network = system.network
if opts.verbose:
print "decomposing network into generations"
gens, coups = network.gens()
# =================================================
### define subsets of modes
# =================================================
### pull apart by number of couplings
if opts.verbose:
print "separating modes by num_k:"
num_k_modeNos = {}
for num_k in opts.num_k:
plots = (opts.res or opts.amp_res) ################################################## # # load data # ################################################## filename1, filename2 = args if opts.verbose: print "loading data from", filename1 t_P1, q1, N_m1 = nmu.load_out(filename1, tmin=opts.tmin, tmax=opts.tmax, downsample=opts.downsample1) if opts.logfilename1: if opts.verbose: print "loading system information from", opts.logfilename1 system1 = nmu.load_log(opts.logfilename1) if opts.verbose: print "loading data from ", filename2 t_P2, q2, N_m2 = nmu.load_out(filename2, tmin=opts.tmin, tmax=opts.tmax, downsample=opts.downsample2) if opts.logfilename2: if opts.verbose: print "loading system information from", opts.logfilename2 system2 = nmu.load_log(opts.logfilename2) ################################################## # # check for network consistency/compatibility # ##################################################
filenames = sorted( [l.strip() for l in cachefile.readlines() ])
cachefile.close()
else:
import os
import glob
if opts.verbose: print "reading in all *.log filenames from "+os.getcwd()
filenames = sorted( glob.glob("./*.log") )
dat = []
### read data from logs
for filename in filenames:
if opts.verbose: print "reading system from "+filename
try:
sys = nmu.load_log(filename)
except:
print "\tERROR when reading "+filename
continue
net = sys.network
### pull out data
gens, coups = net.gens()
this_log = [ filename, (sys.Mprim, sys.Rprim, sys.Mcomp, sys.eccentricity, sys.Porb), (len(net), len(net.to_triples())) ]
these_gens = []
for genNo, gen in enumerate(gens):
max_m = -infty
outfile = open(opts.outfilename, "a")
else:
outfile = open(opts.outfilename, "w")
else:
outfile = sys.stdout
if not opts.logfilename:
opts.logfilename = raw_input("logfilename = ")
if opts.onward_logfilename and (not opts.onward):
opts.onward = raw_input("onward filename = ")
if opts.init_time != "none":
opts.init_time = float(opts.init_time) # stupid fix to get control sequence correct later on
system = nm_u.load_log(opts.logfilename)
Porb = system.Porb
Oorb = 2*np.pi/Porb
N_m = len(system.network) # number of modes
#dimension = np.array(2*(N_m), dtype="i") # dimensionality of the problem
dimension = 2*(N_m) # dimensionality of the problem
if opts.equilib_IC:
if "dxdt" in opts.function:
tcurrent = "x"
elif "dqdt" in opts.function:
tcurrent = "q"
else:
raise ValueError, "could not determine tcurrent from --function=%s"%opts.function
