fig = p.figure(figsize=(3.375, 2.0))
data = []
column_dic = []
for file_name in data_file_list:
param_dict = extract_params_from_file_name(file_name)
full_path = os.path.join(prefix, file_name)
if not os.path.exists(full_path):
print('Path {} doesn\'t exist'.format(full_path))
continue
data.append(
np.loadtxt(full_path,
dtype=np.float64,
delimiter=';',
unpack=True,
skiprows=1))
column_names = read_column_names(full_path)
column_dic.append(dict(zip(column_names, range(len(column_names)))))
p.subplot(1, 2, 1)
handle1, = p.loglog(data[0][column_dic[0]['NAtoms']],
1 - data[0][column_dic[0]['F_CJ_tNoInteraction_one']],
color='#009900',
linestyle=':',
label=r'$t_{\rm b}=1$, $\Lambda$-type',
linewidth=common_line_width)
handle2, = p.loglog(data[1][column_dic[1]['NAtoms']],
1 - data[1][column_dic[1]['F_CJ_tNoInteraction_one']],
'k-.',
label=r'$t_{\rm b}=1$, dual-V',
linewidth=common_line_width)
handle3, = p.loglog(data[0][column_dic[0]['NAtoms']],
data_dualV_file_name = 'grating_dispersion_relation_dualV_N_10000_g1d_0.1_Deltac_-90_Omega_1_OmegaPeriods_5000_seed_12345.txt'
param_dict_dualV = extract_params_from_file_name(data_dualV_file_name)
full_path_dualV = os.path.join(prefix, data_dualV_file_name)
full_path_lambda = os.path.join(prefix, data_lambda_file_name)
data_dualV = p.loadtxt(full_path_dualV,
dtype=p.float64,
delimiter=';',
unpack=True,
skiprows=1)
data_lambda = p.loadtxt(full_path_lambda,
dtype=p.float64,
delimiter=';',
unpack=True,
skiprows=1)
column_names_lambda = read_column_names(full_path_lambda)
column_dic_lambda = dict(
zip(column_names_lambda, range(len(column_names_lambda))))
column_names_dualV = read_column_names(full_path_dualV)
column_dic_dualV = dict(zip(column_names_dualV,
range(len(column_names_dualV))))
g1d = 0.1
Deltac = -90
Omega = 1
usetex()
p.figure(figsize=(3.3, 2.5))
ax = p.subplot(111)
delta_array_log = np.linspace(-6, 0, 500)