Subroutine */
| 33803 | } /* dorghr_ */ |
| 33804 | |
| 33805 | /* Subroutine */ int dorgl2_(integer *m, integer *n, integer *k, doublereal * |
| 33806 | a, integer *lda, doublereal *tau, doublereal *work, integer *info) |
| 33807 | { |
| 33808 | /* System generated locals */ |
| 33809 | integer a_dim1, a_offset, i__1, i__2; |
| 33810 | doublereal d__1; |
| 33811 | |
| 33812 | /* Local variables */ |
| 33813 | static integer i__, j, l; |
| 33814 | extern /* Subroutine */ int dscal_(integer *, doublereal *, doublereal *, |
| 33815 | integer *), dlarf_(char *, integer *, integer *, doublereal *, |
| 33816 | integer *, doublereal *, doublereal *, integer *, doublereal *), xerbla_(char *, integer *); |
| 33817 | |
| 33818 | |
| 33819 | /* |
| 33820 | -- LAPACK routine (version 3.2) -- |
| 33821 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 33822 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 33823 | November 2006 |
| 33824 | |
| 33825 | |
| 33826 | Purpose |
| 33827 | ======= |
| 33828 | |
| 33829 | DORGL2 generates an m by n real matrix Q with orthonormal rows, |
| 33830 | which is defined as the first m rows of a product of k elementary |
| 33831 | reflectors of order n |
| 33832 | |
| 33833 | Q = H(k) . . . H(2) H(1) |
| 33834 | |
| 33835 | as returned by DGELQF. |
| 33836 | |
| 33837 | Arguments |
| 33838 | ========= |
| 33839 | |
| 33840 | M (input) INTEGER |
| 33841 | The number of rows of the matrix Q. M >= 0. |
| 33842 | |
| 33843 | N (input) INTEGER |
| 33844 | The number of columns of the matrix Q. N >= M. |
| 33845 | |
| 33846 | K (input) INTEGER |
| 33847 | The number of elementary reflectors whose product defines the |
| 33848 | matrix Q. M >= K >= 0. |
| 33849 | |
| 33850 | A (input/output) DOUBLE PRECISION array, dimension (LDA,N) |
| 33851 | On entry, the i-th row must contain the vector which defines |
| 33852 | the elementary reflector H(i), for i = 1,2,...,k, as returned |
| 33853 | by DGELQF in the first k rows of its array argument A. |
| 33854 | On exit, the m-by-n matrix Q. |
| 33855 | |
| 33856 | LDA (input) INTEGER |
| 33857 | The first dimension of the array A. LDA >= max(1,M). |
| 33858 | |
| 33859 | TAU (input) DOUBLE PRECISION array, dimension (K) |
| 33860 | TAU(i) must contain the scalar factor of the elementary |
| 33861 | reflector H(i), as returned by DGELQF. |
| 33862 |