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