Re: [Help-gsl] Complex vectors

[Brian Gough]

Cory Burkhardt writes:
 > I am attempting to perform FFT operations on complex vectors.  In the FFT
 > section of the documentation, it makes a comment about using a
 > gsl_complex_vector type, but this type is not described anywhere in the
 > documentation. 

Hello,

Thanks for your message. That's a typo for gsl_vector_complex --
thanks for bringing it to my attention.

 > In the section on Vectors, the gsl_vector's stride member is
 > described, but nowhere does it describe how to modify the size and stride
 > members.  If I want to create a complex vector with 5 complex numbers, would
 > I create a gsl_vector of size 10 and manually modify the stride to be 2?

You could make a vector like that with 

   #include <gsl/gsl_vector_complex.h>
   ...
   gsl_vector_complex * v =  gsl_vector_complex_alloc (5);

However, to use the FFTs you don't need a gsl_vector_complex type, you
can just malloc an array of 2*N doubles and pass it directly -- see
the examples in the FFT chapter for details.  

-- 
Brian Gough

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