71 parameter ( nmax = 3 )
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 zgeqp3( -1, 0, a, 1, ip, tau, w, lw, rw, info )
123 CALL chkxer(
'ZGEQP3', infot, nout, lerr, ok )
125 CALL zgeqp3( 1, -1, a, 1, ip, tau, w, lw, rw, info )
126 CALL chkxer(
'ZGEQP3', infot, nout, lerr, ok )
128 CALL zgeqp3( 2, 3, a, 1, ip, tau, w, lw, rw, info )
129 CALL chkxer(
'ZGEQP3', infot, nout, lerr, ok )
131 CALL zgeqp3( 2, 2, a, 2, ip, tau, w, lw-10, rw, info )
132 CALL chkxer(
'ZGEQP3', infot, nout, lerr, ok )
137 CALL alaesm( path, ok, nout )
logical function lsamen(N, CA, CB)
LSAMEN
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine zgeqp3(M, N, A, LDA, JPVT, TAU, WORK, LWORK, RWORK, INFO)
ZGEQP3