Проблема с подпрограммой FFTW2 в C - PullRequest
Новая
       4

Проблема с подпрограммой FFTW2 в C

0 голосов
/ 11 мая 2011

В настоящее время я работаю с библиотекой FFTW2 (самое быстрое преобразование Фурье на Западе), и после написания успешной процедуры на Фортране я перехожу к C. Однако у меня возникают некоторые проблемы с назначением данных, когда япопробуйте ввести данные для преобразования (т. е. значения sin (x)).В настоящее время мой код: 

#include <stdio.h>
#include <fftw.h>
#include <math.h>
#include <complex.h>

void compute_fft(){
  int i;
  const int n[8];
  fftw_complex in[8];
  fftwnd_plan p;

  in[0]->re = 0;
  in[1]->re = (sqrt(2)/2);
  in[2]->re = 1;
  in[3]->re = (sqrt(2)/2);
  in[4]->re = 0;
  in[5]->re = -(sqrt(2)/2);
  in[6]->re = -1;
  in[7]->re = -(sqrt(2)/2);

  for(i = 0; i < 8; i++){
    (in[i])->im = 0;
  }

  p = fftwnd_create_plan(8, n, FFTW_FORWARD, FFTW_ESIMATE | FFTW_IN_PLACE);
  fftwnd_one(p, &in[0], NULL);
  fftwnd_destroy_plan(p);

  printf("Sin\n");
  for(i = 0; i < 8; i++){
    printf("%d\n", n[i]);
  }
}

Я использовал эту ссылку http://fftw.org/fftw2_doc/ для документации / учебных целей, и в настоящее время моя ошибка - «неверный аргумент типа a-> a (haveafftw_complexa) ", а над символами 'a' по обе стороны от -> и fftw_complex над ними стоит карат.Я использую в основном тот же шаблон, который работал, когда я писал это на Фортране, и я, похоже, следую учебному пособию, здесь на самом деле только это задание, которое я испортил синтаксис.Для записи, я использую компилятор nvcc, если это имеет значение (ранее я пробовал gcc).Если кто-то здесь уже использовал каждый FFTW2 в C, не могли бы вы помочь исправить мою ошибку?Спасибо!

Ответы [ 2 ]

0 голосов
/ 21 июля 2011

Поскольку ffwt_complex является двойным [2], являющимся первым измерением ([0]) для реальных данных и вторым ([1]) для мнимых данных, безопасное решение: 

  in[0][0] = 0;
  in[1][0] = (sqrt(2)/2);
  in[2][0] = 1;
  in[3][0] = (sqrt(2)/2);
  in[4][0] = 0;
  in[5][0] = -(sqrt(2)/2);
  in[6][0] = -1;
  in[7][0] = -(sqrt(2)/2);

 for(i = 0; i < 8; i++){
   in[i][1] = 0;
  }
0 голосов
/ 11 мая 2011

Это может быть связано с тем, что ваш массив "in" является массивом fftw_complex, поэтому вместо использования в [0] -> re = 0 вы должны использовать его как в [0] .re = 0. Если fftw_complex снова не будет typedefined для некоторого массива. 

fftw_complex in[8];
in[0].re = 0;
in[1].re = (sqrt(2)/2);
in[2].re = 1;
in[3].re = (sqrt(2)/2);
in[4].re = 0;
in[5].re = -(sqrt(2)/2);
in[6].re = -1;
in[7].re = -(sqrt(2)/2);
...