BUG: Use lop1p to improve numerical precision

Use log1p(-x) instead of log(1 - x)
Seperate legacy version from current

closes #17020

1 files changed, 8 insertions, 7 deletions

diff --git a/numpy/random/_examples/cython/setup.py b/numpy/random/_examples/cython/setup.py
index 83f06fde8..7e0dd3e05 100644
--- a/numpy/random/_examples/cython/setup.py
+++ b/numpy/random/_examples/cython/setup.py
@@ -4,19 +4,20 @@ Build the Cython demonstrations of low-level access to NumPy random
 
 Usage: python setup.py build_ext -i
 """
+from distutils.core import setup
+from os.path import dirname, join, abspath
 
 import numpy as np
-from distutils.core import setup
 from Cython.Build import cythonize
+from numpy.distutils.misc_util import get_info
 from setuptools.extension import Extension
-from os.path import join, dirname
 
 path = dirname(__file__)
 src_dir = join(dirname(path), '..', 'src')
 defs = [('NPY_NO_DEPRECATED_API', 0)]
 inc_path = np.get_include()
-# not so nice. We need the random/lib library from numpy
-lib_path = join(np.get_include(), '..', '..', 'random', 'lib')
+lib_path = [abspath(join(np.get_include(), '..', '..', 'random', 'lib'))]
+lib_path += get_info('npymath')['library_dirs']
 
 extending = Extension("extending",
                       sources=[join('.', 'extending.pyx')],
@@ -29,10 +30,10 @@ extending = Extension("extending",
 distributions = Extension("extending_distributions",
                           sources=[join('.', 'extending_distributions.pyx')],
                           include_dirs=[inc_path],
-                          library_dirs=[lib_path],
-                          libraries=['npyrandom'],
+                          library_dirs=lib_path,
+                          libraries=['npyrandom', 'npymath'],
                           define_macros=defs,
-                         )
+                          )
 
 extensions = [extending, distributions]