Parallel Triangular Matrix Inversion With Increased Parallelism And Less Use of Shared Memory
نویسنده
چکیده
A triangular matrix is a special kind of square matrix. A square matrix is called lower triangular if all the entries above the main diagonal are zero. Conversely a square matrix is called upper triangular if all the entries below the main diagonal are zero. [6] For a square matrix M , M−1 is the inverse matrix where M ×M−1 = I and I denotes the n× n identity matrix. We say that the problem size of computing the inversion of a triangular Matrix M with dimension n, is n. Matrix inversion has innately some sequential steps, and the size of the problem at each step is remarkably large comaring to the size of the original problem. Many approaches have been proposed to either reduce the size of the sequential parts, increase concurrency [2] or use the memory efficiently[1]. In this project, a parallel algorithm for computing inversion of a triangular is proposed, and implemented in OpenMP. This algorithm relies on extended concurrency and efficient processor communication.
منابع مشابه
High-performance and Parallel Inversion of a Symmetric Positive Definite Matrix
We present families of algorithms for operations related to the computation of the inverse of a Symmetric Positive Definite (SPD) matrix: Cholesky factorization, inversion of a triangular matrix, multiplication of a triangular matrix by its transpose, and one-sweep inversion of an SPD matrix. These algorithms are systematically derived and implemented via the Formal Linear Algebra Methodology E...
متن کاملParallel Communication-Free Algorithm for Triangular Matrix Inversion on Heterogeneous Platforms
We address in this paper the parallelization of a recursive algorithm for triangular matrix inversion (TMI) based on the ‘Divide and Conquer’ (D&C) paradigm. A series of different versions of an original sequential algorithm are first presented. A theoretical performance study permits to establish an accurate comparison between the designed algorithms. Afterwards, we develop an optimal parallel...
متن کاملInverse Matrix using Gauss Elimination Method by OpenMP
OpenMP is an implementation program interface that might be utilized to explicitly immediate multi-threaded and it shared memory parallelism. OpenMP platform for specifications multi-processing via concurrent work between interested parties of hardware and software industry, governments and academia. OpenMP is not needs implemented identically by all vendors and it is not proposed for distribut...
متن کاملTEL-AVIV UNIVERSITY RAYMOND AND BEVERLY SACKLER FACULTY OF EXACT SCIENCES SCHOOL OF COMPUTER SCIENCE Designing Communication-Efficient Matrix Algorithms in Distributed-Memory Cilk
This thesis studies the relationship between parallelism, space and communication in dense matrix algorithms. We study existing matrix multiplication algorithms, specifically those that are designed for shared-memory multiprocessor machines (SMP’s). These machines are rapidly becoming commodity in the computer industry, but exploiting their computing power remains difficult. We improve algorith...
متن کاملA Simple OpenMP Scheme for Parallel Iteration Solvers in Finite Element Analysis
This study develops an OpenMP scheme to parallel the preconditioned conjugate gradient methods (PCG) in shared memory computers. The proposed method is simple and systematic, so a minor change in traditional PCG methods may produce effective parallelism. At first, the global stiffness matrix is renumbered in order to produce a parallel three-line form matrix, and a subroutine only needs to be c...
متن کامل