Subroutine */
| 33675 | } /* sorghr_ */ |
| 33676 | |
| 33677 | /* Subroutine */ int sorgl2_(integer *m, integer *n, integer *k, real *a, |
| 33678 | integer *lda, real *tau, real *work, integer *info) |
| 33679 | { |
| 33680 | /* System generated locals */ |
| 33681 | integer a_dim1, a_offset, i__1, i__2; |
| 33682 | real r__1; |
| 33683 | |
| 33684 | /* Local variables */ |
| 33685 | static integer i__, j, l; |
| 33686 | extern /* Subroutine */ int sscal_(integer *, real *, real *, integer *), |
| 33687 | slarf_(char *, integer *, integer *, real *, integer *, real *, |
| 33688 | real *, integer *, real *), xerbla_(char *, integer *); |
| 33689 | |
| 33690 | |
| 33691 | /* |
| 33692 | -- LAPACK routine (version 3.2) -- |
| 33693 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 33694 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 33695 | November 2006 |
| 33696 | |
| 33697 | |
| 33698 | Purpose |
| 33699 | ======= |
| 33700 | |
| 33701 | SORGL2 generates an m by n real matrix Q with orthonormal rows, |
| 33702 | which is defined as the first m rows of a product of k elementary |
| 33703 | reflectors of order n |
| 33704 | |
| 33705 | Q = H(k) . . . H(2) H(1) |
| 33706 | |
| 33707 | as returned by SGELQF. |
| 33708 | |
| 33709 | Arguments |
| 33710 | ========= |
| 33711 | |
| 33712 | M (input) INTEGER |
| 33713 | The number of rows of the matrix Q. M >= 0. |
| 33714 | |
| 33715 | N (input) INTEGER |
| 33716 | The number of columns of the matrix Q. N >= M. |
| 33717 | |
| 33718 | K (input) INTEGER |
| 33719 | The number of elementary reflectors whose product defines the |
| 33720 | matrix Q. M >= K >= 0. |
| 33721 | |
| 33722 | A (input/output) REAL array, dimension (LDA,N) |
| 33723 | On entry, the i-th row must contain the vector which defines |
| 33724 | the elementary reflector H(i), for i = 1,2,...,k, as returned |
| 33725 | by SGELQF in the first k rows of its array argument A. |
| 33726 | On exit, the m-by-n matrix Q. |
| 33727 | |
| 33728 | LDA (input) INTEGER |
| 33729 | The first dimension of the array A. LDA >= max(1,M). |
| 33730 | |
| 33731 | TAU (input) REAL array, dimension (K) |
| 33732 | TAU(i) must contain the scalar factor of the elementary |
| 33733 | reflector H(i), as returned by SGELQF. |
| 33734 |