[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