After poking around a bit, I found a way to do this:
Code:
import cobayaimport numpy as npfrom cobaya.yaml import yaml_loadfrom cobaya.model import get_modelyaml_file='./XXX.yaml' #your yaml fileinfo = yaml_load(yaml_file)model = get_model(info)point = dict(zip(model.parameterization.sampled_params(),model.prior.sample(ignore_external=True)[0]))logposterior = model.logposterior(point, as_dict=True)Cls = model.provider.get_Cl(ell_factor=True)np.savez('Cls.npz', Cls)Statistics: Posted by Tristan L Smith — April 09 2025
