Re: [Help-gsl] gsl_complex to fftw_complex type

[Chia-Hung Lin]

2007/8/15, Oliver Jennrich <address@hidden>:
> On 8/15/07, Eric Germaneau <address@hidden> wrote:
> > Hello,
> >
> > I wondering whether its possible to convert gsl_complex to fftw_complex
> > type.
>
> Well, gsl_complex is defined as
>
> typedef struct
>   {
>     double dat[2];
>   }
> gsl_complex;
>
> whereas fftw_complex is either directly defined as _complex double or
> as double[2].
>
> In any event, fftw_complex should be bit compatible with gsl_complex,
> so a simple typecasting should do it.
>

So is it safe to use the following method ?

gsl_complex *a;
fftw_complex *b;
a = new gsl_complex[10];
b = reinterpret_cast<fftw_complex*>(a);

I tried this cast and it seemed to work, but I still want to be sure
that it's safe to do so.

Any comment would be very appreciated.

-- 
Chia-Hung Lin
Institute of Photonics and Optoelectronics, National Taiwan University