78 COMPLEX*16 a( nmax, nmax ), af( nmax, nmax ),
b( nmax ),
79 $ w( nmax ),
x( nmax )
91 COMMON / infoc / infot, nout, ok, lerr
92 COMMON / srnamc / srnamt
95 INTRINSIC dble, dcmplx
100 WRITE( nout, fmt = * )
106 a( i,
j ) = dcmplx( 1.d0 / dble( i+
j ),
107 $ -1.d0 / dble( i+
j ) )
108 af( i,
j ) = dcmplx( 1.d0 / dble( i+
j ),
109 $ -1.d0 / dble( i+
j ) )
123 CALL
zgerqf( -1, 0, a, 1,
b, w, 1, info )
124 CALL
chkxer(
'ZGERQF', infot, nout, lerr, ok )
126 CALL
zgerqf( 0, -1, a, 1,
b, w, 1, info )
127 CALL
chkxer(
'ZGERQF', infot, nout, lerr, ok )
129 CALL
zgerqf( 2, 1, a, 1,
b, w, 2, info )
130 CALL
chkxer(
'ZGERQF', infot, nout, lerr, ok )
132 CALL
zgerqf( 2, 1, a, 2,
b, w, 1, info )
133 CALL
chkxer(
'ZGERQF', infot, nout, lerr, ok )
139 CALL
zgerq2( -1, 0, a, 1,
b, w, info )
140 CALL
chkxer(
'ZGERQ2', infot, nout, lerr, ok )
142 CALL
zgerq2( 0, -1, a, 1,
b, w, info )
143 CALL
chkxer(
'ZGERQ2', infot, nout, lerr, ok )
145 CALL
zgerq2( 2, 1, a, 1,
b, w, info )
146 CALL
chkxer(
'ZGERQ2', infot, nout, lerr, ok )
152 CALL
zgerqs( -1, 0, 0, a, 1,
x,
b, 1, w, 1, info )
153 CALL
chkxer(
'ZGERQS', infot, nout, lerr, ok )
155 CALL
zgerqs( 0, -1, 0, a, 1,
x,
b, 1, w, 1, info )
156 CALL
chkxer(
'ZGERQS', infot, nout, lerr, ok )
158 CALL
zgerqs( 2, 1, 0, a, 2,
x,
b, 1, w, 1, info )
159 CALL
chkxer(
'ZGERQS', infot, nout, lerr, ok )
161 CALL
zgerqs( 0, 0, -1, a, 1,
x,
b, 1, w, 1, info )
162 CALL
chkxer(
'ZGERQS', infot, nout, lerr, ok )
164 CALL
zgerqs( 2, 2, 0, a, 1,
x,
b, 2, w, 1, info )
165 CALL
chkxer(
'ZGERQS', infot, nout, lerr, ok )
167 CALL
zgerqs( 2, 2, 0, a, 2,
x,
b, 1, w, 1, info )
168 CALL
chkxer(
'ZGERQS', infot, nout, lerr, ok )
170 CALL
zgerqs( 1, 1, 2, a, 1,
x,
b, 1, w, 1, info )
171 CALL
chkxer(
'ZGERQS', infot, nout, lerr, ok )
177 CALL
zungrq( -1, 0, 0, a, 1,
x, w, 1, info )
178 CALL
chkxer(
'ZUNGRQ', infot, nout, lerr, ok )
180 CALL
zungrq( 0, -1, 0, a, 1,
x, w, 1, info )
181 CALL
chkxer(
'ZUNGRQ', infot, nout, lerr, ok )
183 CALL
zungrq( 2, 1, 0, a, 2,
x, w, 2, info )
184 CALL
chkxer(
'ZUNGRQ', infot, nout, lerr, ok )
186 CALL
zungrq( 0, 0, -1, a, 1,
x, w, 1, info )
187 CALL
chkxer(
'ZUNGRQ', infot, nout, lerr, ok )
189 CALL
zungrq( 1, 2, 2, a, 1,
x, w, 1, info )
190 CALL
chkxer(
'ZUNGRQ', infot, nout, lerr, ok )
192 CALL
zungrq( 2, 2, 0, a, 1,
x, w, 2, info )
193 CALL
chkxer(
'ZUNGRQ', infot, nout, lerr, ok )
195 CALL
zungrq( 2, 2, 0, a, 2,
x, w, 1, info )
196 CALL
chkxer(
'ZUNGRQ', infot, nout, lerr, ok )
202 CALL
zungr2( -1, 0, 0, a, 1,
x, w, info )
203 CALL
chkxer(
'ZUNGR2', infot, nout, lerr, ok )
205 CALL
zungr2( 0, -1, 0, a, 1,
x, w, info )
206 CALL
chkxer(
'ZUNGR2', infot, nout, lerr, ok )
208 CALL
zungr2( 2, 1, 0, a, 2,
x, w, info )
209 CALL
chkxer(
'ZUNGR2', infot, nout, lerr, ok )
211 CALL
zungr2( 0, 0, -1, a, 1,
x, w, info )
212 CALL
chkxer(
'ZUNGR2', infot, nout, lerr, ok )
214 CALL
zungr2( 1, 2, 2, a, 2,
x, w, info )
215 CALL
chkxer(
'ZUNGR2', infot, nout, lerr, ok )
217 CALL
zungr2( 2, 2, 0, a, 1,
x, w, info )
218 CALL
chkxer(
'ZUNGR2', infot, nout, lerr, ok )
224 CALL
zunmrq(
'/',
'N', 0, 0, 0, a, 1,
x, af, 1, w, 1, info )
225 CALL
chkxer(
'ZUNMRQ', infot, nout, lerr, ok )
227 CALL
zunmrq(
'L',
'/', 0, 0, 0, a, 1,
x, af, 1, w, 1, info )
228 CALL
chkxer(
'ZUNMRQ', infot, nout, lerr, ok )
230 CALL
zunmrq(
'L',
'N', -1, 0, 0, a, 1,
x, af, 1, w, 1, info )
231 CALL
chkxer(
'ZUNMRQ', infot, nout, lerr, ok )
233 CALL
zunmrq(
'L',
'N', 0, -1, 0, a, 1,
x, af, 1, w, 1, info )
234 CALL
chkxer(
'ZUNMRQ', infot, nout, lerr, ok )
236 CALL
zunmrq(
'L',
'N', 0, 0, -1, a, 1,
x, af, 1, w, 1, info )
237 CALL
chkxer(
'ZUNMRQ', infot, nout, lerr, ok )
239 CALL
zunmrq(
'L',
'N', 0, 1, 1, a, 1,
x, af, 1, w, 1, info )
240 CALL
chkxer(
'ZUNMRQ', infot, nout, lerr, ok )
242 CALL
zunmrq(
'R',
'N', 1, 0, 1, a, 1,
x, af, 1, w, 1, info )
243 CALL
chkxer(
'ZUNMRQ', infot, nout, lerr, ok )
245 CALL
zunmrq(
'L',
'N', 2, 1, 2, a, 1,
x, af, 2, w, 1, info )
246 CALL
chkxer(
'ZUNMRQ', infot, nout, lerr, ok )
248 CALL
zunmrq(
'R',
'N', 1, 2, 2, a, 1,
x, af, 1, w, 1, info )
249 CALL
chkxer(
'ZUNMRQ', infot, nout, lerr, ok )
251 CALL
zunmrq(
'L',
'N', 2, 1, 0, a, 1,
x, af, 1, w, 1, info )
252 CALL
chkxer(
'ZUNMRQ', infot, nout, lerr, ok )
254 CALL
zunmrq(
'L',
'N', 1, 2, 0, a, 1,
x, af, 1, w, 1, info )
255 CALL
chkxer(
'ZUNMRQ', infot, nout, lerr, ok )
257 CALL
zunmrq(
'R',
'N', 2, 1, 0, a, 1,
x, af, 2, w, 1, info )
258 CALL
chkxer(
'ZUNMRQ', infot, nout, lerr, ok )
264 CALL
zunmr2(
'/',
'N', 0, 0, 0, a, 1,
x, af, 1, w, info )
265 CALL
chkxer(
'ZUNMR2', infot, nout, lerr, ok )
267 CALL
zunmr2(
'L',
'/', 0, 0, 0, a, 1,
x, af, 1, w, info )
268 CALL
chkxer(
'ZUNMR2', infot, nout, lerr, ok )
270 CALL
zunmr2(
'L',
'N', -1, 0, 0, a, 1,
x, af, 1, w, info )
271 CALL
chkxer(
'ZUNMR2', infot, nout, lerr, ok )
273 CALL
zunmr2(
'L',
'N', 0, -1, 0, a, 1,
x, af, 1, w, info )
274 CALL
chkxer(
'ZUNMR2', infot, nout, lerr, ok )
276 CALL
zunmr2(
'L',
'N', 0, 0, -1, a, 1,
x, af, 1, w, info )
277 CALL
chkxer(
'ZUNMR2', infot, nout, lerr, ok )
279 CALL
zunmr2(
'L',
'N', 0, 1, 1, a, 1,
x, af, 1, w, info )
280 CALL
chkxer(
'ZUNMR2', infot, nout, lerr, ok )
282 CALL
zunmr2(
'R',
'N', 1, 0, 1, a, 1,
x, af, 1, w, info )
283 CALL
chkxer(
'ZUNMR2', infot, nout, lerr, ok )
285 CALL
zunmr2(
'L',
'N', 2, 1, 2, a, 1,
x, af, 2, w, info )
286 CALL
chkxer(
'ZUNMR2', infot, nout, lerr, ok )
288 CALL
zunmr2(
'R',
'N', 1, 2, 2, a, 1,
x, af, 1, w, info )
289 CALL
chkxer(
'ZUNMR2', infot, nout, lerr, ok )
291 CALL
zunmr2(
'L',
'N', 2, 1, 0, a, 1,
x, af, 1, w, info )
292 CALL
chkxer(
'ZUNMR2', infot, nout, lerr, ok )
296 CALL
alaesm( path, ok, nout )
subroutine zgerqf(M, N, A, LDA, TAU, WORK, LWORK, INFO)
ZGERQF
subroutine zerrrq(PATH, NUNIT)
ZERRRQ
subroutine zunmrq(SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, WORK, LWORK, INFO)
ZUNMRQ
subroutine zungr2(M, N, K, A, LDA, TAU, WORK, INFO)
ZUNGR2 generates all or part of the unitary matrix Q from an RQ factorization determined by cgerqf (u...
subroutine zunmr2(SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, WORK, INFO)
ZUNMR2 multiplies a general matrix by the unitary matrix from a RQ factorization determined by cgerqf...
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
subroutine zgerqs(M, N, NRHS, A, LDA, TAU, B, LDB, WORK, LWORK, INFO)
ZGERQS
subroutine zungrq(M, N, K, A, LDA, TAU, WORK, LWORK, INFO)
ZUNGRQ
subroutine zgerq2(M, N, A, LDA, TAU, WORK, INFO)
ZGERQ2 computes the RQ factorization of a general rectangular matrix using an unblocked algorithm...
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
subroutine alaesm(PATH, OK, NOUT)
ALAESM