def METHOD_NAME():
    rstate = get_rstate()
    g = Gaussian()
    sampler = dynesty.DynamicNestedSampler(g.loglikelihood_with_blob,
                                           g.prior_transform,
                                           g.ndim,
                                           nlive=nlive,
                                           rstate=rstate,
                                           blob=True)
    sampler.run_nested(print_progress=printing, dlogz_init=1)
    res = sampler.results
    assert res['blob'].shape == (len(res['samples']), 3)
    assert np.all(res['blob'] == res['samples'])