In [19]:
from time import time
from multiprocessing import cpu_count, Pool
import warnings
from itertools import izip
from os import path
from ast import literal_eval
import numpy as np
import emcee as mc
from scipy.linalg import inv
import h5py
from pearce.emulator import OriginalRecipe, ExtraCrispy, SpicyBuffalo, NashvilleHot
In [20]:
lnprob = lambda x: 0
In [21]:
num_params = 4
nwalkers = 50
pos0 = np.random.randn(nwalkers, num_params)
In [22]:
from pearce.inference import run_mcmc_config
In [31]:
fname = '/scratch/users/swmclau2/PearceMCMC/pearce_mcmc_nh_gg.hdf5'
In [40]:
run_mcmc_config(fname)
In [44]:
import h5py
In [53]:
f = h5py.File(fname, 'r')
In [54]:
f.attrs['nwalkers']
Out[54]:
In [36]:
f.attrs['nsteps'] = 5
In [37]:
f.attrs.keys()
Out[37]:
In [43]:
f['chain'].value
Out[43]:
In [55]:
f.close()
In [ ]:
f.attrs