OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
zerred.f
Go to the documentation of this file.
1*> \brief \b ZERRED
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 ZERRED( 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*> ZERRED tests the error exits for the eigenvalue driver routines for
25*> DOUBLE COMPLEX PRECISION matrices:
26*>
27*> PATH driver description
28*> ---- ------ -----------
29*> ZEV ZGEEV find eigenvalues/eigenvectors for nonsymmetric A
30*> ZES ZGEES find eigenvalues/Schur form for nonsymmetric A
31*> ZVX ZGEEVX ZGEEV + balancing and condition estimation
32*> ZSX ZGEESX ZGEES + balancing and condition estimation
33*> ZBD ZGESVD compute SVD of an M-by-N matrix A
34*> ZGESDD compute SVD of an M-by-N matrix A(by divide and
35*> conquer)
36*> ZGEJSV compute SVD of an M-by-N matrix A where M >= N
37*> ZGESVDX compute SVD of an M-by-N matrix A(by bisection
38*> and inverse iteration)
39*> ZGESVDQ compute SVD of an M-by-N matrix A(with a
40*> QR-Preconditioned )
41*> \endverbatim
42*
43* Arguments:
44* ==========
45*
46*> \param[in] PATH
47*> \verbatim
48*> PATH is CHARACTER*3
49*> The LAPACK path name for the routines to be tested.
50*> \endverbatim
51*>
52*> \param[in] NUNIT
53*> \verbatim
54*> NUNIT is INTEGER
55*> The unit number for output.
56*> \endverbatim
57*
58* Authors:
59* ========
60*
61*> \author Univ. of Tennessee
62*> \author Univ. of California Berkeley
63*> \author Univ. of Colorado Denver
64*> \author NAG Ltd.
65*
66*> \ingroup complex16_eig
67*
68* =====================================================================
69 SUBROUTINE zerred( PATH, NUNIT )
70*
71* -- LAPACK test routine --
72* -- LAPACK is a software package provided by Univ. of Tennessee, --
73* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
74*
75* .. Scalar Arguments ..
76 CHARACTER*3 PATH
77 INTEGER NUNIT
78* ..
79*
80* =====================================================================
81*
82* .. Parameters ..
83 INTEGER NMAX, LW
84 parameter( nmax = 4, lw = 5*nmax )
85 DOUBLE PRECISION ONE, ZERO
86 parameter( one = 1.0d0, zero = 0.0d0 )
87* ..
88* .. Local Scalars ..
89 CHARACTER*2 C2
90 INTEGER I, IHI, ILO, INFO, J, NS, NT, SDIM
91 DOUBLE PRECISION ABNRM
92* ..
93* .. Local Arrays ..
94 LOGICAL B( NMAX )
95 INTEGER IW( 4*NMAX )
96 DOUBLE PRECISION R1( NMAX ), R2( NMAX ), RW( LW ), S( NMAX )
97 COMPLEX*16 A( NMAX, NMAX ), U( NMAX, NMAX ),
98 $ VL( NMAX, NMAX ), VR( NMAX, NMAX ),
99 $ VT( NMAX, NMAX ), W( 10*NMAX ), X( NMAX )
100* ..
101* .. External Subroutines ..
102 EXTERNAL chkxer, zgees, zgeesx, zgeev, zgeevx, zgesvj,
103 $ zgesdd, zgesvd, zgesvdx, zgesvq
104* ..
105* .. External Functions ..
106 LOGICAL LSAMEN, ZSLECT
107 EXTERNAL lsamen, zslect
108* ..
109* .. Intrinsic Functions ..
110 INTRINSIC len_trim
111* ..
112* .. Arrays in Common ..
113 LOGICAL SELVAL( 20 )
114 DOUBLE PRECISION SELWI( 20 ), SELWR( 20 )
115* ..
116* .. Scalars in Common ..
117 LOGICAL LERR, OK
118 CHARACTER*32 SRNAMT
119 INTEGER INFOT, NOUT, SELDIM, SELOPT
120* ..
121* .. Common blocks ..
122 COMMON / infoc / infot, nout, ok, lerr
123 COMMON / srnamc / srnamt
124 COMMON / sslct / selopt, seldim, selval, selwr, selwi
125* ..
126* .. Executable Statements ..
127*
128 nout = nunit
129 WRITE( nout, fmt = * )
130 c2 = path( 2: 3 )
131*
132* Initialize A
133*
134 DO 20 j = 1, nmax
135 DO 10 i = 1, nmax
136 a( i, j ) = zero
137 10 CONTINUE
138 20 CONTINUE
139 DO 30 i = 1, nmax
140 a( i, i ) = one
141 30 CONTINUE
142 ok = .true.
143 nt = 0
144*
145 IF( lsamen( 2, c2, 'EV' ) ) THEN
146*
147* Test ZGEEV
148*
149 srnamt = 'ZGEEV '
150 infot = 1
151 CALL zgeev( 'X', 'N', 0, a, 1, x, vl, 1, vr, 1, w, 1, rw,
152 $ info )
153 CALL chkxer( 'ZGEEV ', infot, nout, lerr, ok )
154 infot = 2
155 CALL zgeev( 'N', 'X', 0, a, 1, x, vl, 1, vr, 1, w, 1, rw,
156 $ info )
157 CALL chkxer( 'ZGEEV ', infot, nout, lerr, ok )
158 infot = 3
159 CALL zgeev( 'N', 'N', -1, a, 1, x, vl, 1, vr, 1, w, 1, rw,
160 $ info )
161 CALL chkxer( 'ZGEEV ', infot, nout, lerr, ok )
162 infot = 5
163 CALL zgeev( 'N', 'N', 2, a, 1, x, vl, 1, vr, 1, w, 4, rw,
164 $ info )
165 CALL chkxer( 'ZGEEV ', infot, nout, lerr, ok )
166 infot = 8
167 CALL zgeev( 'V', 'N', 2, a, 2, x, vl, 1, vr, 1, w, 4, rw,
168 $ info )
169 CALL chkxer( 'ZGEEV ', infot, nout, lerr, ok )
170 infot = 10
171 CALL zgeev( 'N', 'V', 2, a, 2, x, vl, 1, vr, 1, w, 4, rw,
172 $ info )
173 CALL chkxer( 'ZGEEV ', infot, nout, lerr, ok )
174 infot = 12
175 CALL zgeev( 'V', 'V', 1, a, 1, x, vl, 1, vr, 1, w, 1, rw,
176 $ info )
177 CALL chkxer( 'ZGEEV ', infot, nout, lerr, ok )
178 nt = nt + 7
179*
180 ELSE IF( lsamen( 2, c2, 'ES' ) ) THEN
181*
182* Test ZGEES
183*
184 srnamt = 'ZGEES '
185 infot = 1
186 CALL zgees( 'X', 'N', zslect, 0, a, 1, sdim, x, vl, 1, w, 1,
187 $ rw, b, info )
188 CALL chkxer( 'ZGEES ', infot, nout, lerr, ok )
189 infot = 2
190 CALL zgees( 'N', 'X', zslect, 0, a, 1, sdim, x, vl, 1, w, 1,
191 $ rw, b, info )
192 CALL chkxer( 'ZGEES ', infot, nout, lerr, ok )
193 infot = 4
194 CALL zgees( 'N', 'S', zslect, -1, a, 1, sdim, x, vl, 1, w, 1,
195 $ rw, b, info )
196 CALL chkxer( 'ZGEES ', infot, nout, lerr, ok )
197 infot = 6
198 CALL zgees( 'N', 'S', zslect, 2, a, 1, sdim, x, vl, 1, w, 4,
199 $ rw, b, info )
200 CALL chkxer( 'ZGEES ', infot, nout, lerr, ok )
201 infot = 10
202 CALL zgees( 'V', 'S', zslect, 2, a, 2, sdim, x, vl, 1, w, 4,
203 $ rw, b, info )
204 CALL chkxer( 'ZGEES ', infot, nout, lerr, ok )
205 infot = 12
206 CALL zgees( 'N', 'S', zslect, 1, a, 1, sdim, x, vl, 1, w, 1,
207 $ rw, b, info )
208 CALL chkxer( 'ZGEES ', infot, nout, lerr, ok )
209 nt = nt + 6
210*
211 ELSE IF( lsamen( 2, c2, 'VX' ) ) THEN
212*
213* Test ZGEEVX
214*
215 srnamt = 'ZGEEVX'
216 infot = 1
217 CALL zgeevx( 'X', 'N', 'N', 'N', 0, a, 1, x, vl, 1, vr, 1, ilo,
218 $ ihi, s, abnrm, r1, r2, w, 1, rw, info )
219 CALL chkxer( 'ZGEEVX', infot, nout, lerr, ok )
220 infot = 2
221 CALL zgeevx( 'N', 'X', 'N', 'N', 0, a, 1, x, vl, 1, vr, 1, ilo,
222 $ ihi, s, abnrm, r1, r2, w, 1, rw, info )
223 CALL chkxer( 'ZGEEVX', infot, nout, lerr, ok )
224 infot = 3
225 CALL zgeevx( 'N', 'N', 'X', 'N', 0, a, 1, x, vl, 1, vr, 1, ilo,
226 $ ihi, s, abnrm, r1, r2, w, 1, rw, info )
227 CALL chkxer( 'ZGEEVX', infot, nout, lerr, ok )
228 infot = 4
229 CALL zgeevx( 'N', 'N', 'N', 'X', 0, a, 1, x, vl, 1, vr, 1, ilo,
230 $ ihi, s, abnrm, r1, r2, w, 1, rw, info )
231 CALL chkxer( 'ZGEEVX', infot, nout, lerr, ok )
232 infot = 5
233 CALL zgeevx( 'N', 'N', 'N', 'N', -1, a, 1, x, vl, 1, vr, 1,
234 $ ilo, ihi, s, abnrm, r1, r2, w, 1, rw, info )
235 CALL chkxer( 'ZGEEVX', infot, nout, lerr, ok )
236 infot = 7
237 CALL zgeevx( 'n', 'n', 'n', 'n', 2, A, 1, X, VL, 1, VR, 1, ILO,
238 $ IHI, S, ABNRM, R1, R2, W, 4, RW, INFO )
239 CALL CHKXER( 'zgeevx', INFOT, NOUT, LERR, OK )
240 INFOT = 10
241 CALL ZGEEVX( 'n', 'v', 'n', 'n', 2, A, 2, X, VL, 1, VR, 1, ILO,
242 $ IHI, S, ABNRM, R1, R2, W, 4, RW, INFO )
243 CALL CHKXER( 'zgeevx', INFOT, NOUT, LERR, OK )
244 INFOT = 12
245 CALL ZGEEVX( 'n', 'n', 'v', 'n', 2, A, 2, X, VL, 1, VR, 1, ILO,
246 $ IHI, S, ABNRM, R1, R2, W, 4, RW, INFO )
247 CALL CHKXER( 'zgeevx', INFOT, NOUT, LERR, OK )
248 INFOT = 20
249 CALL ZGEEVX( 'n', 'n', 'n', 'n', 1, A, 1, X, VL, 1, VR, 1, ILO,
250 $ IHI, S, ABNRM, R1, R2, W, 1, RW, INFO )
251 CALL CHKXER( 'zgeevx', INFOT, NOUT, LERR, OK )
252 INFOT = 20
253 CALL ZGEEVX( 'n', 'n', 'v', 'v', 1, A, 1, X, VL, 1, VR, 1, ILO,
254 $ IHI, S, ABNRM, R1, R2, W, 2, RW, INFO )
255 CALL CHKXER( 'zgeevx', INFOT, NOUT, LERR, OK )
256 NT = NT + 10
257*
258 ELSE IF( LSAMEN( 2, C2, 'sx' ) ) THEN
259*
260* Test ZGEESX
261*
262 SRNAMT = 'zgeesx'
263 INFOT = 1
264 CALL ZGEESX( 'x', 'n', ZSLECT, 'n', 0, A, 1, SDIM, X, VL, 1,
265 $ R1( 1 ), R2( 1 ), W, 1, RW, B, INFO )
266 CALL CHKXER( 'zgeesx', INFOT, NOUT, LERR, OK )
267 INFOT = 2
268 CALL ZGEESX( 'n', 'x', ZSLECT, 'n', 0, A, 1, SDIM, X, VL, 1,
269 $ R1( 1 ), R2( 1 ), W, 1, RW, B, INFO )
270 CALL CHKXER( 'zgeesx', INFOT, NOUT, LERR, OK )
271 INFOT = 4
272 CALL ZGEESX( 'n', 'n', ZSLECT, 'x', 0, A, 1, SDIM, X, VL, 1,
273 $ R1( 1 ), R2( 1 ), W, 1, RW, B, INFO )
274 CALL CHKXER( 'zgeesx', INFOT, NOUT, LERR, OK )
275 INFOT = 5
276 CALL ZGEESX( 'n', 'n', ZSLECT, 'n', -1, A, 1, SDIM, X, VL, 1,
277 $ R1( 1 ), R2( 1 ), W, 1, RW, B, INFO )
278 CALL CHKXER( 'zgeesx', INFOT, NOUT, LERR, OK )
279 INFOT = 7
280 CALL ZGEESX( 'n', 'n', ZSLECT, 'n', 2, A, 1, SDIM, X, VL, 1,
281 $ R1( 1 ), R2( 1 ), W, 4, RW, B, INFO )
282 CALL CHKXER( 'zgeesx', INFOT, NOUT, LERR, OK )
283 INFOT = 11
284 CALL ZGEESX( 'v', 'n', ZSLECT, 'n', 2, A, 2, SDIM, X, VL, 1,
285 $ R1( 1 ), R2( 1 ), W, 4, RW, B, INFO )
286 CALL CHKXER( 'zgeesx', INFOT, NOUT, LERR, OK )
287 INFOT = 15
288 CALL ZGEESX( 'n', 'n', ZSLECT, 'n', 1, A, 1, SDIM, X, VL, 1,
289 $ R1( 1 ), R2( 1 ), W, 1, RW, B, INFO )
290 CALL CHKXER( 'zgeesx', INFOT, NOUT, LERR, OK )
291 NT = NT + 7
292*
293 ELSE IF( LSAMEN( 2, C2, 'bd' ) ) THEN
294*
295* Test ZGESVD
296*
297 SRNAMT = 'zgesvd'
298 INFOT = 1
299 CALL ZGESVD( 'x', 'n', 0, 0, A, 1, S, U, 1, VT, 1, W, 1, RW,
300 $ INFO )
301 CALL CHKXER( 'zgesvd', INFOT, NOUT, LERR, OK )
302 INFOT = 2
303 CALL ZGESVD( 'n', 'x', 0, 0, A, 1, S, U, 1, VT, 1, W, 1, RW,
304 $ INFO )
305 CALL CHKXER( 'zgesvd', INFOT, NOUT, LERR, OK )
306 INFOT = 2
307 CALL ZGESVD( 'o', 'o', 0, 0, A, 1, S, U, 1, VT, 1, W, 1, RW,
308 $ INFO )
309 CALL CHKXER( 'zgesvd', INFOT, NOUT, LERR, OK )
310 INFOT = 3
311 CALL ZGESVD( 'n', 'n', -1, 0, A, 1, S, U, 1, VT, 1, W, 1, RW,
312 $ INFO )
313 CALL CHKXER( 'zgesvd', INFOT, NOUT, LERR, OK )
314 INFOT = 4
315 CALL ZGESVD( 'n', 'n', 0, -1, A, 1, S, U, 1, VT, 1, W, 1, RW,
316 $ INFO )
317 CALL CHKXER( 'zgesvd', INFOT, NOUT, LERR, OK )
318 INFOT = 6
319 CALL ZGESVD( 'n', 'n', 2, 1, A, 1, S, U, 1, VT, 1, W, 5, RW,
320 $ INFO )
321 CALL CHKXER( 'zgesvd', INFOT, NOUT, LERR, OK )
322 INFOT = 9
323 CALL ZGESVD( 'a', 'n', 2, 1, A, 2, S, U, 1, VT, 1, W, 5, RW,
324 $ INFO )
325 CALL CHKXER( 'zgesvd', INFOT, NOUT, LERR, OK )
326 INFOT = 11
327 CALL ZGESVD( 'n', 'a', 1, 2, A, 1, S, U, 1, VT, 1, W, 5, RW,
328 $ INFO )
329 CALL CHKXER( 'zgesvd', INFOT, NOUT, LERR, OK )
330 NT = NT + 8
331 IF( OK ) THEN
332 WRITE( NOUT, FMT = 9999 )SRNAMT( 1:LEN_TRIM( SRNAMT ) ),
333 $ NT
334 ELSE
335 WRITE( NOUT, FMT = 9998 )
336 END IF
337*
338* Test ZGESDD
339*
340 SRNAMT = 'zgesdd'
341 INFOT = 1
342 CALL ZGESDD( 'x', 0, 0, A, 1, S, U, 1, VT, 1, W, 1, RW, IW,
343 $ INFO )
344 CALL CHKXER( 'zgesdd', INFOT, NOUT, LERR, OK )
345 INFOT = 2
346 CALL ZGESDD( 'n', -1, 0, A, 1, S, U, 1, VT, 1, W, 1, RW, IW,
347 $ INFO )
348 CALL CHKXER( 'zgesdd', INFOT, NOUT, LERR, OK )
349 INFOT = 3
350 CALL ZGESDD( 'n', 0, -1, A, 1, S, U, 1, VT, 1, W, 1, RW, IW,
351 $ INFO )
352 CALL CHKXER( 'zgesdd', INFOT, NOUT, LERR, OK )
353 INFOT = 5
354 CALL ZGESDD( 'n', 2, 1, A, 1, S, U, 1, VT, 1, W, 5, RW, IW,
355 $ INFO )
356 CALL CHKXER( 'zgesdd', INFOT, NOUT, LERR, OK )
357 INFOT = 8
358 CALL ZGESDD( 'a', 2, 1, A, 2, S, U, 1, VT, 1, W, 5, RW, IW,
359 $ INFO )
360 CALL CHKXER( 'zgesdd', INFOT, NOUT, LERR, OK )
361 INFOT = 10
362 CALL ZGESDD( 'a', 1, 2, A, 1, S, U, 1, VT, 1, W, 5, RW, IW,
363 $ INFO )
364 CALL CHKXER( 'zgesdd', INFOT, NOUT, LERR, OK )
365 NT = NT - 2
366 IF( OK ) THEN
367 WRITE( NOUT, FMT = 9999 )SRNAMT( 1:LEN_TRIM( SRNAMT ) ),
368 $ NT
369 ELSE
370 WRITE( NOUT, FMT = 9998 )
371 END IF
372*
373* Test ZGEJSV
374*
375 SRNAMT = 'zgejsv'
376 INFOT = 1
377 CALL ZGEJSV( 'x', 'u', 'v', 'r', 'n', 'n',
378 $ 0, 0, A, 1, S, U, 1, VT, 1,
379 $ W, 1, RW, 1, IW, INFO)
380 CALL CHKXER( 'zgejsv', INFOT, NOUT, LERR, OK )
381 INFOT = 2
382 CALL ZGEJSV( 'g', 'x', 'v', 'r', 'n', 'n',
383 $ 0, 0, A, 1, S, U, 1, VT, 1,
384 $ W, 1, RW, 1, IW, INFO)
385 CALL CHKXER( 'zgejsv', INFOT, NOUT, LERR, OK )
386 INFOT = 3
387 CALL ZGEJSV( 'g', 'u', 'x', 'r', 'n', 'n',
388 $ 0, 0, A, 1, S, U, 1, VT, 1,
389 $ W, 1, RW, 1, IW, INFO)
390 CALL CHKXER( 'zgejsv', INFOT, NOUT, LERR, OK )
391 INFOT = 4
392 CALL ZGEJSV( 'g', 'u', 'v', 'x', 'n', 'n',
393 $ 0, 0, A, 1, S, U, 1, VT, 1,
394 $ W, 1, RW, 1, IW, INFO)
395 CALL CHKXER( 'zgejsv', INFOT, NOUT, LERR, OK )
396 INFOT = 5
397 CALL ZGEJSV( 'g', 'u', 'v', 'r', 'x', 'n',
398 $ 0, 0, A, 1, S, U, 1, VT, 1,
399 $ W, 1, RW, 1, IW, INFO)
400 CALL CHKXER( 'zgejsv', INFOT, NOUT, LERR, OK )
401 INFOT = 6
402 CALL ZGEJSV( 'g', 'u', 'v', 'r', 'n', 'x',
403 $ 0, 0, A, 1, S, U, 1, VT, 1,
404 $ W, 1, RW, 1, IW, INFO)
405 CALL CHKXER( 'zgejsv', INFOT, NOUT, LERR, OK )
406 INFOT = 7
407 CALL ZGEJSV( 'g', 'u', 'v', 'r', 'n', 'n',
408 $ -1, 0, A, 1, S, U, 1, VT, 1,
409 $ W, 1, RW, 1, IW, INFO)
410 CALL CHKXER( 'zgejsv', INFOT, NOUT, LERR, OK )
411 INFOT = 8
412 CALL ZGEJSV( 'g', 'u', 'v', 'r', 'n', 'n',
413 $ 0, -1, A, 1, S, U, 1, VT, 1,
414 $ W, 1, RW, 1, IW, INFO)
415 CALL CHKXER( 'zgejsv', INFOT, NOUT, LERR, OK )
416 INFOT = 10
417 CALL ZGEJSV( 'g', 'u', 'v', 'r', 'N', 'N',
418 $ 2, 1, a, 1, s, u, 1, vt, 1,
419 $ w, 1, rw, 1, iw, info)
420 CALL chkxer( 'ZGEJSV', infot, nout, lerr, ok )
421 infot = 13
422 CALL zgejsv( 'G', 'U', 'V', 'R', 'N', 'N',
423 $ 2, 2, a, 2, s, u, 1, vt, 2,
424 $ w, 1, rw, 1, iw, info)
425 CALL chkxer( 'ZGEJSV', infot, nout, lerr, ok )
426 infot = 15
427 CALL zgejsv( 'G', 'U', 'V', 'r', 'n', 'n',
428 $ 2, 2, A, 2, S, U, 2, VT, 1,
429 $ W, 1, RW, 1, IW, INFO)
430 CALL CHKXER( 'zgejsv', INFOT, NOUT, LERR, OK )
431 NT = 11
432 IF( OK ) THEN
433 WRITE( NOUT, FMT = 9999 )SRNAMT( 1:LEN_TRIM( SRNAMT ) ),
434 $ NT
435 ELSE
436 WRITE( NOUT, FMT = 9998 )
437 END IF
438*
439* Test ZGESVDX
440*
441 SRNAMT = 'zgesvdx'
442 INFOT = 1
443 CALL ZGESVDX( 'x', 'n', 'a', 0, 0, A, 1, ZERO, ZERO,
444 $ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
445 CALL CHKXER( 'zgesvdx', INFOT, NOUT, LERR, OK )
446 INFOT = 2
447 CALL ZGESVDX( 'n', 'x', 'a', 0, 0, A, 1, ZERO, ZERO,
448 $ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
449 CALL CHKXER( 'zgesvdx', INFOT, NOUT, LERR, OK )
450 INFOT = 3
451 CALL ZGESVDX( 'n', 'n', 'x', 0, 0, A, 1, ZERO, ZERO,
452 $ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
453 CALL CHKXER( 'zgesvdx', INFOT, NOUT, LERR, OK )
454 INFOT = 4
455 CALL ZGESVDX( 'n', 'n', 'a', -1, 0, A, 1, ZERO, ZERO,
456 $ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
457 CALL CHKXER( 'zgesvdx', INFOT, NOUT, LERR, OK )
458 INFOT = 5
459 CALL ZGESVDX( 'n', 'n', 'a', 0, -1, A, 1, ZERO, ZERO,
460 $ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
461 CALL CHKXER( 'zgesvdx', INFOT, NOUT, LERR, OK )
462 INFOT = 7
463 CALL ZGESVDX( 'n', 'n', 'a', 2, 1, A, 1, ZERO, ZERO,
464 $ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
465 CALL CHKXER( 'zgesvdx', INFOT, NOUT, LERR, OK )
466 INFOT = 8
467 CALL ZGESVDX( 'n', 'n', 'v', 2, 1, A, 2, -ONE, ZERO,
468 $ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
469 CALL CHKXER( 'zgesvdx', INFOT, NOUT, LERR, OK )
470 INFOT = 9
471 CALL ZGESVDX( 'n', 'n', 'v', 2, 1, A, 2, ONE, ZERO,
472 $ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
473 CALL CHKXER( 'zgesvdx', INFOT, NOUT, LERR, OK )
474 INFOT = 10
475 CALL ZGESVDX( 'n', 'n', 'i', 2, 2, A, 2, ZERO, ZERO,
476 $ 0, 1, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
477 CALL CHKXER( 'zgesvdx', INFOT, NOUT, LERR, OK )
478 INFOT = 11
479 CALL ZGESVDX( 'v', 'n', 'i', 2, 2, A, 2, ZERO, ZERO,
480 $ 1, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
481 CALL CHKXER( 'zgesvdx', INFOT, NOUT, LERR, OK )
482 INFOT = 15
483 CALL ZGESVDX( 'v', 'n', 'a', 2, 2, A, 2, ZERO, ZERO,
484 $ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
485 CALL CHKXER( 'zgesvdx', INFOT, NOUT, LERR, OK )
486 INFOT = 17
487 CALL ZGESVDX( 'n', 'v', 'a', 2, 2, A, 2, ZERO, ZERO,
488 $ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
489 CALL CHKXER( 'zgesvdx', INFOT, NOUT, LERR, OK )
490 NT = 12
491 IF( OK ) THEN
492 WRITE( NOUT, FMT = 9999 )SRNAMT( 1:LEN_TRIM( SRNAMT ) ),
493 $ NT
494 ELSE
495 WRITE( NOUT, FMT = 9998 )
496 END IF
497*
498* Test ZGESVDQ
499*
500 SRNAMT = 'zgesvdq'
501 INFOT = 1
502 CALL ZGESVDQ( 'x', 'p', 't', 'a', 'a', 0, 0, A, 1, S, U,
503 $ 0, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
504 CALL CHKXER( 'zgesvdq', INFOT, NOUT, LERR, OK )
505 INFOT = 2
506 CALL ZGESVDQ( 'a', 'x', 't', 'a', 'a', 0, 0, A, 1, S, U,
507 $ 0, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
508 CALL CHKXER( 'zgesvdq', INFOT, NOUT, LERR, OK )
509 INFOT = 3
510 CALL ZGESVDQ( 'a', 'p', 'x', 'a', 'a', 0, 0, A, 1, S, U,
511 $ 0, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
512 CALL CHKXER( 'zgesvdq', INFOT, NOUT, LERR, OK )
513 INFOT = 4
514 CALL ZGESVDQ( 'a', 'p', 't', 'x', 'a', 0, 0, A, 1, S, U,
515 $ 0, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
516 CALL CHKXER( 'zgesvdq', INFOT, NOUT, LERR, OK )
517 INFOT = 5
518 CALL ZGESVDQ( 'a', 'p', 't', 'a', 'x', 0, 0, A, 1, S, U,
519 $ 0, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
520 CALL CHKXER( 'zgesvdq', INFOT, NOUT, LERR, OK )
521 INFOT = 6
522 CALL ZGESVDQ( 'a', 'p', 't', 'a', 'a', -1, 0, A, 1, S, U,
523 $ 0, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
524 CALL CHKXER( 'zgesvdq', INFOT, NOUT, LERR, OK )
525 INFOT = 7
526 CALL ZGESVDQ( 'a', 'p', 't', 'a', 'a', 0, 1, A, 1, S, U,
527 $ 0, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
528 CALL CHKXER( 'zgesvdq', INFOT, NOUT, LERR, OK )
529 INFOT = 9
530 CALL ZGESVDQ( 'a', 'p', 't', 'a', 'a', 1, 1, A, 0, S, U,
531 $ 0, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
532 CALL CHKXER( 'zgesvdq', INFOT, NOUT, LERR, OK )
533 INFOT = 12
534 CALL ZGESVDQ( 'a', 'p', 't', 'a', 'a', 1, 1, A, 1, S, U,
535 $ -1, VT, 0, NS, IW, 1, W, 1, RW, 1, INFO )
536 CALL CHKXER( 'zgesvdq', INFOT, NOUT, LERR, OK )
537 INFOT = 14
538 CALL ZGESVDQ( 'a', 'p', 't', 'a', 'a', 1, 1, A, 1, S, U,
539 $ 1, VT, -1, NS, IW, 1, W, 1, RW, 1, INFO )
540 CALL CHKXER( 'zgesvdq', INFOT, NOUT, LERR, OK )
541 INFOT = 17
542 CALL ZGESVDQ( 'a', 'p', 't', 'a', 'a', 1, 1, A, 1, S, U,
543 $ 1, VT, 1, NS, IW, -5, W, 1, RW, 1, INFO )
544 CALL CHKXER( 'zgesvdq', INFOT, NOUT, LERR, OK )
545 NT = 11
546 IF( OK ) THEN
547 WRITE( NOUT, FMT = 9999 )SRNAMT( 1:LEN_TRIM( SRNAMT ) ),
548 $ NT
549 ELSE
550 WRITE( NOUT, FMT = 9998 )
551 END IF
552 END IF
553*
554* Print a summary line.
555*
556.NOT. IF( LSAMEN( 2, C2, 'bd' ) ) THEN
557 IF( OK ) THEN
558 WRITE( NOUT, FMT = 9999 )SRNAMT( 1:LEN_TRIM( SRNAMT ) ),
559 $ NT
560 ELSE
561 WRITE( NOUT, FMT = 9998 )
562 END IF
563 END IF
564*
565 9999 FORMAT( 1X, A, ' passed the tests of the error exits(', I3,
566 $ ' tests done)' )
567 9998 FORMAT( ' *** ', A, ' failed the tests of the error exits ***' )
568 RETURN
569*
570* End of ZERRED
571*
572 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
zgesvj
subroutine zgesvj(joba, jobu, jobv, m, n, a, lda, sva, mv, v, ldv, cwork, lwork, rwork, lrwork, info)
ZGESVJ
Definition zgesvj.f:351
zgeesx
subroutine zgeesx(jobvs, sort, select, sense, n, a, lda, sdim, w, vs, ldvs, rconde, rcondv, work, lwork, rwork, bwork, info)
ZGEESX computes the eigenvalues, the Schur form, and, optionally, the matrix of Schur vectors for GE ...
Definition zgeesx.f:239
zgeev
subroutine zgeev(jobvl, jobvr, n, a, lda, w, vl, ldvl, vr, ldvr, work, lwork, rwork, info)
ZGEEV computes the eigenvalues and, optionally, the left and/or right eigenvectors for GE matrices
Definition zgeev.f:180
zgees
subroutine zgees(jobvs, sort, select, n, a, lda, sdim, w, vs, ldvs, work, lwork, rwork, bwork, info)
ZGEES computes the eigenvalues, the Schur form, and, optionally, the matrix of Schur vectors for GE m...
Definition zgees.f:197
zgeevx
subroutine zgeevx(balanc, jobvl, jobvr, sense, n, a, lda, w, vl, ldvl, vr, ldvr, ilo, ihi, scale, abnrm, rconde, rcondv, work, lwork, rwork, info)
ZGEEVX computes the eigenvalues and, optionally, the left and/or right eigenvectors for GE matrices
Definition zgeevx.f:288
zgesdd
subroutine zgesdd(jobz, m, n, a, lda, s, u, ldu, vt, ldvt, work, lwork, rwork, iwork, info)
ZGESDD
Definition zgesdd.f:227
zgesvdx
subroutine zgesvdx(jobu, jobvt, range, m, n, a, lda, vl, vu, il, iu, ns, s, u, ldu, vt, ldvt, work, lwork, rwork, iwork, info)
ZGESVDX computes the singular value decomposition (SVD) for GE matrices
Definition zgesvdx.f:270
zgesvd
subroutine zgesvd(jobu, jobvt, m, n, a, lda, s, u, ldu, vt, ldvt, work, lwork, rwork, info)
ZGESVD computes the singular value decomposition (SVD) for GE matrices
Definition zgesvd.f:214
zgejsv
subroutine zgejsv(joba, jobu, jobv, jobr, jobt, jobp, m, n, a, lda, sva, u, ldu, v, ldv, cwork, lwork, rwork, lrwork, iwork, info)
ZGEJSV
Definition zgejsv.f:569
zgesvdq
subroutine zgesvdq(joba, jobp, jobr, jobu, jobv, m, n, a, lda, s, u, ldu, v, ldv, numrank, iwork, liwork, cwork, lcwork, rwork, lrwork, info)
ZGESVDQ computes the singular value decomposition (SVD) with a QR-Preconditioned QR SVD Method for GE...
Definition zgesvdq.f:413
zerred
subroutine zerred(path, nunit)
ZERRED
Definition zerred.f:70