Re: [Help-gsl] FFTs and views

[Brian Gough]

At Tue, 19 Dec 2006 10:50:47 -0600,
lilly canelkatz wrote:
> I get views of the rows of this matrix
>   thisrowz=gsl_matrix_complex_row(thiszmatrix , thisrownumber);
> and would like to use them as input to 1-D FFTs:
>     gsl_fft_complex_radix2_forward(thisrowz, 1, NOCOLS)
>  I can't however because
> gsl_fft_complex_radix2_forward  takes a gsl_complex_packed_array as its
> first argument rather than
>  a complex vector view. Is there an elegant way to handle this please?

Hello,
You can call it like this:

gsl_fft_complex_radix2_forward(thisrowz.vector.data, 
                               thisrowz.vector.stride, 
                               thisrowz.vector.size)

gsl_complex_packed_array is equivalent to double*

-- 
Brian Gough
(GSL Maintainer)

Network Theory Ltd,
Publishing the GSL Manual - http://www.network-theory.co.uk/gsl/manual/