Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

A parallel implementation of the invariant subspace decomposition algorithm for dense symmetric matrices

Conference ·
OSTI ID:54430
;  [1];  [2]
  1. Supercomputing Research Center, Bowie, MD (United States)
  2. CONVEX Computer Corp., Richardson, TX (United States)
We give an overview of the Invariant Subspace Decomposition Algorithm for dense symmetric matrices (SYISDA) by first describing the algorithm, followed by a discussion of a parallel implementation of SYISDA on the Intel Delta. Our implementation utilizes an optimized parallel matrix multiplication implementation we have developed. Load balancing in the costly early stages of the algorithm is accomplished without redistribution of data between stages through the use of the block scattered decomposition. Computation of the invariant subspaces at each stage is done using a new tridiagonalization scheme due to Bischof and Sun.
OSTI ID:
54430
Report Number(s):
DOE/ER/25151--1-Vol.1; CONF-930331--Vol.1; CNN: Contract DM28E04120
Country of Publication:
United States
Language:
English