Mathematics – Numerical Analysis
Scientific paper
2006-12-10
Numer. Math. 108 (2007), no. 1, 59-91
Mathematics
Numerical Analysis
26 pages; final version; to appear in Numerische Mathematik
Scientific paper
10.1007/s00211-007-0114-x
In an earlier paper, we showed that a large class of fast recursive matrix multiplication algorithms is stable in a normwise sense, and that in fact if multiplication of $n$-by-$n$ matrices can be done by any algorithm in $O(n^{\omega + \eta})$ operations for any $\eta > 0$, then it can be done stably in $O(n^{\omega + \eta})$ operations for any $\eta > 0$. Here we extend this result to show that essentially all standard linear algebra operations, including LU decomposition, QR decomposition, linear equation solving, matrix inversion, solving least squares problems, (generalized) eigenvalue problems and the singular value decomposition can also be done stably (in a normwise sense) in $O(n^{\omega + \eta})$ operations.
Demmel James
Dumitriu Ioana
Holtz Olga
No associations
LandOfFree
Fast linear algebra is stable does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with Fast linear algebra is stable, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Fast linear algebra is stable will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-183578