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
dgeqpf( -1, 0, a, 1, ip, tau, w, info )
118 CALL
chkxer(
'DGEQPF', infot, nout, lerr, ok )
120 CALL
dgeqpf( 0, -1, a, 1, ip, tau, w, info )
121 CALL
chkxer(
'DGEQPF', infot, nout, lerr, ok )
123 CALL
dgeqpf( 2, 0, a, 1, ip, tau, w, info )
124 CALL
chkxer(
'DGEQPF', infot, nout, lerr, ok )
130 CALL
dgeqp3( -1, 0, a, 1, ip, tau, w, lw, info )
131 CALL
chkxer(
'DGEQP3', infot, nout, lerr, ok )
133 CALL
dgeqp3( 1, -1, a, 1, ip, tau, w, lw, info )
134 CALL
chkxer(
'DGEQP3', infot, nout, lerr, ok )
136 CALL
dgeqp3( 2, 3, a, 1, ip, tau, w, lw, info )
137 CALL
chkxer(
'DGEQP3', infot, nout, lerr, ok )
139 CALL
dgeqp3( 2, 2, a, 2, ip, tau, w, lw-10, info )
140 CALL
chkxer(
'DGEQP3', infot, nout, lerr, ok )
145 CALL
alaesm( path, ok, nout )
subroutine derrqp(PATH, NUNIT)
DERRQP
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
logical function lsamen(N, CA, CB)
LSAMEN
subroutine dgeqpf(M, N, A, LDA, JPVT, TAU, WORK, INFO)
DGEQPF
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine dgeqp3(M, N, A, LDA, JPVT, TAU, WORK, LWORK, INFO)
DGEQP3