[Numpy-discussion] Complex matrix seg fault

Jonathan Tu jhtu at princeton.edu
Fri Feb 4 09:13:29 EST 2011


Hi,

When I run the following code, I get a seg fault on certain cluster nodes, but not others:

>>> import numpy as N
>>> A = N.random.rand(262144,10)+1j*N.random.rand(262144,10)
>>> x = N.random.rand(262144,)
>>> d, res, rank, s = N.linalg.lstsq( A, x )

My numpy is configured as such:

lapack_opt_info:
    libraries = ['mkl_lapack', 'mkl_intel_lp64', 'mkl_intel_thread', 'mkl_core', 'mkl_def', 'mkl_mc', 'iomp5', 'pthread']
    library_dirs = ['/home/jhtu/local/lib/mkl_links']
    define_macros = [('SCIPY_MKL_H', None)]
    include_dirs = ['/opt/intel/Compiler/11.1/072/mkl/include']
blas_opt_info:
    libraries = ['mkl_intel_lp64', 'mkl_intel_thread', 'mkl_core', 'mkl_def', 'mkl_mc', 'iomp5', 'pthread']
    library_dirs = ['/home/jhtu/local/lib/mkl_links']
    define_macros = [('SCIPY_MKL_H', None)]
    include_dirs = ['/opt/intel/Compiler/11.1/072/mkl/include']
lapack_mkl_info:
    libraries = ['mkl_lapack', 'mkl_intel_lp64', 'mkl_intel_thread', 'mkl_core', 'mkl_def', 'mkl_mc', 'iomp5', 'pthread']
    library_dirs = ['/home/jhtu/local/lib/mkl_links']
    define_macros = [('SCIPY_MKL_H', None)]
    include_dirs = ['/opt/intel/Compiler/11.1/072/mkl/include']
blas_mkl_info:
    libraries = ['mkl_intel_lp64', 'mkl_intel_thread', 'mkl_core', 'mkl_def', 'mkl_mc', 'iomp5', 'pthread']
    library_dirs = ['/home/jhtu/local/lib/mkl_links']
    define_macros = [('SCIPY_MKL_H', None)]
    include_dirs = ['/opt/intel/Compiler/11.1/072/mkl/include']
mkl_info:
    libraries = ['mkl_intel_lp64', 'mkl_intel_thread', 'mkl_core', 'mkl_def', 'mkl_mc', 'iomp5', 'pthread']
    library_dirs = ['/home/jhtu/local/lib/mkl_links']
    define_macros = [('SCIPY_MKL_H', None)]
    include_dirs = ['/opt/intel/Compiler/11.1/072/mkl/include']

All unit tests pass except the following, but I was told (by this mailing list) that this shouldn't be an issue:

======================================================================
FAIL: test_special_values (test_umath_complex.TestClog)
----------------------------------------------------------------------
Traceback (most recent call last):
  File "/home/jhtu/local/lib/python2.7/site-packages/numpy/testing/decorators.py", line 215, in knownfailer
    return f(*args, **kwargs)
  File "/home/jhtu/local/lib/python2.7/site-packages/numpy/core/tests/test_umath_complex.py", line 312, in test_special_values
    assert_almost_equal(np.log(np.conj(xa[i])), np.conj(np.log(xa[i])))
  File "/home/jhtu/local/lib/python2.7/site-packages/numpy/testing/utils.py", line 443, in assert_almost_equal
    raise AssertionError(msg)
AssertionError: 
Arrays are not almost equal
 ACTUAL: array([-inf+3.14159265j])
 DESIRED: array([-inf-3.14159265j])

----------------------------------------------------------------------
Ran 2648 tests in 15.852s

FAILED (KNOWNFAIL=4, SKIP=1, failures=1)
<nose.result.TextTestResult run=2648 errors=0 failures=1>

Anyone have ideas why this would fail on some nodes but not others?  (I mean the code snippet fails on some nodes.  The above unit test always fails.)






Thanks,



Jonathan Tu





More information about the NumPy-Discussion mailing list