A high-performance sampler for Polya-Gamma random variates. - View it on GitHub