cgemqrt - N matrix C with Q*C, C*Q, Q**H*C, or C*Q**H depending on values of SIDE and TRANS
SUBROUTINE CGEMQRT(SIDE, TRANS, M, N, K, NB, V, LDV, T, LDT, C, LDC, WORK, INFO) CHARACTER*1 SIDE, TRANS INTEGER INFO, K, LDV, LDC, M, N, NB, LDT COMPLEX V(LDV,*), C(LDC,*), T(LDT,*), WORK(*) SUBROUTINE CGEMQRT_64(SIDE, TRANS, M, N, K, NB, V, LDV, T, LDT, C, LDC, WORK, INFO) CHARACTER*1 SIDE, TRANS INTEGER*8 INFO, K, LDV, LDC, M, N, NB, LDT COMPLEX V(LDV,*), C(LDC,*), T(LDT,*), WORK(*) F95 INTERFACE SUBROUTINE GEMQRT(SIDE, TRANS, M, N, K, NB, V, LDV, T, LDT, C, LDC, WORK, INFO) INTEGER :: M, N, K, NB, LDV, LDT, LDC, INFO CHARACTER(LEN=1) :: SIDE, TRANS COMPLEX, DIMENSION(:,:) :: V, T, C COMPLEX, DIMENSION(:) :: WORK SUBROUTINE GEMQRT_64(SIDE, TRANS, M, N, K, NB, V, LDV, T, LDT, C, LDC, WORK, INFO) INTEGER(8) :: M, N, K, NB, LDV, LDT, LDC, INFO CHARACTER(LEN=1) :: SIDE, TRANS COMPLEX, DIMENSION(:,:) :: V, T, C COMPLEX, DIMENSION(:) :: WORK C INTERFACE #include <sunperf.h> void cgemqrt (char side, char trans, int m, int n, int k, int nb, floatcomplex *v, int ldv, floatcomplex *t, int ldt, floatcom- plex *c, int ldc, int *info); void cgemqrt_64 (char side, char trans, long m, long n, long k, long nb, floatcomplex *v, long ldv, floatcomplex *t, long ldt, floatcomplex *c, long ldc, long *info);
Oracle Solaris Studio Performance Library cgemqrt(3P) NAME cgemqrt - overwrite the general complex M-by-N matrix C with Q*C, C*Q, Q**H*C, or C*Q**H depending on values of SIDE and TRANS SYNOPSIS SUBROUTINE CGEMQRT(SIDE, TRANS, M, N, K, NB, V, LDV, T, LDT, C, LDC, WORK, INFO) CHARACTER*1 SIDE, TRANS INTEGER INFO, K, LDV, LDC, M, N, NB, LDT COMPLEX V(LDV,*), C(LDC,*), T(LDT,*), WORK(*) SUBROUTINE CGEMQRT_64(SIDE, TRANS, M, N, K, NB, V, LDV, T, LDT, C, LDC, WORK, INFO) CHARACTER*1 SIDE, TRANS INTEGER*8 INFO, K, LDV, LDC, M, N, NB, LDT COMPLEX V(LDV,*), C(LDC,*), T(LDT,*), WORK(*) F95 INTERFACE SUBROUTINE GEMQRT(SIDE, TRANS, M, N, K, NB, V, LDV, T, LDT, C, LDC, WORK, INFO) INTEGER :: M, N, K, NB, LDV, LDT, LDC, INFO CHARACTER(LEN=1) :: SIDE, TRANS COMPLEX, DIMENSION(:,:) :: V, T, C COMPLEX, DIMENSION(:) :: WORK SUBROUTINE GEMQRT_64(SIDE, TRANS, M, N, K, NB, V, LDV, T, LDT, C, LDC, WORK, INFO) INTEGER(8) :: M, N, K, NB, LDV, LDT, LDC, INFO CHARACTER(LEN=1) :: SIDE, TRANS COMPLEX, DIMENSION(:,:) :: V, T, C COMPLEX, DIMENSION(:) :: WORK C INTERFACE #include <sunperf.h> void cgemqrt (char side, char trans, int m, int n, int k, int nb, floatcomplex *v, int ldv, floatcomplex *t, int ldt, floatcom- plex *c, int ldc, int *info); void cgemqrt_64 (char side, char trans, long m, long n, long k, long nb, floatcomplex *v, long ldv, floatcomplex *t, long ldt, floatcomplex *c, long ldc, long *info); PURPOSE cgemqrt overwrites the general complex M-by-N matrix C with SIDE = 'L' SIDE = 'R' TRANS = 'N': Q*C C*Q TRANS = 'C': Q**H*C C*Q**H where Q is a complex orthogonal matrix defined as the product of K ele- mentary reflectors: Q = H(1)*H(2) . . . H(K) = I - V*T*V**H generated using the compact WY representation as returned by CGEQRT. Q is of order M if SIDE = 'L' and of order N if SIDE = 'R'. ARGUMENTS SIDE (input) SIDE is CHARACTER*1 = 'L': apply Q or Q**H from the Left; = 'R': apply Q or Q**H from the Right. TRANS (input) TRANS is CHARACTER*1 = 'N': No transpose, apply Q; = 'C': Transpose, apply Q**H. M (input) M is INTEGER The number of rows of the matrix C. M >= 0. N (input) N is INTEGER The number of columns of the matrix C. N >= 0. K (input) K is INTEGER The number of elementary reflectors whose product defines the matrix Q. If SIDE = 'L', M >= K >= 0; if SIDE = 'R', N >= K >= 0. NB (input) NB is INTEGER The block size used for the storage of T. K >= NB >= 1. This must be the same value of NB used to generate T in CGE- QRT. V (input) V is COMPLEX array, dimension (LDV,K) The i-th column must contain the vector which defines the elementary reflector H(i), for i = 1,2,...,k, as returned by CGEQRT in the first K columns of its array argument A. LDV (input) LDV is INTEGER The leading dimension of the array V. If SIDE = 'L', LDA >= max(1,M); if SIDE = 'R', LDA >= max(1,N). T (input) T is COMPLEX array, dimension (LDT,K) The upper triangular factors of the block reflectors as returned by CGEQRT, stored as a NB-by-N matrix. LDT (input) LDT is INTEGER The leading dimension of the array T. LDT >= NB. C (input/output) C is COMPLEX array, dimension (LDC,N) On entry, the M-by-N matrix C. On exit, C is overwritten by Q*C, Q**H*C, C*Q**H or C*Q. LDC (input) LDC is INTEGER The leading dimension of the array C. LDC >= max(1,M). WORK (output) WORK is COMPLEX array. The dimension of WORK is N*NB if SIDE = 'L', or M*NB if SIDE = 'R'. INFO (output) INFO is INTEGER = 0: successful exit < 0: if INFO = -i, the i-th argument had an illegal value 7 Nov 2015 cgemqrt(3P)