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
