Subroutine */
| 33064 | } /* slauum_ */ |
| 33065 | |
| 33066 | /* Subroutine */ int sorg2r_(integer *m, integer *n, integer *k, real *a, |
| 33067 | integer *lda, real *tau, real *work, integer *info) |
| 33068 | { |
| 33069 | /* System generated locals */ |
| 33070 | integer a_dim1, a_offset, i__1, i__2; |
| 33071 | real r__1; |
| 33072 | |
| 33073 | /* Local variables */ |
| 33074 | static integer i__, j, l; |
| 33075 | extern /* Subroutine */ int sscal_(integer *, real *, real *, integer *), |
| 33076 | slarf_(char *, integer *, integer *, real *, integer *, real *, |
| 33077 | real *, integer *, real *), xerbla_(char *, integer *); |
| 33078 | |
| 33079 | |
| 33080 | /* |
| 33081 | -- LAPACK routine (version 3.2) -- |
| 33082 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 33083 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 33084 | November 2006 |
| 33085 | |
| 33086 | |
| 33087 | Purpose |
| 33088 | ======= |
| 33089 | |
| 33090 | SORG2R generates an m by n real matrix Q with orthonormal columns, |
| 33091 | which is defined as the first n columns of a product of k elementary |
| 33092 | reflectors of order m |
| 33093 | |
| 33094 | Q = H(1) H(2) . . . H(k) |
| 33095 | |
| 33096 | as returned by SGEQRF. |
| 33097 | |
| 33098 | Arguments |
| 33099 | ========= |
| 33100 | |
| 33101 | M (input) INTEGER |
| 33102 | The number of rows of the matrix Q. M >= 0. |
| 33103 | |
| 33104 | N (input) INTEGER |
| 33105 | The number of columns of the matrix Q. M >= N >= 0. |
| 33106 | |
| 33107 | K (input) INTEGER |
| 33108 | The number of elementary reflectors whose product defines the |
| 33109 | matrix Q. N >= K >= 0. |
| 33110 | |
| 33111 | A (input/output) REAL array, dimension (LDA,N) |
| 33112 | On entry, the i-th column must contain the vector which |
| 33113 | defines the elementary reflector H(i), for i = 1,2,...,k, as |
| 33114 | returned by SGEQRF in the first k columns of its array |
| 33115 | argument A. |
| 33116 | On exit, the m-by-n matrix Q. |
| 33117 | |
| 33118 | LDA (input) INTEGER |
| 33119 | The first dimension of the array A. LDA >= max(1,M). |
| 33120 | |
| 33121 | TAU (input) REAL array, dimension (K) |
| 33122 | TAU(i) must contain the scalar factor of the elementary |
| 33123 | reflector H(i), as returned by SGEQRF. |