[devel] atlas is broken again?

Valery Pipin pip на iszf.irk.ru
Вт Апр 5 17:51:13 UTC 2011


05.04.2011 01:20, REAL пишет:
> 05.04.2011 15:02, Valery Pipin пишет:
>>>> Отсюда, http://www.enthought.com/
>>>
>>> А зачем, если есть в сизифе?
>>>
>> сизифовский медленный, и к тому же периодически ломается.
>
> Поломок не замечал, а то, что MKL быстрее - факт. Увы, его в сизиф не 
> положишь :(
>
> [...]> Возвращаемся к первому сообщению. Вот это с  libatlas-3.9.35:
>
>> In [1]: from numpy.random import *
>>
>> In [2]: from numpy.linalg import *
>>
>> In [3]: m=randn(2000,2000)
>>
>> In [4]: %time U,s,V=svd(m)
>> Ошибка сегментирования
>
> Ну нет у меня такого нигде, ни на x86_64, ни на i586.
Вот  сегодня обновил рабочую станцию, там те же глюки.
>
>> Теперь пробуем для массивов меньших размеров:
> [...]
>> In [6]: m=randn(200,200)
>>
>> In [7]: %time U,s,V=svd(m)
>> --------------------------------------------------------------------------- 
>>
>> LinAlgError                               Traceback (most recent call 
>> last)
> [...]
>>     1319     if results['info']>  0:
>> ->  1320         raise LinAlgError, 'SVD did not converge'
>>     1321     s = s.astype(_realType(result_t))
>>     1322     if compute_uv:
>>
>> LinAlgError: SVD did not converge
>
> Что же у Вас такое может быть, что даже на 200x200 вываливается? Увы, 
> я тут бессилен, ибо не могу воспроизвести ситуацию. Может быть, на 
> scipy.org спросите на этот счёт? Или, для начала, всё же дайте вывод gdb.
>
Как это сделать? Я поставил debuginfo пакеты.


Подробная информация о списке рассылки Devel