From a649c0e6c363957bf122ec4ac7d5534eb819d260 Mon Sep 17 00:00:00 2001
From: Rayne Liu <rayne.liu@atlas1>
Date: Tue, 13 Oct 2020 12:28:16 +0000
Subject: [PATCH] Minor updates
---
code/Mock_Interpolate-0001_t_10M_wandt.py | 12 ++++++------
...condor_submit_MockInterpolate0001_t_10M_wandt.sub | 2 +-
2 files changed, 7 insertions(+), 7 deletions(-)
mode change 100644 => 100755 code/Mock_Interpolate-0001_t_10M_wandt.py
mode change 100644 => 100755 code/condor_submit_MockInterpolate0001_t_10M_wandt.sub
diff --git a/code/Mock_Interpolate-0001_t_10M_wandt.py b/code/Mock_Interpolate-0001_t_10M_wandt.py
old mode 100644
new mode 100755
index aa87b59..c6801a9
--- a/code/Mock_Interpolate-0001_t_10M_wandt.py
+++ b/code/Mock_Interpolate-0001_t_10M_wandt.py
@@ -42,12 +42,12 @@ tshift=10
vary_fund = True
#sampler parameters
-npoints = 20#1200
-nwalkers = 10#256
+npoints = 220
+nwalkers = 100
ntemps=12
dim = nmax+1
ndim = 4*dim
-burnin = 10#600 #How many points do you burn before doing the corner plot. You need to watch the convergence of the chain plot a bit.
+burnin = 150 #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 = 32 #corner plot parameter - how many bins you want
datacolor = '#105670' #'#4fa3a7'
@@ -155,8 +155,8 @@ def log_prior(theta):
return 0.0
elif tshift == 10:
- if all([0.5 <= omega0 <= 0.6, 0.5 <= omega1 <= 0.6, 9 <= tau0 <= 15, 0. <= tau1 <= 9., \
- 0. <= xvar0 <= 3., 0. <= xvar1 <= 6., -np.pi <= yvar0 <= np.pi, -np.pi <= yvar1 <= np.pi]):
+ if all([0.52 <= omega0 <= 0.64, 0.43 <= omega1 <= 0.58, 4. <= tau0 <= 18., 0. <= tau1 <= 10., \
+ 0. <= xvar0 <= 1., 0. <= xvar1 <= 1.5, -np.pi <= yvar0 <= np.pi, -np.pi <= yvar1 <= np.pi]):
return 0.0
return -np.inf
@@ -188,7 +188,7 @@ def log_probability(theta):
#pool = choose_pool(1)
#pool.size = 1
np.random.seed(42)
-pos = np.array([random.uniform(0.55,0.57), random.uniform(0.5,0.52), random.uniform(10., 13.7), random.uniform(4.,7.), random.uniform(0.3,0.5), random.uniform(0.3, 0.5), random.uniform(-2., 2.), random.uniform(-2., 2.)])
+pos = np.array([random.uniform(0.57,0.63), random.uniform(0.5,0.54), random.uniform(8., 13.7), random.uniform(4.,7.), random.uniform(0.3,0.5), random.uniform(0.3, 0.5), random.uniform(-1., 1.), random.uniform(-1., 1.)])
pos = list(pos)
pos += 1e-5 * np.random.randn(ntemps, nwalkers, ndim)
with Pool() as pool:
diff --git a/code/condor_submit_MockInterpolate0001_t_10M_wandt.sub b/code/condor_submit_MockInterpolate0001_t_10M_wandt.sub
old mode 100644
new mode 100755
index 59f4edd..15e87f8
--- a/code/condor_submit_MockInterpolate0001_t_10M_wandt.sub
+++ b/code/condor_submit_MockInterpolate0001_t_10M_wandt.sub
@@ -12,7 +12,7 @@ initialdir = .
notify_user = rl746@cornell.edu
notification = Complete
arguments = "-processid $(Process)"
-request_memory = 64GB
+request_memory = 192GB
request_cpus = 1
on_exit_remove = (ExitBySignal == False) || ((ExitBySignal == True) && (ExitSignal != 11))
accounting_group = aei.dev.test_dynesty
--
GitLab