|[eigen] question on numerical stability of non-pivoting householder algos|
[ Thread Index |
| More lists.tuxfamily.org/eigen Archives
- To: eigen <eigen@xxxxxxxxxxxxxxxxxxx>
- Subject: [eigen] question on numerical stability of non-pivoting householder algos
- From: Benoit Jacob <jacob.benoit.1@xxxxxxxxx>
- Date: Sun, 8 Aug 2010 15:57:51 -0400
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:mime-version:received:received:date:message-id :subject:from:to:content-type; bh=I+x8Mvfn+r8meXdHN/EDccqdB8DKfq7vJlyn6VDfu3k=; b=mUbGL3b84xJzYNyiDLx8vk2bd3e6VNvgN/JhhJ8sPpbFFoZMorq0LHj7hgeiFlmMvW o8KeNFPEbk2XpTelbAuPlCKXMIx9D3AyWXhn07Cz0dsRrbHVW8TSZPtSzPbkzvkjAd/v EgyrQ506E/Vv/C/0hXd+ZJLYAMwou5uea6t3U=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:date:message-id:subject:from:to:content-type; b=dsWx1+DlxqYo7GvQhAa0KULYzoKvaj4Mgxbl+Ckmk4NHmZcsVyP9bt2QkThstHzxsW ShR1k1wpBmv3GjeynCm3cKRcsrjasad/BMJfYv5LMvVmSvTGIwJI+a6bwfq01PzAy2gy kyNo5EkjKjWcTvxtRiC1oO61iUSdnSh3o6Qh0=
I'd like to ask a question to the numerics gurus here.
I've always been considering that plain HouseholderQR (no pivoting)
should be considered not-very-stable on singular matrices, since a
tiny perturbation on the matrix could result in a completely different
Is this correct? For sure HouseholderQR doesn't give as precise
results as pivoting variants, in my experiments. But how bad can it
be? Is there some upper bound on the imprecision on ill-conditioned
matrices or can it get arbitrarily bad?
Related question: bidiagonalization and tridiagonalizations are close
cousins of QR and will suffer from exactly the same problem. Worse:
contrary to QR, they can't be made pivoting (at least not in the
direction that matters the most).
Since bidiag and tridiag are the basic ingredient for SVD and
SelfadjointEigensolver respectively, this means that SVD and
SelfadjointEigensolver suffer from the same inaccuracy/instability
issues on near-singular matrices as HouseholderQR does.
I'm asking mostly because we're writing this page,
and need to know what to say about the accuracy/reliability of various algos.