From ea0c03d15e2b590cfd13487093bf7ab4d31029e6 Mon Sep 17 00:00:00 2001
From: Rayne Liu <rayne.liu@atlas1>
Date: Wed, 30 Sep 2020 12:47:30 +0000
Subject: [PATCH] 8 cores
---
code/NR_Interpolate-0001_t_10M.py | 10 +++++-----
code/condor_submit_Interpolate0001_t_10M.sub | 0
2 files changed, 5 insertions(+), 5 deletions(-)
mode change 100644 => 100755 code/NR_Interpolate-0001_t_10M.py
mode change 100644 => 100755 code/condor_submit_Interpolate0001_t_10M.sub
diff --git a/code/NR_Interpolate-0001_t_10M.py b/code/NR_Interpolate-0001_t_10M.py
old mode 100644
new mode 100755
index 6622dd7..72af4e5
--- a/code/NR_Interpolate-0001_t_10M.py
+++ b/code/NR_Interpolate-0001_t_10M.py
@@ -39,12 +39,12 @@ tshift=10
vary_fund = True
#sampler parameters
-npoints = 100
-nwalkers = 50
+npoints = 600
+nwalkers = 300
ntemps=16
dim = nmax+1
ndim = 4*dim
-burnin = 50 #How many points do you burn before doing the corner plot. You need to watch the convergence of the chain plot a bit.
+burnin = 300 #How many points do you burn before doing the corner plot. You need to watch the convergence of the chain plot a bit.
#This is trivial but often forgotten: this cannot be more than npoints! I usually use half the points.
numbins = 42 #corner plot parameter - how many bins you want
datacolor = '#105670' #'#4fa3a7'
@@ -175,7 +175,7 @@ def log_prior(theta):
if all([-0.06 <= alpha0 <= 0.06, -0.32 <= alpha1 <= -0.08, -0.19 <= beta0 <= 1.0, 0. <= beta1 <= 1.5, 0 <= xvar0 <= 1.1, 0 <= xvar1 <= 1.2, -np.pi <= yvar0 <= np.pi, -np.pi <= yvar1 <= np.pi]):
return 0.0
elif tshift == 10:
- if all([-0.08 <= alpha0 <= 0.12, -0.6 <= alpha1 <= 0.4, -0.3 <= beta0 <= 1.5, -1. <= beta1 <= 3.2, 0 <= xvar0 <= 1.1, 0 <= xvar1 <= 1.2, -np.pi <= yvar0 <= np.pi, -np.pi <= yvar1 <= np.pi]):
+ if all([-0.04 <= alpha0 <= 0.04, -0.1 <= alpha1 <= 0.15, -0.4 <= beta0 <= 0.4, -1. <= beta1 <= 2, 0 <= xvar0 <= 1.1, 0 <= xvar1 <= 1.2, -np.pi <= yvar0 <= np.pi, -np.pi <= yvar1 <= np.pi]):
return 0.0
"""
if nmax == 0:
@@ -298,4 +298,4 @@ for yi in range(naxes):
ax.axvline(median[xi], color=mediancolor)
ax.axhline(median[yi], color=mediancolor)
ax.plot(median[xi], median[yi], color = mediancolor, marker = 's')
-figcorn.savefig(rootpath + '/plotsmc/0001_10M_interpolated_cornerplot.png', format='png', bbox_inches='tight', dpi=300)
\ No newline at end of file
+figcorn.savefig(rootpath + '/plotsmc/0001_10M_interpolated_cornerplot.png', format='png', bbox_inches='tight', dpi=300)
diff --git a/code/condor_submit_Interpolate0001_t_10M.sub b/code/condor_submit_Interpolate0001_t_10M.sub
old mode 100644
new mode 100755
--
GitLab