dorgrq - generate an M-by-N real matrix Q with orthonormal rows,
SUBROUTINE DORGRQ( M, N, K, A, LDA, TAU, WORK, LDWORK, INFO) INTEGER M, N, K, LDA, LDWORK, INFO DOUBLE PRECISION A(LDA,*), TAU(*), WORK(*)
SUBROUTINE DORGRQ_64( M, N, K, A, LDA, TAU, WORK, LDWORK, INFO) INTEGER*8 M, N, K, LDA, LDWORK, INFO DOUBLE PRECISION A(LDA,*), TAU(*), WORK(*)
SUBROUTINE ORGRQ( M, [N], [K], A, [LDA], TAU, [WORK], [LDWORK], * [INFO]) INTEGER :: M, N, K, LDA, LDWORK, INFO REAL(8), DIMENSION(:) :: TAU, WORK REAL(8), DIMENSION(:,:) :: A
SUBROUTINE ORGRQ_64( M, [N], [K], A, [LDA], TAU, [WORK], [LDWORK], * [INFO]) INTEGER(8) :: M, N, K, LDA, LDWORK, INFO REAL(8), DIMENSION(:) :: TAU, WORK REAL(8), DIMENSION(:,:) :: A
#include <sunperf.h>
void dorgrq(int m, int n, int k, double *a, int lda, double *tau, int *info);
void dorgrq_64(long m, long n, long k, double *a, long lda, double *tau, long *info);
dorgrq generates an M-by-N real matrix Q with orthonormal rows, which is defined as the last M rows of a product of K elementary reflectors of order N
Q = H(1) H(2) . . . H(k)
as returned by SGERQF.
TAU(i)
must contain the scalar factor of the elementary
reflector H(i), as returned by SGERQF.
WORK(1)
returns the optimal LDWORK.
If LDWORK = -1, then a workspace query is assumed; the routine only calculates the optimal size of the WORK array, returns this value as the first entry of the WORK array, and no error message related to LDWORK is issued by XERBLA.
= 0: successful exit
< 0: if INFO = -i, the i-th argument has an illegal value