Loading…
A Lanczos bidiagonalization algorithm for Hankel matrices
This paper presents a fast algorithm for bidiagonalizing a Hankel matrix. An m × n Hankel matrix is reduced to a real bidiagonal matrix in O ( ( m + n ) n log ( m + n ) ) floating-point operations (flops) using the Lanczos method with modified partial orthogonalization and reset schemes to improve i...
Saved in:
Published in: | Linear algebra and its applications 2009-03, Vol.430 (5), p.1531-1543 |
---|---|
Main Authors: | , , |
Format: | Article |
Language: | English |
Subjects: | |
Citations: | Items that this one cites Items that cite this one |
Online Access: | Get full text |
Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
Summary: | This paper presents a fast algorithm for bidiagonalizing a Hankel matrix. An
m
×
n
Hankel matrix is reduced to a real bidiagonal matrix in
O
(
(
m
+
n
)
n
log
(
m
+
n
)
)
floating-point operations (flops) using the Lanczos method with modified partial orthogonalization and reset schemes to improve its stability. Performance improvement is achieved by exploiting the Hankel structure, as fast Hankel matrix–vector multiplication is used. The accuracy and efficiency of the algorithm are demonstrated by our numerical experiments. |
---|---|
ISSN: | 0024-3795 1873-1856 |
DOI: | 10.1016/j.laa.2008.01.012 |