Subroutine */
| 35941 | } /* dormlq_ */ |
| 35942 | |
| 35943 | /* Subroutine */ int dormql_(char *side, char *trans, integer *m, integer *n, |
| 35944 | integer *k, doublereal *a, integer *lda, doublereal *tau, doublereal * |
| 35945 | c__, integer *ldc, doublereal *work, integer *lwork, integer *info) |
| 35946 | { |
| 35947 | /* System generated locals */ |
| 35948 | address a__1[2]; |
| 35949 | integer a_dim1, a_offset, c_dim1, c_offset, i__1, i__2, i__3[2], i__4, |
| 35950 | i__5; |
| 35951 | char ch__1[2]; |
| 35952 | |
| 35953 | /* Local variables */ |
| 35954 | static integer i__; |
| 35955 | static doublereal t[4160] /* was [65][64] */; |
| 35956 | static integer i1, i2, i3, ib, nb, mi, ni, nq, nw, iws; |
| 35957 | static logical left; |
| 35958 | extern logical lsame_(char *, char *); |
| 35959 | static integer nbmin, iinfo; |
| 35960 | extern /* Subroutine */ int dorm2l_(char *, char *, integer *, integer *, |
| 35961 | integer *, doublereal *, integer *, doublereal *, doublereal *, |
| 35962 | integer *, doublereal *, integer *), dlarfb_(char |
| 35963 | *, char *, char *, char *, integer *, integer *, integer *, |
| 35964 | doublereal *, integer *, doublereal *, integer *, doublereal *, |
| 35965 | integer *, doublereal *, integer *), dlarft_(char *, char *, integer *, integer *, doublereal |
| 35966 | *, integer *, doublereal *, doublereal *, integer *), xerbla_(char *, integer *); |
| 35967 | extern integer ilaenv_(integer *, char *, char *, integer *, integer *, |
| 35968 | integer *, integer *, ftnlen, ftnlen); |
| 35969 | static logical notran; |
| 35970 | static integer ldwork, lwkopt; |
| 35971 | static logical lquery; |
| 35972 | |
| 35973 | |
| 35974 | /* |
| 35975 | -- LAPACK routine (version 3.2) -- |
| 35976 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 35977 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 35978 | November 2006 |
| 35979 | |
| 35980 | |
| 35981 | Purpose |
| 35982 | ======= |
| 35983 | |
| 35984 | DORMQL overwrites the general real M-by-N matrix C with |
| 35985 | |
| 35986 | SIDE = 'L' SIDE = 'R' |
| 35987 | TRANS = 'N': Q * C C * Q |
| 35988 | TRANS = 'T': Q**T * C C * Q**T |
| 35989 | |
| 35990 | where Q is a real orthogonal matrix defined as the product of k |
| 35991 | elementary reflectors |
| 35992 | |
| 35993 | Q = H(k) . . . H(2) H(1) |
| 35994 | |
| 35995 | as returned by DGEQLF. Q is of order M if SIDE = 'L' and of order N |
| 35996 | if SIDE = 'R'. |
| 35997 | |
| 35998 | Arguments |
| 35999 | ========= |
| 36000 |