gsl-ocaml / lib / mlgsl_fft.c

Diff from to

lib/mlgsl_fft.c

 
 CAMLprim value ml_gsl_fft_complex_rad2_forward(value dif, value stride, value data)
 {
-  size_t N = Double_array_length(data);
+  size_t N = Double_array_length(data) / 2;
   size_t c_stride = Opt_arg(stride, Int_val, 1);
   int c_dif = Opt_arg(dif, Bool_val, 0);
 
 CAMLprim value ml_gsl_fft_complex_rad2_transform(value dif, value stride, 
 					value data, value sign)
 {
-  size_t N = Double_array_length(data);
+  size_t N = Double_array_length(data) / 2;
   size_t c_stride = Opt_arg(stride, Int_val, 1);
   int c_dif = Opt_arg(dif, Bool_val, 0);
   gsl_fft_direction c_sign =
 
 CAMLprim value ml_gsl_fft_complex_rad2_backward(value dif, value stride, value data)
 {
-  size_t N = Double_array_length(data);
+  size_t N = Double_array_length(data) / 2;
   size_t c_stride = Opt_arg(stride, Int_val, 1);
   int c_dif = Opt_arg(dif, Bool_val, 0);
 
 
 CAMLprim value ml_gsl_fft_complex_rad2_inverse(value dif, value stride, value data)
 {
-  size_t N = Double_array_length(data);
+  size_t N = Double_array_length(data) / 2;
   size_t c_stride = Opt_arg(stride, Int_val, 1);
   int c_dif = Opt_arg(dif, Bool_val, 0);
 
Tip: Filter by directory path e.g. /media app.js to search for public/media/app.js.
Tip: Use camelCasing e.g. ProjME to search for ProjectModifiedEvent.java.
Tip: Filter by extension type e.g. /repo .js to search for all .js files in the /repo directory.
Tip: Separate your search with spaces e.g. /ssh pom.xml to search for src/ssh/pom.xml.
Tip: Use ↑ and ↓ arrow keys to navigate and return to view the file.
Tip: You can also navigate files with Ctrl+j (next) and Ctrl+k (previous) and view the file with Ctrl+o.
Tip: You can also navigate files with Alt+j (next) and Alt+k (previous) and view the file with Alt+o.