FFTW 软件名称: the Faster Fourier Transform in the West
FFTW 软件简介:FFTW是一个快速计算离散傅里叶变换的标准C语言程序集,可计算一维或多维实和复数据以及任意规模的DFT。
FFTW 所需环境:Intel Compiler / GNU Compiler + MPI Compiler (如果编译mpi版的fftw)
FFTW 下载网址:http://www.fftw.org/download.html
安装步骤:
1、解压下载的安装包:tar -zxvf fftw-3.3.10.tar.gz
为了更好的发挥FFTW的功能,采用分别优化编译安装单精度版本和双精度版本的策略。
2、单精度版本的优化编译安装:(假设预期的安装路径为/home/seismic/Software/fftw3)
cd fftw-3.3.10
./configure --prefix=/home/seismic/Software/fftw3 CC=gcc F77=gfortran --enable-shared --enable-static --enable-float --enable-sse --enable-sse2 --enable-avx --enable-avx2 --enable-fma --enable-mpi
make
make install
3、双精度版本的优化编译安装:(两次安装过程类似,且相互独立)
./configure --prefix=/home/seismic/Software/fftw3 CC=gcc F77=gfortran --enable-shared --enable-static --enable-sse2 --enable-avx --enable-avx2 --enable-fma --enable-mpi
make
make install
安装参考:fftw编译安装教程-超算小站
4、安装后的代码测试
新建一个C++文本,写入附录代码,编译运行,正确输出结果即安装成功。
touch testcode.cpp
vi testcode.cpp
g++ testcode.cpp -o testcode -lm -lstdc++ -lfftw3
./testcode
(编译过程中,-lstdc++ 加上 g++ 的部分库,-lm 加上数学库 ,-lfftw3 加上fftw的库)
附录代码如下:
#include<fftw3.h>
#include<stdio.h>
#define N 8
int main()
{
int i;
fftw_complex *in,*out;
fftw_plan p;
in = (fftw_complex*) fftw_malloc(sizeof(fftw_complex) * N);
out = (fftw_complex*) fftw_malloc(sizeof(fftw_complex) * N);
if((in==NULL)||(out==NULL))
{
printf("Error:insufficient available memory\n");
}
else
{
for(i=0; i<N; i++)
{
in[i][0] = i+1;
in[i][1] = 0;
}
}
p = fftw_plan_dft_1d(N, in, out, FFTW_FORWARD,FFTW_ESTIMATE);
fftw_execute(p); /* repeat as needed */
fftw_destroy_plan(p);
fftw_cleanup();
for(i=0;i<N;i++)/*OUTPUT*/
{
printf("%f,%fi\n",in[i][0],in[i][1]);
}
printf("\n");
for(i=0;i<N;i++)/*OUTPUT*/
{
printf("%f,%fi\n",out[i][0],out[i][1]);
}
if(in!=NULL) fftw_free(in);
if(out!=NULL) fftw_free(out);
return 0;
}