Title
Vectorizing matrix operations arising from PDE discretization on 9-point stencils
Abstract
When solving a system of PDEs, discretized on 9-point stencils over a nonrectangular domain, the linear systems that arise will have matrices with an irregular block structure. In this paper we discuss the vectorization of the matrix-vector multiply and of the Incomplete LU factorization and backsolve for these types of matrices. The performance of the matrix-vector multiply is already optimal for a small number of grid points (one result per clock cycle). For the ILU factorization and backsolve the vector performance is not as satisfying, partly because the resulting vector length is generally small and partly because of the heavy use of indirect addressing. A comparison with the general-purpose routines from the SLAP library shows a significant gain in computational time.
Year
DOI
Venue
1994
10.1007/BF01666907
The Journal of Supercomputing
Keywords
Field
DocType
Vectorization,nonsymmetric sparse linear systems,nonrectangular domain,matrix-vector multiplication,ILU preconditioning,hyperplane method
Discretization,Euclidean vector,Computer science,Matrix (mathematics),Parallel computing,Algorithm,Vectorization (mathematics),Incomplete LU factorization,Factorization,Matrix multiplication,Sparse matrix
Journal
Volume
Issue
ISSN
8
1
0920-8542
Citations 
PageRank 
References 
2
1.25
3
Authors
2
Name
Order
Citations
PageRank
Joke G Blom116123.43
J. G. Verwer213139.71