[devel] atlas is broken again?

[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.

-- 

REAL aka Евгений Ростовцев, программист ЦНИТ КемГУ