The Highly Parallel Incomplete Gram-Schmidt Preconditioner

نویسندگان

  • Laurence T. Yang
  • Hai-Xiang Lin
چکیده

In this paper we study the parallel aspects of IMGS, Incomplete Modiied Gram-Schmidt preconditioner which can be used for ef-ciently solving sparse and large linear systems and least squares problems on massively parallel distributed memory computers. The performance of this preconditioning technique on this kind of architecture is always limited because of the global communication required for the inner products, even for ParIMGS, a parallel version of IMGS where we create some possibilities such that the computation can be overlapped with the communication. We will describe a more eecient alternative, namely Improved ParIMGS (IParIMGS) which avoids the global communication of inner products and only requires local communications. Therefore, the cost of communication can be signiicantly reduced. Several numerical experiments carried out on Parsytec GC/PowerPlus are presented as well.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Computing a block incomplete LU preconditioner as the by-product of block left-looking A-biconjugation process

In this paper, we present a block version of incomplete LU preconditioner which is computed as the by-product of block A-biconjugation process. The pivot entries of this block preconditioner are one by one or two by two blocks. The L and U factors of this block preconditioner are computed separately. The block pivot selection of this preconditioner is inherited from one of the block versions of...

متن کامل

Preconditioners for the conjugate gradient algorithm using Gram-Schmidt and least squares methods

This paper is devoted to the study of some preconditioners for the conjugate gradient algorithm used to solve large sparse linear and symmetric positive definite systems. The construction of a preconditioner based on the Gram–Schmidt orthogonalization process and the least squares method is presented. Some results on the condition number of the preconditioned system are provided. Finally, numer...

متن کامل

Solving Sparse Least Squares Problems on Massively Distributed Memory Computers

In this paper we study the parallel aspects of PCGLS, a basic iterative method whose main idea is to organize the computation of conjugate gradient method with preconditioner applied to normal equations, and Incomplete Modified Gram-Schmidt (IMGS) preconditioner for solving sparse least squares problems on massively parallel distributed memory computers. The performance of these methods on this...

متن کامل

New IMGS-based Preconditioners for Least Squares Problems

Convergence acceleration by preconditioning is usually essential when solving the standard least squares problems by an iterative method. IMGS, is an incomplete modiied version of Gram-Schmidt orthogonalization to obtain an incomplete orthogonal factorization pre-conditioner M = R, where A = Q R + E is an approximation of a QR factorization, Q is an orthogonal matrix and R is upper triangular m...

متن کامل

CIMGS: An Incomplete Orthogonal FactorizationPreconditioner

A new preconditioner for symmetric positive definite systems is proposed, analyzed, and tested. The preconditioner, compressed incomplete modified Gram–Schmidt (CIMGS), is based on an incomplete orthogonal factorization. CIMGS is robust both theoretically and empirically, existing (in exact arithmetic) for any full rank matrix. Numerically it is more robust than an incomplete Cholesky factoriza...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:

دوره   شماره 

صفحات  -

تاریخ انتشار 1997