Path: csiph.com!x330-a1.tempe.blueboxinc.net!usenet.pasdenom.info!gegeweb.org!de-l.enfer-du-nord.net!feeder1.enfer-du-nord.net!tudelft.nl!txtfeed1.tudelft.nl!multikabel.net!newsfeed20.multikabel.net!news2.euro.net!newsgate.cistron.nl!newsgate.news.xs4all.nl!post.news.xs4all.nl!not-for-mail Return-Path: X-Original-To: python-list@python.org Delivered-To: python-list@mail.python.org X-Spam-Status: OK 0.008 X-Spam-Evidence: '*H*': 0.98; '*S*': 0.00; 'numerical': 0.05; 'linear': 0.07; 'squares': 0.07; 'python': 0.08; 'received:80.91': 0.09; 'received:80.91.229': 0.09; 'received:gmane.org': 0.09; 'received:list': 0.09; 'underlying': 0.09; 'am,': 0.12; 'enigma': 0.16; 'kern': 0.16; 'numpy': 0.16; 'overwhelmed': 0.16; 'wrote:': 0.18; 'header:In-Reply-To:1': 0.22; 'ken': 0.23; 'though.': 0.23; 'url:doc': 0.25; 'performing': 0.26; 'bit': 0.28; 'interpret': 0.29; 'seem': 0.29; 'least': 0.30; 'actually': 0.31; "i've": 0.32; 'requests': 0.32; "can't": 0.33; 'header:User-Agent:1': 0.33; 'header:X-Complaints-To:1': 0.34; 'to:addr:python-list': 0.35; 'received:org': 0.36; 'but': 0.37; 'using': 0.37; 'url:org': 0.39; 'except': 0.39; 'user': 0.40; 'to:addr:python.org': 0.40; 'your': 0.61; 'world': 0.61; 'received:86': 0.63; 'our': 0.63; 'believe': 0.65; 'brand': 0.76; 'algebra': 0.84; 'roberts': 0.84; 'url:reference': 0.84; 'eco': 0.91 X-Injected-Via-Gmane: http://gmane.org/ To: python-list@python.org From: Robert Kern Subject: Re: Numerical Linear Algebra in arbitrary precision Date: Fri, 17 Feb 2012 10:26:56 +0000 References: Mime-Version: 1.0 Content-Type: text/plain; charset=ISO-8859-1; format=flowed Content-Transfer-Encoding: 7bit X-Gmane-NNTP-Posting-Host: host86-160-161-180.range86-160.btcentralplus.com User-Agent: Mozilla/5.0 (Macintosh; Intel Mac OS X 10.7; rv:10.0.1) Gecko/20120208 Thunderbird/10.0.1 In-Reply-To: X-BeenThere: python-list@python.org X-Mailman-Version: 2.1.12 Precedence: list List-Id: General discussion list for the Python programming language List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , Newsgroups: comp.lang.python Message-ID: Lines: 27 NNTP-Posting-Host: 2001:888:2000:d::a6 X-Trace: 1329474426 news.xs4all.nl 6947 [2001:888:2000:d::a6]:53634 X-Complaints-To: abuse@xs4all.nl Xref: x330-a1.tempe.blueboxinc.net comp.lang.python:20547 On 2/17/12 6:09 AM, Tim Roberts wrote: > Ken wrote: >> >> Brand new Python user and a bit overwhelmed with the variety of >> packages available. Any recommendation for performing numerical >> linear algebra (specifically least squares and generalized least >> squares using QR or SVD) in arbitrary precision? I've been looking at >> mpmath but can't seem to find much info on built in functions except >> for LU decomposition/solve. > > It is been my experience that numpy is the best place to start with > requests like this, although I don't know whether it will actually solve > your specific tasks: > > http://docs.scipy.org/doc/numpy/reference/routines.linalg.html This will not do arbitrary-precision, though. We use the double- and single-precision routines from LAPACK. -- Robert Kern "I have come to believe that the whole world is an enigma, a harmless enigma that is made terrible by our own mad attempt to interpret it as though it had an underlying truth." -- Umberto Eco