按位交换的并行FFT.pptVIP

按位交换的并行FFT 超立方体上的按位交换并行FFT 二维网格上的按位交换并行FFT 并行FFT中的额外计算 超立方体上的按位交换并行FFT 设要在 p=2d个进程上并行计算n=2r个点的FFT，将这n个点平均分为p段，使得每段含有连续的n/p个点 对0~n-1中的任何一个点l，记其二进制形式为(b0b1…br-1)，则将CT_FFT中的Rl与Sl分配到进程(b0b1…bd-1)上 只有前d=logp次外迭代需要进行通信，而后r-d次外迭代无需任何通信 超立方体上的按位交换并行FFT(续) 超立方体上的按位交换并行FFT(续) 并行执行时间为 二维网格上的按位交换并行FFT 设n=2r，p=2d，r与d都是偶数，且处理器组织成物理上的2d/2×2d/2二维网格 对任意一个l＝0~n-1，如果其二进制形式为(b0b1…br-1)，则将其分布到第(b0b1…bd/2-1, bd/2bd/2+1…bd-1)个进程 一次消息传递可能需要经过网格中多条物理链路；同一时候不同进程对之间的消息传递可能需要经过同一物理链路 二维网格上的按位交换并行FFT(续) 二维网格上的按位交换并行FFT(续) 第m logp/2次外迭代，通信只出现在同行进程间，通信所需要的时间为s+b(n/p)2logp/2-1-m 之后的log p / 2次外迭代的总通信时间与前log p / 2次外迭代的总通信时间相同 并行执行时间为 Tp = c n log n / p + s log p + 2 b (n / p) (p1/2 – 1) 对应于通信起步时间s的等效率函数为 W=?(p log p) 对应于数据传输b的等效率函数为 并行FFT中的额外计算 在串行FFT中，旋转因子一共只有n个，其计算复杂性为O(n)，对总计算复杂性影响不大 由于计算一个旋转因子的时间远小于传输一个旋转因子的时间，所以对有的旋转因子，可能需要在不同的进程上分别进行计算，引入了额外计算 8点FFT的例子 并行FFT中的额外计算(续) 设n=2r，p=2d，对任何一个点(b0b1…br-1)，在数据分布时，将分布到进程(b0b1…bd-1)上 从CT_FFT知，在第m次外迭代，点(b0b1…br-1)的值更新只需要用到幂次为(bmbm-1…b00…0)的旋转因子 当m=0,…,d-1时，第m次外迭代上每个进程只用到一个旋转因子 对m=d,…,r-1时，第m次外迭代上每个进程需要用到2m-d+1个旋转因子 并行FFT中的额外计算(续) 每个进程最多需要用到一共logp + 21 +…+2r-d即2(n/p-1)+logp个旋转因子 假设每个旋转因子的计算时间为c’，则在超立方体上的并行执行时间应该修改为 * 国防科学技术大学并行与分布处理重点实验室 * 国防科学技术大学并行与分布处理重点实验室 当p固定不变，而n不断增大时，并行效率必然趋于1 等效率函数，记k＝Ep/(1- Ep)，则 当kb/c小于1时为W=?(plogp) 当kb/c大于1时为W=?(pkb/clogp) 24 25 26 27 16 17 18 19 8 9 10 11 0 1 2 3 4 5 6 7 12 13 14 15 20 21 22 23 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 等效率函数在量阶上并没有变化