From: Z F <mail4me9999@yahoo.com>
To: Tommy Nordgren <tommy.nordgren@comhem.se>
Cc: gsl-discuss@sources.redhat.com
Subject: Re: Toeplitz solver
Date: Fri, 02 Nov 2007 01:40:00 -0000 [thread overview]
Message-ID: <438054.31674.qm@web56615.mail.re3.yahoo.com> (raw)
In-Reply-To: <52AA1BF5-4A7B-4574-A0BF-28F247BEECF0@comhem.se>
--- Tommy Nordgren <tommy.nordgren@comhem.se> wrote:
>
> 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.
THanks for the comment, I knew about that.
I also know what the problem was. The formulae are correct, but the
description of the iteration process is wrong. I have changed it and
it is going.. (If you want to know, NR says first find G and H and then
find X, where is one has to first find X and then G and H. See the
section on that in NR)
The question remains of there is a reference for a clear description
of N log N method to solve Toeplitz equations.
Thanks..
ZF
__________________________________________________
Do You Yahoo!?
Tired of spam? Yahoo! Mail has the best spam protection around
http://mail.yahoo.com
prev parent reply other threads:[~2007-11-02 1:40 UTC|newest]
Thread overview: 3+ messages / expand[flat|nested] mbox.gz Atom feed top
2007-11-01 5:06 Z F
2007-11-01 21:16 ` Tommy Nordgren
2007-11-02 1:40 ` Z F [this message]
Reply instructions:
You may reply publicly to this message via plain-text email
using any one of the following methods:
* Save the following mbox file, import it into your mail client,
and reply-to-all from there: mbox
Avoid top-posting and favor interleaved quoting:
https://en.wikipedia.org/wiki/Posting_style#Interleaved_style
* Reply using the --to, --cc, and --in-reply-to
switches of git-send-email(1):
git send-email \
--in-reply-to=438054.31674.qm@web56615.mail.re3.yahoo.com \
--to=mail4me9999@yahoo.com \
--cc=gsl-discuss@sources.redhat.com \
--cc=tommy.nordgren@comhem.se \
/path/to/YOUR_REPLY
https://kernel.org/pub/software/scm/git/docs/git-send-email.html
* If your mail client supports setting the In-Reply-To header
via mailto: links, try the mailto: link
Be sure your reply has a Subject: header at the top and a blank line
before the message body.
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).