Subroutine */
| 34082 | } /* sorglq_ */ |
| 34083 | |
| 34084 | /* Subroutine */ int sorgqr_(integer *m, integer *n, integer *k, real *a, |
| 34085 | integer *lda, real *tau, real *work, integer *lwork, integer *info) |
| 34086 | { |
| 34087 | /* System generated locals */ |
| 34088 | integer a_dim1, a_offset, i__1, i__2, i__3; |
| 34089 | |
| 34090 | /* Local variables */ |
| 34091 | static integer i__, j, l, ib, nb, ki, kk, nx, iws, nbmin, iinfo; |
| 34092 | extern /* Subroutine */ int sorg2r_(integer *, integer *, integer *, real |
| 34093 | *, integer *, real *, real *, integer *), slarfb_(char *, char *, |
| 34094 | char *, char *, integer *, integer *, integer *, real *, integer * |
| 34095 | , real *, integer *, real *, integer *, real *, integer *), xerbla_(char *, integer *); |
| 34096 | extern integer ilaenv_(integer *, char *, char *, integer *, integer *, |
| 34097 | integer *, integer *, ftnlen, ftnlen); |
| 34098 | extern /* Subroutine */ int slarft_(char *, char *, integer *, integer *, |
| 34099 | real *, integer *, real *, real *, integer *); |
| 34100 | static integer ldwork, lwkopt; |
| 34101 | static logical lquery; |
| 34102 | |
| 34103 | |
| 34104 | /* |
| 34105 | -- LAPACK routine (version 3.2) -- |
| 34106 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 34107 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 34108 | November 2006 |
| 34109 | |
| 34110 | |
| 34111 | Purpose |
| 34112 | ======= |
| 34113 | |
| 34114 | SORGQR generates an M-by-N real matrix Q with orthonormal columns, |
| 34115 | which is defined as the first N columns of a product of K elementary |
| 34116 | reflectors of order M |
| 34117 | |
| 34118 | Q = H(1) H(2) . . . H(k) |
| 34119 | |
| 34120 | as returned by SGEQRF. |
| 34121 | |
| 34122 | Arguments |
| 34123 | ========= |
| 34124 | |
| 34125 | M (input) INTEGER |
| 34126 | The number of rows of the matrix Q. M >= 0. |
| 34127 | |
| 34128 | N (input) INTEGER |
| 34129 | The number of columns of the matrix Q. M >= N >= 0. |
| 34130 | |
| 34131 | K (input) INTEGER |
| 34132 | The number of elementary reflectors whose product defines the |
| 34133 | matrix Q. N >= K >= 0. |
| 34134 | |
| 34135 | A (input/output) REAL array, dimension (LDA,N) |
| 34136 | On entry, the i-th column must contain the vector which |
| 34137 | defines the elementary reflector H(i), for i = 1,2,...,k, as |
| 34138 | returned by SGEQRF in the first k columns of its array |
| 34139 | argument A. |
| 34140 | On exit, the M-by-N matrix Q. |
| 34141 |