71 parameter ( nmax = 3 )
79 DOUBLE PRECISION a( nmax, nmax ), tau( nmax ), w( 3*nmax+1 )
94 COMMON / infoc / infot, nout, ok, lerr
95 COMMON / srnamc / srnamt
100 WRITE( nout, fmt = * )
109 IF(
lsamen( 2, c2,
'QP' ) )
THEN
117 CALL dgeqp3( -1, 0, a, 1, ip, tau, w, lw, info )
118 CALL chkxer(
'DGEQP3', infot, nout, lerr, ok )
120 CALL dgeqp3( 1, -1, a, 1, ip, tau, w, lw, info )
121 CALL chkxer(
'DGEQP3', infot, nout, lerr, ok )
123 CALL dgeqp3( 2, 3, a, 1, ip, tau, w, lw, info )
124 CALL chkxer(
'DGEQP3', infot, nout, lerr, ok )
126 CALL dgeqp3( 2, 2, a, 2, ip, tau, w, lw-10, info )
127 CALL chkxer(
'DGEQP3', infot, nout, lerr, ok )
132 CALL alaesm( path, ok, nout )
logical function lsamen(N, CA, CB)
LSAMEN
subroutine dgeqp3(M, N, A, LDA, JPVT, TAU, WORK, LWORK, INFO)
DGEQP3
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)