MPI程序设计PC11.pptVIP

现代密码学理论与实践之五 并 行 计 算 中国科学技术大学计算机科学与技术系 国家高性能计算中心(合肥) 2003年9月 第三篇 并行数值算法 第八章 基本通讯操作 第九章 稠密矩阵运算 第十章 线性方程组的求解 第十一章 快速傅里叶变换  第十一章 快速傅里叶变换 11.1 快速傅里叶变换 11.2 并行FFT算法    11.1 快速傅里叶变换(FFT) 11.1.1 离散傅里叶变换(DFT) 11.1.2 DFT的顺序代码 11.1.3 串行FFT递归算法 11.1.4 串行FFT非递归算法 离散傅里叶变换(DFT) 定义 给定向量A=(a0,a1,…,an-1)T,DFT将A变换为B=(b0,b1,…,bn-1)T 11.1 快速傅里叶变换(FFT) 11.1.1 离散傅里叶变换(DFT) 11.1.2 DFT的顺序代码 11.1.3 串行FFT递归算法 11.1.4 串行FFT非递归算法 DFT的顺序代码 代码1 for j=0 to n-1 do b[j]=0 for k=0 to n-1 do b[j]=b[j]+ωk*ja[k] end for end for 注：代码1需要计算ωk*j 代码2的复杂度为O(n2) 11.1 快速傅里叶变换(FFT) 11.1.1 离散傅里叶变换(DFT) 11.1.2 DFT的顺序代码 11.1.3 串行FFT递归算法 11.1.4 串行FFT非递归算法 串行FFT递归算法 蝶式递归计算原理 令 为n/2次单位元根，则有 . 将b向量的偶数项 和奇数项 分别记为 和 注意推导中反复使用 串行FFT递归算法 串行FFT递归算法 串行FFT递归算法 FFT的蝶式递归计算图(由计算原理推出) 串行FFT递归算法 特别地，n=8的FFT蝶式计算图(展开的) 串行FFT递归算法 SISD上的FFT分治递归算法 输入: a=(a0,a1,…,an-1); 输出: B=(b0,b1,…,bn-1) Procedure RFFT(a,b) begin if n=1 then b0=a0 else (1)RFFT(a0,a2,…,an-2, u0,u1,…,un/2-1) (2)RFFT(a1,a3,…,an-1, v0,v1,…,vn/2-1) (3)z=1 (4)for j=0 to n-1 do (4.1)bj=uj mod n/2+zvj mod n/2 (4.2)z=zω endfor 注: (1)算法时间复杂度t(n)=2t(n/2)+O(n) endif 解得 t(n)=O(nlogn) end (2)算法原理？ 11.1 快速傅里叶变换(FFT) 11.1.1 离散傅里叶变换(DFT) 11.1.2 DFT的顺序代码 11.1.3 串行FFT递归算法 11.1.4 串行FFT非递归算法 串行FFT非递归算法 蝶式计算示例 串行FFT非递归算法 蝶式计算流图 串行FFT非递归算法 行0 行1 行2 行3 行4 行5 行6 行7 如: b6=[(a0+a4)-(a2+a6)]-[(a1+a5)-(a3+a7)]ω2 注: ①下行线结点处的权因子的确定问题； ②bi的下标确定:取行号的位序反。如，行3: 3=(011)2 ==(110)2=6, == 行3的输出为b6 第十一章 快速傅里叶变换 1