From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: X-Spam-Checker-Version: SpamAssassin 3.2.5 (2008-06-10) on sa.int.altlinux.org X-Spam-Level: * X-Spam-Status: No, score=1.1 required=5.0 tests=BAYES_50,DNS_FROM_OPENWHOIS autolearn=no version=3.2.5 Message-ID: <4D9AD0BB.4080109@mmedia2.kemsu.ru> Date: Tue, 05 Apr 2011 15:20:11 +0700 From: REAL User-Agent: Mozilla/5.0 (X11; U; Linux i686; en-US; rv:1.9.2.16pre) Gecko/20110308 Thunderbird/3.1.10pre MIME-Version: 1.0 To: ALT Linux Team development discussions References: <201104042044.11351.pip@iszf.irk.ru> <201104050034.12244.pip@iszf.irk.ru> <4D9AC7A7.1030306@mmedia2.kemsu.ru> <201104050102.18009.pip@iszf.irk.ru> In-Reply-To: <201104050102.18009.pip@iszf.irk.ru> Content-Type: text/plain; charset=KOI8-R; format=flowed Content-Transfer-Encoding: 8bit Subject: Re: [devel] atlas is broken again? X-BeenThere: devel@lists.altlinux.org X-Mailman-Version: 2.1.12 Precedence: list Reply-To: ALT Linux Team development discussions List-Id: ALT Linux Team development discussions List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Tue, 05 Apr 2011 08:15:27 -0000 Archived-At: List-Archive: List-Post: 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 Евгений Ростовцев, программист ЦНИТ КемГУ