Hi, The documantation for gsl_fft_real_radix2_transform says that the output can be sorted using gsl_fft_halfcomplex_unpack. The documantation of gsl_fft_halfcomplex_unpack says it can also be used to sort the output of gsl_fft_real_transform, but the output of gsl_fft_real_radix2_transform and gsl_fft_real_transform don't have the same storage. Actually I verified than gsl_fft_real_radix2_transform + gsl_fft_halfcomplex_unpack doesn't sort the data in the right way: the result is different from what you get converting the data to a complex array and transforming it with gsl_fft_complex_radix2_forward, but the result is correct if the output of gsl_fft_real_radix2_transform is sorted by hand.
Thank you. -- Cristian Toraci DIST - University of Genova Viale Francesco Causa, 13 16145 Genova, Italy tel. +39 010 3532789 fax. +39 010 3532948 [EMAIL PROTECTED] www.bio.dist.unige.it _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
