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 )