# sky sampler: cos(dec) uniform in [-1, 1), adaptive sampling # dec_sampler = mcsampler.dec_samp_vector sampler.add_parameter("declination", \ pdf = dec_sampler, \ cdf_inv = None, \ left_limit = param_limits["declination"][0], \ right_limit = param_limits["declination"][1], \ prior_pdf = mcsampler.uniform_samp_dec, \ adaptive_sampling = not opts.no_adapt) # # Determine pinned and non-pinned parameters # pinned_params = common_cl.get_pinned_params(opts) unpinned_params = common_cl.get_unpinned_params(opts, sampler.params) print "{0:<25s} {1:>5s} {2:>5s} {3:>20s} {4:<10s}".format("parameter", "lower limit", "upper limit", "pinned?", "pin value") plen = len(sorted(sampler.params, key=lambda p: len(p))[-1]) for p in sampler.params: if pinned_params.has_key(p): pinned, value = True, "%1.3g" % pinned_params[p] else: pinned, value = False, "" if isinstance(p, tuple): for subp, subl, subr in zip(p, sampler.llim[p], sampler.rlim[p]): subp = subp + " "*min(0, plen-len(subp)) print "|{0:<25s} {1:>1.3g} {2:>1.3g} {3:>20s} {4:<10s}".format(subp, subl, subr, str(False), "") else: p = p + " "*min(0, plen-len(p))
# sky sampler: cos(dec) uniform in [-1, 1), adaptive sampling # dec_sampler = mcsampler.dec_samp_vector sampler.add_parameter("declination", \ pdf = dec_sampler, \ cdf_inv = None, \ left_limit = param_limits["declination"][0], \ right_limit = param_limits["declination"][1], \ prior_pdf = mcsampler.uniform_samp_dec, \ adaptive_sampling = not opts.no_adapt) # # Determine pinned and non-pinned parameters # pinned_params = common_cl.get_pinned_params(opts) unpinned_params = common_cl.get_unpinned_params(opts, sampler.params) print "{0:<25s} {1:>5s} {2:>5s} {3:>20s} {4:<10s}".format( "parameter", "lower limit", "upper limit", "pinned?", "pin value") plen = len(sorted(sampler.params, key=lambda p: len(p))[-1]) for p in sampler.params: if pinned_params.has_key(p): pinned, value = True, "%1.3g" % pinned_params[p] else: pinned, value = False, "" if isinstance(p, tuple): for subp, subl, subr in zip(p, sampler.llim[p], sampler.rlim[p]): subp = subp + " " * min(0, plen - len(subp)) print "|{0:<25s} {1:>1.3g} {2:>1.3g} {3:>20s} {4:<10s}".format( subp, subl, subr, str(False), "")