Toeplitz solver

Tommy Nordgren tommy.nordgren@comhem.se
Thu Nov 1 21:16:00 GMT 2007


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





More information about the Gsl-discuss mailing list