From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 7055 invoked by alias); 1 Nov 2007 21:16:24 -0000 Received: (qmail 7045 invoked by uid 22791); 1 Nov 2007 21:16:24 -0000 X-Spam-Check-By: sourceware.org Received: from ch-smtp01.sth.basefarm.net (HELO ch-smtp01.sth.basefarm.net) (80.76.149.212) by sourceware.org (qpsmtpd/0.31) with ESMTP; Thu, 01 Nov 2007 21:16:19 +0000 Received: from c213-89-160-115.bredband.comhem.se ([213.89.160.115]:33159 helo=[192.168.0.2]) by ch-smtp01.sth.basefarm.net with esmtp (Exim 4.68) (envelope-from ) id 1InhOi-0001c1-5H; Thu, 01 Nov 2007 22:16:16 +0100 In-Reply-To: <383921.47800.qm@web56606.mail.re3.yahoo.com> References: <383921.47800.qm@web56606.mail.re3.yahoo.com> Mime-Version: 1.0 (Apple Message framework v752.2) Content-Type: text/plain; charset=US-ASCII; delsp=yes; format=flowed Message-Id: <52AA1BF5-4A7B-4574-A0BF-28F247BEECF0@comhem.se> Cc: gsl-discuss@sources.redhat.com Content-Transfer-Encoding: 7bit From: Tommy Nordgren Subject: Re: Toeplitz solver Date: Thu, 01 Nov 2007 21:16:00 -0000 To: Z F X-Mailer: Apple Mail (2.752.2) X-Scan-Result: No virus found in message 1InhOi-0001c1-5H. Mailing-List: contact gsl-discuss-help@sourceware.org; run by ezmlm Precedence: bulk List-Id: List-Subscribe: List-Archive: List-Post: List-Help: , Sender: gsl-discuss-owner@sourceware.org X-SW-Source: 2007-q4/txt/msg00013.txt.bz2 On 1 nov 2007, at 06.06, Z F wrote: > Hello everybody, > > Could someone suggest a (comprehendible) reference on Toeplitz > equations > solver for a general Toeplitz matrix. I am interested in fast methods > :) > > Numerical recepies has a mistake in the ALGORITHM, and in the code, of > course. (The third edition has this old bug also). It addresses vector > indecies which are out of range. (look at eq 2.8.23 and 2.824, terms > with R_{m+1} and R_{-m-1} are out of range of the vector size) > > And as to myself, I can not seem to trace were the problem happens. > > Thanks a lot, > > ZF > Numerical Recipes adds an offset to pointers returned from allocation routines, so that the useful lower index can be other than zero. ------------------------------------- This sig is dedicated to the advancement of Nuclear Power Tommy Nordgren tommy.nordgren@comhem.se