40 Years of Computing at Newcastle

Department Technical Report Series No. 401

Parallel Block LU Factorisation in the Global Element Method Part 2: Local Memory Implementation

C. Phillips

University of Newcastle upon Tyne. 1992

Abstract

The Global Element Methods is a (spectral) numerical technique for solving elliptic and parabolic differential equations in two space dimensions. The solution, defined as a set of polynomial coefficients, is obtained from a linear system of equations in which the coefficient matrix is block sparse. We describe a parallel LU decomposition of the matrix targeted at a local memory multiprocessor.