This is the mail archive of the gsl-discuss@sourceware.org mailing list for the GSL project.


Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]
Other format: [Raw text]

Re: Toeplitz solver



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





Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]