%0 Journal Article %T Parallel LSQR algorithms used in seismic tomography
地震层析成像LSQR算法的并行化 %A LIU Jin-Song %A LIU Fu-Tian %A LIU Jun %A HAO Tian-Yao %A
刘劲松 %A 刘福田 %A 刘俊 %A 郝天珧 %J 地球物理学报 %D 2006 %I %X We discuss the LSQR algorithms used in earthquake travel time tomography. We keep the epicenter terms in the equation for regional events, and then use the orthogonal projection method to eliminate the epicenter terms. For tele-events, the classic smoothing process is used. The number of non-zero elements in the partial derivative matrix is increased by several times because of the orthogonal projection and smoothing processes. For a large scale inversion problem, the amount of non-zero elements can be dozens of Gigabytes or hundreds of Gigabytes. The huge amount of memory requirement becomes the bottle neck of LSQR algorithms. matrix, designed an efficient data structure for the sparse matrix, used a distributed memory and computation scheme for matrix computation, and implemented it on a multi-processor super-computer. We have derived an estimation formula of parallel efficiency and tested two real tomography models. %K Seismic tomography %K Traveltime inversion %K LSQR algorithms %K Parallel computation
地震层析成像 %K 走时反演 %K LSQR算法 %K 并行计算 %U http://www.alljournals.cn/get_abstract_url.aspx?pcid=E62459D214FD64A3C8082E4ED1ABABED5711027BBBDDD35B&cid=1E44AE713D8A6DE0&jid=14DC41C59CBF6770055A7D610D53AE46&aid=868C69416A6D2347&yid=37904DC365DD7266&vid=2A3781E88AB1776F&iid=0B39A22176CE99FB&sid=DFBC046213B3DD86&eid=6313C162FF75889A&journal_id=0001-5733&journal_name=地球物理学报&referenced_num=2&reference_num=17