Re: [Bug-gsl] Bug in gsl_matrix_complex_view_array and gsl_vector_complex_view_array

[Brian Gough]

At Fri, 11 Sep 2009 17:09:30 +0200,
Igor Gotlibovych wrote:
> I have discovered a bug with gsl_matrix_complex_view_array, seen
> from the definition
> 
> gsl_matrix_complex_view_array (double * base,
>                              const size_t n1, const size_t n2);
>  
> I expected base to be of type gsl_complex*, to be consistent with
> corresponding functions for other types. The same problem exists for
> gsl_vector_complex_view_array.

Hello,

For arrays of complex types we recommend using double *, or its
equivalent gsl_complex_packed_array defined in gsl_complex.h.

We don't use arrays of type gsl_complex * as these may differ in size,
due to the possibility of padding at the end of the gsl_complex
struct.

-- 
Brian Gough

GNU Scientific Library -
http://www.gnu.org/software/gsl/