From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 8670 invoked by alias); 1 Nov 2007 05:06:59 -0000 Received: (qmail 8659 invoked by uid 22791); 1 Nov 2007 05:06:58 -0000 X-Spam-Check-By: sourceware.org Received: from web56606.mail.re3.yahoo.com (HELO web56606.mail.re3.yahoo.com) (66.196.97.50) by sourceware.org (qpsmtpd/0.31) with SMTP; Thu, 01 Nov 2007 05:06:56 +0000 Received: (qmail 49698 invoked by uid 60001); 1 Nov 2007 05:06:54 -0000 X-YMail-OSG: YbbAr.EVM1mewsGEGz.R_cnQ9Qc.xogGWOJ8xVBauF2BQSqS9BbKIjavhiOotJFWvKvZPrKDgUjFIXGeLIh6qK2rD4ZAZLZNVxbY.ZdQW46pndEsuGGcOZ8CESxv6jz8ynpt1M1yGTUbORHVuxpiBhVnzA-- Received: from [24.185.43.112] by web56606.mail.re3.yahoo.com via HTTP; Wed, 31 Oct 2007 22:06:54 PDT Date: Thu, 01 Nov 2007 05:06:00 -0000 From: Z F Subject: Toeplitz solver To: gsl-discuss@sources.redhat.com MIME-Version: 1.0 Content-Type: text/plain; charset=iso-8859-1 Content-Transfer-Encoding: 8bit Message-ID: <383921.47800.qm@web56606.mail.re3.yahoo.com> 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/msg00012.txt.bz2 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 __________________________________________________ Do You Yahoo!? Tired of spam? Yahoo! Mail has the best spam protection around http://mail.yahoo.com