FFTW¶
简介¶
FFTW是由麻省理工学院计算机科学实验室超级计算技术组开发的一套离散傅立叶变换(DFT)的计算库,开源、高效和标准C语言编写的代码使其得到了非常广泛的应用。Intel的数学库和Scilib(类似于Matlab的科学计算软件)都使用FFTW做FFT计算。FFTW是计算离散Fourier变换(DFT)的快速C程序的一个完整集合。
FFTW使用说明¶
思源一号上的FFTW¶
创建fftwtest目录并进入该目录:
mkdir fftwtest
cd fftwtest
在该目录下编写如下myfftw.c文件:
#include <stdio.h>
#include "fftw3.h"
int main()
{
fftw_complex *in, *out;
fftw_plan p;
int N= 8;
int i;
int j;
in = (fftw_complex*) fftw_malloc(sizeof(fftw_complex) * N);
out = (fftw_complex*) fftw_malloc(sizeof(fftw_complex) * N);
for( i=0; i < N; i++)
{
in[i][0] = 1.0;
in[i][1] = 0.0;
printf("%6.2f ",in[i][0]);
}
printf("\n");
p=fftw_plan_dft_1d(N,in,out, FFTW_FORWARD, FFTW_ESTIMATE);
fftw_execute(p); /* repeat as needed*/
for(j = 0;j < N;j++)
{
printf("%6.2f ",out[j][0]);
}
printf("\n");
fftw_destroy_plan(p);
fftw_free(in);
fftw_free(out);
return 0;
}
在该目录下编写如下fftwtest.slurm脚本:
#!/bin/bash
#SBATCH --job-name=fftwtest
#SBATCH --partition=64c512g
#SBATCH --ntasks-per-node=1
#SBATCH -n 1
#SBATCH --output=%j.out
#SBATCH --error=%j.err
ulimit -s unlimited
ulimit -l unlimited
module load gcc
module load fftw/3.3.10-gcc-8.3.1-openmpi
gcc myfftw.c -o myfftw -lfftw3 -lm
./myfftw
使用如下命令提交作业:
sbatch fftwtest.slurm
作业完成后在.out文件中可看到如下结果:
1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00
8.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
pi2.0上的FFTW¶
此步骤和上文完全相同;
此步骤和上文完全相同;
在该目录下编写如下fftwtest.slurm脚本:
#!/bin/bash
#SBATCH --job-name=fftwtest
#SBATCH --partition=small
#SBATCH --ntasks-per-node=1
#SBATCH -n 1
#SBATCH --output=%j.out
#SBATCH --error=%j.err
ulimit -s unlimited
ulimit -l unlimited
module load gcc
module load fftw/3.3.8-gcc-4.8.5
gcc myfftw.c -o myfftw -lfftw3 -lm
./myfftw
使用如下命令提交作业:
sbatch fftwtest.slurm
作业完成后在.out文件中可看到如下结果:
1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00
8.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
参考资料¶
最后更新:
2024 年 11 月 14 日