MT Occam并行反演方案及性能分析
|
摘要
针对MT Occam反演速度慢的问题,分析了实现其并行计算的可行性和途径,提出了基于频点和拉格朗日乘子(μ-1)扫描点划分的大粒度并行反演方案和具体实现步骤。对方案的加速性能进行了理论分析,得出了理论加速比曲线,明确了方案的适用性。采用基于消息传递的PVM并行平台,以主-从编程模式,实现了偏导数矩阵和μ值求取的并行计算。算例表明,综合并行效率较高,反演速度提高显著。
This paper analyzed the feasibility and ways for parallel implementation of the MT(magneto telluric) Occam inversion in order to increase the inversion speed.A coarse-grained parallel inversion scheme based on frequency point and Lagrange multiplier(μ-1) scanning point division was proposed and a concrete parallel procedure was put forward.Theoretical performance of the scheme was analyzed and the speed-up curve obtained,and the applicability of the scheme was clarified.Parallel computations for the partial derivative matrix and the μ seeking were realized by using the master-slave programming model under the message-passing PVM parallel platform.The experiment result indicated that the all-around parallel efficiency was ideal and the inversion speed increased obviously.
This paper analyzed the feasibility and ways for parallel implementation of the MT(magneto telluric) Occam inversion in order to increase the inversion speed.A coarse-grained parallel inversion scheme based on frequency point and Lagrange multiplier(μ-1) scanning point division was proposed and a concrete parallel procedure was put forward.Theoretical performance of the scheme was analyzed and the speed-up curve obtained,and the applicability of the scheme was clarified.Parallel computations for the partial derivative matrix and the μ seeking were realized by using the master-slave programming model under the message-passing PVM parallel platform.The experiment result indicated that the all-around parallel efficiency was ideal and the inversion speed increased obviously.
引文
[1]吴小平.大地电磁数据的Occam反演改进[J].地球物理学报,1998,41(4):547-553.
[2]陈小斌.大地电磁正反演新算法研究及资料处理与解释的可视化集成系统开发[D].北京:中国地震局地质研究所,2003.
[3]Constable S C,Parker R L,Constable C G.Occam’s Inversion:A Practical Algorithm for Generating Smooth Models from EMSounding Data[J].Geophysics,1987,52(1):289-300.
[4]Harry F,Jordan Gita Alaghband.并行处理基本原理[M].迟利华,刘杰,译.北京:清华大学出版社,2004.
[5]陈国良.并行计算——结构.算法.编程[M].北京:高等教育出版社,1999.
[6]饶文碧,张莉,李卓球,等.网络分布式环境下的并行计算方法[J].武汉工业大学学报,2000,22(1):54-56.
[2]陈小斌.大地电磁正反演新算法研究及资料处理与解释的可视化集成系统开发[D].北京:中国地震局地质研究所,2003.
[3]Constable S C,Parker R L,Constable C G.Occam’s Inversion:A Practical Algorithm for Generating Smooth Models from EMSounding Data[J].Geophysics,1987,52(1):289-300.
[4]Harry F,Jordan Gita Alaghband.并行处理基本原理[M].迟利华,刘杰,译.北京:清华大学出版社,2004.
[5]陈国良.并行计算——结构.算法.编程[M].北京:高等教育出版社,1999.
[6]饶文碧,张莉,李卓球,等.网络分布式环境下的并行计算方法[J].武汉工业大学学报,2000,22(1):54-56.
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心