for seed in all_seeds:
print("seed = " + str(seed))
local_path = f"/regular/{system}/{task}/N_{N}/alpha_{alpha:0.1f}/seed_{seed}"
data_path = get_root() + local_path
config_list = []
config_list.append('[lind_ham]')
config_list.append('gen_type = ' + str(gen_type))
config_list.append('N = ' + str(N))
config_list.append('alpha = ' + str(alpha))
config_list.append('seed = ' + str(seed))
config_list.append('num_seeds = ' + str(num_seeds))
config_list += get_regular_global_config(system,
task,
save_rho='true',
name_precision=4)
pathlib.Path(data_path).mkdir(parents=True, exist_ok=True)
file_config = open(data_path + '/config.ini', 'w')
file_config.write('\n'.join(config_list))
fn_suffix = f"gen({gen_type})_N({N})_alpha({alpha:0.4f})_seed({seed})"
fn_test = data_path + '/lindbladian_evals_' + fn_suffix + '.txt'
if not os.path.isfile(fn_test):
print('file does not exist')
if segment == 'short':
os.system('sbatch run_short.sh ' + data_path)
data_path = get_root() + local_path
config_list = []
config_list.append('[mbl]')
config_list.append('num_spins = ' + str(num_spins))
config_list.append('seed = ' + str(seed))
config_list.append('num_seeds = ' + str(num_seeds))
config_list.append('diss_type = ' + str(diss_type))
config_list.append('diss_phase = ' + str(diss_phase))
config_list.append('diss_gamma = ' + str(diss_gamma))
config_list.append('W = ' + str(W))
config_list.append('U = ' + str(U))
config_list.append('J = ' + str(J))
config_list += get_regular_global_config(system,
task,
save_rho='false')
if task == 'odeint':
config_list += get_odeint_config(step, total_num_periods,
current_num_periods,
current_num_time_points,
is_continue, data_path + '/')
elif task == 'smallest_eigen_vector':
config_list += get_smallest_eigen_vector_config(
max_num_iterations, tolerance)
elif task == 'all_evals':
config_list += get_all_evals_config(max_num_iterations, tolerance)
pathlib.Path(data_path).mkdir(parents=True, exist_ok=True)
file_config = open(data_path + '/config.ini', 'w')
data_path = get_root() + local_path
config_list = []
config_list.append('[dimer]')
config_list.append('num_particles = ' + str(N))
config_list.append('diss_type = ' + str(diss_type))
config_list.append('diss_gamma = ' + str(diss_gamma))
config_list.append('E = ' + str(E))
config_list.append('U = ' + str(U))
config_list.append('J = ' + str(J))
config_list.append('drv_type = ' + str(drv_type))
config_list.append('drv_ampl = ' + str(drv_ampl))
config_list.append('drv_freq = ' + str(drv_freq))
config_list.append('drv_phase = ' + str(drv_phase))
config_list += get_regular_global_config(system, task)
if task == 'odeint':
config_list += get_odeint_config(step, total_num_periods,
current_num_periods,
current_num_time_points,
is_continue, data_path + '/')
elif task == 'smallest_eigen_vector':
config_list += get_smallest_eigen_vector_config(
max_num_iterations, tolerance)
elif task == 'all_evals':
config_list += get_all_evals_config(max_num_iterations, tolerance)
pathlib.Path(data_path).mkdir(parents=True, exist_ok=True)
file_config = open(data_path + '/config.ini', 'w')
file_config.write('\n'.join(config_list))