天然地震走时反演矩阵顺序三角化算法的优化和并行化

被引:3
作者
刘劲松 [1 ]
刘福田 [1 ]
刘俊 [2 ]
郝天珧 [1 ]
机构
[1] 中国科学院地质与地球物理研究所
[2] 中国科学院网络中心高性能计算中心
关键词
地震层析成像; 走时反演; QR分解; Givens变换; 顺序正交三角化; 并行计算;
D O I
暂无
中图分类号
P631.4 [地震勘探];
学科分类号
摘要
本文讨论了用于天然地震层析成像的矩阵顺序三角化算法,揭示了偏导数矩阵中非零元素的分布规律,在分析算法特点的基础上提出了对程序的几点优化措施,包括:采用合理的存储结构和计算次序以节省内存寻址时间;以一个一维数组存储上三角方阵以节省出一半的内存;对关键的Givens变换子程序进行优化和矢量化等,最终使程序的运行效率提高了数十倍.同时结合算法特点研究了矩阵顺序三角化的并行化方法,提出了两种并行化方案,在联想深腾6800超级计算机上实现了这两种并行算法;对计算量均衡等问题进行了讨论.
引用
收藏
页码:911 / 915
页数:5
相关论文
共 6 条
[1]   波动方程叠前深度偏移及并行计算在南黄海地区的应用 [J].
朱振宇 ;
刘洪 ;
李幼铭 .
地球物理学进展, 2003, (02) :302-305
[2]   高精度混合法叠前深度偏移及其并行实现 [J].
张文生 ;
张关泉 ;
宋海斌 .
地球物理学报, 2001, (04) :542-551
[3]   细胞自动机地震波模拟的并行化算法 [J].
王真理 ;
李幼铭 .
地球物理学报, 1999, (03) :410-415
[4]   用于速度图象重建的层析成象法 [J].
刘福田 ;
李强 ;
吴华 ;
胡戈 ;
刘建华 .
地球物理学报, 1989, (01) :46-61
[5]  
矩阵计算[M]. 科学出版社 , (美)G.H.戈卢布(GeneH.Golub), 2001
[6]  
Sparse QR factorization on a massively parallel computer[J] . Steven G. Kratzer.The Journal of Supercomputing . 1992 (3)