79 DOUBLE PRECISION a( nmax, nmax ), t( nmax, nmax ), w( nmax ),
92 COMMON / infoc / infot, nout, ok, lerr
93 COMMON / srnamc / srnamt
101 WRITE( nout, fmt = * )
107 a( i,
j ) = 1.d0 / dble( i+
j )
108 c( i,
j ) = 1.d0 / dble( i+
j )
109 t( i,
j ) = 1.d0 / dble( i+
j )
121 CALL
dgeqrt( -1, 0, 1, a, 1, t, 1, w, info )
122 CALL
chkxer(
'DGEQRT', infot, nout, lerr, ok )
124 CALL
dgeqrt( 0, -1, 1, a, 1, t, 1, w, info )
125 CALL
chkxer(
'DGEQRT', infot, nout, lerr, ok )
127 CALL
dgeqrt( 0, 0, 0, a, 1, t, 1, w, info )
128 CALL
chkxer(
'DGEQRT', infot, nout, lerr, ok )
130 CALL
dgeqrt( 2, 1, 1, a, 1, t, 1, w, info )
131 CALL
chkxer(
'DGEQRT', infot, nout, lerr, ok )
133 CALL
dgeqrt( 2, 2, 2, a, 2, t, 1, w, info )
134 CALL
chkxer(
'DGEQRT', infot, nout, lerr, ok )
140 CALL
dgeqrt2( -1, 0, a, 1, t, 1, info )
141 CALL
chkxer(
'DGEQRT2', infot, nout, lerr, ok )
143 CALL
dgeqrt2( 0, -1, a, 1, t, 1, info )
144 CALL
chkxer(
'DGEQRT2', infot, nout, lerr, ok )
146 CALL
dgeqrt2( 2, 1, a, 1, t, 1, info )
147 CALL
chkxer(
'DGEQRT2', infot, nout, lerr, ok )
149 CALL
dgeqrt2( 2, 2, a, 2, t, 1, info )
150 CALL
chkxer(
'DGEQRT2', infot, nout, lerr, ok )
156 CALL
dgeqrt3( -1, 0, a, 1, t, 1, info )
157 CALL
chkxer(
'DGEQRT3', infot, nout, lerr, ok )
159 CALL
dgeqrt3( 0, -1, a, 1, t, 1, info )
160 CALL
chkxer(
'DGEQRT3', infot, nout, lerr, ok )
162 CALL
dgeqrt3( 2, 1, a, 1, t, 1, info )
163 CALL
chkxer(
'DGEQRT3', infot, nout, lerr, ok )
165 CALL
dgeqrt3( 2, 2, a, 2, t, 1, info )
166 CALL
chkxer(
'DGEQRT3', infot, nout, lerr, ok )
172 CALL
dgemqrt(
'/',
'N', 0, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
173 CALL
chkxer(
'DGEMQRT', infot, nout, lerr, ok )
175 CALL
dgemqrt(
'L',
'/', 0, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
176 CALL
chkxer(
'DGEMQRT', infot, nout, lerr, ok )
178 CALL
dgemqrt(
'L',
'N', -1, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
179 CALL
chkxer(
'DGEMQRT', infot, nout, lerr, ok )
181 CALL
dgemqrt(
'L',
'N', 0, -1, 0, 1, a, 1, t, 1, c, 1, w, info )
182 CALL
chkxer(
'DGEMQRT', infot, nout, lerr, ok )
184 CALL
dgemqrt(
'L',
'N', 0, 0, -1, 1, a, 1, t, 1, c, 1, w, info )
185 CALL
chkxer(
'DGEMQRT', infot, nout, lerr, ok )
187 CALL
dgemqrt(
'R',
'N', 0, 0, -1, 1, a, 1, t, 1, c, 1, w, info )
188 CALL
chkxer(
'DGEMQRT', infot, nout, lerr, ok )
190 CALL
dgemqrt(
'L',
'N', 0, 0, 0, 0, a, 1, t, 1, c, 1, w, info )
191 CALL
chkxer(
'DGEMQRT', infot, nout, lerr, ok )
193 CALL
dgemqrt(
'R',
'N', 1, 2, 1, 1, a, 1, t, 1, c, 1, w, info )
194 CALL
chkxer(
'DGEMQRT', infot, nout, lerr, ok )
196 CALL
dgemqrt(
'L',
'N', 2, 1, 1, 1, a, 1, t, 1, c, 1, w, info )
197 CALL
chkxer(
'DGEMQRT', infot, nout, lerr, ok )
199 CALL
dgemqrt(
'R',
'N', 1, 1, 1, 1, a, 1, t, 0, c, 1, w, info )
200 CALL
chkxer(
'DGEMQRT', infot, nout, lerr, ok )
202 CALL
dgemqrt(
'L',
'N', 1, 1, 1, 1, a, 1, t, 1, c, 0, w, info )
203 CALL
chkxer(
'DGEMQRT', infot, nout, lerr, ok )
207 CALL
alaesm( path, ok, nout )
subroutine dgeqrt2(M, N, A, LDA, T, LDT, INFO)
DGEQRT2 computes a QR factorization of a general real or complex matrix using the compact WY represen...
subroutine derrqrt(PATH, NUNIT)
DERRQRT
subroutine dgemqrt(SIDE, TRANS, M, N, K, NB, V, LDV, T, LDT, C, LDC, WORK, INFO)
DGEMQRT
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
recursive subroutine dgeqrt3(M, N, A, LDA, T, LDT, INFO)
DGEQRT3 recursively computes a QR factorization of a general real or complex matrix using the compact...
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine dgeqrt(M, N, NB, A, LDA, T, LDT, WORK, INFO)
DGEQRT