sorgqr - generate an M-by-N real matrix Q with orthonormal columns,
SUBROUTINE SORGQR( M, N, K, A, LDA, TAU, WORK, LDWORK, INFO) INTEGER M, N, K, LDA, LDWORK, INFO REAL A(LDA,*), TAU(*), WORK(*)
SUBROUTINE SORGQR_64( M, N, K, A, LDA, TAU, WORK, LDWORK, INFO) INTEGER*8 M, N, K, LDA, LDWORK, INFO REAL A(LDA,*), TAU(*), WORK(*)
SUBROUTINE ORGQR( M, [N], [K], A, [LDA], TAU, [WORK], [LDWORK], * [INFO]) INTEGER :: M, N, K, LDA, LDWORK, INFO REAL, DIMENSION(:) :: TAU, WORK REAL, DIMENSION(:,:) :: A
SUBROUTINE ORGQR_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 sorgqr(int m, int n, int k, float *a, int lda, float *tau, int *info);
void sorgqr_64(long m, long n, long k, float *a, long lda, float *tau, long *info);
sorgqr generates an M-by-N real matrix Q with orthonormal columns, which is defined as the first N columns of a product of K elementary reflectors of order M
Q = H(1) H(2) . . . H(k)
as returned by SGEQRF.
TAU(i)
must contain the scalar factor of the elementary
reflector H(i), as returned by SGEQRF.
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