Skip to content
Snippets Groups Projects
Commit 08ad986a authored by Colm Talbot's avatar Colm Talbot
Browse files

Merge branch 'master' into 'project_restructuring'

# Conflicts:
#   examples/open_data_examples/GW150914.py
parents 63745698 cda72bcb
No related branches found
No related tags found
1 merge request!55Project restructuring
Pipeline #
examples/logo/a.jpg

12.5 KiB

examples/logo/image.jpg

1.73 KiB

examples/logo/k.jpg

12.5 KiB

examples/logo/p.jpg

9.83 KiB

""" Script used to generate the samples for the tupak logo """
import tupak
import numpy as np
import scipy.interpolate as si
from skimage import io
class Likelihood(tupak.Likelihood):
def __init__(self, interp):
self.interp = interp
self.parameters = dict(x=None, y=None)
def log_likelihood(self):
return -1/(self.interp(self.parameters['x'], self.parameters['y'])[0])
for letter in ['t', 'u', 'p', 'a', 'k']:
img = 1-io.imread('{}.jpg'.format(letter), as_grey=True)[::-1, :]
x = np.arange(img.shape[0])
y = np.arange(img.shape[1])
interp = si.interpolate.interp2d(x, y, img.T)
likelihood = Likelihood(interp)
priors = {}
priors['x'] = tupak.prior.Uniform(0, max(x), 'x')
priors['y'] = tupak.prior.Uniform(0, max(y), 'y')
result = tupak.run_sampler(
likelihood=likelihood, priors=priors, sampler='nestle', npoints=5000,
label=letter)
fig = result.plot_corner(quantiles=None, smooth1d=4)
examples/logo/t.jpg

11 KiB

examples/logo/u.jpg

11.7 KiB

......@@ -14,3 +14,4 @@ nestle
deepdish
ptemcee
mock
lalsuite
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment