40 Years of Computing at Newcastle

Department Technical Report Series No. 385

Parallel Block LU Factorisation in the Global Element Method Part 1: Share Memory Implementation

C. Phillips

University of Newcastle upon Tyne. 1992

Abstract

The Global Element Method 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 shared memory multiprocessor.