polyagamma | flexible sampler of the Pólya-Gamma distribution | GPU library

import array
import numpy as np
from polyagamma import random_polyagamma

# generate a PG(1, 0) sample
o = random_polyagamma()

# Get a 5 by 1 array of PG(1, 2) variates.
o = random_polyagamma(z=2, size=5)

# We can pass sequences as input. Numpy's b  

from cpython.pycapsule cimport PyCapsule_GetPointer
from polyagamma cimport random_polyagamma_fill, DEVROYE
from numpy.random cimport bitgen_t
import numpy as np

# assuming there exists an instance of the Generator class called `rng`.
bitgenerator =  

$ pip install --pre -U polyagamma

$ conda install -c conda-forge polyagamma

$ git clone https://github.com/zoj613/polyagamma.git
$ cd polya-gamma/
# install dependencies (make sure pip is up to date before hand)
$ poetry install --no-root
$ make in