factor de escala para el CUFFT

Question

Escribí un programa de prueba simple, donde estaba haciendo complejos a complejos FT y acabo de generar algunos datos 1..50 y los metí en la parte real e imaginaria para cada índice de la matriz.

Cuando hago una operación como ésta IFFT (FFT (A)) = A

para probarlos, estoy consiguiendo resultados diferentes para cada biblioteca.

FFTW, tengo que dividir la salida por Len (A) para volver al original Un

Sin embargo, hacer esto adelante y luego revertir FFT con CUFFT, parece que tengo que dividir por (sqrt (2) * 50) para volver a los datos originales.

¿De dónde viene este factor de raíz cuadrada extra?

de acuerdo con la documentación de CUFFT: CUFFT realiza FFT no normalizadas; es decir, realizar una FFT hacia delante en un conjunto de datos de entrada seguido de una FFT inversa en el conjunto resultante arroja datos que son iguales a la entrada escalada por el número de elementos . El tamaño de la transformación por el recíproco del tamaño del conjunto de datos se deja para que el usuario lo realice como lo considere oportuno.

Gracias de antemano

Answer 1

Esto terminó siendo un problema con la forma en que se calculaban los valores absolutos del número complejo. en la biblioteca std :: complex, estaba calculando la distancia del vector.

Answer 2

CUFFT tiene el mismo comportamiento que FFTW, se calcula la FFT no normalizadas. IFFT(FFT(A))=n A donde n es la longitud del vector. La longitud n está en número de muestras (no flotantes ni bytes). Hay algunas diferencias de relleno entre FFTW y CUFFT con C2R y R2C que pueden arruinar una comparación simple, pero no para C2C. Verificaría dos veces la configuración de datos y los cálculos de longitud, y verificaría su plan tanto en FFTW como en CUFFT.