Results 1 
6 of
6
Computing RankRevealing QR Factorizations of Dense Matrices
 Argonne Preprint ANLMCSP5590196, Argonne National Laboratory
, 1996
"... this paper, and we give only a brief synopsis here. For details, the reader is referred to the code. Test matrices 1 through 5 were designed to exercise column pivoting. Matrix 6 was designed to test the behavior of the condition estimation in the presence of clusters for the smallest singular value ..."
Abstract

Cited by 26 (1 self)
 Add to MetaCart
this paper, and we give only a brief synopsis here. For details, the reader is referred to the code. Test matrices 1 through 5 were designed to exercise column pivoting. Matrix 6 was designed to test the behavior of the condition estimation in the presence of clusters for the smallest singular value. For the other cases, we employed the LAPACK matrix generator xLATMS, which generates random symmetric matrices by multiplying a diagonal matrix with prescribed singular values by random orthogonal matrices from the left and right. For the break1 distribution, all singular values are 1.0 except for one. In the arithmetic and geometric distributions, they decay from 1.0 to a specified smallest singular value in an arithmetic and geometric fashion, respectively. In the "reversed" distributions, the order of the diagonal entries was reversed. For test cases 7 though 12, we used xLATMS to generate a matrix of order
A BLAS3 version of the QR factorization with column pivoting
 SIAM J. SCI. COMPUT
, 1995
"... The QR factorization with column pivoting (QRP), originally suggested by Golub and Businger in 1965, is a popular approach to computing rankrevealing factorizations. Using BLAS Level 1, it was implemented in LINPACK, and, using BLAS Level 2, in LAPACK. While the BLAS Level2version delivers, in gen ..."
Abstract

Cited by 17 (2 self)
 Add to MetaCart
The QR factorization with column pivoting (QRP), originally suggested by Golub and Businger in 1965, is a popular approach to computing rankrevealing factorizations. Using BLAS Level 1, it was implemented in LINPACK, and, using BLAS Level 2, in LAPACK. While the BLAS Level2version delivers, in general, superior performance, it may result in worse performance for large matrix sizes due to cache e ects. We introduce a modi cation of the QRP algorithm which allows the use of BLAS Level 3 kernels while maintaining the numerical behavior of the LINPACK and LAPACK implementations. Experimental comparisons of this approach with the LINPACK and LAPACK implementations on IBM RS/6000, SGI R8000, and DEC Alpha platforms show considerable performance improvements.
On Iterative Algorithms for Linear Least Squares Problems With Bound Constraints
, 1995
"... Three new iterative methods for the solution of the linear least squares problem with bound constraints are presented and their performance analyzed. The first is a modification of a method proposed by Lotstedt, while the two others are characterized by a technique allowing for fast active set chang ..."
Abstract

Cited by 15 (2 self)
 Add to MetaCart
Three new iterative methods for the solution of the linear least squares problem with bound constraints are presented and their performance analyzed. The first is a modification of a method proposed by Lotstedt, while the two others are characterized by a technique allowing for fast active set changes resulting in noticeable improvements on the speed at which constraints active at the solution are identified. The numerical efficiency of these algorithms is experimentally studied, with particular emphasis on dependence of the starting point choice and the use of preconditioning for illconditioned problems.
Mathematical Models for Transportation Demand Analysis
, 1996
"... this paper, we will concentrate on the overspeci#cation arising from the ASCs and will not consider other possible errors sources. ..."
Abstract

Cited by 4 (1 self)
 Add to MetaCart
this paper, we will concentrate on the overspeci#cation arising from the ASCs and will not consider other possible errors sources.
BlockPartitioned Algorithms for Solving the Linear Least Squares Problem
"... The linear least squares problem arises in many areas of sciences and engineerings. When the coe cient matrix has full rank, the solution can be obtained in a fast way by using QR factorization with BLAS3. In contrast, when the matrix is rankdeficient, or the rank is unknown, other slower methods ..."
Abstract
 Add to MetaCart
The linear least squares problem arises in many areas of sciences and engineerings. When the coe cient matrix has full rank, the solution can be obtained in a fast way by using QR factorization with BLAS3. In contrast, when the matrix is rankdeficient, or the rank is unknown, other slower methods should be applied: the SVD or the complete orthogonal decompositions. The SVD gives more reliable determination of rank but is computationally more expensive. On the other hand, the complete orthogonal decomposition is faster and in practice works well. We present several new implementations for solving the linear least squares problem by means of the complete orthogonal decomposition that are faster than the algorithms currently included in LAPACK. Experimental comparison of our methods with the LAPACK implementations on a wide range of platforms (such as IBM RS/6000370, SUN HyperSPARC, SGI R8000, DEC Alpha/AXP, HP 9000/715, etc.) show considerable performance improvements. Some of the new code has been already included in the latest release of LAPACK (3.0). In addition, for fullrank matrices the performances of the new methods are very close to the performance of the fast method based on QR factorization with BLAS3, thus providing a valuable general tool for fullrank matrices and rankdeficient matrices, as well as those matrices with unknown rank.
par
"... Résolution de problèmes de moindres carrés linéaires denses de grande taille sur des calculateurs parallèles distribués. Application au calcul de champ de gravité terrestre. Solving large dense linear least squares problems on parallel distributed computers. Application to the Earth’s gravity field ..."
Abstract
 Add to MetaCart
Résolution de problèmes de moindres carrés linéaires denses de grande taille sur des calculateurs parallèles distribués. Application au calcul de champ de gravité terrestre. Solving large dense linear least squares problems on parallel distributed computers. Application to the Earth’s gravity field computation. Thèse présentée le 21 Mars 2006 à Toulouse devant le jury composé de: