BUG: use scipy.special.expm1 instead of np.expm1