Bug 24088 - numpy.linalg.svd fails
Summary: numpy.linalg.svd fails
Status: CLOSED WORKSFORME
Alias: None
Product: Sisyphus
Classification: Development
Component: python-module-numpy (show other bugs)
Version: unstable
Hardware: x86 Linux
: P3 normal
Assignee: Grigory Ustinov
QA Contact: qa-sisyphus
URL:
Keywords:
Depends on:
Blocks:
 
Reported: 2010-09-15 23:26 MSD by Valery Pipin
Modified: 2010-09-20 08:44 MSD (History)
20 users (show)

See Also:


Attachments

Note You need to log in before you can comment on or make changes to this bug.
Description Valery Pipin 2010-09-15 23:26:46 MSD
Делаем следующий тест.

from numpy.random import *
from numpy.linalg import *
m=randn(1000,1000)

%time U,s,V=svd(m)
## it should work

m=randn(2000,2000)
%time U,s,V=svd(m)

Получаем:

ERROR: An unexpected error occurred while tokenizing input
The following traceback may be corrupted or invalid
The error message is: ('EOF in multi-line statement', (236, 0))

---------------------------------------------------------------------------
MemoryError                               Traceback (most recent call last)

/home/va/<ipython console> in <module>()

/usr/lib/python2.6/site-packages/IPython/iplib.pyc in ipmagic(self, arg_s)
   1180         else:
   1181             magic_args = self.var_expand(magic_args,1)
-> 1182             return fn(magic_args)
   1183 
   1184     def ipalias(self,arg_s):

/usr/lib/python2.6/site-packages/IPython/Magic.pyc in magic_time(self, parameter_s)                                                                             
   1969         else:
   1970             st = clk()
-> 1971             exec code in glob
   1972             end = clk()
   1973             out = None

/home/va/<timed exec> in <module>()

/usr/lib/python2.6/site-packages/numpy/linalg/linalg.pyc in svd(a, full_matrices, compute_uv)                                                                   
   1314                                  work, -1, iwork, 0)
   1315         lwork = int(work[0])
-> 1316         work = zeros((lwork,), t)
   1317         results = lapack_routine(option, m, n, a, m, s, u, m, vt, nvt,
   1318                                  work, lwork, iwork, 0)

MemoryError: 

Модуль linalg из scipy с этой задачей справился.
Все это на 32бит сизифе из-под виртуалбокс (ему доступно 512Mb RAM) . 
python-module-numpy-2.0.0-alt1.svn20100607.1
python-module-scipy-0.9.0-alt1.svn20100607
Comment 1 Andrey Rahmatullin 2010-09-15 23:31:23 MSD
У меня работает на последнем 32бит сизифе.
Comment 2 Andrey Rahmatullin 2010-09-15 23:31:55 MSD
А, MemoryError.
Памяти 4 Гб.
Comment 3 Valery Pipin 2010-09-15 23:38:06 MSD
(В ответ на комментарий №2)
> А, MemoryError.
> Памяти 4 Гб.

Да если памяти добавить исчезает. Но ведь модуль из scipy справился и на 512-ти
Comment 4 real@altlinux.org 2010-09-16 06:18:14 MSD
Мне негде проверить на такой конфигурации. Лично у меня всё работает. Так что лучше подождать подтверждения от других пользователей, прежде чем обращаться в апстрим. Или всё же прямо так им настучать? :)
Comment 5 Valery Pipin 2010-09-16 06:27:17 MSD
(In reply to comment #4)
> Мне негде проверить на такой конфигурации. Лично у меня всё работает. Так что
> лучше подождать подтверждения от других пользователей, прежде чем обращаться в
> апстрим. Или всё же прямо так им настучать? :)
Не надо. Можно вообще закрыть пока.
Но я вами восхищаюсь. Как на апстриме можно работать?
Comment 6 real@altlinux.org 2010-09-16 06:35:57 MSD
"Можно вообще закрыть пока."

Я всё-таки подожду, вдруг кто-нибудь подтвердит такое падение в схожих условиях...

"Как на апстриме можно работать?"

Не понял...
Comment 7 Valery Pipin 2010-09-16 06:52:57 MSD
(In reply to comment #6)
> "Можно вообще закрыть пока."
> 
> Я всё-таки подожду, вдруг кто-нибудь подтвердит такое падение в схожих
> условиях...
> 
> "Как на апстриме можно работать?"
> 
> Не понял...
Ну так 
Scipy и numpy очень активно развиваются вероятность частных поломок очень велика. Сизиф пожалуй единственный репозитарий где можно найти практически апстримную версию этих программ.
Comment 8 real@altlinux.org 2010-09-16 07:05:58 MSD
На то он и сизиф ;)
Comment 9 real@altlinux.org 2010-09-20 08:44:34 MSD
Так, ну раз пока никто ничего не подтвердил, закрою, если что, переоткрывайте.