sorgrq - generate an M-by-N real matrix Q with orthonormal rows,
SUBROUTINE SORGRQ( M, N, K, A, LDA, TAU, WORK, LDWORK, INFO) INTEGER M, N, K, LDA, LDWORK, INFO REAL A(LDA,*), TAU(*), WORK(*)
SUBROUTINE SORGRQ_64( M, N, K, A, LDA, TAU, WORK, LDWORK, INFO) INTEGER*8 M, N, K, LDA, LDWORK, INFO REAL A(LDA,*), TAU(*), WORK(*)
SUBROUTINE ORGRQ( M, [N], [K], A, [LDA], TAU, [WORK], [LDWORK], * [INFO]) INTEGER :: M, N, K, LDA, LDWORK, INFO REAL, DIMENSION(:) :: TAU, WORK REAL, DIMENSION(:,:) :: A
SUBROUTINE ORGRQ_64( M, [N], [K], A, [LDA], TAU, [WORK], [LDWORK], * [INFO]) INTEGER(8) :: M, N, K, LDA, LDWORK, INFO REAL, DIMENSION(:) :: TAU, WORK REAL, DIMENSION(:,:) :: A
#include <sunperf.h>
void sorgrq(int m, int n, int k, float *a, int lda, float *tau, int *info);
void sorgrq_64(long m, long n, long k, float *a, long lda, float *tau, long *info);
sorgrq 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