elif args.params == "emcee":
eparams = np.median(np.load("eparams_emcee.npy"), axis=0)
print("Using emcee median")
else:
import sys
sys.exit()
# Print out the emulator parameters in an easily-readable format
lambda_xi = eparams[0]
hparams = eparams[1:].reshape((my_pca.m, -1))
print("Emulator parameters are:")
print("lambda_xi", lambda_xi)
for row in hparams:
print(row)
emulator = Emulator(my_pca, eparams)
# We will want to produce interpolated plots spanning each parameter dimension,
# for each eigenspectrum.
# Create a list of parameter blocks.
# Go through each parameter, and create a list of all parameter combination of
# the other two parameters.
unique_points = [
np.unique(my_pca.gparams[:, i]) for i in range(len(Starfish.parname))
]
blocks = []
for ipar, pname in enumerate(Starfish.parname):
upars = unique_points.copy()
dim = upars.pop(ipar)
ndim = len(dim)
