78 COMPLEX a( nmax, nmax ), af( nmax, nmax ),
b( nmax ),
79 $ w( nmax ),
x( nmax )
91 COMMON / infoc / infot, nout, ok, lerr
92 COMMON / srnamc / srnamt
100 WRITE( nout, fmt = * )
106 a( i,
j ) = cmplx( 1. /
REAL( I+J ), -1. /
REAL( I+J ) )
107 af( i,
j ) = cmplx( 1. /
REAL( I+J ), -1. /
REAL( I+J ) )
121 CALL
cgerqf( -1, 0, a, 1,
b, w, 1, info )
122 CALL
chkxer(
'CGERQF', infot, nout, lerr, ok )
124 CALL
cgerqf( 0, -1, a, 1,
b, w, 1, info )
125 CALL
chkxer(
'CGERQF', infot, nout, lerr, ok )
127 CALL
cgerqf( 2, 1, a, 1,
b, w, 2, info )
128 CALL
chkxer(
'CGERQF', infot, nout, lerr, ok )
130 CALL
cgerqf( 2, 1, a, 2,
b, w, 1, info )
131 CALL
chkxer(
'CGERQF', infot, nout, lerr, ok )
137 CALL
cgerq2( -1, 0, a, 1,
b, w, info )
138 CALL
chkxer(
'CGERQ2', infot, nout, lerr, ok )
140 CALL
cgerq2( 0, -1, a, 1,
b, w, info )
141 CALL
chkxer(
'CGERQ2', infot, nout, lerr, ok )
143 CALL
cgerq2( 2, 1, a, 1,
b, w, info )
144 CALL
chkxer(
'CGERQ2', infot, nout, lerr, ok )
150 CALL
cgerqs( -1, 0, 0, a, 1,
x,
b, 1, w, 1, info )
151 CALL
chkxer(
'CGERQS', infot, nout, lerr, ok )
153 CALL
cgerqs( 0, -1, 0, a, 1,
x,
b, 1, w, 1, info )
154 CALL
chkxer(
'CGERQS', infot, nout, lerr, ok )
156 CALL
cgerqs( 2, 1, 0, a, 2,
x,
b, 1, w, 1, info )
157 CALL
chkxer(
'CGERQS', infot, nout, lerr, ok )
159 CALL
cgerqs( 0, 0, -1, a, 1,
x,
b, 1, w, 1, info )
160 CALL
chkxer(
'CGERQS', infot, nout, lerr, ok )
162 CALL
cgerqs( 2, 2, 0, a, 1,
x,
b, 2, w, 1, info )
163 CALL
chkxer(
'CGERQS', infot, nout, lerr, ok )
165 CALL
cgerqs( 2, 2, 0, a, 2,
x,
b, 1, w, 1, info )
166 CALL
chkxer(
'CGERQS', infot, nout, lerr, ok )
168 CALL
cgerqs( 1, 1, 2, a, 1,
x,
b, 1, w, 1, info )
169 CALL
chkxer(
'CGERQS', infot, nout, lerr, ok )
175 CALL
cungrq( -1, 0, 0, a, 1,
x, w, 1, info )
176 CALL
chkxer(
'CUNGRQ', infot, nout, lerr, ok )
178 CALL
cungrq( 0, -1, 0, a, 1,
x, w, 1, info )
179 CALL
chkxer(
'CUNGRQ', infot, nout, lerr, ok )
181 CALL
cungrq( 2, 1, 0, a, 2,
x, w, 2, info )
182 CALL
chkxer(
'CUNGRQ', infot, nout, lerr, ok )
184 CALL
cungrq( 0, 0, -1, a, 1,
x, w, 1, info )
185 CALL
chkxer(
'CUNGRQ', infot, nout, lerr, ok )
187 CALL
cungrq( 1, 2, 2, a, 1,
x, w, 1, info )
188 CALL
chkxer(
'CUNGRQ', infot, nout, lerr, ok )
190 CALL
cungrq( 2, 2, 0, a, 1,
x, w, 2, info )
191 CALL
chkxer(
'CUNGRQ', infot, nout, lerr, ok )
193 CALL
cungrq( 2, 2, 0, a, 2,
x, w, 1, info )
194 CALL
chkxer(
'CUNGRQ', infot, nout, lerr, ok )
200 CALL
cungr2( -1, 0, 0, a, 1,
x, w, info )
201 CALL
chkxer(
'CUNGR2', infot, nout, lerr, ok )
203 CALL
cungr2( 0, -1, 0, a, 1,
x, w, info )
204 CALL
chkxer(
'CUNGR2', infot, nout, lerr, ok )
206 CALL
cungr2( 2, 1, 0, a, 2,
x, w, info )
207 CALL
chkxer(
'CUNGR2', infot, nout, lerr, ok )
209 CALL
cungr2( 0, 0, -1, a, 1,
x, w, info )
210 CALL
chkxer(
'CUNGR2', infot, nout, lerr, ok )
212 CALL
cungr2( 1, 2, 2, a, 2,
x, w, info )
213 CALL
chkxer(
'CUNGR2', infot, nout, lerr, ok )
215 CALL
cungr2( 2, 2, 0, a, 1,
x, w, info )
216 CALL
chkxer(
'CUNGR2', infot, nout, lerr, ok )
222 CALL
cunmrq(
'/',
'N', 0, 0, 0, a, 1,
x, af, 1, w, 1, info )
223 CALL
chkxer(
'CUNMRQ', infot, nout, lerr, ok )
225 CALL
cunmrq(
'L',
'/', 0, 0, 0, a, 1,
x, af, 1, w, 1, info )
226 CALL
chkxer(
'CUNMRQ', infot, nout, lerr, ok )
228 CALL
cunmrq(
'L',
'N', -1, 0, 0, a, 1,
x, af, 1, w, 1, info )
229 CALL
chkxer(
'CUNMRQ', infot, nout, lerr, ok )
231 CALL
cunmrq(
'L',
'N', 0, -1, 0, a, 1,
x, af, 1, w, 1, info )
232 CALL
chkxer(
'CUNMRQ', infot, nout, lerr, ok )
234 CALL
cunmrq(
'L',
'N', 0, 0, -1, a, 1,
x, af, 1, w, 1, info )
235 CALL
chkxer(
'CUNMRQ', infot, nout, lerr, ok )
237 CALL
cunmrq(
'L',
'N', 0, 1, 1, a, 1,
x, af, 1, w, 1, info )
238 CALL
chkxer(
'CUNMRQ', infot, nout, lerr, ok )
240 CALL
cunmrq(
'R',
'N', 1, 0, 1, a, 1,
x, af, 1, w, 1, info )
241 CALL
chkxer(
'CUNMRQ', infot, nout, lerr, ok )
243 CALL
cunmrq(
'L',
'N', 2, 1, 2, a, 1,
x, af, 2, w, 1, info )
244 CALL
chkxer(
'CUNMRQ', infot, nout, lerr, ok )
246 CALL
cunmrq(
'R',
'N', 1, 2, 2, a, 1,
x, af, 1, w, 1, info )
247 CALL
chkxer(
'CUNMRQ', infot, nout, lerr, ok )
249 CALL
cunmrq(
'L',
'N', 2, 1, 0, a, 1,
x, af, 1, w, 1, info )
250 CALL
chkxer(
'CUNMRQ', infot, nout, lerr, ok )
252 CALL
cunmrq(
'L',
'N', 1, 2, 0, a, 1,
x, af, 1, w, 1, info )
253 CALL
chkxer(
'CUNMRQ', infot, nout, lerr, ok )
255 CALL
cunmrq(
'R',
'N', 2, 1, 0, a, 1,
x, af, 2, w, 1, info )
256 CALL
chkxer(
'CUNMRQ', infot, nout, lerr, ok )
262 CALL
cunmr2(
'/',
'N', 0, 0, 0, a, 1,
x, af, 1, w, info )
263 CALL
chkxer(
'CUNMR2', infot, nout, lerr, ok )
265 CALL
cunmr2(
'L',
'/', 0, 0, 0, a, 1,
x, af, 1, w, info )
266 CALL
chkxer(
'CUNMR2', infot, nout, lerr, ok )
268 CALL
cunmr2(
'L',
'N', -1, 0, 0, a, 1,
x, af, 1, w, info )
269 CALL
chkxer(
'CUNMR2', infot, nout, lerr, ok )
271 CALL
cunmr2(
'L',
'N', 0, -1, 0, a, 1,
x, af, 1, w, info )
272 CALL
chkxer(
'CUNMR2', infot, nout, lerr, ok )
274 CALL
cunmr2(
'L',
'N', 0, 0, -1, a, 1,
x, af, 1, w, info )
275 CALL
chkxer(
'CUNMR2', infot, nout, lerr, ok )
277 CALL
cunmr2(
'L',
'N', 0, 1, 1, a, 1,
x, af, 1, w, info )
278 CALL
chkxer(
'CUNMR2', infot, nout, lerr, ok )
280 CALL
cunmr2(
'R',
'N', 1, 0, 1, a, 1,
x, af, 1, w, info )
281 CALL
chkxer(
'CUNMR2', infot, nout, lerr, ok )
283 CALL
cunmr2(
'L',
'N', 2, 1, 2, a, 1,
x, af, 2, w, info )
284 CALL
chkxer(
'CUNMR2', infot, nout, lerr, ok )
286 CALL
cunmr2(
'R',
'N', 1, 2, 2, a, 1,
x, af, 1, w, info )
287 CALL
chkxer(
'CUNMR2', infot, nout, lerr, ok )
289 CALL
cunmr2(
'L',
'N', 2, 1, 0, a, 1,
x, af, 1, w, info )
290 CALL
chkxer(
'CUNMR2', infot, nout, lerr, ok )
294 CALL
alaesm( path, ok, nout )
subroutine cunmrq(SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, WORK, LWORK, INFO)
CUNMRQ
subroutine cunmr2(SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, WORK, INFO)
CUNMR2 multiplies a general matrix by the unitary matrix from a RQ factorization determined by cgerqf...
subroutine cgerq2(M, N, A, LDA, TAU, WORK, INFO)
CGERQ2 computes the RQ factorization of a general rectangular matrix using an unblocked algorithm...
subroutine cerrrq(PATH, NUNIT)
CERRRQ
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine cungr2(M, N, K, A, LDA, TAU, WORK, INFO)
CUNGR2 generates all or part of the unitary matrix Q from an RQ factorization determined by cgerqf (u...
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
subroutine cungrq(M, N, K, A, LDA, TAU, WORK, LWORK, INFO)
CUNGRQ
subroutine cgerqf(M, N, A, LDA, TAU, WORK, LWORK, INFO)
CGERQF
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
subroutine cgerqs(M, N, NRHS, A, LDA, TAU, B, LDB, WORK, LWORK, INFO)
CGERQS
subroutine alaesm(PATH, OK, NOUT)
ALAESM