Efficiently generate samples from the Polya-Gamma distribution using a NumPy/SciPy compatible interface.