def main(): ### get the folder structure analysistype = "Energy_spectrum" analysisdatabase = '/usr/users/iff_th2/duman/Cells_in_LAMMPS/DATA/' analysisdatabase += analysistype + '/' datafolderbase = '/local/duman/SIMULATIONS/Cells_in_LAMMPS/density_0.8/' savebase = "/usr/users/iff_th2/duman/Cells_in_LAMMPS/PLOTS/" savefolder = "Energy_spectrum" ### make the parameter choice # motility fp = [1.0, 5.0, 10.0] eps = 1.0 areak = 10.0 param = fp param_choice = 'fp' # compressibility # fp = 5.0 # eps = 5.0 # areak = [1.0, 10.0, 100.0] # param = areak # param_choice = 'areak' # adhesion # fp = 3.0 # areak = 10.0 # eps = [0.05, 1.0, 10.0] # param = eps # param_choice = 'eps' data = {} # carries the data per parameter set sims = {} # carries the simulation information per parameter set for p in param: if param_choice == 'areak': datafolder, analysisfile = read_write.gen_folders( eps, fp, p, savefolder, datafolderbase, analysisdatabase) elif param_choice == 'eps': datafolder, analysisfile = read_write.gen_folders( p, fp, areak, savefolder, datafolderbase, analysisdatabase) elif param_choice == 'fp': datafolder, analysisfile = read_write.gen_folders( eps, p, areak, savefolder, datafolderbase, analysisdatabase) sims[p] = read_write.read_sim_info(datafolder) x, y = read_write.read_2d_analysis_data(analysisfile) data[p] = y ### plot the data as a function of the parameter plot_data(x, data, param_choice, sims, savebase, savefolder) return
def main(): ### get the data folder parser = argparse.ArgumentParser() parser.add_argument("-e", "--eps", type=float, help="Strength of LJ potential") parser.add_argument("-f", "--fp", type=float, help="Propulsion force") parser.add_argument("-a", "--areak", type=float, help="Strength of area constraint potential") parser.add_argument("-k", "--kappa", type=float, help="Bending rigidity") parser.add_argument("-fl", "--folder", nargs="?", \ const='/local/duman/SIMULATIONS/Cells_in_LAMMPS/density_0.8/', \ help="Folder containing data, as in /local/duman/SIMULATIONS/Cells_in_LAMMPS/density_0.8/") parser.add_argument("-sb", "--savebase", nargs="?", \ const = "/usr/users/iff_th2/duman/Cells_in_LAMMPS/PLOTS/", \ help="Folder to save the data, as in /usr/users/iff_th2/duman/Cells_in_LAMMPS/PLOTS/") parser.add_argument("-sf", "--savefolder", nargs="?", const="Energy_spectrum", \ help="Specific folder for saving, as in Energy_spectrum") parser.add_argument("-s", "--save_eps", action="store_true", help="Decide whether to save in eps or not") args = parser.parse_args() ### read general information from the folder path1 = 'eps_' + str(args.eps) + '/fp_' + str(args.fp) + \ '/areak_' + str(args.areak) + '/kappa_' + str(args.kappa) + '/' datafolder = args.folder + path1 sim = read_write.read_sim_info(datafolder) ### read the data # path2 = 'eps_' + str(args.eps) + '_fp_' + str(args.fp) + \ # '_areak_' + str(args.areak) + '_kappa_' + str(args.kappa) path2 = 'eps_' + str(args.eps) + '_fp_' + str(args.fp) + \ '_areak_' + str(args.areak) analysisdatabase = '/usr/users/iff_th2/duman/Cells_in_LAMMPS/DATA/' analysisdatabase += args.savefolder + '/' analysisdata = analysisdatabase + args.savefolder + '_' + path2 + '.txt' x, y = read_write.read_2d_analysis_data(analysisdata) ### plot the data in the given time window plot_data(x, y, sim, args.savebase, args.savefolder, args.save_eps) return
def main(): ### get the data folder parser = argparse.ArgumentParser() parser.add_argument("-d", "--density", type=float, help="Density") parser.add_argument("-ka", "--kappa", type=float, help="Bending rigidity") parser.add_argument("-km", "--km", type=float, help="Motor spring strength") parser.add_argument("-p", "--pa", type=float, help="Probability when filaments are antialigned") parser.add_argument("-fl", "--folder", nargs="?", \ const='/local/duman/SIMULATIONS/Motor/Data/', \ help="Folder containing data, as in /local/duman/SIMULATIONS/Motor/Data/") parser.add_argument("-sb", "--savebase", nargs="?", \ const = "/usr/users/iff_th2/duman/Motor/PLOTS/", \ help="Folder to save the data, as in /usr/users/iff_th2/duman/Motor/PLOTS/") parser.add_argument("-sf", "--savefolder", nargs="?", const="Overlap_fnc", \ help="Specific folder for saving, as in Overlap_fnc") parser.add_argument("-s", "--save_eps", action="store_true", help="Decide whether to save in eps or not") args = parser.parse_args() ### read general information from the folder path1 = 'density_' + str(args.density) + '/kappa_' + str(args.kappa) + \ '/km_' + str(args.km) + '/panti_' + str(args.pa) + '/' datafolder = args.folder + path1 sim = read_write.read_sim_info(datafolder) ### read the data path2 = '_density_' + str(args.density) + '_kappa_' + str(args.kappa) + \ '_km_' + str(args.km) + '_panti_' + str(args.pa) analysisdatabase = '/usr/users/iff_th2/duman/Motor/DATA/' analysisdatabase += args.savefolder + '/' analysisdata = analysisdatabase + args.savefolder + path2 + '.txt' x, y = read_write.read_2d_analysis_data(analysisdata) ### plot the data in the given time window plot_data(x, y, sim, args.savebase, args.savefolder, args.save_eps) return
def main(): ### get the data folder parser = argparse.ArgumentParser() parser.add_argument("-e", "--eps", type=float, help="Strength of LJ potential") parser.add_argument("-f", "--fp", type=float, help="Propulsion force") parser.add_argument("-fl", "--folder", nargs="?", \ const='/local/duman/SIMULATIONS/Cells_in_LAMMPS/density_0.8/', \ help="Folder containing data, as in /local/duman/SIMULATIONS/Cells_in_LAMMPS/density_0.8/") parser.add_argument("-sb", "--savebase", nargs="?", \ const = "/usr/users/iff_th2/duman/Cells_in_LAMMPS/PLOTS/", \ help="Folder to save the data, as in /usr/users/iff_th2/duman/Cells_in_LAMMPS/PLOTS/") parser.add_argument("-sf", "--savefolder", nargs="?", const="Static_struct", \ help="Specific folder for saving, as in Static_struct") parser.add_argument("-s", "--save_eps", action="store_true", help="Decide whether to save in eps or not") args = parser.parse_args() ### accummulate data as a function of areak -with a dictionary- analysisdatabase = '/usr/users/iff_th2/duman/Cells_in_LAMMPS/DATA/' analysisdatabase += args.savefolder + '/' areak = [1.0, 10.0, 100.0] data = {} # carries the data per parameter set sims = {} # carries the simulation information per parameter set for a in areak: datafolder, analysisfile = read_write.gen_folders( args.eps, args.fp, a, args.savefolder, args.folder, analysisdatabase) sims[a] = read_write.read_sim_info(datafolder) x, y = read_write.read_2d_analysis_data(analysisfile) data[a] = y ### plot the data as a function of the parameter plot_data(x, data, sims, args.savebase, args.savefolder, args.save_eps) return
def main(): parser = argparse.ArgumentParser() parser.add_argument("-e", "--eps", type=float, \ help="Strength of LJ potential") parser.add_argument("-f", "--fp", type=float, \ help="Propulsion force") parser.add_argument("-a", "--areak", type=float, \ help="Area constraint potential strength") parser.add_argument("-k", "--kappa", type=float, \ help="Bending rigidity") parser.add_argument("-fl", "--folder", nargs="?", \ const='/local/duman/SIMULATIONS/Cells_in_LAMMPS/density_0.8/', \ help="Folder containing data, as in /local/duman/SIMULATIONS/Cells_in_LAMMPS/density_0.8/") parser.add_argument("-sb", "--savebase", nargs="?", \ const = "/usr/users/iff_th2/duman/Cells_in_LAMMPS/DATA/", \ help="Folder to save the data, as in /usr/users/iff_th2/duman/Cells_in_LAMMPS/DATA/") parser.add_argument("-sf", "--savefolder", nargs="?", \ const="Corr_length", \ help="Specific folder for saving, as in Corr_length") parser.add_argument("-s", "--save", action="store_true", help="Save for the phase diagram or not") args = parser.parse_args() ### read the data and general information from the folder print "Calculating correlation length for the following parameters : ", \ args.eps, ", ", args.fp, ", ", args.areak, ", ", args.kappa ### get the simulation data folder = args.folder + "eps_" + str(args.eps) + \ "/fp_" + str(args.fp) + "/areak_" + str(args.areak) + \ "/kappa_" + str(args.kappa) + "/" sim = read_write.read_sim_info(folder, True) ### calculate correlation length corr_len_folder = "/usr/users/iff_th2/duman/Cells_in_LAMMPS/DATA/" aname = "Sp_velocity_corr_subt" corr_len_folder += aname corr_len_file = corr_len_folder + "/" + aname + "_eps_" + str(args.eps) + \ "_fp_" + str(args.fp) + "_areak_" + str(args.areak) + \ "_kappa_" + str(args.kappa) + ".txt" x, y = read_write.read_2d_analysis_data(corr_len_file) corr_len = calc_corr_len(y, sim) ### write the data savebase = args.savebase + args.savefolder + "/" os.system("mkdir -p " + savebase) savefile = savebase + args.savefolder + "_eps_" + str(args.eps) + \ "_fp_" + str(args.fp) + "_areak_" + str(args.areak) + \ "_kappa_" + str(args.kappa) + ".txt" read_write.write_single_analysis_data(corr_len, savefile) ### collect the data for a phase diagram if args.save: aname = "PHASE_DIAGRAM" savebase = args.savebase + aname + "/" os.system("mkdir -p " + savebase) aname_type = aname + '_from_corr_length' savefile = savebase + aname_type + ".txt" fl = open(savefile, "a") fl.write(str(sim.eps) + "\t" + str(sim.fp) + "\t" + str(sim.areak) \ + "\t" + str(sim.kappa) + "\t" + str(corr_len) + "\n") fl.close() return
def main(): ### get the data folder parser = argparse.ArgumentParser() parser.add_argument("-e", "--eps", type=float, nargs="?", const=-1, \ help="Strength of LJ potential") parser.add_argument("-f", "--fp", type=float, nargs="?", const=-1, \ help="Propulsion force") parser.add_argument("-a", "--areak", type=float, nargs="?", const=-1, \ help="Area constraint potential strength") parser.add_argument("-fl", "--folder", nargs="?", \ const='/local/duman/SIMULATIONS/Cells_in_LAMMPS/density_0.8/', \ help="Folder containing data, as in /local/duman/SIMULATIONS/Cells_in_LAMMPS/density_0.8/") parser.add_argument("-sb", "--savebase", nargs="?", \ const = "/usr/users/iff_th2/duman/Cells_in_LAMMPS/PLOTS/", \ help="Folder to save the data, as in /usr/users/iff_th2/duman/Cells_in_LAMMPS/PLOTS/") parser.add_argument("-sf", "--savefolder", nargs="?", const="Vel_corr_length", \ help="Specific folder for saving, as in Vel_corr_length") parser.add_argument("-s", "--save_eps", action="store_true", help="Decide whether to save in eps or not") args = parser.parse_args() ### accummulate data as a function of areak -with a dictionary- analysisdatabase = '/usr/users/iff_th2/duman/Cells_in_LAMMPS/DATA/' analysisdatabase += args.savefolder + '/' ### detect the parameter choice param = [] param_choice = '' if args.eps == -1: param_choice = 'eps' param = [0.05, 0.5, 1.0, 5.0, 10.0, 20.0] if args.fp == -1: param_choice = 'fp' param = [0.0, 0.5, 1.0, 3.0, 5.0, 10.0] if args.areak == -1: param_choice = 'areak' param = [1.0, 10.0, 100.0] data = {} # carries the data per parameter set sims = {} # carries the simulation information per parameter set for p in param: if param_choice == 'areak': datafolder, analysisfile = read_write.gen_folders( args.eps, args.fp, p, args.savefolder, args.folder, analysisdatabase) elif param_choice == 'eps': datafolder, analysisfile = read_write.gen_folders( p, args.fp, args.areak, args.savefolder, args.folder, analysisdatabase) elif param_choice == 'fp': datafolder, analysisfile = read_write.gen_folders( args.eps, p, args.areak, args.savefolder, args.folder, analysisdatabase) sims[p] = read_write.read_sim_info(datafolder) x, y = read_write.read_2d_analysis_data(analysisfile) data[p] = y ### plot the data as a function of the parameter plot_data(x, data, param_choice, sims, args.savebase, args.savefolder, args.save_eps) return