Maintenance will be performed on git.ligo.org, chat.ligo.org, containers.ligo.org, and docs.ligo.org on Tuesday 26 May 2020 starting at approximately 10am CDT. It is expected to take around 30 minutes and will involve a short period of downtime, around 5 minutes, towards the end of the maintenance period. Please address any questions, comments, or concerns to uwm-help@cgca.uwm.edu.

Commit ca5ada8a authored by Richard O'Shaughnessy's avatar Richard O'Shaughnessy

ILE-no-op: change width of toy model mass distribution, fix output bug

parent 9c2ee3d4
......@@ -271,7 +271,7 @@ Sigma_list = []
rv_list = []
sigma_1d = 0.1
x_0 = numpy.array([35,30,0.0,0.0,0.0,0,0,0])
Sigma= numpy.diag( [1,1, 0.01, 0.01, 0.01, 0.01,0.01,0.01])
Sigma= numpy.diag( [16,16, 0.01, 0.01, 0.01, 0.01,0.01,0.01])
print x_0.shape, Sigma.shape
Sigma_list.append(Sigma)
rv_list.append( scipy.stats.multivariate_normal(mean=x_0, cov=Sigma)) # one-dimensional in one parameter
......@@ -306,8 +306,8 @@ for param in numpy.arange(n_gauss):
# print rv_list[param].pdf(x_now)
Lval+= rv_list[param].pdf(x_now) * wt_array[param]
res = numpy.log(Lval) + 250 # manual offset
var = 0.1
res = Lval*numpy.exp(250) # manual offset
var = 0.1*res*res
neff = 100
ntotal=1e5
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment