zungqr - generate an M-by-N complex matrix Q with orthonormal columns,
SUBROUTINE ZUNGQR( M, N, K, A, LDA, TAU, WORKIN, LWORKIN, INFO) DOUBLE COMPLEX A(LDA,*), TAU(*), WORKIN(*) INTEGER M, N, K, LDA, LWORKIN, INFO
SUBROUTINE ZUNGQR_64( M, N, K, A, LDA, TAU, WORKIN, LWORKIN, INFO) DOUBLE COMPLEX A(LDA,*), TAU(*), WORKIN(*) INTEGER*8 M, N, K, LDA, LWORKIN, INFO
SUBROUTINE UNGQR( M, [N], [K], A, [LDA], TAU, [WORKIN], [LWORKIN], * [INFO]) COMPLEX(8), DIMENSION(:) :: TAU, WORKIN COMPLEX(8), DIMENSION(:,:) :: A INTEGER :: M, N, K, LDA, LWORKIN, INFO
SUBROUTINE UNGQR_64( M, [N], [K], A, [LDA], TAU, [WORKIN], [LWORKIN], * [INFO]) COMPLEX(8), DIMENSION(:) :: TAU, WORKIN COMPLEX(8), DIMENSION(:,:) :: A INTEGER(8) :: M, N, K, LDA, LWORKIN, INFO
#include <sunperf.h>
void zungqr(int m, int n, int k, doublecomplex *a, int lda, doublecomplex *tau, int *info);
void zungqr_64(long m, long n, long k, doublecomplex *a, long lda, doublecomplex *tau, long *info);
zungqr generates an M-by-N complex 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 CGEQRF.
TAU(i)
must contain the scalar factor of the elementary
reflector H(i), as returned by CGEQRF.
WORKIN(1)
returns the optimal LWORKIN.
If LWORKIN = -1, then a workspace query is assumed; the routine only calculates the optimal size of the WORKIN array, returns this value as the first entry of the WORKIN array, and no error message related to LWORKIN is issued by XERBLA.
= 0: successful exit
< 0: if INFO = -i, the i-th argument has an illegal value