OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
zerrgg.f
Go to the documentation of this file.
1*> \brief \b ZERRGG
2*
3* =========== DOCUMENTATION ===========
4*
5* Online html documentation available at
6* http://www.netlib.org/lapack/explore-html/
7*
8* Definition:
9* ===========
10*
11* SUBROUTINE ZERRGG( PATH, NUNIT )
12*
13* .. Scalar Arguments ..
14* CHARACTER*3 PATH
15* INTEGER NUNIT
16* ..
17*
18*
19*> \par Purpose:
20* =============
21*>
22*> \verbatim
23*>
24*> ZERRGG tests the error exits for ZGGES, ZGGESX, ZGGEV, ZGGEVX,
25*> ZGGES3, ZGGEV3, ZGGGLM, ZGGHRD, ZGGLSE, ZGGQRF, ZGGRQF,
26*> ZGGSVD3, ZGGSVP3, ZHGEQZ, ZTGEVC, ZTGEXC, ZTGSEN, ZTGSJA,
27*> ZTGSNA, ZTGSYL, and ZUNCSD.
28*> \endverbatim
29*
30* Arguments:
31* ==========
32*
33*> \param[in] PATH
34*> \verbatim
35*> PATH is CHARACTER*3
36*> The LAPACK path name for the routines to be tested.
37*> \endverbatim
38*>
39*> \param[in] NUNIT
40*> \verbatim
41*> NUNIT is INTEGER
42*> The unit number for output.
43*> \endverbatim
44*
45* Authors:
46* ========
47*
48*> \author Univ. of Tennessee
49*> \author Univ. of California Berkeley
50*> \author Univ. of Colorado Denver
51*> \author NAG Ltd.
52*
53*> \ingroup complex16_eig
54*
55* =====================================================================
56 SUBROUTINE zerrgg( PATH, NUNIT )
57*
58* -- LAPACK test routine --
59* -- LAPACK is a software package provided by Univ. of Tennessee, --
60* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
61*
62* .. Scalar Arguments ..
63 CHARACTER*3 PATH
64 INTEGER NUNIT
65* ..
66*
67* =====================================================================
68*
69* .. Parameters ..
70 INTEGER NMAX, LW
71 parameter( nmax = 3, lw = 6*nmax )
72 DOUBLE PRECISION ONE, ZERO
73 parameter( one = 1.0d+0, zero = 0.0d+0 )
74* ..
75* .. Local Scalars ..
76 CHARACTER*2 C2
77 INTEGER DUMMYK, DUMMYL, I, IFST, IHI, ILO, ILST, INFO,
78 $ J, M, NCYCLE, NT, SDIM, LWORK
79 DOUBLE PRECISION ANRM, BNRM, DIF, SCALE, TOLA, TOLB
80* ..
81* .. Local Arrays ..
82 LOGICAL BW( NMAX ), SEL( NMAX )
83 INTEGER IW( LW ), IDUM(NMAX)
84 DOUBLE PRECISION LS( NMAX ), R1( NMAX ), R2( NMAX ),
85 $ RCE( NMAX ), RCV( NMAX ), RS( NMAX ), RW( LW )
86 COMPLEX*16 A( NMAX, NMAX ), ALPHA( NMAX ),
87 $ B( NMAX, NMAX ), BETA( NMAX ), Q( NMAX, NMAX ),
88 $ TAU( NMAX ), U( NMAX, NMAX ), V( NMAX, NMAX ),
89 $ W( LW ), Z( NMAX, NMAX )
90* ..
91* .. External Functions ..
92 LOGICAL LSAMEN, ZLCTES, ZLCTSX
93 EXTERNAL lsamen, zlctes, zlctsx
94* ..
95* .. External Subroutines ..
96 EXTERNAL chkxer, zgges, zggesx, zggev, zggevx, zggglm,
97 $ zgghrd, zgglse, zggqrf, zggrqf,
98 $ zhgeqz, ztgevc, ztgexc, ztgsen, ztgsja, ztgsna,
99 $ ztgsyl, zuncsd, zgges3, zggev3, zgghd3,
100 $ zggsvd3, zggsvp3, xlaenv
101* ..
102* .. Scalars in Common ..
103 LOGICAL LERR, OK
104 CHARACTER*32 SRNAMT
105 INTEGER INFOT, NOUT
106* ..
107* .. Common blocks ..
108 COMMON / infoc / infot, nout, ok, lerr
109 COMMON / srnamc / srnamt
110* ..
111* .. Executable Statements ..
112*
113 nout = nunit
114 WRITE( nout, fmt = * )
115 c2 = path( 2: 3 )
116*
117* Set the variables to innocuous values.
118*
119 DO 20 j = 1, nmax
120 sel( j ) = .true.
121 DO 10 i = 1, nmax
122 a( i, j ) = zero
123 b( i, j ) = zero
124 10 CONTINUE
125 20 CONTINUE
126 DO 30 i = 1, nmax
127 a( i, i ) = one
128 b( i, i ) = one
129 30 CONTINUE
130 ok = .true.
131 tola = 1.0d0
132 tolb = 1.0d0
133 ifst = 1
134 ilst = 1
135 nt = 0
136 lwork = 1
137*
138* Call XLAENV to set the parameters used in CLAQZ0
139*
140 CALL xlaenv( 12, 10 )
141 CALL xlaenv( 13, 12 )
142 CALL xlaenv( 14, 13 )
143 CALL xlaenv( 15, 2 )
144 CALL xlaenv( 17, 10 )
145*
146* Test error exits for the GG path.
147*
148 IF( lsamen( 2, c2, 'GG' ) ) THEN
149*
150* ZGGHRD
151*
152 srnamt = 'ZGGHRD'
153 infot = 1
154 CALL zgghrd( '/', 'N', 0, 1, 0, a, 1, b, 1, q, 1, z, 1, info )
155 CALL chkxer( 'ZGGHRD', infot, nout, lerr, ok )
156 infot = 2
157 CALL zgghrd( 'N', '/', 0, 1, 0, a, 1, b, 1, q, 1, z, 1, info )
158 CALL chkxer( 'ZGGHRD', infot, nout, lerr, ok )
159 infot = 3
160 CALL zgghrd( 'N', 'N', -1, 0, 0, a, 1, b, 1, q, 1, z, 1, info )
161 CALL chkxer( 'ZGGHRD', infot, nout, lerr, ok )
162 infot = 4
163 CALL zgghrd( 'N', 'N', 0, 0, 0, a, 1, b, 1, q, 1, z, 1, info )
164 CALL chkxer( 'zgghrd', INFOT, NOUT, LERR, OK )
165 INFOT = 5
166 CALL ZGGHRD( 'n', 'n', 0, 1, 1, a, 1, b, 1, q, 1, z, 1, info )
167 CALL chkxer( 'ZGGHRD', infot, nout, lerr, ok )
168 infot = 7
169 CALL zgghrd( 'N', 'N', 2, 1, 1, a, 1, b, 2, q, 1, z, 1, info )
170 CALL chkxer( 'ZGGHRD', infot, nout, lerr, ok )
171 infot = 9
172 CALL zgghrd( 'N', 'N', 2, 1, 1, a, 2, b, 1, q, 1, z, 1, info )
173 CALL chkxer( 'ZGGHRD', infot, nout, lerr, ok )
174 infot = 11
175 CALL zgghrd( 'V', 'N', 2, 1, 1, a, 2, b, 2, q, 1, z, 1, info )
176 CALL chkxer( 'ZGGHRD', infot, nout, lerr, ok )
177 infot = 13
178 CALL zgghrd( 'N', 'V', 2, 1, 1, a, 2, b, 2, q, 1, z, 1, info )
179 CALL chkxer( 'ZGGHRD', infot, nout, lerr, ok )
180 nt = nt + 9
181*
182* ZGGHD3
183*
184 srnamt = 'ZGGHD3'
185 infot = 1
186 CALL zgghd3( '/', 'N', 0, 1, 0, a, 1, b, 1, q, 1, z, 1, w, lw,
187 $ info )
188 CALL chkxer( 'ZGGHD3', infot, nout, lerr, ok )
189 infot = 2
190 CALL zgghd3( 'N', '/', 0, 1, 0, a, 1, b, 1, q, 1, z, 1, w, lw,
191 $ info )
192 CALL chkxer( 'ZGGHD3', infot, nout, lerr, ok )
193 infot = 3
194 CALL zgghd3( 'N', 'N', -1, 0, 0, a, 1, b, 1, q, 1, z, 1, w, lw,
195 $ info )
196 CALL chkxer( 'ZGGHD3', infot, nout, lerr, ok )
197 infot = 4
198 CALL zgghd3( 'N', 'N', 0, 0, 0, a, 1, b, 1, q, 1, z, 1, w, lw,
199 $ info )
200 CALL chkxer( 'ZGGHD3', infot, nout, lerr, ok )
201 infot = 5
202 CALL zgghd3( 'N', 'N', 0, 1, 1, a, 1, b, 1, q, 1, z, 1, w, lw,
203 $ info )
204 CALL chkxer( 'ZGGHD3', infot, nout, lerr, ok )
205 infot = 7
206 CALL zgghd3( 'N', 'N', 2, 1, 1, a, 1, b, 2, q, 1, z, 1, w, lw,
207 $ info )
208 CALL chkxer( 'ZGGHD3', infot, nout, lerr, ok )
209 infot = 9
210 CALL zgghd3( 'N', 'N', 2, 1, 1, a, 2, b, 1, q, 1, z, 1, w, lw,
211 $ info )
212 CALL chkxer( 'ZGGHD3', infot, nout, lerr, ok )
213 infot = 11
214 CALL zgghd3( 'V', 'N', 2, 1, 1, a, 2, b, 2, q, 1, z, 1, w, lw,
215 $ info )
216 CALL chkxer( 'ZGGHD3', infot, nout, lerr, ok )
217 infot = 13
218 CALL zgghd3( 'N', 'V', 2, 1, 1, a, 2, b, 2, q, 1, z, 1, w, lw,
219 $ info )
220 CALL chkxer( 'ZGGHD3', infot, nout, lerr, ok )
221 nt = nt + 9
222*
223* ZHGEQZ
224*
225 srnamt = 'ZHGEQZ'
226 infot = 1
227 CALL zhgeqz( '/', 'N', 'N', 0, 1, 0, a, 1, b, 1, alpha, beta,
228 $ q, 1, z, 1, w, 1, rw, info )
229 CALL chkxer( 'ZHGEQZ', infot, nout, lerr, ok )
230 infot = 2
231 CALL zhgeqz( 'E', '/', 'N', 0, 1, 0, a, 1, b, 1, alpha, beta,
232 $ q, 1, z, 1, w, 1, rw, info )
233 CALL chkxer( 'ZHGEQZ', infot, nout, lerr, ok )
234 infot = 3
235 CALL zhgeqz( 'E', 'N', '/', 0, 1, 0, a, 1, b, 1, alpha, beta,
236 $ q, 1, z, 1, w, 1, rw, info )
237 CALL chkxer( 'ZHGEQZ', infot, nout, lerr, ok )
238 infot = 4
239 CALL zhgeqz( 'e', 'n', 'n', -1, 0, 0, A, 1, B, 1, ALPHA, BETA,
240 $ Q, 1, Z, 1, W, 1, RW, INFO )
241 CALL CHKXER( 'zhgeqz', INFOT, NOUT, LERR, OK )
242 INFOT = 5
243 CALL ZHGEQZ( 'e', 'n', 'n', 0, 0, 0, A, 1, B, 1, ALPHA, BETA,
244 $ Q, 1, Z, 1, W, 1, RW, INFO )
245 CALL CHKXER( 'zhgeqz', INFOT, NOUT, LERR, OK )
246 INFOT = 6
247 CALL ZHGEQZ( 'e', 'n', 'n', 0, 1, 1, A, 1, B, 1, ALPHA, BETA,
248 $ Q, 1, Z, 1, W, 1, RW, INFO )
249 CALL CHKXER( 'zhgeqz', INFOT, NOUT, LERR, OK )
250 INFOT = 8
251 CALL ZHGEQZ( 'e', 'n', 'n', 2, 1, 1, A, 1, B, 2, ALPHA, BETA,
252 $ Q, 1, Z, 1, W, 1, RW, INFO )
253 CALL CHKXER( 'zhgeqz', INFOT, NOUT, LERR, OK )
254 INFOT = 10
255 CALL ZHGEQZ( 'e', 'n', 'n', 2, 1, 1, A, 2, B, 1, ALPHA, BETA,
256 $ Q, 1, Z, 1, W, 1, RW, INFO )
257 CALL CHKXER( 'zhgeqz', INFOT, NOUT, LERR, OK )
258 INFOT = 14
259 CALL ZHGEQZ( 'e', 'v', 'n', 2, 1, 1, A, 2, B, 2, ALPHA, BETA,
260 $ Q, 1, Z, 1, W, 1, RW, INFO )
261 CALL CHKXER( 'zhgeqz', INFOT, NOUT, LERR, OK )
262 INFOT = 16
263 CALL ZHGEQZ( 'e', 'n', 'v', 2, 1, 1, A, 2, B, 2, ALPHA, BETA,
264 $ Q, 1, Z, 1, W, 1, RW, INFO )
265 CALL CHKXER( 'zhgeqz', INFOT, NOUT, LERR, OK )
266 NT = NT + 10
267*
268* ZTGEVC
269*
270 SRNAMT = 'ztgevc'
271 INFOT = 1
272 CALL ZTGEVC( '/', 'a', SEL, 0, A, 1, B, 1, Q, 1, Z, 1, 0, M, W,
273 $ RW, INFO )
274 CALL CHKXER( 'ztgevc', INFOT, NOUT, LERR, OK )
275 INFOT = 2
276 CALL ZTGEVC( 'r', '/', SEL, 0, A, 1, B, 1, Q, 1, Z, 1, 0, M, W,
277 $ RW, INFO )
278 CALL CHKXER( 'ztgevc', INFOT, NOUT, LERR, OK )
279 INFOT = 4
280 CALL ZTGEVC( 'r', 'a', SEL, -1, A, 1, B, 1, Q, 1, Z, 1, 0, M,
281 $ W, RW, INFO )
282 CALL CHKXER( 'ztgevc', INFOT, NOUT, LERR, OK )
283 INFOT = 6
284 CALL ZTGEVC( 'r', 'a', SEL, 2, A, 1, B, 2, Q, 1, Z, 2, 0, M, W,
285 $ RW, INFO )
286 CALL CHKXER( 'ztgevc', INFOT, NOUT, LERR, OK )
287 INFOT = 8
288 CALL ZTGEVC( 'r', 'a', SEL, 2, A, 2, B, 1, Q, 1, Z, 2, 0, M, W,
289 $ RW, INFO )
290 CALL CHKXER( 'ztgevc', INFOT, NOUT, LERR, OK )
291 INFOT = 10
292 CALL ZTGEVC( 'l', 'a', SEL, 2, A, 2, B, 2, Q, 1, Z, 1, 0, M, W,
293 $ RW, INFO )
294 CALL CHKXER( 'ztgevc', INFOT, NOUT, LERR, OK )
295 INFOT = 12
296 CALL ZTGEVC( 'r', 'a', SEL, 2, A, 2, B, 2, Q, 1, Z, 1, 0, M, W,
297 $ RW, INFO )
298 CALL CHKXER( 'ztgevc', INFOT, NOUT, LERR, OK )
299 INFOT = 13
300 CALL ZTGEVC( 'r', 'a', SEL, 2, A, 2, B, 2, Q, 1, Z, 2, 1, M, W,
301 $ RW, INFO )
302 CALL CHKXER( 'ztgevc', INFOT, NOUT, LERR, OK )
303 NT = NT + 8
304*
305* Test error exits for the GSV path.
306*
307 ELSE IF( LSAMEN( 3, PATH, 'gsv' ) ) THEN
308*
309* ZGGSVD3
310*
311 SRNAMT = 'zggsvd3'
312 INFOT = 1
313 CALL ZGGSVD3( '/', 'n', 'n', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B,
314 $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM,
315 $ INFO )
316 CALL CHKXER( 'zggsvd3', INFOT, NOUT, LERR, OK )
317 INFOT = 2
318 CALL ZGGSVD3( 'n', '/', 'n', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B,
319 $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM,
320 $ INFO )
321 CALL CHKXER( 'zggsvd3', INFOT, NOUT, LERR, OK )
322 INFOT = 3
323 CALL ZGGSVD3( 'n', 'n', '/', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B,
324 $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM,
325 $ INFO )
326 CALL CHKXER( 'zggsvd3', INFOT, NOUT, LERR, OK )
327 INFOT = 4
328 CALL ZGGSVD3( 'n', 'n', 'n', -1, 0, 0, DUMMYK, DUMMYL, A, 1, B,
329 $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM,
330 $ INFO )
331 CALL CHKXER( 'zggsvd3', INFOT, NOUT, LERR, OK )
332 INFOT = 5
333 CALL ZGGSVD3( 'n', 'n', 'n', 0, -1, 0, DUMMYK, DUMMYL, A, 1, B,
334 $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM,
335 $ INFO )
336 CALL CHKXER( 'zggsvd3', INFOT, NOUT, LERR, OK )
337 INFOT = 6
338 CALL ZGGSVD3( 'n', 'n', 'n', 0, 0, -1, DUMMYK, DUMMYL, A, 1, B,
339 $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM,
340 $ INFO )
341 CALL CHKXER( 'zggsvd3', INFOT, NOUT, LERR, OK )
342 INFOT = 10
343 CALL ZGGSVD3( 'n', 'n', 'n', 2, 1, 1, DUMMYK, DUMMYL, A, 1, B,
344 $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM,
345 $ INFO )
346 CALL CHKXER( 'zggsvd3', INFOT, NOUT, LERR, OK )
347 INFOT = 12
348 CALL ZGGSVD3( 'n', 'n', 'n', 1, 1, 2, DUMMYK, DUMMYL, A, 1, B,
349 $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM,
350 $ INFO )
351 CALL CHKXER( 'zggsvd3', INFOT, NOUT, LERR, OK )
352 INFOT = 16
353 CALL ZGGSVD3( 'u', 'n', 'n', 2, 2, 2, DUMMYK, DUMMYL, A, 2, B,
354 $ 2, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM,
355 $ INFO )
356 CALL CHKXER( 'zggsvd3', INFOT, NOUT, LERR, OK )
357 INFOT = 18
358 CALL ZGGSVD3( 'n', 'v', 'n', 2, 2, 2, DUMMYK, DUMMYL, A, 2, B,
359 $ 2, R1, R2, U, 2, V, 1, Q, 1, W, LWORK, RW, IDUM,
360 $ INFO )
361 CALL CHKXER( 'zggsvd3', INFOT, NOUT, LERR, OK )
362 INFOT = 20
363 CALL ZGGSVD3( 'n', 'n', 'q', 2, 2, 2, DUMMYK, DUMMYL, A, 2, B,
364 $ 2, R1, R2, U, 2, V, 2, Q, 1, W, LWORK, RW, IDUM,
365 $ INFO )
366 CALL CHKXER( 'zggsvd3', INFOT, NOUT, LERR, OK )
367 NT = NT + 11
368*
369* ZGGSVP3
370*
371 SRNAMT = 'zggsvp3'
372 INFOT = 1
373 CALL ZGGSVP3( '/', 'n', 'n', 0, 0, 0, A, 1, B, 1, TOLA, TOLB,
374 $ DUMMYK, DUMMYL, U, 1, V, 1, Q, 1, IW, RW, TAU, W,
375 $ LWORK, INFO )
376 CALL CHKXER( 'zggsvp3', INFOT, NOUT, LERR, OK )
377 INFOT = 2
378 CALL ZGGSVP3( 'n', '/', 'n', 0, 0, 0, A, 1, B, 1, TOLA, TOLB,
379 $ DUMMYK, DUMMYL, U, 1, V, 1, Q, 1, IW, RW, TAU, W,
380 $ LWORK, INFO )
381 CALL CHKXER( 'zggsvp3', INFOT, NOUT, LERR, OK )
382 INFOT = 3
383 CALL ZGGSVP3( 'n', 'n', '/', 0, 0, 0, A, 1, B, 1, TOLA, TOLB,
384 $ DUMMYK, DUMMYL, U, 1, V, 1, Q, 1, IW, RW, TAU, W,
385 $ LWORK, INFO )
386 CALL CHKXER( 'zggsvp3', INFOT, NOUT, LERR, OK )
387 INFOT = 4
388 CALL ZGGSVP3( 'n', 'n', 'n', -1, 0, 0, A, 1, B, 1, TOLA, TOLB,
389 $ DUMMYK, DUMMYL, U, 1, V, 1, Q, 1, IW, RW, TAU, W,
390 $ LWORK, INFO )
391 CALL CHKXER( 'zggsvp3', INFOT, NOUT, LERR, OK )
392 INFOT = 5
393 CALL ZGGSVP3( 'n', 'n', 'n', 0, -1, 0, A, 1, B, 1, TOLA, TOLB,
394 $ DUMMYK, DUMMYL, U, 1, V, 1, Q, 1, IW, RW, TAU, W,
395 $ LWORK, INFO )
396 CALL CHKXER( 'zggsvp3', INFOT, NOUT, LERR, OK )
397 INFOT = 6
398 CALL ZGGSVP3( 'n', 'n', 'n', 0, 0, -1, A, 1, B, 1, TOLA, TOLB,
399 $ DUMMYK, DUMMYL, U, 1, V, 1, Q, 1, IW, RW, TAU, W,
400 $ LWORK, INFO )
401 CALL CHKXER( 'zggsvp3', INFOT, NOUT, LERR, OK )
402 INFOT = 8
403 CALL ZGGSVP3( 'n', 'n', 'n', 2, 1, 1, A, 1, B, 1, TOLA, TOLB,
404 $ DUMMYK, DUMMYL, U, 1, V, 1, Q, 1, IW, RW, TAU, W,
405 $ LWORK, INFO )
406 CALL CHKXER( 'zggsvp3', INFOT, NOUT, LERR, OK )
407 INFOT = 10
408 CALL ZGGSVP3( 'n', 'n', 'n', 1, 2, 1, A, 1, B, 1, TOLA, TOLB,
409 $ DUMMYK, DUMMYL, U, 1, V, 1, Q, 1, IW, RW, TAU, W,
410 $ LWORK, INFO )
411 CALL CHKXER( 'zggsvp3', INFOT, NOUT, LERR, OK )
412 INFOT = 16
413 CALL ZGGSVP3( 'u', 'n', 'n', 2, 2, 2, A, 2, B, 2, TOLA, TOLB,
414 $ DUMMYK, DUMMYL, U, 1, V, 1, Q, 1, IW, RW, TAU, W,
415 $ LWORK, INFO )
416 CALL CHKXER( 'zggsvp3', INFOT, NOUT, LERR, OK )
417 INFOT = 18
418 CALL ZGGSVP3( 'n', 'v', 'n', 2, 2, 2, A, 2, B, 2, TOLA, TOLB,
419 $ DUMMYK, DUMMYL, U, 2, V, 1, Q, 1, IW, RW, TAU, W,
420 $ LWORK, INFO )
421 CALL CHKXER( 'zggsvp3', INFOT, NOUT, LERR, OK )
422 INFOT = 20
423 CALL ZGGSVP3( 'n', 'n', 'q', 2, 2, 2, A, 2, B, 2, TOLA, TOLB,
424 $ DUMMYK, DUMMYL, U, 2, V, 2, Q, 1, IW, RW, TAU, W,
425 $ LWORK, INFO )
426 CALL CHKXER( 'zggsvp3', INFOT, NOUT, LERR, OK )
427 NT = NT + 11
428*
429* ZTGSJA
430*
431 SRNAMT = 'ztgsja'
432 INFOT = 1
433 CALL ZTGSJA( '/', 'n', 'n', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B,
434 $ 1, TOLA, TOLB, R1, R2, U, 1, V, 1, Q, 1, W,
435 $ NCYCLE, INFO )
436 CALL CHKXER( 'ztgsja', INFOT, NOUT, LERR, OK )
437 INFOT = 2
438 CALL ZTGSJA( 'n', '/', 'n', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B,
439 $ 1, TOLA, TOLB, R1, R2, U, 1, V, 1, Q, 1, W,
440 $ NCYCLE, INFO )
441 CALL CHKXER( 'ztgsja', INFOT, NOUT, LERR, OK )
442 INFOT = 3
443 CALL ZTGSJA( 'n', 'n', '/', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B,
444 $ 1, TOLA, TOLB, R1, R2, U, 1, V, 1, Q, 1, W,
445 $ NCYCLE, INFO )
446 CALL CHKXER( 'ztgsja', INFOT, NOUT, LERR, OK )
447 INFOT = 4
448 CALL ZTGSJA( 'n', 'n', 'n', -1, 0, 0, DUMMYK, DUMMYL, A, 1, B,
449 $ 1, TOLA, TOLB, R1, R2, U, 1, V, 1, Q, 1, W,
450 $ NCYCLE, INFO )
451 CALL CHKXER( 'ztgsja', INFOT, NOUT, LERR, OK )
452 INFOT = 5
453 CALL ZTGSJA( 'n', 'n', 'n', 0, -1, 0, DUMMYK, DUMMYL, A, 1, B,
454 $ 1, TOLA, TOLB, R1, R2, U, 1, V, 1, Q, 1, W,
455 $ NCYCLE, INFO )
456 CALL CHKXER( 'ztgsja', INFOT, NOUT, LERR, OK )
457 INFOT = 6
458 CALL ZTGSJA( 'n', 'n', 'n', 0, 0, -1, DUMMYK, DUMMYL, A, 1, B,
459 $ 1, TOLA, TOLB, R1, R2, U, 1, V, 1, Q, 1, W,
460 $ NCYCLE, INFO )
461 CALL CHKXER( 'ztgsja', INFOT, NOUT, LERR, OK )
462 INFOT = 10
463 CALL ZTGSJA( 'n', 'n', 'n', 0, 0, 0, DUMMYK, DUMMYL, A, 0, B,
464 $ 1, TOLA, TOLB, R1, R2, U, 1, V, 1, Q, 1, W,
465 $ NCYCLE, INFO )
466 CALL CHKXER( 'ztgsja', INFOT, NOUT, LERR, OK )
467 INFOT = 12
468 CALL ZTGSJA( 'n', 'n', 'n', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B,
469 $ 0, TOLA, TOLB, R1, R2, U, 1, V, 1, Q, 1, W,
470 $ NCYCLE, INFO )
471 CALL CHKXER( 'ztgsja', INFOT, NOUT, LERR, OK )
472 INFOT = 18
473 CALL ZTGSJA( 'u', 'n', 'n', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B,
474 $ 1, TOLA, TOLB, R1, R2, U, 0, V, 1, Q, 1, W,
475 $ NCYCLE, INFO )
476 CALL CHKXER( 'ztgsja', INFOT, NOUT, LERR, OK )
477 INFOT = 20
478 CALL ZTGSJA( 'n', 'v', 'n', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B,
479 $ 1, TOLA, TOLB, R1, R2, U, 1, V, 0, Q, 1, W,
480 $ NCYCLE, INFO )
481 CALL CHKXER( 'ztgsja', INFOT, NOUT, LERR, OK )
482 INFOT = 22
483 CALL ZTGSJA( 'n', 'n', 'q', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B,
484 $ 1, TOLA, TOLB, R1, R2, U, 1, V, 1, Q, 0, W,
485 $ NCYCLE, INFO )
486 CALL CHKXER( 'ztgsja', INFOT, NOUT, LERR, OK )
487 NT = NT + 11
488*
489* Test error exits for the GLM path.
490*
491 ELSE IF( LSAMEN( 3, PATH, 'glm' ) ) THEN
492*
493* ZGGGLM
494*
495 SRNAMT = 'zggglm'
496 INFOT = 1
497 CALL ZGGGLM( -1, 0, 0, A, 1, B, 1, TAU, ALPHA, BETA, W, LW,
498 $ INFO )
499 CALL CHKXER( 'zggglm', INFOT, NOUT, LERR, OK )
500 INFOT = 2
501 CALL ZGGGLM( 0, -1, 0, A, 1, B, 1, TAU, ALPHA, BETA, W, LW,
502 $ INFO )
503 CALL CHKXER( 'zggglm', INFOT, NOUT, LERR, OK )
504 INFOT = 2
505 CALL ZGGGLM( 0, 1, 0, A, 1, B, 1, TAU, ALPHA, BETA, W, LW,
506 $ INFO )
507 CALL CHKXER( 'zggglm', INFOT, NOUT, LERR, OK )
508 INFOT = 3
509 CALL ZGGGLM( 0, 0, -1, A, 1, B, 1, TAU, ALPHA, BETA, W, LW,
510 $ INFO )
511 CALL CHKXER( 'zggglm', INFOT, NOUT, LERR, OK )
512 INFOT = 3
513 CALL ZGGGLM( 1, 0, 0, A, 1, B, 1, TAU, ALPHA, BETA, W, LW,
514 $ INFO )
515 CALL CHKXER( 'zggglm', INFOT, NOUT, LERR, OK )
516 INFOT = 5
517 CALL ZGGGLM( 0, 0, 0, A, 0, B, 1, TAU, ALPHA, BETA, W, LW,
518 $ INFO )
519 CALL CHKXER( 'zggglm', INFOT, NOUT, LERR, OK )
520 INFOT = 7
521 CALL ZGGGLM( 0, 0, 0, A, 1, B, 0, TAU, ALPHA, BETA, W, LW,
522 $ INFO )
523 CALL CHKXER( 'zggglm', infot, nout, lerr, ok )
524 infot = 12
525 CALL zggglm( 1, 1, 1, a, 1, b, 1, tau, alpha, beta, w, 1,
526 $ info )
527 CALL chkxer( 'ZGGGLM', infot, nout, lerr, ok )
528 nt = nt + 8
529*
530* Test error exits for the LSE path.
531*
532 ELSE IF( lsamen( 3, path, 'LSE' ) ) THEN
533*
534* ZGGLSE
535*
536 srnamt = 'ZGGLSE'
537 infot = 1
538 CALL zgglse( -1, 0, 0, a, 1, b, 1, tau, alpha, beta, w, lw,
539 $ info )
540 CALL chkxer( 'ZGGLSE', infot, nout, lerr, ok )
541 infot = 2
542 CALL zgglse( 0, -1, 0, a, 1, b, 1, tau, alpha, beta, w, lw,
543 $ info )
544 CALL chkxer( 'ZGGLSE', infot, nout, lerr, ok )
545 infot = 3
546 CALL zgglse( 0, 0, -1, a, 1, b, 1, tau, alpha, beta, w, lw,
547 $ info )
548 CALL chkxer( 'ZGGLSE', infot, nout, lerr, ok )
549 infot = 3
550 CALL zgglse( 0, 0, 1, a, 1, b, 1, tau, alpha, beta, w, lw,
551 $ info )
552 CALL chkxer( 'ZGGLSE', infot, nout, lerr, ok )
553 infot = 3
554 CALL zgglse( 0, 1, 0, a, 1, b, 1, tau, alpha, beta, w, lw,
555 $ info )
556 CALL chkxer( 'ZGGLSE', infot, nout, lerr, ok )
557 infot = 5
558 CALL zgglse( 0, 0, 0, a, 0, b, 1, tau, alpha, beta, w, lw,
559 $ info )
560 CALL chkxer( 'ZGGLSE', infot, nout, lerr, ok )
561 infot = 7
562 CALL zgglse( 0, 0, 0, a, 1, b, 0, tau, alpha, beta, w, lw,
563 $ info )
564 CALL chkxer( 'ZGGLSE', infot, nout, lerr, ok )
565 infot = 12
566 CALL zgglse( 1, 1, 1, a, 1, b, 1, tau, alpha, beta, w, 1,
567 $ info )
568 CALL chkxer( 'ZGGLSE', infot, nout, lerr, ok )
569 nt = nt + 8
570*
571* Test error exits for the CSD path.
572*
573 ELSE IF( lsamen( 3, path, 'CSD' ) ) THEN
574*
575* ZUNCSD
576*
577 srnamt = 'ZUNCSD'
578 infot = 7
579 CALL zuncsd( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
580 $ -1, 0, 0, a, 1, a,
581 $ 1, a, 1, a, 1, rs,
582 $ a, 1, a, 1, a, 1, a,
583 $ 1, w, lw, rw, lw, iw, info )
584 CALL chkxer( 'ZUNCSD', infot, nout, lerr, ok )
585 infot = 8
586 CALL zuncsd( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
587 $ 1, -1, 0, a, 1, a,
588 $ 1, a, 1, a, 1, rs,
589 $ a, 1, a, 1, a, 1, a,
590 $ 1, w, lw, rw, lw, iw, info )
591 CALL chkxer( 'zuncsd', INFOT, NOUT, LERR, OK )
592 INFOT = 9
593 CALL ZUNCSD( 'y', 'y', 'y', 'y', 'n', 'n',
594 $ 1, 1, -1, A, 1, A,
595 $ 1, A, 1, A, 1, RS,
596 $ A, 1, A, 1, A, 1, A,
597 $ 1, W, LW, RW, LW, IW, INFO )
598 CALL CHKXER( 'zuncsd', INFOT, NOUT, LERR, OK )
599 INFOT = 11
600 CALL ZUNCSD( 'y', 'y', 'y', 'y', 'n', 'n',
601 $ 1, 1, 1, A, -1, A,
602 $ 1, A, 1, A, 1, RS,
603 $ A, 1, A, 1, A, 1, A,
604 $ 1, W, LW, RW, LW, IW, INFO )
605 CALL CHKXER( 'zuncsd', INFOT, NOUT, LERR, OK )
606 INFOT = 20
607 CALL ZUNCSD( 'y', 'y', 'y', 'y', 'n', 'N',
608 $ 1, 1, 1, a, 1, a,
609 $ 1, a, 1, a, 1, rs,
610 $ a, -1, a, 1, a, 1, a,
611 $ 1, w, lw, rw, lw, iw, info )
612 CALL chkxer( 'ZUNCSD', infot, nout, lerr, ok )
613 infot = 22
614 CALL zuncsd( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
615 $ 1, 1, 1, a, 1, a,
616 $ 1, a, 1, a, 1, rs,
617 $ a, 1, a, -1, a, 1, a,
618 $ 1, w, lw, rw, lw, iw, info )
619 CALL chkxer( 'ZUNCSD', infot, nout, lerr, ok )
620 infot = 24
621 CALL zuncsd( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
622 $ 1, 1, 1, a, 1, a,
623 $ 1, a, 1, a, 1, rs,
624 $ a, 1, a, 1, a, -1, a,
625 $ 1, w, lw, rw, lw, iw, info )
626 CALL chkxer( 'ZUNCSD', infot, nout, lerr, ok )
627 infot = 26
628 CALL zuncsd( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
629 $ 1, 1, 1, a, 1, a,
630 $ 1, a, 1, a, 1, rs,
631 $ a, 1, a, 1, a, 1, a,
632 $ -1, w, lw, rw, lw, iw, info )
633 CALL chkxer( 'ZUNCSD', infot, nout, lerr, ok )
634 nt = nt + 8
635*
636* Test error exits for the GQR path.
637*
638 ELSE IF( lsamen( 3, path, 'gqr' ) ) THEN
639*
640* ZGGQRF
641*
642 SRNAMT = 'zggqrf'
643 INFOT = 1
644 CALL ZGGQRF( -1, 0, 0, A, 1, ALPHA, B, 1, BETA, W, LW, INFO )
645 CALL CHKXER( 'zggqrf', INFOT, NOUT, LERR, OK )
646 INFOT = 2
647 CALL ZGGQRF( 0, -1, 0, A, 1, ALPHA, B, 1, BETA, W, LW, INFO )
648 CALL CHKXER( 'zggqrf', INFOT, NOUT, LERR, OK )
649 INFOT = 3
650 CALL ZGGQRF( 0, 0, -1, A, 1, ALPHA, B, 1, BETA, W, LW, INFO )
651 CALL CHKXER( 'zggqrf', INFOT, NOUT, LERR, OK )
652 INFOT = 5
653 CALL ZGGQRF( 0, 0, 0, A, 0, ALPHA, B, 1, BETA, W, LW, INFO )
654 CALL CHKXER( 'zggqrf', INFOT, NOUT, LERR, OK )
655 INFOT = 8
656 CALL ZGGQRF( 0, 0, 0, A, 1, ALPHA, B, 0, BETA, W, LW, INFO )
657 CALL CHKXER( 'zggqrf', INFOT, NOUT, LERR, OK )
658 INFOT = 11
659 CALL ZGGQRF( 1, 1, 2, A, 1, ALPHA, B, 1, BETA, W, 1, INFO )
660 CALL CHKXER( 'zggqrf', INFOT, NOUT, LERR, OK )
661 NT = NT + 6
662*
663* ZGGRQF
664*
665 SRNAMT = 'zggrqf'
666 INFOT = 1
667 CALL ZGGRQF( -1, 0, 0, A, 1, ALPHA, B, 1, BETA, W, LW, INFO )
668 CALL CHKXER( 'zggrqf', INFOT, NOUT, LERR, OK )
669 INFOT = 2
670 CALL ZGGRQF( 0, -1, 0, A, 1, ALPHA, B, 1, BETA, W, LW, INFO )
671 CALL CHKXER( 'zggrqf', INFOT, NOUT, LERR, OK )
672 INFOT = 3
673 CALL ZGGRQF( 0, 0, -1, A, 1, ALPHA, B, 1, BETA, W, LW, INFO )
674 CALL CHKXER( 'zggrqf', INFOT, NOUT, LERR, OK )
675 INFOT = 5
676 CALL ZGGRQF( 0, 0, 0, A, 0, ALPHA, B, 1, BETA, W, LW, INFO )
677 CALL CHKXER( 'zggrqf', INFOT, NOUT, LERR, OK )
678 INFOT = 8
679 CALL ZGGRQF( 0, 0, 0, A, 1, ALPHA, B, 0, BETA, W, LW, INFO )
680 CALL CHKXER( 'zggrqf', INFOT, NOUT, LERR, OK )
681 INFOT = 11
682 CALL ZGGRQF( 1, 1, 2, A, 1, ALPHA, B, 1, BETA, W, 1, INFO )
683 CALL CHKXER( 'zggrqf', INFOT, NOUT, LERR, OK )
684 NT = NT + 6
685*
686* Test error exits for the ZGS, ZGV, ZGX, and ZXV paths.
687*
688 ELSE IF( LSAMEN( 3, PATH, 'zgs.OR.' )
689 $ LSAMEN( 3, PATH, 'zgv.OR.' )
690 $ LSAMEN( 3, PATH, 'zgx.OR.' ) LSAMEN( 3, PATH, 'zxv' ) )
691 $ THEN
692*
693* ZGGES
694*
695 SRNAMT = 'zgges '
696 INFOT = 1
697 CALL ZGGES( '/', 'n', 's', ZLCTES, 1, A, 1, B, 1, SDIM, ALPHA,
698 $ BETA, Q, 1, U, 1, W, 1, RW, BW, INFO )
699 CALL CHKXER( 'zgges ', INFOT, NOUT, LERR, OK )
700 INFOT = 2
701 CALL ZGGES( 'n', '/', 's', ZLCTES, 1, A, 1, B, 1, SDIM, ALPHA,
702 $ BETA, Q, 1, U, 1, W, 1, RW, BW, INFO )
703 CALL CHKXER( 'zgges ', INFOT, NOUT, LERR, OK )
704 INFOT = 3
705 CALL ZGGES( 'n', 'v', '/', zlctes, 1, a, 1, b, 1, sdim, alpha,
706 $ beta, q, 1, u, 1, w, 1, rw, bw, info )
707 CALL chkxer( 'ZGGES ', infot, nout, lerr, ok )
708 infot = 5
709 CALL zgges( 'N', 'V', 'S', zlctes, -1, a, 1, b, 1, sdim, alpha,
710 $ beta, q, 1, u, 1, w, 1, rw, bw, info )
711 CALL chkxer( 'ZGGES ', infot, nout, lerr, ok )
712 infot = 7
713 CALL zgges( 'N', 'V', 'S', zlctes, 1, a, 0, b, 1, sdim, alpha,
714 $ beta, q, 1, u, 1, w, 1, rw, bw, info )
715 CALL chkxer( 'ZGGES ', infot, nout, lerr, ok )
716 infot = 9
717 CALL zgges( 'N', 'V', 'S', zlctes, 1, a, 1, b, 0, sdim, alpha,
718 $ beta, q, 1, u, 1, w, 1, rw, bw, info )
719 CALL chkxer( 'ZGGES ', infot, nout, lerr, ok )
720 infot = 14
721 CALL zgges( 'N', 'V', 'S', zlctes, 1, a, 1, b, 1, sdim, alpha,
722 $ beta, q, 0, u, 1, w, 1, rw, bw, info )
723 CALL chkxer( 'ZGGES ', infot, nout, lerr, ok )
724 infot = 14
725 CALL zgges( 'V', 'V', 'S', zlctes, 2, a, 2, b, 2, sdim, alpha,
726 $ beta, q, 1, u, 2, w, 1, rw, bw, info )
727 CALL chkxer( 'ZGGES ', infot, nout, lerr, ok )
728 infot = 16
729 CALL zgges( 'N', 'V', 'S', zlctes, 1, a, 1, b, 1, sdim, alpha,
730 $ beta, q, 1, u, 0, w, 1, rw, bw, info )
731 CALL chkxer( 'ZGGES ', infot, nout, lerr, ok )
732 infot = 16
733 CALL zgges( 'V', 'V', 'S', zlctes, 2, a, 2, b, 2, sdim, alpha,
734 $ beta, q, 2, u, 1, w, 1, rw, bw, info )
735 CALL chkxer( 'ZGGES ', infot, nout, lerr, ok )
736 infot = 18
737 CALL zgges( 'V', 'V', 'S', zlctes, 2, a, 2, b, 2, sdim, alpha,
738 $ beta, q, 2, u, 2, w, 1, rw, bw, info )
739 CALL chkxer( 'ZGGES ', infot, nout, lerr, ok )
740 nt = nt + 11
741*
742* ZGGES3
743*
744 srnamt = 'ZGGES3'
745 infot = 1
746 CALL zgges3( '/', 'N', 'S', zlctes, 1, a, 1, b, 1, sdim, alpha,
747 $ beta, q, 1, u, 1, w, 1, rw, bw, info )
748 CALL chkxer( 'ZGGES3', infot, nout, lerr, ok )
749 infot = 2
750 CALL zgges3( 'n', '/', 's', ZLCTES, 1, A, 1, B, 1, SDIM, ALPHA,
751 $ BETA, Q, 1, U, 1, W, 1, RW, BW, INFO )
752 CALL CHKXER( 'zgges3', INFOT, NOUT, LERR, OK )
753 INFOT = 3
754 CALL ZGGES3( 'n', 'v', '/', ZLCTES, 1, A, 1, B, 1, SDIM, ALPHA,
755 $ BETA, Q, 1, U, 1, W, 1, RW, BW, INFO )
756 CALL CHKXER( 'zgges3', INFOT, NOUT, LERR, OK )
757 INFOT = 5
758 CALL ZGGES3( 'n', 'v', 's', ZLCTES, -1, A, 1, B, 1, SDIM,
759 $ ALPHA, BETA, Q, 1, U, 1, W, 1, RW, BW, INFO )
760 CALL CHKXER( 'zgges3', INFOT, NOUT, LERR, OK )
761 INFOT = 7
762 CALL ZGGES3( 'n', 'v', 's', ZLCTES, 1, A, 0, B, 1, SDIM, ALPHA,
763 $ BETA, Q, 1, U, 1, W, 1, RW, BW, INFO )
764 CALL CHKXER( 'zgges3', INFOT, NOUT, LERR, OK )
765 INFOT = 9
766 CALL ZGGES3( 'n', 'v', 's', ZLCTES, 1, A, 1, B, 0, SDIM, ALPHA,
767 $ BETA, Q, 1, U, 1, W, 1, RW, BW, INFO )
768 CALL CHKXER( 'zgges3', INFOT, NOUT, LERR, OK )
769 INFOT = 14
770 CALL ZGGES3( 'n', 'v', 's', ZLCTES, 1, A, 1, B, 1, SDIM, ALPHA,
771 $ BETA, Q, 0, U, 1, W, 1, RW, BW, INFO )
772 CALL CHKXER( 'zgges3', INFOT, NOUT, LERR, OK )
773 INFOT = 14
774 CALL ZGGES3( 'v', 'v', 's', ZLCTES, 2, A, 2, B, 2, SDIM, ALPHA,
775 $ BETA, Q, 1, U, 2, W, 1, RW, BW, INFO )
776 CALL CHKXER( 'zgges3', INFOT, NOUT, LERR, OK )
777 INFOT = 16
778 CALL ZGGES3( 'n', 'v', 's', ZLCTES, 1, A, 1, B, 1, SDIM, ALPHA,
779 $ BETA, Q, 1, U, 0, W, 1, RW, BW, INFO )
780 CALL CHKXER( 'zgges3', INFOT, NOUT, LERR, OK )
781 INFOT = 16
782 CALL ZGGES3( 'v', 'v', 's', ZLCTES, 2, A, 2, B, 2, SDIM, ALPHA,
783 $ BETA, Q, 2, U, 1, W, 1, RW, BW, INFO )
784 CALL CHKXER( 'zgges3', INFOT, NOUT, LERR, OK )
785 INFOT = 18
786 CALL ZGGES3( 'v', 'v', 's', ZLCTES, 2, A, 2, B, 2, SDIM, ALPHA,
787 $ BETA, Q, 2, U, 2, W, 1, RW, BW, INFO )
788 CALL CHKXER( 'zgges3', INFOT, NOUT, LERR, OK )
789 NT = NT + 11
790*
791* ZGGESX
792*
793 SRNAMT = 'zggesx'
794 INFOT = 1
795 CALL ZGGESX( '/', 'n', 's', ZLCTSX, 'n', 1, A, 1, B, 1, SDIM,
796 $ ALPHA, BETA, Q, 1, U, 1, RCE, RCV, W, 1, RW, IW,
797 $ 1, BW, INFO )
798 CALL CHKXER( 'zggesx', INFOT, NOUT, LERR, OK )
799 INFOT = 2
800 CALL ZGGESX( 'n', '/', 's', ZLCTSX, 'n', 1, A, 1, B, 1, SDIM,
801 $ ALPHA, BETA, Q, 1, U, 1, RCE, RCV, W, 1, RW, IW,
802 $ 1, BW, INFO )
803 CALL CHKXER( 'zggesx', INFOT, NOUT, LERR, OK )
804 INFOT = 3
805 CALL ZGGESX( 'v', 'v', '/', ZLCTSX, 'n', 1, A, 1, B, 1, SDIM,
806 $ ALPHA, BETA, Q, 1, U, 1, RCE, RCV, W, 1, RW, IW,
807 $ 1, BW, INFO )
808 CALL CHKXER( 'zggesx', INFOT, NOUT, LERR, OK )
809 INFOT = 5
810 CALL ZGGESX( 'v', 'v', 's', ZLCTSX, '/', 1, A, 1, B, 1, SDIM,
811 $ ALPHA, BETA, Q, 1, U, 1, RCE, RCV, W, 1, RW, IW,
812 $ 1, BW, INFO )
813 CALL CHKXER( 'zggesx', INFOT, NOUT, LERR, OK )
814 INFOT = 6
815 CALL ZGGESX( 'v', 'v', 's', ZLCTSX, 'b', -1, A, 1, B, 1, SDIM,
816 $ ALPHA, BETA, Q, 1, U, 1, RCE, RCV, W, 1, RW, IW,
817 $ 1, BW, INFO )
818 CALL CHKXER( 'zggesx', INFOT, NOUT, LERR, OK )
819 INFOT = 8
820 CALL ZGGESX( 'v', 'v', 's', ZLCTSX, 'b', 1, A, 0, B, 1, SDIM,
821 $ ALPHA, BETA, Q, 1, U, 1, RCE, RCV, W, 1, RW, IW,
822 $ 1, BW, INFO )
823 CALL CHKXER( 'zggesx', INFOT, NOUT, LERR, OK )
824 INFOT = 10
825 CALL ZGGESX( 'v', 'v', 's', ZLCTSX, 'b', 1, A, 1, B, 0, SDIM,
826 $ ALPHA, BETA, Q, 1, U, 1, RCE, RCV, W, 1, RW, IW,
827 $ 1, BW, INFO )
828 CALL CHKXER( 'zggesx', INFOT, NOUT, LERR, OK )
829 INFOT = 15
830 CALL ZGGESX( 'v', 'v', 's', ZLCTSX, 'b', 1, A, 1, B, 1, SDIM,
831 $ ALPHA, BETA, Q, 0, U, 1, RCE, RCV, W, 1, RW, IW,
832 $ 1, BW, INFO )
833 CALL CHKXER( 'zggesx', INFOT, NOUT, LERR, OK )
834 INFOT = 15
835 CALL ZGGESX( 'v', 'v', 's', ZLCTSX, 'b', 2, A, 2, B, 2, SDIM,
836 $ ALPHA, BETA, Q, 1, U, 1, RCE, RCV, W, 1, RW, IW,
837 $ 1, BW, INFO )
838 CALL CHKXER( 'zggesx', INFOT, NOUT, LERR, OK )
839 INFOT = 17
840 CALL ZGGESX( 'v', 'v', 's', ZLCTSX, 'b', 1, A, 1, B, 1, SDIM,
841 $ ALPHA, BETA, Q, 1, U, 0, RCE, RCV, W, 1, RW, IW,
842 $ 1, BW, INFO )
843 CALL CHKXER( 'zggesx', INFOT, NOUT, LERR, OK )
844 INFOT = 17
845 CALL ZGGESX( 'v', 'v', 's', ZLCTSX, 'b', 2, A, 2, B, 2, SDIM,
846 $ ALPHA, BETA, Q, 2, U, 1, RCE, RCV, W, 1, RW, IW,
847 $ 1, BW, INFO )
848 CALL CHKXER( 'zggesx', INFOT, NOUT, LERR, OK )
849 INFOT = 21
850 CALL ZGGESX( 'v', 'v', 's', ZLCTSX, 'b', 2, A, 2, B, 2, SDIM,
851 $ ALPHA, BETA, Q, 2, U, 2, RCE, RCV, W, 1, RW, IW,
852 $ 1, BW, INFO )
853 CALL CHKXER( 'zggesx', INFOT, NOUT, LERR, OK )
854 INFOT = 24
855 CALL ZGGESX( 'v', 'v', 's', ZLCTSX, 'v', 1, A, 1, B, 1, SDIM,
856 $ ALPHA, BETA, Q, 1, U, 1, RCE, RCV, W, 32, RW, IW,
857 $ 0, BW, INFO )
858 CALL CHKXER( 'zggesx', INFOT, NOUT, LERR, OK )
859 NT = NT + 13
860*
861* ZGGEV
862*
863 SRNAMT = 'zggev '
864 INFOT = 1
865 CALL ZGGEV( '/', 'n', 1, A, 1, B, 1, ALPHA, BETA, Q, 1, U, 1,
866 $ W, 1, RW, INFO )
867 CALL CHKXER( 'zggev ', INFOT, NOUT, LERR, OK )
868 INFOT = 2
869 CALL ZGGEV( 'n', '/', 1, A, 1, B, 1, ALPHA, BETA, Q, 1, U, 1,
870 $ W, 1, RW, INFO )
871 CALL CHKXER( 'zggev ', INFOT, NOUT, LERR, OK )
872 INFOT = 3
873 CALL ZGGEV( 'v', 'v', -1, A, 1, B, 1, ALPHA, BETA, Q, 1, U, 1,
874 $ W, 1, RW, INFO )
875 CALL CHKXER( 'zggev ', INFOT, NOUT, LERR, OK )
876 INFOT = 5
877 CALL ZGGEV( 'v', 'v', 1, A, 0, B, 1, ALPHA, BETA, Q, 1, U, 1,
878 $ W, 1, RW, INFO )
879 CALL CHKXER( 'zggev ', INFOT, NOUT, LERR, OK )
880 INFOT = 7
881 CALL ZGGEV( 'v', 'v', 1, A, 1, B, 0, ALPHA, BETA, Q, 1, U, 1,
882 $ W, 1, RW, INFO )
883 CALL CHKXER( 'zggev ', INFOT, NOUT, LERR, OK )
884 INFOT = 11
885 CALL ZGGEV( 'n', 'v', 1, A, 1, B, 1, ALPHA, BETA, Q, 0, U, 1,
886 $ W, 1, RW, INFO )
887 CALL CHKXER( 'zggev ', INFOT, NOUT, LERR, OK )
888 INFOT = 11
889 CALL ZGGEV( 'v', 'v', 2, A, 2, B, 2, ALPHA, BETA, Q, 1, U, 2,
890 $ W, 1, RW, INFO )
891 CALL CHKXER( 'zggev ', INFOT, NOUT, LERR, OK )
892 INFOT = 13
893 CALL ZGGEV( 'v', 'n', 2, A, 2, B, 2, ALPHA, BETA, Q, 2, U, 0,
894 $ W, 1, RW, INFO )
895 CALL CHKXER( 'zggev ', INFOT, NOUT, LERR, OK )
896 INFOT = 13
897 CALL ZGGEV( 'v', 'v', 2, A, 2, B, 2, ALPHA, BETA, Q, 2, U, 1,
898 $ W, 1, RW, INFO )
899 CALL CHKXER( 'zggev ', INFOT, NOUT, LERR, OK )
900 INFOT = 15
901 CALL ZGGEV( 'v', 'v', 1, A, 1, B, 1, ALPHA, BETA, Q, 1, U, 1,
902 $ W, 1, RW, INFO )
903 CALL CHKXER( 'zggev ', INFOT, NOUT, LERR, OK )
904 NT = NT + 10
905*
906* ZGGEV3
907*
908 SRNAMT = 'zggev3'
909 INFOT = 1
910 CALL ZGGEV3( '/', 'n', 1, A, 1, B, 1, ALPHA, BETA, Q, 1, U, 1,
911 $ W, 1, RW, INFO )
912 CALL CHKXER( 'zggev3', INFOT, NOUT, LERR, OK )
913 INFOT = 2
914 CALL ZGGEV3( 'n', '/', 1, A, 1, B, 1, ALPHA, BETA, Q, 1, U, 1,
915 $ W, 1, RW, INFO )
916 CALL CHKXER( 'zggev3', INFOT, NOUT, LERR, OK )
917 INFOT = 3
918 CALL ZGGEV3( 'v', 'v', -1, A, 1, B, 1, ALPHA, BETA, Q, 1, U, 1,
919 $ W, 1, RW, INFO )
920 CALL CHKXER( 'zggev3', INFOT, NOUT, LERR, OK )
921 INFOT = 5
922 CALL ZGGEV3( 'v', 'v', 1, A, 0, B, 1, ALPHA, BETA, Q, 1, U, 1,
923 $ W, 1, RW, INFO )
924 CALL CHKXER( 'zggev3', INFOT, NOUT, LERR, OK )
925 INFOT = 7
926 CALL ZGGEV3( 'v', 'v', 1, A, 1, B, 0, ALPHA, BETA, Q, 1, U, 1,
927 $ W, 1, RW, INFO )
928 CALL CHKXER( 'zggev3', INFOT, NOUT, LERR, OK )
929 INFOT = 11
930 CALL ZGGEV3( 'n', 'v', 1, A, 1, B, 1, ALPHA, BETA, Q, 0, U, 1,
931 $ W, 1, RW, INFO )
932 CALL CHKXER( 'zggev3', INFOT, NOUT, LERR, OK )
933 INFOT = 11
934 CALL ZGGEV3( 'v', 'v', 2, A, 2, B, 2, ALPHA, BETA, Q, 1, U, 2,
935 $ W, 1, RW, INFO )
936 CALL CHKXER( 'zggev3', INFOT, NOUT, LERR, OK )
937 INFOT = 13
938 CALL ZGGEV3( 'v', 'n', 2, A, 2, B, 2, ALPHA, BETA, Q, 2, U, 0,
939 $ W, 1, RW, INFO )
940 CALL CHKXER( 'zggev3', INFOT, NOUT, LERR, OK )
941 INFOT = 13
942 CALL ZGGEV3( 'v', 'v', 2, A, 2, B, 2, ALPHA, BETA, Q, 2, U, 1,
943 $ W, 1, RW, INFO )
944 CALL CHKXER( 'zggev3', INFOT, NOUT, LERR, OK )
945 INFOT = 15
946 CALL ZGGEV3( 'v', 'v', 1, A, 1, B, 1, ALPHA, BETA, Q, 1, U, 1,
947 $ W, 1, RW, INFO )
948 CALL CHKXER( 'zggev3', INFOT, NOUT, LERR, OK )
949 NT = NT + 10
950*
951* ZGGEVX
952*
953 SRNAMT = 'zggevx'
954 INFOT = 1
955 CALL ZGGEVX( '/', 'n', 'n', 'n', 1, A, 1, B, 1, ALPHA, BETA, Q,
956 $ 1, U, 1, ILO, IHI, LS, RS, ANRM, BNRM, RCE, RCV,
957 $ W, 1, RW, IW, BW, INFO )
958 CALL CHKXER( 'zggevx', INFOT, NOUT, LERR, OK )
959 INFOT = 2
960 CALL ZGGEVX( 'n', '/', 'n', 'n', 1, A, 1, B, 1, ALPHA, BETA, Q,
961 $ 1, U, 1, ILO, IHI, LS, RS, ANRM, BNRM, RCE, RCV,
962 $ W, 1, RW, IW, BW, INFO )
963 CALL CHKXER( 'zggevx', INFOT, NOUT, LERR, OK )
964 INFOT = 3
965 CALL ZGGEVX( 'n', 'n', '/', 'n', 1, A, 1, B, 1, ALPHA, BETA, Q,
966 $ 1, U, 1, ILO, IHI, LS, RS, ANRM, BNRM, RCE, RCV,
967 $ W, 1, RW, IW, BW, INFO )
968 CALL CHKXER( 'zggevx', INFOT, NOUT, LERR, OK )
969 INFOT = 4
970 CALL ZGGEVX( 'n', 'n', 'n', '/', 1, A, 1, B, 1, ALPHA, BETA, Q,
971 $ 1, U, 1, ILO, IHI, LS, RS, ANRM, BNRM, RCE, RCV,
972 $ W, 1, RW, IW, BW, INFO )
973 CALL CHKXER( 'zggevx', INFOT, NOUT, LERR, OK )
974 INFOT = 5
975 CALL ZGGEVX( 'n', 'n', 'n', 'n', -1, A, 1, B, 1, ALPHA, BETA,
976 $ Q, 1, U, 1, ILO, IHI, LS, RS, ANRM, BNRM, RCE,
977 $ RCV, W, 1, RW, IW, BW, INFO )
978 CALL CHKXER( 'zggevx', INFOT, NOUT, LERR, OK )
979 INFOT = 7
980 CALL ZGGEVX( 'n', 'n', 'n', 'n', 1, A, 0, B, 1, ALPHA, BETA, Q,
981 $ 1, U, 1, ILO, IHI, LS, RS, ANRM, BNRM, RCE, RCV,
982 $ W, 1, RW, IW, BW, INFO )
983 CALL CHKXER( 'zggevx', INFOT, NOUT, LERR, OK )
984 INFOT = 9
985 CALL ZGGEVX( 'n', 'n', 'n', 'n', 1, A, 1, B, 0, ALPHA, BETA, Q,
986 $ 1, U, 1, ILO, IHI, LS, RS, ANRM, BNRM, RCE, RCV,
987 $ W, 1, RW, IW, BW, INFO )
988 CALL CHKXER( 'zggevx', INFOT, NOUT, LERR, OK )
989 INFOT = 13
990 CALL ZGGEVX( 'n', 'n', 'n', 'n', 1, A, 1, B, 1, ALPHA, BETA, Q,
991 $ 0, U, 1, ILO, IHI, LS, RS, ANRM, BNRM, RCE, RCV,
992 $ W, 1, RW, IW, BW, INFO )
993 CALL CHKXER( 'zggevx', INFOT, NOUT, LERR, OK )
994 INFOT = 13
995 CALL ZGGEVX( 'n', 'v', 'n', 'n', 2, A, 2, B, 2, ALPHA, BETA, Q,
996 $ 1, U, 2, ILO, IHI, LS, RS, ANRM, BNRM, RCE, RCV,
997 $ W, 1, RW, IW, BW, INFO )
998 CALL CHKXER( 'zggevx', INFOT, NOUT, LERR, OK )
999 INFOT = 15
1000 CALL ZGGEVX( 'n', 'n', 'n', 'n', 1, A, 1, B, 1, ALPHA, BETA, Q,
1001 $ 1, U, 0, ILO, IHI, LS, RS, ANRM, BNRM, RCE, RCV,
1002 $ W, 1, RW, IW, BW, INFO )
1003 CALL CHKXER( 'zggevx', INFOT, NOUT, LERR, OK )
1004 INFOT = 15
1005 CALL ZGGEVX( 'n', 'n', 'v', 'n', 2, A, 2, B, 2, ALPHA, BETA, Q,
1006 $ 2, U, 1, ILO, IHI, LS, RS, ANRM, BNRM, RCE, RCV,
1007 $ W, 1, RW, IW, BW, INFO )
1008 CALL CHKXER( 'zggevx', INFOT, NOUT, LERR, OK )
1009 INFOT = 25
1010 CALL ZGGEVX( 'n', 'n', 'v', 'n', 2, a, 2, b, 2, alpha, beta, q,
1011 $ 2, u, 2, ilo, ihi, ls, rs, anrm, bnrm, rce, rcv,
1012 $ w, 0, rw, iw, bw, info )
1013 CALL chkxer( 'ZGGEVX', infot, nout, lerr, ok )
1014 nt = nt + 12
1015*
1016* ZTGEXC
1017*
1018 srnamt = 'ZTGEXC'
1019 infot = 3
1020 CALL ztgexc( .true., .true., -1, a, 1, b, 1, q, 1, z, 1, ifst,
1021 $ ilst, info )
1022 CALL chkxer( 'ZTGEXC', infot, nout, lerr, ok )
1023 infot = 5
1024 CALL ztgexc( .true., .true., 1, a, 0, b, 1, q, 1, z, 1, ifst,
1025 $ ilst, info )
1026 CALL chkxer( 'ZTGEXC', infot, nout, lerr, ok )
1027 infot = 7
1028 CALL ztgexc( .true., .true., 1, a, 1, b, 0, q, 1, z, 1, ifst,
1029 $ ilst, info )
1030 CALL chkxer( 'ZTGEXC', infot, nout, lerr, ok )
1031 infot = 9
1032 CALL ztgexc( .false., .true., 1, a, 1, b, 1, q, 0, z, 1, ifst,
1033 $ ilst, info )
1034 CALL chkxer( 'ZTGEXC', infot, nout, lerr, ok )
1035 infot = 9
1036 CALL ztgexc( .true., .true., 1, a, 1, b, 1, q, 0, z, 1, ifst,
1037 $ ilst, info )
1038 CALL chkxer( 'ZTGEXC', infot, nout, lerr, ok )
1039 infot = 11
1040 CALL ztgexc( .true., .false., 1, a, 1, b, 1, q, 1, z, 0, ifst,
1041 $ ilst, info )
1042 CALL chkxer( 'ZTGEXC', infot, nout, lerr, ok )
1043 infot = 11
1044 CALL ztgexc( .true., .true., 1, a, 1, b, 1, q, 1, z, 0, ifst,
1045 $ ilst, info )
1046 CALL chkxer( 'ZTGEXC', infot, nout, lerr, ok )
1047 nt = nt + 7
1048*
1049* ZTGSEN
1050*
1051 srnamt = 'ZTGSEN'
1052 infot = 1
1053 CALL ztgsen( -1, .true., .true., sel, 1, a, 1, b, 1, alpha,
1054 $ beta, q, 1, z, 1, m, tola, tolb, rcv, w, 1, iw, 1,
1055 $ info )
1056 CALL chkxer( 'ZTGSEN', infot, nout, lerr, ok )
1057 infot = 5
1058 CALL ztgsen( 1, .true., .true., sel, -1, a, 1, b, 1, alpha,
1059 $ beta, q, 1, z, 1, m, tola, tolb, rcv, w, 1, iw, 1,
1060 $ info )
1061 CALL chkxer( 'ZTGSEN', infot, nout, lerr, ok )
1062 infot = 7
1063 CALL ztgsen( 1, .true., .true., sel, 1, a, 0, b, 1, alpha,
1064 $ beta, q, 1, z, 1, m, tola, tolb, rcv, w, 1, iw, 1,
1065 $ info )
1066 CALL chkxer( 'ZTGSEN', infot, nout, lerr, ok )
1067 infot = 9
1068 CALL ztgsen( 1, .true., .true., sel, 1, a, 1, b, 0, alpha,
1069 $ beta, q, 1, z, 1, m, tola, tolb, rcv, w, 1, iw, 1,
1070 $ info )
1071 CALL chkxer( 'ZTGSEN', infot, nout, lerr, ok )
1072 infot = 13
1073 CALL ztgsen( 1, .true., .true., sel, 1, a, 1, b, 1, alpha,
1074 $ beta, q, 0, z, 1, m, tola, tolb, rcv, w, 1, iw, 1,
1075 $ info )
1076 CALL chkxer( 'ZTGSEN', infot, nout, lerr, ok )
1077 infot = 15
1078 CALL ztgsen( 1, .true., .true., sel, 1, a, 1, b, 1, alpha,
1079 $ beta, q, 1, z, 0, m, tola, tolb, rcv, w, 1, iw, 1,
1080 $ info )
1081 CALL chkxer( 'ZTGSEN', infot, nout, lerr, ok )
1082 infot = 21
1083 CALL ztgsen( 3, .true., .true., sel, 1, a, 1, b, 1, alpha,
1084 $ beta, q, 1, z, 1, m, tola, tolb, rcv, w, -5, iw,
1085 $ 1, info )
1086 CALL chkxer( 'ZTGSEN', infot, nout, lerr, ok )
1087 infot = 23
1088 CALL ztgsen( 0, .true., .true., sel, 1, a, 1, b, 1, alpha,
1089 $ beta, q, 1, z, 1, m, tola, tolb, rcv, w, 20, iw,
1090 $ 0, info )
1091 CALL chkxer( 'ZTGSEN', infot, nout, lerr, ok )
1092 infot = 23
1093 CALL ztgsen( 1, .true., .true., sel, 1, a, 1, b, 1, alpha,
1094 $ beta, q, 1, z, 1, m, tola, tolb, rcv, w, 20, iw,
1095 $ 0, info )
1096 CALL chkxer( 'ZTGSEN', infot, nout, lerr, ok )
1097 infot = 23
1098 CALL ztgsen( 5, .true., .true., sel, 1, a, 1, b, 1, alpha,
1099 $ beta, q, 1, z, 1, m, tola, tolb, rcv, w, 20, iw,
1100 $ 1, info )
1101 CALL chkxer( 'ZTGSEN', infot, nout, lerr, ok )
1102 nt = nt + 11
1103*
1104* ZTGSNA
1105*
1106 srnamt = 'ZTGSNA'
1107 infot = 1
1108 CALL ztgsna( '/', 'A', sel, 1, a, 1, b, 1, q, 1, u, 1, r1, r2,
1109 $ 1, m, w, 1, iw, info )
1110 CALL chkxer( 'ZTGSNA', infot, nout, lerr, ok )
1111 infot = 2
1112 CALL ztgsna( 'B', '/', sel, 1, a, 1, b, 1, q, 1, u, 1, r1, r2,
1113 $ 1, m, w, 1, iw, info )
1114 CALL chkxer( 'ZTGSNA', infot, nout, lerr, ok )
1115 infot = 4
1116 CALL ztgsna( 'B', 'A', sel, -1, a, 1, b, 1, q, 1, u, 1, r1, r2,
1117 $ 1, m, w, 1, iw, info )
1118 CALL chkxer( 'ZTGSNA', infot, nout, lerr, ok )
1119 infot = 6
1120 CALL ztgsna( 'B', 'A', sel, 1, a, 0, b, 1, q, 1, u, 1, r1, r2,
1121 $ 1, m, w, 1, iw, info )
1122 CALL chkxer( 'ZTGSNA', infot, nout, lerr, ok )
1123 infot = 8
1124 CALL ztgsna( 'B', 'A', sel, 1, a, 1, b, 0, q, 1, u, 1, r1, r2,
1125 $ 1, m, w, 1, iw, info )
1126 CALL chkxer( 'ZTGSNA', infot, nout, lerr, ok )
1127 infot = 10
1128 CALL ztgsna( 'E', 'A', sel, 1, a, 1, b, 1, q, 0, u, 1, r1, r2,
1129 $ 1, m, w, 1, iw, info )
1130 CALL chkxer( 'ZTGSNA', infot, nout, lerr, ok )
1131 infot = 12
1132 CALL ztgsna( 'E', 'A', sel, 1, a, 1, b, 1, q, 1, u, 0, r1, r2,
1133 $ 1, m, w, 1, iw, info )
1134 CALL chkxer( 'ZTGSNA', infot, nout, lerr, ok )
1135 infot = 15
1136 CALL ztgsna( 'E', 'A', sel, 1, a, 1, b, 1, q, 1, u, 1, r1, r2,
1137 $ 0, m, w, 1, iw, info )
1138 CALL chkxer( 'ZTGSNA', infot, nout, lerr, ok )
1139 infot = 18
1140 CALL ztgsna( 'E', 'A', sel, 1, a, 1, b, 1, q, 1, u, 1, r1, r2,
1141 $ 1, m, w, 0, iw, info )
1142 CALL chkxer( 'ZTGSNA', infot, nout, lerr, ok )
1143 nt = nt + 9
1144*
1145* ZTGSYL
1146*
1147 srnamt = 'ZTGSYL'
1148 infot = 1
1149 CALL ztgsyl( '/', 0, 1, 1, a, 1, b, 1, q, 1, u, 1, v, 1, z, 1,
1150 $ scale, dif, w, 1, iw, info )
1151 CALL chkxer( 'ZTGSYL', infot, nout, lerr, ok )
1152 infot = 2
1153 CALL ztgsyl( 'N', -1, 1, 1, a, 1, b, 1, q, 1, u, 1, v, 1, z, 1,
1154 $ scale, dif, w, 1, iw, info )
1155 CALL chkxer( 'ZTGSYL', infot, nout, lerr, ok )
1156 infot = 3
1157 CALL ztgsyl( 'N', 0, 0, 1, a, 1, b, 1, q, 1, u, 1, v, 1, z, 1,
1158 $ scale, dif, w, 1, iw, info )
1159 CALL chkxer( 'ZTGSYL', infot, nout, lerr, ok )
1160 infot = 4
1161 CALL ztgsyl( 'N', 0, 1, 0, a, 1, b, 1, q, 1, u, 1, v, 1, z, 1,
1162 $ scale, dif, w, 1, iw, info )
1163 CALL chkxer( 'ZTGSYL', infot, nout, lerr, ok )
1164 infot = 6
1165 CALL ztgsyl( 'N', 0, 1, 1, a, 0, b, 1, q, 1, u, 1, v, 1, z, 1,
1166 $ scale, dif, w, 1, iw, info )
1167 CALL chkxer( 'ZTGSYL', infot, nout, lerr, ok )
1168 infot = 8
1169 CALL ztgsyl( 'n', 0, 1, 1, A, 1, B, 0, Q, 1, U, 1, V, 1, Z, 1,
1170 $ SCALE, DIF, W, 1, IW, INFO )
1171 CALL CHKXER( 'ztgsyl', INFOT, NOUT, LERR, OK )
1172 INFOT = 10
1173 CALL ZTGSYL( 'n', 0, 1, 1, A, 1, B, 1, Q, 0, U, 1, V, 1, Z, 1,
1174 $ SCALE, DIF, W, 1, IW, INFO )
1175 CALL CHKXER( 'ztgsyl', INFOT, NOUT, LERR, OK )
1176 INFOT = 12
1177 CALL ZTGSYL( 'n', 0, 1, 1, A, 1, B, 1, Q, 1, U, 0, V, 1, Z, 1,
1178 $ SCALE, DIF, W, 1, IW, INFO )
1179 CALL CHKXER( 'ztgsyl', INFOT, NOUT, LERR, OK )
1180 INFOT = 14
1181 CALL ZTGSYL( 'n', 0, 1, 1, A, 1, B, 1, Q, 1, U, 1, V, 0, Z, 1,
1182 $ SCALE, DIF, W, 1, IW, INFO )
1183 CALL CHKXER( 'ztgsyl', INFOT, NOUT, LERR, OK )
1184 INFOT = 16
1185 CALL ZTGSYL( 'n', 0, 1, 1, A, 1, B, 1, Q, 1, U, 1, V, 1, Z, 0,
1186 $ SCALE, DIF, W, 1, IW, INFO )
1187 CALL CHKXER( 'ztgsyl', INFOT, NOUT, LERR, OK )
1188 INFOT = 20
1189 CALL ZTGSYL( 'n', 1, 1, 1, A, 1, B, 1, Q, 1, U, 1, V, 1, Z, 1,
1190 $ SCALE, DIF, W, 1, IW, INFO )
1191 CALL CHKXER( 'ztgsyl', INFOT, NOUT, LERR, OK )
1192 INFOT = 20
1193 CALL ZTGSYL( 'n', 2, 1, 1, A, 1, B, 1, Q, 1, U, 1, V, 1, Z, 1,
1194 $ SCALE, DIF, W, 1, IW, INFO )
1195 CALL CHKXER( 'ztgsyl', INFOT, NOUT, LERR, OK )
1196 NT = NT + 12
1197 END IF
1198*
1199* Print a summary line.
1200*
1201 IF( OK ) THEN
1202 WRITE( NOUT, FMT = 9999 )PATH, NT
1203 ELSE
1204 WRITE( NOUT, FMT = 9998 )PATH
1205 END IF
1206*
1207 9999 FORMAT( 1X, A3, ' routines passed the tests of the error exits(',
1208 $ I3, ' tests done)' )
1209 9998 FORMAT( ' *** ', A3, ' routines failed the tests of the error ',
1210 $ 'exits ***' )
1211*
1212 RETURN
1213*
1214* End of ZERRGG
1215*
1216 END
chkxer
subroutine chkxer(srnamt, infot, nout, lerr, ok)
Definition cblat2.f:3196
the
end diagonal values have been computed in the(sparse) matrix id.SOL
xlaenv
subroutine xlaenv(ispec, nvalue)
XLAENV
Definition xlaenv.f:81
ztgevc
subroutine ztgevc(side, howmny, select, n, s, lds, p, ldp, vl, ldvl, vr, ldvr, mm, m, work, rwork, info)
ZTGEVC
Definition ztgevc.f:219
ztgexc
subroutine ztgexc(wantq, wantz, n, a, lda, b, ldb, q, ldq, z, ldz, ifst, ilst, info)
ZTGEXC
Definition ztgexc.f:200
zhgeqz
subroutine zhgeqz(job, compq, compz, n, ilo, ihi, h, ldh, t, ldt, alpha, beta, q, ldq, z, ldz, work, lwork, rwork, info)
ZHGEQZ
Definition zhgeqz.f:284
zgges3
subroutine zgges3(jobvsl, jobvsr, sort, selctg, n, a, lda, b, ldb, sdim, alpha, beta, vsl, ldvsl, vsr, ldvsr, work, lwork, rwork, bwork, info)
ZGGES3 computes the eigenvalues, the Schur form, and, optionally, the matrix of Schur vectors for GE ...
Definition zgges3.f:269
zggesx
subroutine zggesx(jobvsl, jobvsr, sort, selctg, sense, n, a, lda, b, ldb, sdim, alpha, beta, vsl, ldvsl, vsr, ldvsr, rconde, rcondv, work, lwork, rwork, iwork, liwork, bwork, info)
ZGGESX computes the eigenvalues, the Schur form, and, optionally, the matrix of Schur vectors for GE ...
Definition zggesx.f:330
zgges
subroutine zgges(jobvsl, jobvsr, sort, selctg, n, a, lda, b, ldb, sdim, alpha, beta, vsl, ldvsl, vsr, ldvsr, work, lwork, rwork, bwork, info)
ZGGES computes the eigenvalues, the Schur form, and, optionally, the matrix of Schur vectors for GE m...
Definition zgges.f:270
zggev
subroutine zggev(jobvl, jobvr, n, a, lda, b, ldb, alpha, beta, vl, ldvl, vr, ldvr, work, lwork, rwork, info)
ZGGEV computes the eigenvalues and, optionally, the left and/or right eigenvectors for GE matrices
Definition zggev.f:217
zggev3
subroutine zggev3(jobvl, jobvr, n, a, lda, b, ldb, alpha, beta, vl, ldvl, vr, ldvr, work, lwork, rwork, info)
ZGGEV3 computes the eigenvalues and, optionally, the left and/or right eigenvectors for GE matrices (...
Definition zggev3.f:216
zggevx
subroutine zggevx(balanc, jobvl, jobvr, sense, n, a, lda, b, ldb, alpha, beta, vl, ldvl, vr, ldvr, ilo, ihi, lscale, rscale, abnrm, bbnrm, rconde, rcondv, work, lwork, rwork, iwork, bwork, info)
ZGGEVX computes the eigenvalues and, optionally, the left and/or right eigenvectors for GE matrices
Definition zggevx.f:374
zggsvd3
subroutine zggsvd3(jobu, jobv, jobq, m, n, p, k, l, a, lda, b, ldb, alpha, beta, u, ldu, v, ldv, q, ldq, work, lwork, rwork, iwork, info)
ZGGSVD3 computes the singular value decomposition (SVD) for OTHER matrices
Definition zggsvd3.f:353
zgghd3
subroutine zgghd3(compq, compz, n, ilo, ihi, a, lda, b, ldb, q, ldq, z, ldz, work, lwork, info)
ZGGHD3
Definition zgghd3.f:227
zgghrd
subroutine zgghrd(compq, compz, n, ilo, ihi, a, lda, b, ldb, q, ldq, z, ldz, info)
ZGGHRD
Definition zgghrd.f:204
zggsvp3
subroutine zggsvp3(jobu, jobv, jobq, m, p, n, a, lda, b, ldb, tola, tolb, k, l, u, ldu, v, ldv, q, ldq, iwork, rwork, tau, work, lwork, info)
ZGGSVP3
Definition zggsvp3.f:278
zuncsd
recursive subroutine zuncsd(jobu1, jobu2, jobv1t, jobv2t, trans, signs, m, p, q, x11, ldx11, x12, ldx12, x21, ldx21, x22, ldx22, theta, u1, ldu1, u2, ldu2, v1t, ldv1t, v2t, ldv2t, work, lwork, rwork, lrwork, iwork, info)
ZUNCSD
Definition zuncsd.f:320
ztgsna
subroutine ztgsna(job, howmny, select, n, a, lda, b, ldb, vl, ldvl, vr, ldvr, s, dif, mm, m, work, lwork, iwork, info)
ZTGSNA
Definition ztgsna.f:311
zggqrf
subroutine zggqrf(n, m, p, a, lda, taua, b, ldb, taub, work, lwork, info)
ZGGQRF
Definition zggqrf.f:215
ztgsja
subroutine ztgsja(jobu, jobv, jobq, m, p, n, k, l, a, lda, b, ldb, tola, tolb, alpha, beta, u, ldu, v, ldv, q, ldq, work, ncycle, info)
ZTGSJA
Definition ztgsja.f:379
ztgsen
subroutine ztgsen(ijob, wantq, wantz, select, n, a, lda, b, ldb, alpha, beta, q, ldq, z, ldz, m, pl, pr, dif, work, lwork, iwork, liwork, info)
ZTGSEN
Definition ztgsen.f:433
zggrqf
subroutine zggrqf(m, p, n, a, lda, taua, b, ldb, taub, work, lwork, info)
ZGGRQF
Definition zggrqf.f:214
zggglm
subroutine zggglm(n, m, p, a, lda, b, ldb, d, x, y, work, lwork, info)
ZGGGLM
Definition zggglm.f:185
zgglse
subroutine zgglse(m, n, p, a, lda, b, ldb, c, d, x, work, lwork, info)
ZGGLSE solves overdetermined or underdetermined systems for OTHER matrices
Definition zgglse.f:180
ztgsyl
subroutine ztgsyl(trans, ijob, m, n, a, lda, b, ldb, c, ldc, d, ldd, e, lde, f, ldf, scale, dif, work, lwork, iwork, info)
ZTGSYL
Definition ztgsyl.f:295
zerrgg
subroutine zerrgg(path, nunit)
ZERRGG
Definition zerrgg.f:57