PolyaGamma sampling wrapper for Python.  R/C++ code from Jesse Windle