An efficient and flexible sampler of the Pólya-Gamma distribution with a NumPy/SciPy compatible interface.