PARALLEL COMPUTATION OF ENTRIES OF A-1

被引:14
作者
Amestoy, Patrick R. [1 ]
Duff, Iain S. [2 ,3 ]
L'Excellent, Jean-Yves [4 ,5 ]
Rouet, Francois-Henry [1 ,6 ]
机构
[1] Univ Toulouse, INPT ENSEEIHT IRIT, F-31071 Toulouse, France
[2] CERFACS, F-31057 Toulouse, France
[3] Rutherford Appleton Lab, Didcot OX11 0QX, Oxon, England
[4] Univ Lyon, Inria, F-69364 Lyon 07, France
[5] Univ Lyon, Lab LIP UMR CNRS ENS Lyon Inria 5668, F-69364 Lyon 07, France
[6] Univ Calif Berkeley, Lawrence Berkeley Natl Lab, Berkeley, CA 94720 USA
基金
英国工程与自然科学研究理事会;
关键词
sparse matrices; direct methods for linear system and matrix inversion; parallel algorithms; ALGORITHM; SYSTEMS; MATRIX;
D O I
10.1137/120902616
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
In this paper, we consider the computation in parallel of several entries of the inverse of a large sparse matrix. We assume that the matrix has already been factorized by a direct method and that the factors are distributed. Entries are efficiently computed by exploiting sparsity of the right-hand sides and the solution vectors in the triangular solution phase. We demonstrate that in this setting, parallelism and computational efficiency are two contrasting objectives. We develop an efficient approach and show its efficiency on a general purpose parallel multifrontal solver.
引用
收藏
页码:C268 / C284
页数:17
相关论文
共 17 条
[1]   ON COMPUTING INVERSE ENTRIES OF A SPARSE MATRIX IN AN OUT-OF-CORE ENVIRONMENT [J].
Amestoy, Patrick R. ;
Duff, Iain S. ;
L'Excellent, Jean-Yves ;
Robert, Yves ;
Rouet, Francois-Henry ;
Ucar, Bora .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2012, 34 (04) :A1975-A1999
[2]   Hybrid scheduling for the parallel solution of linear systems [J].
Amestoy, PR ;
Guermouche, A ;
L'Excellent, JY ;
Pralet, S .
PARALLEL COMPUTING, 2006, 32 (02) :136-156
[3]   A fully asynchronous multifrontal solver using distributed dynamic scheduling [J].
Amestoy, PR ;
Duff, IS ;
L'Excellent, JY ;
Koster, J .
SIAM JOURNAL ON MATRIX ANALYSIS AND APPLICATIONS, 2001, 23 (01) :15-41
[4]  
Bjorck A, 1996, NUMERICAL METHODS L
[5]   Integral SPI observation of the galactic central radian:: Contribution of discrete sources and implication for the diffuse emission [J].
Bouchet, L ;
Roques, JP ;
Mandrou, P ;
Strong, A ;
Diehl, R ;
Lebrun, F ;
Terrier, R .
ASTROPHYSICAL JOURNAL, 2005, 635 (02) :1103-1115
[6]  
Campbell Y., 1995, TR95021 U FLOR COMP
[7]   A scalable distributed method for quantum-scale device simulation [J].
Cauley, Stephen ;
Jain, Jitesh ;
Koh, Cheng-Kok ;
Balakrishnan, Venkataramanan .
JOURNAL OF APPLIED PHYSICS, 2007, 101 (12)
[8]   TASK-SCHEDULING FOR PARALLEL SPARSE CHOLESKY FACTORIZATION [J].
GEIST, GA ;
NG, E .
INTERNATIONAL JOURNAL OF PARALLEL PROGRAMMING, 1989, 18 (04) :291-314
[9]   A FAST PARALLEL ALGORITHM FOR SELECTED INVERSION OF STRUCTURED SPARSE MATRICES WITH APPLICATION TO 2D ELECTRONIC STRUCTURE CALCULATIONS [J].
Lin, Lin ;
Yang, Chao ;
Lu, Jianfeng ;
Ying, Lexing ;
E, Weinan .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2011, 33 (03) :1329-1351
[10]  
Lin L, 2009, COMMUN MATH SCI, V7, P755