Self-scaling fast rotations for stiff and equality-constrained linear least squares problems

被引:14
作者
Anda, AA
Park, H
机构
[1] Computer Science Department, University of Minnesota, Minneapolis
关键词
D O I
10.1016/0024-3795(94)00092-1
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
We present algorithms which apply self-scaling fast plane rotations to the QR decomposition for stiff least squares problems. We show that both fast and standard Givens rotation-based algorithms produce accurate results, regardless of row sorting and even with extremely large weights, when equality-constrained least squares problems are solved by the weighting method. Numerical test results show that the Householder QR decomposition algorithm is sensitive to row sorting and produces less accurate results when the weights are large, and that the modified Gram-Schmidt algorithm is less sensitive to row sorting. This makes the fast plane rotation a method of choice for the QR decomposition of stiff matrices, since it is also competitive in computational complexity. Based on the above results, an efficient algorithm is also derived for the application where the least squares solutions are required for various constrained matrices for each fixed data matrix.
引用
收藏
页码:137 / 161
页数:25
相关论文
共 27 条
[1]  
ANDA A, UNPUB IMPLEMENTATION
[2]   FAST PLANE ROTATIONS WITH DYNAMIC SCALING [J].
ANDA, AA ;
PARK, HS .
SIAM JOURNAL ON MATRIX ANALYSIS AND APPLICATIONS, 1994, 15 (01) :162-174
[3]   SCALED GIVENS ROTATIONS FOR THE SOLUTION OF LINEAR LEAST-SQUARES PROBLEMS ON SYSTOLIC ARRAYS [J].
BARLOW, JL ;
IPSEN, ICF .
SIAM JOURNAL ON SCIENTIFIC AND STATISTICAL COMPUTING, 1987, 8 (05) :716-733
[4]   ITERATIVE METHODS FOR EQUALITY-CONSTRAINED LEAST-SQUARES PROBLEMS [J].
BARLOW, JL ;
NICHOLS, NK ;
PLEMMONS, RJ .
SIAM JOURNAL ON SCIENTIFIC AND STATISTICAL COMPUTING, 1988, 9 (05) :892-906
[6]  
BJORCK A, 1984, SIAM J SCI STAT COMP, V5, P394, DOI 10.1137/0905029
[7]  
BJORCK A, 1992, SIAM J MATRIX ANAL A, V13, P176
[8]  
BJORCK A, 1991, ALGORITHMS LINEAR LE
[9]  
Bjorck A., 1990, Handbook of Numerical Analysis, V1, P465, DOI DOI 10.1016/S1570-8659(05)80036-5
[10]   A ONE-SIDED JACOBI ALGORITHM FOR COMPUTING THE SINGULAR VALUE DECOMPOSITION ON A VECTOR COMPUTER [J].
DERIJK, PPM .
SIAM JOURNAL ON SCIENTIFIC AND STATISTICAL COMPUTING, 1989, 10 (02) :359-371