def _pdf(self, x, b):
        # We use (x**2 + b**2)/2 = ((x-b)**2)/2 + xb.
        # The factor of exp(-xb) is then included in the i0e function
        # in place of the modified Bessel function, i0, improving
        # numerical stability for large values of xb.
        return x * exp(-(x-b)*(x-b)/2.0) * special.i0e(x*b)
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/rice_gen._pdf/original