-
Notifications
You must be signed in to change notification settings - Fork 1
/
run_fast_euler_one_sys.py
127 lines (107 loc) · 5.54 KB
/
run_fast_euler_one_sys.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
import os
import sys
import numpy as np
# Variables to run jobs
json_spec_dir="/scratch/users/tabakg/qsd_output/json_spec/"
output_dir="/scratch/users/tabakg/qsd_output/fast_out/"
dev_dir="/scratch/users/tabakg/qsd_dev"
fast_sim_dir="/scratch/users/tabakg/qsd_dev/fast_euler"
sys.path.append(dev_dir)
# Variables for each job
memory = 16000
# Create subdirectories for job, error, and output files
job_dir = os.path.join(json_spec_dir, ".job")
out_dir = os.path.join(json_spec_dir, ".out")
for new_dir in [output_dir,job_dir,out_dir]:
if not os.path.exists(new_dir):
os.mkdir(new_dir)
OVERWRITE=False
REGIME = 'kerr_bistableA21.75'
delta_ts = [1e-6]
num_systems_arr=[1]
NUM_SEEDS=32
ntraj=1
Nfock_a=30
Nfock_j=2
duration=50.0
downsample_dict = {1e-6:1000}
partition = 'normal'
partition_dict = {1e-6: False}
noise_amp=1.0
trans_phase=1.0
Rs=[0.0]
min_EPSs = [0. if r ==0. else (1-np.sqrt(1-r**2))/r for r in Rs]
EPSs = {R: [E] for R, E in zip(Rs, min_EPSs)}
lambds_dict={R : [0.0] for R in Rs}
SEEDs=range(1, NUM_SEEDS + 1)
DRIVES=False, ##Driving second system?
sdeint_method_name = 'itoImplicitEuler'
for R in Rs:
for EPS in EPSs[R]:
for lambd in lambds_dict[R]:
for num_systems in num_systems_arr:
for delta_t in delta_ts:
for seed in SEEDs:
for drive_second_system in DRIVES:
downsample=downsample_dict[delta_t]
long_run = partition_dict[delta_t]
args = (REGIME,
seed,
ntraj,
delta_t,
Nfock_a,
Nfock_j,
duration,
downsample,
sdeint_method_name,
num_systems,
R,
EPS,
noise_amp,
lambd,
trans_phase,
drive_second_system)
param_str = make_params_string(args)
json_spec_name = "json_spec_" + param_str + ".json"
json_spec_loc = os.path.join(json_spec_dir, json_spec_name)
json_spec_exists = os.path.isfile(json_spec_loc)
psis_name = "psis_" + param_str + ".json"
expects_name = "expects_" + param_str + ".json"
psis_loc = os.path.join(output_dir, psis_name)
expects_loc = os.path.join(output_dir, expects_name)
print("OVERWRITE is %s and file %s existence is %s" %(OVERWRITE, json_spec_name, json_spec_exists))
print("If overwriting or file does not exist, going to process new seed.")
if OVERWRITE or not json_spec_exists:
print ("Processing seed %s" %(seed))
## Write job to file
filey_loc = os.path.join(job_dir, "%s.job" %(seed))
filey = open(filey_loc, "w")
filey.writelines("#!/bin/bash\n")
filey.writelines("#SBATCH --job-name=QSD\n")
filey.writelines("#SBATCH --output=%s/qsd_%s.out\n" %(out_dir,seed))
filey.writelines("#SBATCH --error=%s/qsd_%s.err\n" %(out_dir,seed))
if long_run:
filey.writelines("#SBATCH --time=7-00:00\n")
else:
filey.writelines("#SBATCH --time=2-00:00\n")
filey.writelines("#SBATCH --mem=%s\n" %(memory))
script_name = os.path.join(dev_dir, "generate_num_model.py")
fast_sim_name = os.path.join(fast_sim_dir, "fast_sim")
## Python script
to_file = ("python %s --output_dir '%s' --Nfock_a %s --ntraj %s "
"--seed %s --regime '%s' --num_systems %s "
"--delta_t %s --duration %s --downsample %s --sdeint_method_name '%s' "
"--R %s --eps %s --noise_amp %s --lambda %s"
%(script_name, json_spec_dir, Nfock_a, ntraj,
seed, REGIME, num_systems,
delta_t, duration, downsample, sdeint_method_name,
R, EPS, noise_amp, lambd))
if drive_second_system:
to_file = to_file + " --drive_second_system %s" % drive_second_system
to_file += "\n"
filey.writelines(to_file)
## C++ code
filey.writelines("%s %s %s %s" %(fast_sim_name, json_spec_loc, psis_loc, expects_loc))
filey.close()
## run batch script
os.system("sbatch -p %s %s" %(partition, filey_loc))