[SciPy-Dev] Gmres with Sparse Complex Matrices gives Segmentation fault (or NULL result)

[Martin Campos Pinto]

Hi, thank you for the quick answer.

The run `scipy.test("full", verbose=2)' fails with the following error 
message (I'm only pasting the last few lines, but I can send the whole 
output if necessary...)


======================================================================
ERROR: Failure: ImportError 
(dlopen(/Library/Frameworks/EPD64.framework/Versions/7.2/lib/python2.7/site-packages/scipy/linalg/atlas_version.so, 
2): Symbol not found: _ATL_buildinfo
   Referenced from: 
/Library/Frameworks/EPD64.framework/Versions/7.2/lib/python2.7/site-packages/scipy/linalg/atlas_version.so
   Expected in: dynamic lookup
)
----------------------------------------------------------------------
Traceback (most recent call last):
   File 
"/Library/Frameworks/EPD64.framework/Versions/7.2/lib/python2.7/site-packages/nose/loader.py", 
line 390, in loadTestsFromName
     addr.filename, addr.module)
   File 
"/Library/Frameworks/EPD64.framework/Versions/7.2/lib/python2.7/site-packages/nose/importer.py", 
line 39, in importFromPath
     return self.importFromDir(dir_path, fqname)
   File 
"/Library/Frameworks/EPD64.framework/Versions/7.2/lib/python2.7/site-packages/nose/importer.py", 
line 86, in importFromDir
     mod = load_module(part_fqname, fh, filename, desc)
   File 
"/Library/Frameworks/EPD64.framework/Versions/7.2/lib/python2.7/site-packages/scipy/linalg/tests/test_atlas_version.py", 
line 6, in <module>
     import scipy.linalg.atlas_version
ImportError: 
dlopen(/Library/Frameworks/EPD64.framework/Versions/7.2/lib/python2.7/site-packages/scipy/linalg/atlas_version.so, 
2): Symbol not found: _ATL_buildinfo
   Referenced from: 
/Library/Frameworks/EPD64.framework/Versions/7.2/lib/python2.7/site-packages/scipy/linalg/atlas_version.so
   Expected in: dynamic lookup


======================================================================
FAIL: test_dot (test_blas.TestFBLAS1Simple)
----------------------------------------------------------------------
Traceback (most recent call last):
   File 
"/Library/Frameworks/EPD64.framework/Versions/7.2/lib/python2.7/site-packages/scipy/lib/blas/tests/test_blas.py", 
line 71, in test_dot
     assert_almost_equal(f([3j,-4,3-4j],[2,3,1]),-9+2j)
   File 
"/Library/Frameworks/EPD64.framework/Versions/7.2/lib/python2.7/site-packages/numpy/testing/utils.py", 
line 448, in assert_almost_equal
     raise AssertionError(msg)
AssertionError:
Arrays are not almost equal to 7 decimals
  ACTUAL: 0j
  DESIRED: (-9+2j)

======================================================================
FAIL: test_complex_dotc (test_blas.TestFBLAS1Simple)
----------------------------------------------------------------------
Traceback (most recent call last):
   File 
"/Library/Frameworks/EPD64.framework/Versions/7.2/lib/python2.7/site-packages/scipy/linalg/tests/test_blas.py", 
line 121, in test_complex_dotc
     assert_almost_equal(f([3j,-4,3-4j],[2,3j,1]),3-14j)
   File 
"/Library/Frameworks/EPD64.framework/Versions/7.2/lib/python2.7/site-packages/numpy/testing/utils.py", 
line 448, in assert_almost_equal
     raise AssertionError(msg)
AssertionError:
Arrays are not almost equal to 7 decimals
  ACTUAL: 0j
  DESIRED: (3-14j)

======================================================================
FAIL: test_complex_dotu (test_blas.TestFBLAS1Simple)
----------------------------------------------------------------------
Traceback (most recent call last):
   File 
"/Library/Frameworks/EPD64.framework/Versions/7.2/lib/python2.7/site-packages/scipy/linalg/tests/test_blas.py", 
line 115, in test_complex_dotu
     assert_almost_equal(f([3j,-4,3-4j],[2,3,1]),-9+2j)
   File 
"/Library/Frameworks/EPD64.framework/Versions/7.2/lib/python2.7/site-packages/numpy/testing/utils.py", 
line 448, in assert_almost_equal
     raise AssertionError(msg)
AssertionError:
Arrays are not almost equal to 7 decimals
  ACTUAL: 0j
  DESIRED: (-9+2j)

----------------------------------------------------------------------
Ran 5045 tests in 899.955s

FAILED (KNOWNFAIL=14, SKIP=29, errors=1, failures=3)





Le 11/24/12 10:40 AM, Pauli Virtanen a écrit :
> 24.11.2012 07:50, Martin Campos Pinto kirjoitti:
> [clip]
>> I have looked on the web for solutions but haven't found any. Some very old posts
>> indicate similar errors but they don't come with an answer,
>> and I imagine that if those were due to bugs, they should have been fixed by now...
>>
>> Am I doing something stupid here, or is that a real problem ? Is somebody aware of a solution ?
>> (I am using scipy version 0.10.1)
> I don't get a crash with this code example. Valgrind also does not show
> anything, so I suppose this is something occurring only on OSX.
>
> There are some known issues on OSX with Accelerate, but I'm not sure if
> they are in play here.
>
> Can you run `scipy.test("full", verbose=2)` without seeing problems?
>

-- 
Martin Campos Pinto
http://www-irma.u-strasbg.fr/~campos