# 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