powered by
rejection sampler to update local shrinkage parameters lambda.
rejection_sampler(eps, a, b)