R_avg0s = np.linspace(5, 12, num=8, endpoint=True)
sigmas = np.array([np.array([0, 0.001, 0.01, 0.1, 1.0], float)[i_s]],
float)
#fig, axarr = plt.subplots(2,4)
#fig.set_size_inches(width,height)
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
Rs = np.linspace(4, 13, num=100, endpoint=True)
ax.plot(Rs, 0 * np.ones(len(Rs)), chi_0(Rs), 'k-')
scan = {}
markers = ["o-", "v-", "^-", "<-", ">-"]
for i, R_avg0 in enumerate(R_avg0s):
scan['R_avg0'] = str(R_avg0)
scan['chi_0'] = chi_0(R_avg0)
for j, sigma in enumerate(sigmas):
scan['sigma_R0'] = str(sigma)
sigmas = np.array([0,0.001,0.01,0.1,1.0],float)
fig, axarr = plt.subplots(2,4)
fig.set_size_inches(width,height)
scan = {}
markers = ["o-","v-","^-","<-",">-"]
for i,R_avg0 in enumerate(R_avg0s):
scan['R_avg0'] = str(R_avg0)
scan['chi_0'] = chi_0(R_avg0)
for j,sigma in enumerate(sigmas):
scan['sigma_R0'] = str(sigma)
rp = ReadParams(scan=scan,datfile=datfile)
ts = rp.list_of_t_vals()
Ns = np.empty([len(ts)],float)
for i_t,t in enumerate(ts):
datfile = data_path + "/data/input.dat"
R_avg0s = np.linspace(5, 12, num=8, endpoint=True)
sigmas = np.array([0, 0.001, 0.01, 0.1, 1.0], float)
scan = {}
markers = ["o-", "v-", "^-", "<-", ">-"]
for i, R_avg0 in enumerate(R_avg0s):
scan['R_avg0'] = str(R_avg0)
scan['chi_0'] = str(chi_0(R_avg0))
for j, sigma in enumerate(sigmas):
scan['sigma_R0'] = str(sigma)
rp = ReadParams(scan=scan, datfile=datfile)
ts = rp.list_of_t_vals()
print(len(ts))
fig, axarr = plt.subplots(3, len(ts) // 3)
fig.set_size_inches(width, height)