Parallel Implementation of a Block Algorithm for Matrix 1-Norm Estimation

Authors:
Sheung Hun Cheng;Nicholas J. Higham
Affiliations:
-;-
Venue:
Euro-Par '01 Proceedings of the 7th International Euro-Par Conference Manchester on Parallel Processing
Year:
2001

Citing 6
Cited 0

Experience with a matrix norm estimator

SIAM Journal on Scientific and Statistical Computing
FORTRAN codes for estimating the one-norm of a real or complex matrix, with applications to condition estimation

ACM Transactions on Mathematical Software (TOMS)
ScaLAPACK user's guide

ScaLAPACK user's guide
LAPACK Users' guide (third ed.)

LAPACK Users' guide (third ed.)
A Block Algorithm for Matrix 1-Norm Estimation, with an Application to 1-Norm Pseudospectra

SIAM Journal on Matrix Analysis and Applications
Accuracy and Stability of Numerical Algorithms

Accuracy and Stability of Numerical Algorithms

Quantified Score

Hi-index	0.00

Visualization

Abstract

We describe a parallel Fortran 77 implementation, in ScaLA-PACK style, of a block matrix 1-norm estimator of Higham and Tisseur. This estimator differs from that underlying the existing ScaLAPACK code, PxLACON, in that it iterates with a matrix with t columns, where t ≥ 1 is a parameter, rather than with a vector, ands o the basic computational kernel is level 3 BLAS operations. Our experiments on an SGI Origin2000 show that with t = 2 or 4 the new code offers better estimates than PDLACON with a similar execution time. Moreover, with t 4, estimates exact over 90% of the time are achieved with execution time growing much slower than t.