55 SUBROUTINE zerrqp( PATH, NUNIT )
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)
subroutine alaesm(PATH, OK, NOUT)
ALAESM