Approximate algorithm for horseshoe prior
