79 DOUBLE PRECISION rw( 2*nmax )
80 COMPLEX*16 a( nmax, nmax ), tau( nmax ),
96 COMMON / infoc / infot, nout, ok, lerr
97 COMMON / srnamc / srnamt
107 a( 1, 1 ) = dcmplx( 1.0d+0, -1.0d+0 )
108 a( 1, 2 ) = dcmplx( 2.0d+0, -2.0d+0 )
109 a( 2, 2 ) = dcmplx( 3.0d+0, -3.0d+0 )
110 a( 2, 1 ) = dcmplx( 4.0d+0, -4.0d+0 )
112 WRITE( nout, fmt = * )
116 IF(
lsamen( 2, c2,
'QP' ) )
THEN
122 CALL
zgeqpf( -1, 0, a, 1, ip, tau, w, rw, info )
123 CALL
chkxer(
'ZGEQPF', infot, nout, lerr, ok )
125 CALL
zgeqpf( 0, -1, a, 1, ip, tau, w, rw, info )
126 CALL
chkxer(
'ZGEQPF', infot, nout, lerr, ok )
128 CALL
zgeqpf( 2, 0, a, 1, ip, tau, w, rw, info )
129 CALL
chkxer(
'ZGEQPF', infot, nout, lerr, ok )
135 CALL
zgeqp3( -1, 0, a, 1, ip, tau, w, lw, rw, info )
136 CALL
chkxer(
'ZGEQP3', infot, nout, lerr, ok )
138 CALL
zgeqp3( 1, -1, a, 1, ip, tau, w, lw, rw, info )
139 CALL
chkxer(
'ZGEQP3', infot, nout, lerr, ok )
141 CALL
zgeqp3( 2, 3, a, 1, ip, tau, w, lw, rw, info )
142 CALL
chkxer(
'ZGEQP3', infot, nout, lerr, ok )
144 CALL
zgeqp3( 2, 2, a, 2, ip, tau, w, lw-10, rw, info )
145 CALL
chkxer(
'ZGEQP3', infot, nout, lerr, ok )
150 CALL
alaesm( path, ok, nout )
subroutine zerrqp(PATH, NUNIT)
ZERRQP
subroutine zgeqp3(M, N, A, LDA, JPVT, TAU, WORK, LWORK, RWORK, INFO)
ZGEQP3
subroutine zgeqpf(M, N, A, LDA, JPVT, TAU, WORK, RWORK, INFO)
ZGEQPF
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
LOGICAL function lsamen(N, CA, CB)
LSAMEN
subroutine alaesm(PATH, OK, NOUT)
ALAESM