subroutine clarft(DIRECT, STOREV, N, K, V, LDV, TAU, T, LDT)
CLARFT forms the triangular factor T of a block reflector H = I - vtvH
subroutine clarfb(SIDE, TRANS, DIRECT, STOREV, M, N, K, V, LDV, T, LDT, C, LDC, WORK, LDWORK)
CLARFB applies a block reflector or its conjugate-transpose to a general rectangular matrix.
subroutine cung2r(M, N, K, A, LDA, TAU, WORK, INFO)
CUNG2R
integer function ilaenv(ISPEC, NAME, OPTS, N1, N2, N3, N4)
ILAENV