OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
serrst.f
Go to the documentation of this file.
1*> \brief \b SERRST
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 SERRST( 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*> SERRST tests the error exits for SSYTRD, SORGTR, SORMTR, SSPTRD,
25*> SOPGTR, SOPMTR, SSTEQR, SSTERF, SSTEBZ, SSTEIN, SPTEQR, SSBTRD,
26*> SSYEV, SSYEVX, SSYEVD, SSBEV, SSBEVX, SSBEVD,
27*> SSPEV, SSPEVX, SSPEVD, SSTEV, SSTEVX, SSTEVD, and SSTEDC.
28*> SSYEVD_2STAGE, SSYEVR_2STAGE, SSYEVX_2STAGE,
29*> SSYEV_2STAGE, SSBEV_2STAGE, SSBEVD_2STAGE,
30*> SSBEVX_2STAGE, SSYTRD_2STAGE, SSYTRD_SY2SB,
31*> SSYTRD_SB2ST
32*> \endverbatim
33*
34* Arguments:
35* ==========
36*
37*> \param[in] PATH
38*> \verbatim
39*> PATH is CHARACTER*3
40*> The LAPACK path name for the routines to be tested.
41*> \endverbatim
42*>
43*> \param[in] NUNIT
44*> \verbatim
45*> NUNIT is INTEGER
46*> The unit number for output.
47*> \endverbatim
48*
49* Authors:
50* ========
51*
52*> \author Univ. of Tennessee
53*> \author Univ. of California Berkeley
54*> \author Univ. of Colorado Denver
55*> \author NAG Ltd.
56*
57*> \ingroup single_eig
58*
59* =====================================================================
60 SUBROUTINE serrst( PATH, NUNIT )
61*
62* -- LAPACK test routine --
63* -- LAPACK is a software package provided by Univ. of Tennessee, --
64* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
65*
66* .. Scalar Arguments ..
67 CHARACTER*3 PATH
68 INTEGER NUNIT
69* ..
70*
71* =====================================================================
72*
73* NMAX has to be at least 3 or LIW may be too small
74* .. Parameters ..
75 INTEGER NMAX, LIW, LW
76 parameter( nmax = 3, liw = 12*nmax, lw = 20*nmax )
77* ..
78* .. Local Scalars ..
79 CHARACTER*2 C2
80 INTEGER I, INFO, J, M, N, NSPLIT, NT
81* ..
82* .. Local Arrays ..
83 INTEGER I1( NMAX ), I2( NMAX ), I3( NMAX ), IW( LIW )
84 REAL A( NMAX, NMAX ), C( NMAX, NMAX ), D( NMAX ),
85 $ E( NMAX ), Q( NMAX, NMAX ), R( NMAX ),
86 $ TAU( NMAX ), W( LW ), X( NMAX ),
87 $ Z( NMAX, NMAX )
88* ..
89* .. External Functions ..
90 LOGICAL LSAMEN
91 EXTERNAL lsamen
92* ..
93* .. External Subroutines ..
94 EXTERNAL chkxer, sopgtr, sopmtr, sorgtr, sormtr, spteqr,
95 $ ssbev, ssbevd, ssbevx, ssbtrd, sspev, sspevd,
96 $ sspevx, ssptrd, sstebz, sstedc, sstein, ssteqr,
97 $ ssterf, sstev, sstevd, sstevr, sstevx, ssyev,
98 $ ssyevd, ssyevr, ssyevx, ssytrd,
99 $ ssyevd_2stage, ssyevr_2stage, ssyevx_2stage,
100 $ ssyev_2stage, ssbev_2stage, ssbevd_2stage,
101 $ ssbevx_2stage, ssytrd_2stage, ssytrd_sy2sb,
102 $ ssytrd_sb2st
103* ..
104* .. Scalars in Common ..
105 LOGICAL LERR, OK
106 CHARACTER*32 SRNAMT
107 INTEGER INFOT, NOUT
108* ..
109* .. Common blocks ..
110 COMMON / infoc / infot, nout, ok, lerr
111 COMMON / srnamc / srnamt
112* ..
113* .. Intrinsic Functions ..
114 INTRINSIC real
115* ..
116* .. Executable Statements ..
117*
118 nout = nunit
119 WRITE( nout, fmt = * )
120 c2 = path( 2: 3 )
121*
122* Set the variables to innocuous values.
123*
124 DO 20 j = 1, nmax
125 DO 10 i = 1, nmax
126 a( i, j ) = 1. / real( i+j )
127 10 CONTINUE
128 20 CONTINUE
129 DO 30 j = 1, nmax
130 d( j ) = real( j )
131 e( j ) = 0.0
132 i1( j ) = j
133 i2( j ) = j
134 tau( j ) = 1.
135 30 CONTINUE
136 ok = .true.
137 nt = 0
138*
139* Test error exits for the ST path.
140*
141 IF( lsamen( 2, c2, 'ST' ) ) THEN
142*
143* SSYTRD
144*
145 srnamt = 'SSYTRD'
146 infot = 1
147 CALL ssytrd( '/', 0, a, 1, d, e, tau, w, 1, info )
148 CALL chkxer( 'SSYTRD', infot, nout, lerr, ok )
149 infot = 2
150 CALL ssytrd( 'U', -1, a, 1, d, e, tau, w, 1, info )
151 CALL chkxer( 'SSYTRD', infot, nout, lerr, ok )
152 infot = 4
153 CALL ssytrd( 'U', 2, a, 1, d, e, tau, w, 1, info )
154 CALL chkxer( 'SSYTRD', infot, nout, lerr, ok )
155 infot = 9
156 CALL ssytrd( 'U', 0, a, 1, d, e, tau, w, 0, info )
157 CALL chkxer( 'SSYTRD', infot, nout, lerr, ok )
158 nt = nt + 4
159*
160* SSYTRD_2STAGE
161*
162 srnamt = 'SSYTRD_2STAGE'
163 infot = 1
164 CALL ssytrd_2stage( '/', 'U', 0, a, 1, d, e, tau,
165 $ c, 1, w, 1, info )
166 CALL chkxer( 'SSYTRD_2STAGE', infot, nout, lerr, ok )
167 infot = 1
168 CALL ssytrd_2stage( 'H', 'U', 0, a, 1, d, e, tau,
169 $ c, 1, w, 1, info )
170 CALL chkxer( 'SSYTRD_2STAGE', infot, nout, lerr, ok )
171 infot = 2
172 CALL ssytrd_2stage( 'N', '/', 0, a, 1, d, e, tau,
173 $ c, 1, w, 1, info )
174 CALL chkxer( 'SSYTRD_2STAGE', infot, nout, lerr, ok )
175 infot = 3
176 CALL ssytrd_2stage( 'N', 'U', -1, a, 1, d, e, tau,
177 $ c, 1, w, 1, info )
178 CALL chkxer( 'SSYTRD_2STAGE', infot, nout, lerr, ok )
179 infot = 5
180 CALL ssytrd_2stage( 'N', 'U', 2, a, 1, d, e, tau,
181 $ c, 1, w, 1, info )
182 CALL chkxer( 'SSYTRD_2STAGE', infot, nout, lerr, ok )
183 infot = 10
184 CALL ssytrd_2stage( 'N', 'U', 0, a, 1, d, e, tau,
185 $ c, 0, w, 1, info )
186 CALL chkxer( 'SSYTRD_2STAGE', infot, nout, lerr, ok )
187 infot = 12
188 CALL ssytrd_2stage( 'N', 'U', 0, a, 1, d, e, tau,
189 $ c, 1, w, 0, info )
190 CALL chkxer( 'SSYTRD_2STAGE', infot, nout, lerr, ok )
191 nt = nt + 7
192*
193* SSYTRD_SY2SB
194*
195 srnamt = 'SSYTRD_SY2SB'
196 infot = 1
197 CALL ssytrd_sy2sb( '/', 0, 0, a, 1, c, 1, tau, w, 1, info )
198 CALL chkxer( 'SSYTRD_SY2SB', infot, nout, lerr, ok )
199 infot = 2
200 CALL ssytrd_sy2sb( 'U', -1, 0, a, 1, c, 1, tau, w, 1, info )
201 CALL chkxer( 'SSYTRD_SY2SB', infot, nout, lerr, ok )
202 infot = 3
203 CALL ssytrd_sy2sb( 'U', 0, -1, a, 1, c, 1, tau, w, 1, info )
204 CALL chkxer( 'SSYTRD_SY2SB', infot, nout, lerr, ok )
205 infot = 5
206 CALL ssytrd_sy2sb( 'U', 2, 0, a, 1, c, 1, tau, w, 1, info )
207 CALL chkxer( 'SSYTRD_SY2SB', infot, nout, lerr, ok )
208 infot = 7
209 CALL ssytrd_sy2sb( 'U', 0, 2, a, 1, c, 1, tau, w, 1, info )
210 CALL chkxer( 'SSYTRD_SY2SB', infot, nout, lerr, ok )
211 infot = 10
212 CALL ssytrd_sy2sb( 'U', 0, 0, a, 1, c, 1, tau, w, 0, info )
213 CALL chkxer( 'SSYTRD_SY2SB', infot, nout, lerr, ok )
214 nt = nt + 6
215*
216* SSYTRD_SB2ST
217*
218 srnamt = 'SSYTRD_SB2ST'
219 infot = 1
220 CALL ssytrd_sb2st( '/', 'N', 'U', 0, 0, a, 1, d, e,
221 $ c, 1, w, 1, info )
222 CALL chkxer( 'SSYTRD_SB2ST', infot, nout, lerr, ok )
223 infot = 2
224 CALL ssytrd_sb2st( 'Y', '/', 'U', 0, 0, a, 1, d, e,
225 $ c, 1, w, 1, info )
226 CALL chkxer( 'SSYTRD_SB2ST', infot, nout, lerr, ok )
227 infot = 2
228 CALL ssytrd_sb2st( 'Y', 'H', 'U', 0, 0, a, 1, d, e,
229 $ c, 1, w, 1, info )
230 CALL chkxer( 'SSYTRD_SB2ST', infot, nout, lerr, ok )
231 infot = 3
232 CALL ssytrd_sb2st( 'Y', 'N', '/', 0, 0, a, 1, d, e,
233 $ c, 1, w, 1, info )
234 CALL chkxer( 'SSYTRD_SB2ST', infot, nout, lerr, ok )
235 infot = 4
236 CALL ssytrd_sb2st( 'Y', 'N', 'U', -1, 0, a, 1, d, e,
237 $ c, 1, w, 1, info )
238 CALL chkxer( 'SSYTRD_SB2ST', infot, nout, lerr, ok )
239 infot = 5
240 CALL ssytrd_sb2st( 'Y', 'N', 'U', 0, -1, a, 1, d, e,
241 $ c, 1, w, 1, info )
242 CALL chkxer( 'SSYTRD_SB2ST', infot, nout, lerr, ok )
243 infot = 7
244 CALL ssytrd_sb2st( 'Y', 'N', 'U', 0, 1, a, 1, d, e,
245 $ c, 1, w, 1, info )
246 CALL chkxer( 'SSYTRD_SB2ST', infot, nout, lerr, ok )
247 infot = 11
248 CALL ssytrd_sb2st( 'Y', 'N', 'U', 0, 0, a, 1, d, e,
249 $ c, 0, w, 1, info )
250 CALL chkxer( 'SSYTRD_SB2ST', infot, nout, lerr, ok )
251 infot = 13
252 CALL ssytrd_sb2st( 'Y', 'N', 'U', 0, 0, a, 1, d, e,
253 $ c, 1, w, 0, info )
254 CALL chkxer( 'SSYTRD_SB2ST', infot, nout, lerr, ok )
255 nt = nt + 9
256*
257* SORGTR
258*
259 srnamt = 'sorgtr'
260 INFOT = 1
261 CALL SORGTR( '/', 0, A, 1, TAU, W, 1, INFO )
262 CALL CHKXER( 'sorgtr', INFOT, NOUT, LERR, OK )
263 INFOT = 2
264 CALL SORGTR( 'u', -1, A, 1, TAU, W, 1, INFO )
265 CALL CHKXER( 'sorgtr', INFOT, NOUT, LERR, OK )
266 INFOT = 4
267 CALL SORGTR( 'u', 2, A, 1, TAU, W, 1, INFO )
268 CALL CHKXER( 'sorgtr', INFOT, NOUT, LERR, OK )
269 INFOT = 7
270 CALL SORGTR( 'u', 3, A, 3, TAU, W, 1, INFO )
271 CALL CHKXER( 'sorgtr', INFOT, NOUT, LERR, OK )
272 NT = NT + 4
273*
274* SORMTR
275*
276 SRNAMT = 'sormtr'
277 INFOT = 1
278 CALL SORMTR( '/', 'u', 'n', 0, 0, A, 1, TAU, C, 1, W, 1, INFO )
279 CALL CHKXER( 'sormtr', INFOT, NOUT, LERR, OK )
280 INFOT = 2
281 CALL SORMTR( 'l', '/', 'n', 0, 0, A, 1, TAU, C, 1, W, 1, INFO )
282 CALL CHKXER( 'sormtr', INFOT, NOUT, LERR, OK )
283 INFOT = 3
284 CALL SORMTR( 'l', 'u', '/', 0, 0, A, 1, TAU, C, 1, W, 1, INFO )
285 CALL CHKXER( 'sormtr', INFOT, NOUT, LERR, OK )
286 INFOT = 4
287 CALL SORMTR( 'l', 'u', 'n', -1, 0, a, 1, tau, c, 1, w, 1,
288 $ info )
289 CALL chkxer( 'SORMTR', infot, nout, lerr, ok )
290 infot = 5
291 CALL sormtr( 'L', 'U', 'N', 0, -1, a, 1, tau, c, 1, w, 1,
292 $ info )
293 CALL chkxer( 'SORMTR', infot, nout, lerr, ok )
294 infot = 7
295 CALL sormtr( 'L', 'U', 'N', 2, 0, a, 1, tau, c, 2, w, 1, info )
296 CALL chkxer( 'SORMTR', infot, nout, lerr, ok )
297 infot = 7
298 CALL sormtr( 'R', 'U', 'N', 0, 2, a, 1, tau, c, 1, w, 1, info )
299 CALL chkxer( 'SORMTR', infot, nout, lerr, ok )
300 infot = 10
301 CALL sormtr( 'L', 'U', 'N', 2, 0, a, 2, tau, c, 1, w, 1, info )
302 CALL chkxer( 'SORMTR', infot, nout, lerr, ok )
303 infot = 12
304 CALL sormtr( 'L', 'U', 'N', 0, 2, a, 1, tau, c, 1, w, 1, info )
305 CALL chkxer( 'SORMTR', infot, nout, lerr, ok )
306 infot = 12
307 CALL sormtr( 'R', 'U', 'N', 2, 0, a, 1, tau, c, 2, w, 1, info )
308 CALL chkxer( 'SORMTR', infot, nout, lerr, ok )
309 nt = nt + 10
310*
311* SSPTRD
312*
313 srnamt = 'SSPTRD'
314 infot = 1
315 CALL ssptrd( '/', 0, a, d, e, tau, info )
316 CALL chkxer( 'SSPTRD', infot, nout, lerr, ok )
317 infot = 2
318 CALL ssptrd( 'U', -1, a, d, e, tau, info )
319 CALL chkxer( 'SSPTRD', infot, nout, lerr, ok )
320 nt = nt + 2
321*
322* SOPGTR
323*
324 srnamt = 'SOPGTR'
325 infot = 1
326 CALL sopgtr( '/', 0, A, TAU, Z, 1, W, INFO )
327 CALL CHKXER( 'sopgtr', INFOT, NOUT, LERR, OK )
328 INFOT = 2
329 CALL SOPGTR( 'u', -1, A, TAU, Z, 1, W, INFO )
330 CALL CHKXER( 'sopgtr', INFOT, NOUT, LERR, OK )
331 INFOT = 6
332 CALL SOPGTR( 'u', 2, A, TAU, Z, 1, W, INFO )
333 CALL CHKXER( 'sopgtr', INFOT, NOUT, LERR, OK )
334 NT = NT + 3
335*
336* SOPMTR
337*
338 SRNAMT = 'sopmtr'
339 INFOT = 1
340 CALL SOPMTR( '/', 'u', 'n', 0, 0, A, TAU, C, 1, W, INFO )
341 CALL CHKXER( 'sopmtr', INFOT, NOUT, LERR, OK )
342 INFOT = 2
343 CALL SOPMTR( 'l', '/', 'n', 0, 0, A, TAU, C, 1, W, INFO )
344 CALL CHKXER( 'sopmtr', INFOT, NOUT, LERR, OK )
345 INFOT = 3
346 CALL SOPMTR( 'l', 'u', '/', 0, 0, A, TAU, C, 1, W, INFO )
347 CALL CHKXER( 'sopmtr', INFOT, NOUT, LERR, OK )
348 INFOT = 4
349 CALL SOPMTR( 'l', 'u', 'n', -1, 0, A, TAU, C, 1, W, INFO )
350 CALL CHKXER( 'sopmtr', INFOT, NOUT, LERR, OK )
351 INFOT = 5
352 CALL SOPMTR( 'l', 'u', 'n', 0, -1, A, TAU, C, 1, W, INFO )
353 CALL CHKXER( 'sopmtr', INFOT, NOUT, LERR, OK )
354 INFOT = 9
355 CALL SOPMTR( 'l', 'u', 'n', 2, 0, A, TAU, C, 1, W, INFO )
356 CALL CHKXER( 'sopmtr', INFOT, NOUT, LERR, OK )
357 NT = NT + 6
358*
359* SPTEQR
360*
361 SRNAMT = 'spteqr'
362 INFOT = 1
363 CALL SPTEQR( '/', 0, D, E, Z, 1, W, INFO )
364 CALL CHKXER( 'spteqr', INFOT, NOUT, LERR, OK )
365 INFOT = 2
366 CALL SPTEQR( 'n', -1, D, E, Z, 1, W, INFO )
367 CALL CHKXER( 'spteqr', INFOT, NOUT, LERR, OK )
368 INFOT = 6
369 CALL SPTEQR( 'v', 2, D, E, Z, 1, W, INFO )
370 CALL CHKXER( 'spteqr', INFOT, NOUT, LERR, OK )
371 NT = NT + 3
372*
373* SSTEBZ
374*
375 SRNAMT = 'sstebz'
376 INFOT = 1
377 CALL SSTEBZ( '/', 'e', 0, 0.0, 1.0, 1, 0, 0.0, D, E, M, NSPLIT,
378 $ X, I1, I2, W, IW, INFO )
379 CALL CHKXER( 'sstebz', INFOT, NOUT, LERR, OK )
380 INFOT = 2
381 CALL SSTEBZ( 'a', '/', 0, 0.0, 0.0, 0, 0, 0.0, D, E, M, NSPLIT,
382 $ X, I1, I2, W, IW, INFO )
383 CALL CHKXER( 'sstebz', INFOT, NOUT, LERR, OK )
384 INFOT = 3
385 CALL SSTEBZ( 'a', 'e', -1, 0.0, 0.0, 0, 0, 0.0, D, E, M,
386 $ NSPLIT, X, I1, I2, W, IW, INFO )
387 CALL CHKXER( 'sstebz', INFOT, NOUT, LERR, OK )
388 INFOT = 5
389 CALL SSTEBZ( 'v', 'e', 0, 0.0, 0.0, 0, 0, 0.0, D, E, M, NSPLIT,
390 $ X, I1, I2, W, IW, INFO )
391 CALL CHKXER( 'sstebz', INFOT, NOUT, LERR, OK )
392 INFOT = 6
393 CALL SSTEBZ( 'i', 'e', 0, 0.0, 0.0, 0, 0, 0.0, D, E, M, NSPLIT,
394 $ X, I1, I2, W, IW, INFO )
395 CALL CHKXER( 'sstebz', INFOT, NOUT, LERR, OK )
396 INFOT = 6
397 CALL SSTEBZ( 'i', 'e', 1, 0.0, 0.0, 2, 1, 0.0, D, E, M, NSPLIT,
398 $ X, I1, I2, W, IW, INFO )
399 CALL CHKXER( 'sstebz', INFOT, NOUT, LERR, OK )
400 INFOT = 7
401 CALL SSTEBZ( 'i', 'e', 1, 0.0, 0.0, 1, 0, 0.0, D, E, M, NSPLIT,
402 $ X, I1, I2, W, IW, INFO )
403 CALL CHKXER( 'sstebz', INFOT, NOUT, LERR, OK )
404 INFOT = 7
405 CALL SSTEBZ( 'i', 'e', 1, 0.0, 0.0, 1, 2, 0.0, D, E, M, NSPLIT,
406 $ X, I1, I2, W, IW, INFO )
407 CALL CHKXER( 'sstebz', INFOT, NOUT, LERR, OK )
408 NT = NT + 8
409*
410* SSTEIN
411*
412 SRNAMT = 'sstein'
413 INFOT = 1
414 CALL SSTEIN( -1, D, E, 0, X, I1, I2, Z, 1, W, IW, I3, INFO )
415 CALL CHKXER( 'sstein', INFOT, NOUT, LERR, OK )
416 INFOT = 4
417 CALL SSTEIN( 0, D, E, -1, X, I1, I2, Z, 1, W, IW, I3, INFO )
418 CALL CHKXER( 'sstein', INFOT, NOUT, LERR, OK )
419 INFOT = 4
420 CALL SSTEIN( 0, D, E, 1, X, I1, I2, Z, 1, W, IW, I3, INFO )
421 CALL CHKXER( 'sstein', INFOT, NOUT, LERR, OK )
422 INFOT = 9
423 CALL SSTEIN( 2, D, E, 0, X, I1, I2, Z, 1, W, IW, I3, INFO )
424 CALL CHKXER( 'sstein', INFOT, NOUT, LERR, OK )
425 NT = NT + 4
426*
427* SSTEQR
428*
429 SRNAMT = 'ssteqr'
430 INFOT = 1
431 CALL SSTEQR( '/', 0, D, E, Z, 1, W, INFO )
432 CALL CHKXER( 'ssteqr', INFOT, NOUT, LERR, OK )
433 INFOT = 2
434 CALL SSTEQR( 'n', -1, D, E, Z, 1, W, INFO )
435 CALL CHKXER( 'ssteqr', INFOT, NOUT, LERR, OK )
436 INFOT = 6
437 CALL SSTEQR( 'v', 2, D, E, Z, 1, W, INFO )
438 CALL CHKXER( 'ssteqr', INFOT, NOUT, LERR, OK )
439 NT = NT + 3
440*
441* SSTERF
442*
443 SRNAMT = 'ssterf'
444 INFOT = 1
445 CALL SSTERF( -1, D, E, INFO )
446 CALL CHKXER( 'ssterf', INFOT, NOUT, LERR, OK )
447 NT = NT + 1
448*
449* SSTEDC
450*
451 SRNAMT = 'sstedc'
452 INFOT = 1
453 CALL SSTEDC( '/', 0, D, E, Z, 1, W, 1, IW, 1, INFO )
454 CALL CHKXER( 'sstedc', INFOT, NOUT, LERR, OK )
455 INFOT = 2
456 CALL SSTEDC( 'n', -1, D, E, Z, 1, W, 1, IW, 1, INFO )
457 CALL CHKXER( 'sstedc', INFOT, NOUT, LERR, OK )
458 INFOT = 6
459 CALL SSTEDC( 'v', 2, D, E, Z, 1, W, 23, IW, 28, INFO )
460 CALL CHKXER( 'sstedc', INFOT, NOUT, LERR, OK )
461 INFOT = 8
462 CALL SSTEDC( 'n', 1, D, E, Z, 1, W, 0, IW, 1, INFO )
463 CALL CHKXER( 'sstedc', INFOT, NOUT, LERR, OK )
464 INFOT = 8
465 CALL SSTEDC( 'i', 2, D, E, Z, 2, W, 0, IW, 12, INFO )
466 CALL CHKXER( 'sstedc', INFOT, NOUT, LERR, OK )
467 INFOT = 8
468 CALL SSTEDC( 'v', 2, D, E, Z, 2, W, 0, IW, 28, INFO )
469 CALL CHKXER( 'sstedc', INFOT, NOUT, LERR, OK )
470 INFOT = 10
471 CALL SSTEDC( 'n', 1, D, E, Z, 1, W, 1, IW, 0, INFO )
472 CALL CHKXER( 'sstedc', INFOT, NOUT, LERR, OK )
473 INFOT = 10
474 CALL SSTEDC( 'i', 2, D, E, Z, 2, W, 19, IW, 0, INFO )
475 CALL CHKXER( 'sstedc', INFOT, NOUT, LERR, OK )
476 INFOT = 10
477 CALL SSTEDC( 'v', 2, D, E, Z, 2, W, 23, IW, 0, INFO )
478 CALL CHKXER( 'sstedc', INFOT, NOUT, LERR, OK )
479 NT = NT + 9
480*
481* SSTEVD
482*
483 SRNAMT = 'sstevd'
484 INFOT = 1
485 CALL SSTEVD( '/', 0, D, E, Z, 1, W, 1, IW, 1, INFO )
486 CALL CHKXER( 'sstevd', INFOT, NOUT, LERR, OK )
487 INFOT = 2
488 CALL SSTEVD( 'n', -1, D, E, Z, 1, W, 1, IW, 1, INFO )
489 CALL CHKXER( 'sstevd', INFOT, NOUT, LERR, OK )
490 INFOT = 6
491 CALL SSTEVD( 'v', 2, D, E, Z, 1, W, 19, IW, 12, INFO )
492 CALL CHKXER( 'sstevd', INFOT, NOUT, LERR, OK )
493 INFOT = 8
494 CALL SSTEVD( 'n', 1, D, E, Z, 1, W, 0, IW, 1, INFO )
495 CALL CHKXER( 'sstevd', INFOT, NOUT, LERR, OK )
496 INFOT = 8
497 CALL SSTEVD( 'v', 2, D, E, Z, 2, W, 12, IW, 12, INFO )
498 CALL CHKXER( 'sstevd', INFOT, NOUT, LERR, OK )
499 INFOT = 10
500 CALL SSTEVD( 'n', 0, D, E, Z, 1, W, 1, IW, 0, INFO )
501 CALL CHKXER( 'sstevd', INFOT, NOUT, LERR, OK )
502 INFOT = 10
503 CALL SSTEVD( 'v', 2, D, E, Z, 2, W, 19, IW, 11, INFO )
504 CALL CHKXER( 'sstevd', INFOT, NOUT, LERR, OK )
505 NT = NT + 7
506*
507* SSTEV
508*
509 SRNAMT = 'sstev '
510 INFOT = 1
511 CALL SSTEV( '/', 0, d, e, z, 1, w, info )
512 CALL chkxer( 'SSTEV ', infot, nout, lerr, ok )
513 infot = 2
514 CALL sstev( 'N', -1, d, e, z, 1, w, info )
515 CALL chkxer( 'SSTEV ', infot, nout, lerr, ok )
516 infot = 6
517 CALL sstev( 'V', 2, d, e, z, 1, w, info )
518 CALL chkxer( 'sstev ', INFOT, NOUT, LERR, OK )
519 NT = NT + 3
520*
521* SSTEVX
522*
523 SRNAMT = 'sstevx'
524 infot = 1
525 CALL sstevx( '/', 'A', 0, d, e, 0.0, 0.0, 0, 0, 0.0, m, x, z,
526 $ 1, w, iw, i3, info )
527 CALL chkxer( 'SSTEVX', infot, nout, lerr, ok )
528 infot = 2
529 CALL sstevx( 'N', '/', 0, d, e, 0.0, 1.0, 1, 0, 0.0, m, x, z,
530 $ 1, w, iw, i3, info )
531 CALL chkxer( 'SSTEVX', infot, nout, lerr, ok )
532 infot = 3
533 CALL sstevx( 'N', 'A', -1, d, e, 0.0, 0.0, 0, 0, 0.0, m, x, z,
534 $ 1, w, iw, i3, info )
535 CALL chkxer( 'SSTEVX', infot, nout, lerr, ok )
536 infot = 7
537 CALL sstevx( 'N', 'V', 1, d, e, 0.0, 0.0, 0, 0, 0.0, m, x, z,
538 $ 1, w, iw, i3, info )
539 CALL chkxer( 'SSTEVX', infot, nout, lerr, ok )
540 infot = 8
541 CALL sstevx( 'N', 'I', 1, d, e, 0.0, 0.0, 0, 0, 0.0, m, x, z,
542 $ 1, w, iw, i3, info )
543 CALL chkxer( 'SSTEVX', infot, nout, lerr, ok )
544 infot = 8
545 CALL sstevx( 'N', 'I', 1, d, e, 0.0, 0.0, 2, 1, 0.0, m, x, z,
546 $ 1, w, iw, i3, info )
547 CALL chkxer( 'SSTEVX', infot, nout, lerr, ok )
548 infot = 9
549 CALL sstevx( 'N', 'I', 2, d, e, 0.0, 0.0, 2, 1, 0.0, m, x, z,
550 $ 1, w, iw, i3, info )
551 CALL chkxer( 'sstevx', INFOT, NOUT, LERR, OK )
552 INFOT = 9
553 CALL SSTEVX( 'n', 'i', 1, D, E, 0.0, 0.0, 1, 2, 0.0, M, X, Z,
554 $ 1, W, IW, I3, INFO )
555 CALL CHKXER( 'sstevx', INFOT, NOUT, LERR, OK )
556 INFOT = 14
557 CALL SSTEVX( 'v', 'a', 2, D, E, 0.0, 0.0, 0, 0, 0.0, M, X, Z,
558 $ 1, W, IW, I3, INFO )
559 CALL CHKXER( 'sstevx', INFOT, NOUT, LERR, OK )
560 NT = NT + 9
561*
562* SSTEVR
563*
564 N = 1
565 SRNAMT = 'sstevr'
566 INFOT = 1
567 CALL SSTEVR( '/', 'a', 0, D, E, 0.0, 0.0, 1, 1, 0.0, M, R, Z,
568 $ 1, IW, X, 20*N, IW( 2*N+1 ), 10*N, INFO )
569 CALL CHKXER( 'sstevr', INFOT, NOUT, LERR, OK )
570 INFOT = 2
571 CALL SSTEVR( 'v', '/', 0, D, E, 0.0, 0.0, 1, 1, 0.0, M, R, Z,
572 $ 1, IW, X, 20*N, IW( 2*N+1 ), 10*N, INFO )
573 CALL CHKXER( 'sstevr', INFOT, NOUT, LERR, OK )
574 INFOT = 3
575 CALL SSTEVR( 'v', 'a', -1, D, E, 0.0, 0.0, 1, 1, 0.0, M, R, Z,
576 $ 1, IW, X, 20*N, IW( 2*N+1 ), 10*N, INFO )
577 CALL CHKXER( 'sstevr', INFOT, NOUT, LERR, OK )
578 INFOT = 7
579 CALL SSTEVR( 'v', 'v', 1, D, E, 0.0, 0.0, 1, 1, 0.0, M, R, Z,
580 $ 1, IW, X, 20*N, IW( 2*N+1 ), 10*N, INFO )
581 CALL CHKXER( 'sstevr', INFOT, NOUT, LERR, OK )
582 INFOT = 8
583 CALL SSTEVR( 'v', 'i', 1, D, E, 0.0, 0.0, 0, 1, 0.0, M, W, Z,
584 $ 1, IW, X, 20*N, IW( 2*N+1 ), 10*N, INFO )
585 CALL CHKXER( 'sstevr', INFOT, NOUT, LERR, OK )
586 INFOT = 9
587 N = 2
588 CALL SSTEVR( 'v', 'i', 2, D, E, 0.0, 0.0, 2, 1, 0.0, M, W, Z,
589 $ 1, IW, X, 20*N, IW( 2*N+1 ), 10*N, INFO )
590 CALL CHKXER( 'sstevr', infot, nout, lerr, ok )
591 infot = 14
592 n = 1
593 CALL sstevr( 'V', 'I', 1, d, e, 0.0, 0.0, 1, 1, 0.0, m, w, z,
594 $ 0, iw, x, 20*n, iw( 2*n+1 ), 10*n, info )
595 CALL chkxer( 'SSTEVR', infot, nout, lerr, ok )
596 infot = 17
597 CALL sstevr( 'V', 'I', 1, d, e, 0.0, 0.0, 1, 1, 0.0, m, w, z,
598 $ 1, iw, x, 20*n-1, iw( 2*n+1 ), 10*n, info )
599 CALL chkxer( 'SSTEVR', infot, nout, lerr, ok )
600 infot = 19
601 CALL sstevr( 'V', 'I', 1, d, e, 0.0, 0.0, 1, 1, 0.0, m, w, z,
602 $ 1, iw, x, 20*n, iw( 2*n+1 ), 10*n-1, info )
603 CALL chkxer( 'SSTEVR', infot, nout, lerr, ok )
604 nt = nt + 9
605*
606* SSYEVD
607*
608 srnamt = 'SSYEVD'
609 infot = 1
610 CALL ssyevd( '/', 'U', 0, a, 1, x, w, 1, iw, 1, info )
611 CALL chkxer( 'SSYEVD', infot, nout, lerr, ok )
612 infot = 2
613 CALL ssyevd( 'N', '/', 0, a, 1, x, w, 1, iw, 1, info )
614 CALL chkxer( 'SSYEVD', infot, nout, lerr, ok )
615 infot = 3
616 CALL ssyevd( 'N', 'U', -1, a, 1, x, w, 1, iw, 1, info )
617 CALL chkxer( 'SSYEVD', infot, nout, lerr, ok )
618 infot = 5
619 CALL ssyevd( 'N', 'U', 2, a, 1, x, w, 3, iw, 1, info )
620 CALL chkxer( 'SSYEVD', infot, nout, lerr, ok )
621 infot = 8
622 CALL ssyevd( 'N', 'U', 1, a, 1, x, w, 0, iw, 1, info )
623 CALL chkxer( 'SSYEVD', infot, nout, lerr, ok )
624 infot = 8
625 CALL ssyevd( 'N', 'U', 2, a, 2, x, w, 4, iw, 1, info )
626 CALL chkxer( 'SSYEVD', infot, nout, lerr, ok )
627 infot = 8
628 CALL ssyevd( 'V', 'U', 2, a, 2, x, w, 20, iw, 12, info )
629 CALL chkxer( 'SSYEVD', infot, nout, lerr, ok )
630 infot = 10
631 CALL ssyevd( 'N', 'U', 1, a, 1, x, w, 1, iw, 0, info )
632 CALL chkxer( 'SSYEVD', infot, nout, lerr, ok )
633 infot = 10
634 CALL ssyevd( 'N', 'U', 2, a, 2, x, w, 5, iw, 0, info )
635 CALL chkxer( 'SSYEVD', infot, nout, lerr, ok )
636 infot = 10
637 CALL ssyevd( 'V', 'U', 2, a, 2, x, w, 27, iw, 11, info )
638 CALL chkxer( 'SSYEVD', infot, nout, lerr, ok )
639 nt = nt + 10
640*
641* SSYEVD_2STAGE
642*
643 srnamt = 'SSYEVD_2STAGE'
644 infot = 1
645 CALL ssyevd_2stage( '/', 'U', 0, a, 1, x, w, 1, iw, 1, info )
646 CALL chkxer( 'SSYEVD_2STAGE', infot, nout, lerr, ok )
647 infot = 1
648 CALL ssyevd_2stage( 'v', 'u', 0, A, 1, X, W, 1, IW, 1, INFO )
649 CALL CHKXER( 'ssyevd_2stage', INFOT, NOUT, LERR, OK )
650 INFOT = 2
651 CALL SSYEVD_2STAGE( 'n', '/', 0, A, 1, X, W, 1, IW, 1, INFO )
652 CALL CHKXER( 'ssyevd_2stage', INFOT, NOUT, LERR, OK )
653 INFOT = 3
654 CALL SSYEVD_2STAGE( 'n', 'u', -1, A, 1, X, W, 1, IW, 1, INFO )
655 CALL CHKXER( 'ssyevd_2stage', INFOT, NOUT, LERR, OK )
656 INFOT = 5
657 CALL SSYEVD_2STAGE( 'n', 'u', 2, A, 1, X, W, 3, IW, 1, INFO )
658 CALL CHKXER( 'ssyevd_2stage', INFOT, NOUT, LERR, OK )
659 INFOT = 8
660 CALL SSYEVD_2STAGE( 'n', 'u', 1, A, 1, X, W, 0, IW, 1, INFO )
661 CALL CHKXER( 'ssyevd_2stage', INFOT, NOUT, LERR, OK )
662 INFOT = 8
663 CALL SSYEVD_2STAGE( 'n', 'u', 2, A, 2, X, W, 4, IW, 1, INFO )
664 CALL CHKXER( 'ssyevd_2stage', INFOT, NOUT, LERR, OK )
665* INFOT = 8
666* CALL SSYEVD_2STAGE( 'V', 'U', 2, A, 2, X, W, 20, IW, 12, INFO )
667* CALL CHKXER( 'SSYEVD_2STAGE', INFOT, NOUT, LERR, OK )
668 INFOT = 10
669 CALL SSYEVD_2STAGE( 'n', 'u', 1, A, 1, X, W, 1, IW, 0, INFO )
670 CALL CHKXER( 'ssyevd_2stage', INFOT, NOUT, LERR, OK )
671 INFOT = 10
672 CALL SSYEVD_2STAGE( 'n', 'u', 2, A, 2, X, W, 25, IW, 0, INFO )
673 CALL CHKXER( 'ssyevd_2stage', INFOT, NOUT, LERR, OK )
674* INFOT = 10
675* CALL SSYEVD_2STAGE( 'V', 'U', 2, A, 2, X, W, 27, IW, 11, INFO )
676* CALL CHKXER( 'SSYEVD_2STAGE', INFOT, NOUT, LERR, OK )
677 NT = NT + 9
678*
679* SSYEVR
680*
681 SRNAMT = 'ssyevr'
682 N = 1
683 INFOT = 1
684 CALL SSYEVR( '/', 'a', 'u', 0, A, 1, 0.0, 0.0, 1, 1, 0.0, M, R,
685 $ Z, 1, IW, Q, 26*N, IW( 2*N+1 ), 10*N, INFO )
686 CALL CHKXER( 'ssyevr', INFOT, NOUT, LERR, OK )
687 INFOT = 2
688 CALL SSYEVR( 'v', '/', 'u', 0, A, 1, 0.0, 0.0, 1, 1, 0.0, M, R,
689 $ Z, 1, IW, Q, 26*N, IW( 2*N+1 ), 10*N, INFO )
690 CALL CHKXER( 'ssyevr', INFOT, NOUT, LERR, OK )
691 INFOT = 3
692 CALL SSYEVR( 'v', 'a', '/', -1, A, 1, 0.0, 0.0, 1, 1, 0.0, M,
693 $ R, Z, 1, IW, Q, 26*N, IW( 2*N+1 ), 10*N, INFO )
694 CALL CHKXER( 'ssyevr', INFOT, NOUT, LERR, OK )
695 INFOT = 4
696 CALL SSYEVR( 'v', 'a', 'u', -1, A, 1, 0.0, 0.0, 1, 1, 0.0, M,
697 $ R, Z, 1, IW, Q, 26*N, IW( 2*N+1 ), 10*N, INFO )
698 CALL CHKXER( 'ssyevr', INFOT, NOUT, LERR, OK )
699 INFOT = 6
700 CALL SSYEVR( 'v', 'a', 'u', 2, A, 1, 0.0, 0.0, 1, 1, 0.0, M, R,
701 $ Z, 1, IW, Q, 26*N, IW( 2*N+1 ), 10*N, INFO )
702 CALL CHKXER( 'ssyevr', INFOT, NOUT, LERR, OK )
703 INFOT = 8
704 CALL SSYEVR( 'v', 'v', 'u', 1, A, 1, 0.0E0, 0.0E0, 1, 1, 0.0,
705 $ M, R, Z, 1, IW, Q, 26*N, IW( 2*N+1 ), 10*N, INFO )
706 CALL CHKXER( 'ssyevr', INFOT, NOUT, LERR, OK )
707 INFOT = 9
708 CALL SSYEVR( 'v', 'i', 'u', 1, A, 1, 0.0E0, 0.0E0, 0, 1, 0.0,
709 $ M, R, Z, 1, IW, Q, 26*N, IW( 2*N+1 ), 10*N, INFO )
710 CALL CHKXER( 'ssyevr', INFOT, NOUT, LERR, OK )
711 INFOT = 10
712*
713 CALL SSYEVR( 'v', 'i', 'u', 2, A, 2, 0.0E0, 0.0E0, 2, 1, 0.0,
714 $ M, R, Z, 1, IW, Q, 26*N, IW( 2*N+1 ), 10*N, INFO )
715 CALL CHKXER( 'ssyevr', INFOT, NOUT, LERR, OK )
716 INFOT = 15
717 CALL SSYEVR( 'v', 'i', 'u', 1, A, 1, 0.0E0, 0.0E0, 1, 1, 0.0,
718 $ M, R, Z, 0, IW, Q, 26*N, IW( 2*N+1 ), 10*N, INFO )
719 CALL CHKXER( 'ssyevr', INFOT, NOUT, LERR, OK )
720 INFOT = 18
721 CALL SSYEVR( 'v', 'i', 'u', 1, A, 1, 0.0E0, 0.0E0, 1, 1, 0.0,
722 $ M, R, Z, 1, IW, Q, 26*N-1, IW( 2*N+1 ), 10*N,
723 $ INFO )
724 CALL CHKXER( 'ssyevr', INFOT, NOUT, LERR, OK )
725 INFOT = 20
726 CALL SSYEVR( 'v', 'i', 'u', 1, A, 1, 0.0E0, 0.0E0, 1, 1, 0.0,
727 $ M, R, Z, 1, IW, Q, 26*N, IW( 2*N+1 ), 10*N-1,
728 $ INFO )
729 CALL CHKXER( 'ssyevr', INFOT, NOUT, LERR, OK )
730 NT = NT + 11
731*
732* SSYEVR_2STAGE
733*
734 SRNAMT = 'ssyevr_2stage'
735 N = 1
736 INFOT = 1
737 CALL SSYEVR_2STAGE( '/', 'a', 'u', 0, A, 1,
738 $ 0.0E0, 0.0E0, 1, 1, 0.0E0,
739 $ M, R, Z, 1, IW, Q, 26*N, IW( 2*N+1 ), 10*N, INFO )
740 CALL CHKXER( 'ssyevr_2stage', INFOT, NOUT, LERR, OK )
741 INFOT = 1
742 CALL SSYEVR_2STAGE( 'v', 'a', 'u', 0, A, 1,
743 $ 0.0E0, 0.0E0, 1, 1, 0.0E0,
744 $ M, R, Z, 1, IW, Q, 26*N, IW( 2*N+1 ), 10*N, INFO )
745 CALL CHKXER( 'ssyevr_2stage', INFOT, NOUT, LERR, OK )
746 INFOT = 2
747 CALL SSYEVR_2STAGE( 'n', '/', 'u', 0, A, 1,
748 $ 0.0E0, 0.0E0, 1, 1, 0.0E0,
749 $ M, R, Z, 1, IW, Q, 26*N, IW( 2*N+1 ), 10*N, INFO )
750 CALL CHKXER( 'ssyevr_2stage', INFOT, NOUT, LERR, OK )
751 INFOT = 3
752 CALL SSYEVR_2STAGE( 'n', 'a', '/', -1, A, 1,
753 $ 0.0E0, 0.0E0, 1, 1, 0.0E0,
754 $ M, R, Z, 1, IW, Q, 26*N, IW( 2*N+1 ), 10*N, INFO )
755 CALL CHKXER( 'ssyevr_2stage', INFOT, NOUT, LERR, OK )
756 INFOT = 4
757 CALL SSYEVR_2STAGE( 'n', 'a', 'u', -1, A, 1,
758 $ 0.0E0, 0.0E0, 1, 1, 0.0E0,
759 $ M, R, Z, 1, IW, Q, 26*N, IW( 2*N+1 ), 10*N, INFO )
760 CALL CHKXER( 'ssyevr_2stage', INFOT, NOUT, LERR, OK )
761 INFOT = 6
762 CALL SSYEVR_2STAGE( 'n', 'a', 'u', 2, A, 1,
763 $ 0.0E0, 0.0E0, 1, 1, 0.0E0,
764 $ M, R, Z, 1, IW, Q, 26*N, IW( 2*N+1 ), 10*N, INFO )
765 CALL CHKXER( 'ssyevr_2stage', INFOT, NOUT, LERR, OK )
766 INFOT = 8
767 CALL SSYEVR_2STAGE( 'n', 'v', 'u', 1, A, 1,
768 $ 0.0E0, 0.0E0, 1, 1, 0.0E0,
769 $ M, R, Z, 1, IW, Q, 26*N, IW( 2*N+1 ), 10*N, INFO )
770 CALL CHKXER( 'ssyevr_2stage', INFOT, NOUT, LERR, OK )
771 INFOT = 9
772 CALL SSYEVR_2STAGE( 'n', 'i', 'u', 1, A, 1,
773 $ 0.0E0, 0.0E0, 0, 1, 0.0E0,
774 $ M, R, Z, 1, IW, Q, 26*N, IW( 2*N+1 ), 10*N, INFO )
775 CALL CHKXER( 'ssyevr_2stage', INFOT, NOUT, LERR, OK )
776 INFOT = 10
777 CALL SSYEVR_2STAGE( 'n', 'i', 'u', 2, A, 2,
778 $ 0.0E0, 0.0E0, 2, 1, 0.0E0,
779 $ M, R, Z, 1, IW, Q, 26*N, IW( 2*N+1 ), 10*N, INFO )
780 CALL CHKXER( 'ssyevr_2stage', INFOT, NOUT, LERR, OK )
781 INFOT = 15
782 CALL SSYEVR_2STAGE( 'n', 'i', 'u', 1, A, 1,
783 $ 0.0E0, 0.0E0, 1, 1, 0.0E0,
784 $ M, R, Z, 0, IW, Q, 26*N, IW( 2*N+1 ), 10*N, INFO )
785 CALL CHKXER( 'ssyevr_2stage', INFOT, NOUT, LERR, OK )
786 INFOT = 18
787 CALL SSYEVR_2STAGE( 'n', 'i', 'u', 1, A, 1,
788 $ 0.0E0, 0.0E0, 1, 1, 0.0E0,
789 $ M, R, Z, 1, IW, Q, 0, IW( 2*N+1 ), 10*N,
790 $ INFO )
791 CALL CHKXER( 'ssyevr_2stage', INFOT, NOUT, LERR, OK )
792 INFOT = 20
793 CALL SSYEVR_2STAGE( 'n', 'i', 'u', 1, A, 1,
794 $ 0.0E0, 0.0E0, 1, 1, 0.0E0,
795 $ M, R, Z, 1, IW, Q, 26*N, IW( 2*N+1 ), 0,
796 $ INFO )
797 CALL CHKXER( 'ssyevr_2stage', INFOT, NOUT, LERR, OK )
798 NT = NT + 12
799*
800* SSYEV
801*
802 SRNAMT = 'ssyev '
803 INFOT = 1
804 CALL SSYEV( '/', 'u', 0, A, 1, X, W, 1, INFO )
805 CALL CHKXER( 'ssyev ', INFOT, NOUT, LERR, OK )
806 INFOT = 2
807 CALL SSYEV( 'n', '/', 0, A, 1, X, W, 1, INFO )
808 CALL CHKXER( 'ssyev ', INFOT, NOUT, LERR, OK )
809 INFOT = 3
810 CALL SSYEV( 'n', 'u', -1, A, 1, X, W, 1, INFO )
811 CALL CHKXER( 'ssyev ', INFOT, NOUT, LERR, OK )
812 INFOT = 5
813 CALL SSYEV( 'n', 'u', 2, A, 1, X, W, 3, INFO )
814 CALL CHKXER( 'ssyev ', INFOT, NOUT, LERR, OK )
815 INFOT = 8
816 CALL SSYEV( 'n', 'u', 1, A, 1, X, W, 1, INFO )
817 CALL CHKXER( 'ssyev ', INFOT, NOUT, LERR, OK )
818 NT = NT + 5
819*
820* SSYEV_2STAGE
821*
822 SRNAMT = 'ssyev_2stage '
823 INFOT = 1
824 CALL SSYEV_2STAGE( '/', 'u', 0, A, 1, X, W, 1, INFO )
825 CALL CHKXER( 'ssyev_2stage ', INFOT, NOUT, LERR, OK )
826 INFOT = 1
827 CALL SSYEV_2STAGE( 'v', 'u', 0, A, 1, X, W, 1, INFO )
828 CALL CHKXER( 'ssyev_2stage ', INFOT, NOUT, LERR, OK )
829 INFOT = 2
830 CALL SSYEV_2STAGE( 'n', '/', 0, A, 1, X, W, 1, INFO )
831 CALL CHKXER( 'ssyev_2stage ', INFOT, NOUT, LERR, OK )
832 INFOT = 3
833 CALL SSYEV_2STAGE( 'n', 'u', -1, A, 1, X, W, 1, INFO )
834 CALL CHKXER( 'ssyev_2stage ', INFOT, NOUT, LERR, OK )
835 INFOT = 5
836 CALL SSYEV_2STAGE( 'n', 'u', 2, A, 1, X, W, 3, INFO )
837 CALL CHKXER( 'ssyev_2stage ', INFOT, NOUT, LERR, OK )
838 INFOT = 8
839 CALL SSYEV_2STAGE( 'n', 'u', 1, A, 1, X, W, 1, INFO )
840 CALL CHKXER( 'ssyev_2stage ', INFOT, NOUT, LERR, OK )
841 NT = NT + 6
842*
843* SSYEVX
844*
845 SRNAMT = 'ssyevx'
846 INFOT = 1
847 CALL SSYEVX( '/', 'a', 'u', 0, A, 1, 0.0, 0.0, 0, 0, 0.0, M, X,
848 $ Z, 1, W, 1, IW, I3, INFO )
849 CALL CHKXER( 'ssyevx', INFOT, NOUT, LERR, OK )
850 INFOT = 2
851 CALL SSYEVX( 'n', '/', 'u', 0, A, 1, 0.0, 1.0, 1, 0, 0.0, M, X,
852 $ Z, 1, W, 1, IW, I3, INFO )
853 CALL CHKXER( 'ssyevx', INFOT, NOUT, LERR, OK )
854 INFOT = 3
855 CALL SSYEVX( 'n', 'a', '/', 0, A, 1, 0.0, 0.0, 0, 0, 0.0, M, X,
856 $ Z, 1, W, 1, IW, I3, INFO )
857 INFOT = 4
858 CALL SSYEVX( 'n', 'a', 'u', -1, A, 1, 0.0, 0.0, 0, 0, 0.0, M,
859 $ X, Z, 1, W, 1, IW, I3, INFO )
860 CALL CHKXER( 'ssyevx', INFOT, NOUT, LERR, OK )
861 INFOT = 6
862 CALL SSYEVX( 'n', 'a', 'u', 2, A, 1, 0.0, 0.0, 0, 0, 0.0, M, X,
863 $ Z, 1, W, 16, IW, I3, INFO )
864 CALL CHKXER( 'ssyevx', INFOT, NOUT, LERR, OK )
865 INFOT = 8
866 CALL SSYEVX( 'n', 'v', 'u', 1, A, 1, 0.0, 0.0, 0, 0, 0.0, M, X,
867 $ Z, 1, W, 8, IW, I3, INFO )
868 CALL CHKXER( 'ssyevx', INFOT, NOUT, LERR, OK )
869 INFOT = 9
870 CALL SSYEVX( 'n', 'i', 'u', 1, A, 1, 0.0, 0.0, 0, 0, 0.0, M, X,
871 $ Z, 1, W, 8, IW, I3, INFO )
872 CALL CHKXER( 'ssyevx', INFOT, NOUT, LERR, OK )
873 INFOT = 9
874 CALL SSYEVX( 'n', 'i', 'u', 1, A, 1, 0.0, 0.0, 2, 1, 0.0, M, X,
875 $ Z, 1, W, 8, IW, I3, INFO )
876 CALL CHKXER( 'ssyevx', INFOT, NOUT, LERR, OK )
877 INFOT = 10
878 CALL SSYEVX( 'n', 'i', 'u', 2, A, 2, 0.0, 0.0, 2, 1, 0.0, M, X,
879 $ Z, 1, W, 16, IW, I3, INFO )
880 CALL CHKXER( 'ssyevx', INFOT, NOUT, LERR, OK )
881 INFOT = 10
882 CALL SSYEVX( 'n', 'i', 'u', 1, A, 1, 0.0, 0.0, 1, 2, 0.0, M, X,
883 $ Z, 1, W, 8, IW, I3, INFO )
884 CALL CHKXER( 'ssyevx', INFOT, NOUT, LERR, OK )
885 INFOT = 15
886 CALL SSYEVX( 'v', 'a', 'u', 2, A, 2, 0.0, 0.0, 0, 0, 0.0, M, X,
887 $ Z, 1, W, 16, IW, I3, INFO )
888 CALL CHKXER( 'ssyevx', INFOT, NOUT, LERR, OK )
889 INFOT = 17
890 CALL SSYEVX( 'v', 'a', 'u', 1, A, 1, 0.0, 0.0, 0, 0, 0.0, M, X,
891 $ Z, 1, W, 0, IW, I3, INFO )
892 CALL CHKXER( 'ssyevx', INFOT, NOUT, LERR, OK )
893 NT = NT + 12
894*
895* SSYEVX_2STAGE
896*
897 SRNAMT = 'ssyevx_2stage'
898 INFOT = 1
899 CALL SSYEVX_2STAGE( '/', 'a', 'u', 0, A, 1,
900 $ 0.0E0, 0.0E0, 0, 0, 0.0E0,
901 $ M, X, Z, 1, W, 1, IW, I3, INFO )
902 CALL CHKXER( 'ssyevx_2stage', INFOT, NOUT, LERR, OK )
903 INFOT = 1
904 CALL SSYEVX_2STAGE( 'v', 'a', 'u', 0, A, 1,
905 $ 0.0E0, 0.0E0, 0, 0, 0.0E0,
906 $ M, X, Z, 1, W, 1, IW, I3, INFO )
907 CALL CHKXER( 'ssyevx_2stage', INFOT, NOUT, LERR, OK )
908 INFOT = 2
909 CALL SSYEVX_2STAGE( 'n', '/', 'u', 0, A, 1,
910 $ 0.0E0, 1.0E0, 1, 0, 0.0E0,
911 $ M, X, Z, 1, W, 1, IW, I3, INFO )
912 CALL CHKXER( 'ssyevx_2stage', INFOT, NOUT, LERR, OK )
913 INFOT = 3
914 CALL SSYEVX_2STAGE( 'n', 'a', '/', 0, a, 1,
915 $ 0.0e0, 0.0e0, 0, 0, 0.0e0,
916 $ m, x, z, 1, w, 1, iw, i3, info )
917 infot = 4
918 CALL ssyevx_2stage( 'N', 'A', 'U', -1, a, 1,
919 $ 0.0e0, 0.0e0, 0, 0, 0.0e0,
920 $ m, x, z, 1, w, 1, iw, i3, info )
921 CALL chkxer( 'SSYEVX_2STAGE', infot, nout, lerr, ok )
922 infot = 6
923 CALL ssyevx_2stage( 'N', 'A', 'U', 2, a, 1,
924 $ 0.0e0, 0.0e0, 0, 0, 0.0e0,
925 $ m, x, z, 1, w, 16, iw, i3, info )
926 CALL chkxer( 'SSYEVX_2STAGE', infot, nout, lerr, ok )
927 infot = 8
928 CALL ssyevx_2stage( 'N', 'V', 'U', 1, a, 1,
929 $ 0.0e0, 0.0e0, 0, 0, 0.0e0,
930 $ m, x, z, 1, w, 8, iw, i3, info )
931 CALL chkxer( 'SSYEVX_2STAGE', infot, nout, lerr, ok )
932 infot = 9
933 CALL ssyevx_2stage( 'N', 'I', 'U', 1, a, 1,
934 $ 0.0e0, 0.0e0, 0, 0, 0.0e0,
935 $ m, x, z, 1, w, 8, iw, i3, info )
936 CALL chkxer( 'SSYEVX_2STAGE', infot, nout, lerr, ok )
937 infot = 9
938 CALL ssyevx_2stage( 'N', 'I', 'U', 1, a, 1,
939 $ 0.0e0, 0.0e0, 2, 1, 0.0e0,
940 $ m, x, z, 1, w, 8, iw, i3, info )
941 CALL chkxer( 'SSYEVX_2STAGE', infot, nout, lerr, ok )
942 infot = 10
943 CALL ssyevx_2stage( 'n', 'i', 'u', 2, A, 2,
944 $ 0.0E0, 0.0E0, 2, 1, 0.0E0,
945 $ M, X, Z, 1, W, 16, IW, I3, INFO )
946 CALL CHKXER( 'ssyevx_2stage', INFOT, NOUT, LERR, OK )
947 INFOT = 10
948 CALL SSYEVX_2STAGE( 'n', 'i', 'u', 1, A, 1,
949 $ 0.0E0, 0.0E0, 1, 2, 0.0E0,
950 $ M, X, Z, 1, W, 8, IW, I3, INFO )
951 CALL CHKXER( 'ssyevx_2stage', INFOT, NOUT, LERR, OK )
952 INFOT = 15
953 CALL SSYEVX_2STAGE( 'n', 'a', 'u', 2, A, 2,
954 $ 0.0E0, 0.0E0, 0, 0, 0.0E0,
955 $ M, X, Z, 0, W, 16, IW, I3, INFO )
956 CALL CHKXER( 'ssyevx_2stage', INFOT, NOUT, LERR, OK )
957 INFOT = 17
958 CALL SSYEVX_2STAGE( 'n', 'a', 'u', 1, A, 1,
959 $ 0.0E0, 0.0E0, 0, 0, 0.0E0,
960 $ M, X, Z, 1, W, 0, IW, I3, INFO )
961 CALL CHKXER( 'ssyevx_2stage', INFOT, NOUT, LERR, OK )
962 NT = NT + 13
963*
964* SSPEVD
965*
966 SRNAMT = 'sspevd'
967 INFOT = 1
968 CALL SSPEVD( '/', 'u', 0, A, X, Z, 1, W, 1, IW, 1, INFO )
969 CALL CHKXER( 'sspevd', INFOT, NOUT, LERR, OK )
970 INFOT = 2
971 CALL SSPEVD( 'n', '/', 0, A, X, Z, 1, W, 1, IW, 1, INFO )
972 CALL CHKXER( 'sspevd', INFOT, NOUT, LERR, OK )
973 INFOT = 3
974 CALL SSPEVD( 'n', 'u', -1, A, X, Z, 1, W, 1, IW, 1, INFO )
975 CALL CHKXER( 'sspevd', INFOT, NOUT, LERR, OK )
976 INFOT = 7
977 CALL SSPEVD( 'v', 'u', 2, A, X, Z, 1, W, 23, IW, 12, INFO )
978 CALL CHKXER( 'sspevd', INFOT, NOUT, LERR, OK )
979 INFOT = 9
980 CALL SSPEVD( 'n', 'u', 1, A, X, Z, 1, W, 0, IW, 1, INFO )
981 CALL CHKXER( 'sspevd', INFOT, NOUT, LERR, OK )
982 INFOT = 9
983 CALL SSPEVD( 'n', 'u', 2, A, X, Z, 1, W, 3, IW, 1, INFO )
984 CALL CHKXER( 'sspevd', INFOT, NOUT, LERR, OK )
985 INFOT = 9
986 CALL SSPEVD( 'v', 'u', 2, A, X, Z, 2, W, 16, IW, 12, INFO )
987 CALL CHKXER( 'sspevd', INFOT, NOUT, LERR, OK )
988 INFOT = 11
989 CALL SSPEVD( 'n', 'u', 1, A, X, Z, 1, W, 1, IW, 0, INFO )
990 CALL CHKXER( 'sspevd', INFOT, NOUT, LERR, OK )
991 INFOT = 11
992 CALL SSPEVD( 'n', 'u', 2, A, X, Z, 1, W, 4, IW, 0, INFO )
993 CALL CHKXER( 'sspevd', INFOT, NOUT, LERR, OK )
994 INFOT = 11
995 CALL SSPEVD( 'v', 'u', 2, A, X, Z, 2, W, 23, IW, 11, INFO )
996 CALL CHKXER( 'sspevd', INFOT, NOUT, LERR, OK )
997 NT = NT + 10
998*
999* SSPEV
1000*
1001 SRNAMT = 'sspev '
1002 INFOT = 1
1003 CALL SSPEV( '/', 'u', 0, A, W, Z, 1, X, INFO )
1004 CALL CHKXER( 'sspev ', INFOT, NOUT, LERR, OK )
1005 INFOT = 2
1006 CALL SSPEV( 'n', '/', 0, A, W, Z, 1, X, INFO )
1007 CALL CHKXER( 'sspev ', INFOT, NOUT, LERR, OK )
1008 INFOT = 3
1009 CALL SSPEV( 'n', 'u', -1, A, W, Z, 1, X, INFO )
1010 CALL CHKXER( 'sspev ', INFOT, NOUT, LERR, OK )
1011 INFOT = 7
1012 CALL SSPEV( 'v', 'u', 2, A, W, Z, 1, X, INFO )
1013 CALL CHKXER( 'sspev ', INFOT, NOUT, LERR, OK )
1014 NT = NT + 4
1015*
1016* SSPEVX
1017*
1018 SRNAMT = 'sspevx'
1019 INFOT = 1
1020 CALL SSPEVX( '/', 'a', 'u', 0, A, 0.0, 0.0, 0, 0, 0.0, M, X, Z,
1021 $ 1, W, IW, I3, INFO )
1022 CALL CHKXER( 'sspevx', INFOT, NOUT, LERR, OK )
1023 INFOT = 2
1024 CALL SSPEVX( 'n', '/', 'u', 0, A, 0.0, 0.0, 0, 0, 0.0, M, X, Z,
1025 $ 1, W, IW, I3, INFO )
1026 CALL CHKXER( 'sspevx', INFOT, NOUT, LERR, OK )
1027 INFOT = 3
1028 CALL SSPEVX( 'n', 'a', '/', 0, A, 0.0, 0.0, 0, 0, 0.0, M, X, Z,
1029 $ 1, W, IW, I3, INFO )
1030 INFOT = 4
1031 CALL SSPEVX( 'n', 'a', 'u', -1, A, 0.0, 0.0, 0, 0, 0.0, M, X,
1032 $ Z, 1, W, IW, I3, INFO )
1033 CALL CHKXER( 'sspevx', INFOT, NOUT, LERR, OK )
1034 INFOT = 7
1035 CALL SSPEVX( 'n', 'v', 'u', 1, A, 0.0, 0.0, 0, 0, 0.0, M, X, Z,
1036 $ 1, W, IW, I3, INFO )
1037 CALL CHKXER( 'sspevx', INFOT, NOUT, LERR, OK )
1038 INFOT = 8
1039 CALL SSPEVX( 'n', 'i', 'u', 1, A, 0.0, 0.0, 0, 0, 0.0, M, X, Z,
1040 $ 1, W, IW, I3, INFO )
1041 CALL CHKXER( 'sspevx', INFOT, NOUT, LERR, OK )
1042 INFOT = 8
1043 CALL SSPEVX( 'n', 'i', 'u', 1, A, 0.0, 0.0, 2, 1, 0.0, M, X, Z,
1044 $ 1, W, IW, I3, INFO )
1045 CALL CHKXER( 'sspevx', INFOT, NOUT, LERR, OK )
1046 INFOT = 9
1047 CALL SSPEVX( 'n', 'i', 'u', 2, A, 0.0, 0.0, 2, 1, 0.0, M, X, Z,
1048 $ 1, W, IW, I3, INFO )
1049 CALL CHKXER( 'sspevx', INFOT, NOUT, LERR, OK )
1050 INFOT = 9
1051 CALL SSPEVX( 'n', 'i', 'u', 1, A, 0.0, 0.0, 1, 2, 0.0, M, X, Z,
1052 $ 1, W, IW, I3, INFO )
1053 CALL CHKXER( 'sspevx', INFOT, NOUT, LERR, OK )
1054 INFOT = 14
1055 CALL SSPEVX( 'v', 'a', 'u', 2, A, 0.0, 0.0, 0, 0, 0.0, M, X, Z,
1056 $ 1, W, IW, I3, INFO )
1057 CALL CHKXER( 'sspevx', INFOT, NOUT, LERR, OK )
1058 NT = NT + 10
1059*
1060* Test error exits for the SB path.
1061*
1062 ELSE IF( LSAMEN( 2, C2, 'sb' ) ) THEN
1063*
1064* SSBTRD
1065*
1066 SRNAMT = 'ssbtrd'
1067 INFOT = 1
1068 CALL SSBTRD( '/', 'u', 0, 0, A, 1, D, E, Z, 1, W, INFO )
1069 CALL CHKXER( 'ssbtrd', INFOT, NOUT, LERR, OK )
1070 INFOT = 2
1071 CALL SSBTRD( 'n', '/', 0, 0, A, 1, D, E, Z, 1, W, INFO )
1072 CALL CHKXER( 'ssbtrd', INFOT, NOUT, LERR, OK )
1073 INFOT = 3
1074 CALL SSBTRD( 'n', 'u', -1, 0, A, 1, D, E, Z, 1, W, INFO )
1075 CALL CHKXER( 'ssbtrd', INFOT, NOUT, LERR, OK )
1076 INFOT = 4
1077 CALL SSBTRD( 'n', 'u', 0, -1, A, 1, D, E, Z, 1, W, INFO )
1078 CALL CHKXER( 'ssbtrd', INFOT, NOUT, LERR, OK )
1079 INFOT = 6
1080 CALL SSBTRD( 'n', 'u', 1, 1, A, 1, D, E, Z, 1, W, INFO )
1081 CALL CHKXER( 'ssbtrd', INFOT, NOUT, LERR, OK )
1082 INFOT = 10
1083 CALL SSBTRD( 'v', 'u', 2, 0, A, 1, D, E, Z, 1, W, INFO )
1084 CALL CHKXER( 'ssbtrd', INFOT, NOUT, LERR, OK )
1085 NT = NT + 6
1086*
1087* SSYTRD_SB2ST
1088*
1089 SRNAMT = 'ssytrd_sb2st'
1090 INFOT = 1
1091 CALL SSYTRD_SB2ST( '/', 'n', 'u', 0, 0, A, 1, D, E,
1092 $ C, 1, W, 1, INFO )
1093 CALL CHKXER( 'ssytrd_sb2st', INFOT, NOUT, LERR, OK )
1094 INFOT = 2
1095 CALL SSYTRD_SB2ST( 'n', '/', 'u', 0, 0, A, 1, D, E,
1096 $ C, 1, W, 1, INFO )
1097 CALL CHKXER( 'ssytrd_sb2st', INFOT, NOUT, LERR, OK )
1098 INFOT = 2
1099 CALL SSYTRD_SB2ST( 'n', 'h', 'u', 0, 0, A, 1, D, E,
1100 $ C, 1, W, 1, INFO )
1101 CALL CHKXER( 'ssytrd_sb2st', INFOT, NOUT, LERR, OK )
1102 INFOT = 3
1103 CALL SSYTRD_SB2ST( 'n', 'n', '/', 0, 0, A, 1, D, E,
1104 $ C, 1, W, 1, INFO )
1105 CALL CHKXER( 'ssytrd_sb2st', INFOT, NOUT, LERR, OK )
1106 INFOT = 4
1107 CALL SSYTRD_SB2ST( 'n', 'n', 'u', -1, 0, A, 1, D, E,
1108 $ C, 1, W, 1, INFO )
1109 CALL CHKXER( 'ssytrd_sb2st', INFOT, NOUT, LERR, OK )
1110 INFOT = 5
1111 CALL SSYTRD_SB2ST( 'n', 'n', 'u', 0, -1, a, 1, d, e,
1112 $ c, 1, w, 1, info )
1113 CALL chkxer( 'SSYTRD_SB2ST', infot, nout, lerr, ok )
1114 infot = 7
1115 CALL ssytrd_sb2st( 'N', 'N', 'U', 0, 1, a, 1, d, e,
1116 $ c, 1, w, 1, info )
1117 CALL chkxer( 'SSYTRD_SB2ST', infot, nout, lerr, ok )
1118 infot = 11
1119 CALL ssytrd_sb2st( 'N', 'N', 'U', 0, 0, a, 1, d, e,
1120 $ c, 0, w, 1, info )
1121 CALL chkxer( 'SSYTRD_SB2ST', infot, nout, lerr, ok )
1122 infot = 13
1123 CALL ssytrd_sb2st( 'N', 'N', 'U', 0, 0, a, 1, d, e,
1124 $ c, 1, w, 0, info )
1125 CALL chkxer( 'SSYTRD_SB2ST', infot, nout, lerr, ok )
1126 nt = nt + 9
1127*
1128* SSBEVD
1129*
1130 srnamt = 'SSBEVD'
1131 infot = 1
1132 CALL ssbevd( '/', 'U', 0, 0, a, 1, x, z, 1, w, 1, iw, 1, info )
1133 CALL chkxer( 'SSBEVD', infot, nout, lerr, ok )
1134 infot = 2
1135 CALL ssbevd( 'N', '/', 0, 0, a, 1, x, z, 1, w, 1, iw, 1, info )
1136 CALL chkxer( 'SSBEVD', infot, nout, lerr, ok )
1137 infot = 3
1138 CALL ssbevd( 'N', 'U', -1, 0, a, 1, x, z, 1, w, 1, iw, 1,
1139 $ info )
1140 CALL chkxer( 'SSBEVD', infot, nout, lerr, ok )
1141 infot = 4
1142 CALL ssbevd( 'N', 'U', 0, -1, a, 1, x, z, 1, w, 1, iw, 1,
1143 $ info )
1144 CALL chkxer( 'SSBEVD', infot, nout, lerr, ok )
1145 infot = 6
1146 CALL ssbevd( 'N', 'U', 2, 1, a, 1, x, z, 1, w, 4, iw, 1, info )
1147 CALL chkxer( 'SSBEVD', infot, nout, lerr, ok )
1148 infot = 9
1149 CALL ssbevd( 'V', 'U', 2, 1, a, 2, x, z, 1, w, 25, iw, 12,
1150 $ info )
1151 CALL chkxer( 'SSBEVD', infot, nout, lerr, ok )
1152 infot = 11
1153 CALL ssbevd( 'N', 'U', 1, 0, a, 1, x, z, 1, w, 0, iw, 1, info )
1154 CALL chkxer( 'SSBEVD', infot, nout, lerr, ok )
1155 infot = 11
1156 CALL ssbevd( 'N', 'U', 2, 0, a, 1, x, z, 1, w, 3, iw, 1, info )
1157 CALL chkxer( 'SSBEVD', infot, nout, lerr, ok )
1158 infot = 11
1159 CALL ssbevd( 'V', 'U', 2, 0, a, 1, x, z, 2, w, 18, iw, 12,
1160 $ info )
1161 CALL chkxer( 'SSBEVD', infot, nout, lerr, ok )
1162 infot = 13
1163 CALL ssbevd( 'N', 'U', 1, 0, a, 1, x, z, 1, w, 1, iw, 0, info )
1164 CALL chkxer( 'SSBEVD', infot, nout, lerr, ok )
1165 infot = 13
1166 CALL ssbevd( 'V', 'U', 2, 0, a, 1, x, z, 2, w, 25, iw, 11,
1167 $ info )
1168 CALL chkxer( 'SSBEVD', infot, nout, lerr, ok )
1169 nt = nt + 11
1170*
1171* SSBEVD_2STAGE
1172*
1173 srnamt = 'SSBEVD_2STAGE'
1174 infot = 1
1175 CALL ssbevd_2stage( '/', 'U', 0, 0, a, 1, x, z, 1, w,
1176 $ 1, iw, 1, info )
1177 CALL chkxer( 'SSBEVD_2STAGE', infot, nout, lerr, ok )
1178 infot = 1
1179 CALL ssbevd_2stage( 'V', 'U', 0, 0, a, 1, x, z, 1, w,
1180 $ 1, iw, 1, info )
1181 CALL chkxer( 'SSBEVD_2STAGE', infot, nout, lerr, ok )
1182 infot = 2
1183 CALL ssbevd_2stage( 'N', '/', 0, 0, a, 1, x, z, 1, w,
1184 $ 1, iw, 1, info )
1185 CALL chkxer( 'SSBEVD_2STAGE', infot, nout, lerr, ok )
1186 infot = 3
1187 CALL ssbevd_2stage( 'N', 'U', -1, 0, a, 1, x, z, 1, w,
1188 $ 1, iw, 1, info )
1189 CALL chkxer( 'SSBEVD_2STAGE', infot, nout, lerr, ok )
1190 infot = 4
1191 CALL ssbevd_2stage( 'N', 'U', 0, -1, a, 1, x, z, 1, w,
1192 $ 1, iw, 1, info )
1193 CALL chkxer( 'SSBEVD_2STAGE', infot, nout, lerr, ok )
1194 infot = 6
1195 CALL ssbevd_2stage( 'N', 'U', 2, 1, a, 1, x, z, 1, w,
1196 $ 4, iw, 1, info )
1197 CALL chkxer( 'SSBEVD_2STAGE', infot, nout, lerr, ok )
1198* INFOT = 9
1199* CALL SSBEVD_2STAGE( 'V', 'U', 2, 1, A, 2, X, Z, 1, W,
1200* $ 25, IW, 12, INFO )
1201* CALL CHKXER( 'SSBEVD_2STAGE', INFOT, NOUT, LERR, OK )
1202 infot = 11
1203 CALL ssbevd_2stage( 'N', 'U', 1, 0, a, 1, x, z, 1, w,
1204 $ 0, iw, 1, info )
1205 CALL chkxer( 'SSBEVD_2STAGE', infot, nout, lerr, ok )
1206 infot = 11
1207 CALL ssbevd_2stage( 'N', 'U', 2, 0, a, 1, x, z, 1, w,
1208 $ 3, iw, 1, info )
1209 CALL chkxer( 'SSBEVD_2STAGE', infot, nout, lerr, ok )
1210* INFOT = 11
1211* CALL SSBEVD_2STAGE( 'V', 'U', 2, 0, A, 1, X, Z, 2, W,
1212* $ 18, IW, 12, INFO )
1213* CALL CHKXER( 'SSBEVD_2STAGE', INFOT, NOUT, LERR, OK )
1214 infot = 13
1215 CALL ssbevd_2stage( 'N', 'U', 1, 0, a, 1, x, z, 1, w,
1216 $ 1, iw, 0, info )
1217 CALL chkxer( 'SSBEVD_2STAGE', infot, nout, lerr, ok )
1218* INFOT = 13
1219* CALL SSBEVD_2STAGE( 'V', 'U', 2, 0, A, 1, X, Z, 2, W,
1220* $ 25, IW, 11, INFO )
1221* CALL CHKXER( 'SSBEVD_2STAGE', INFOT, NOUT, LERR, OK )
1222* NT = NT + 12
1223 nt = nt + 9
1224*
1225* SSBEV
1226*
1227 srnamt = 'SSBEV '
1228 infot = 1
1229 CALL ssbev( '/', 'U', 0, 0, a, 1, x, z, 1, w, info )
1230 CALL chkxer( 'SSBEV ', infot, nout, lerr, ok )
1231 infot = 2
1232 CALL ssbev( 'N', '/', 0, 0, a, 1, x, z, 1, w, info )
1233 CALL chkxer( 'SSBEV ', infot, nout, lerr, ok )
1234 infot = 3
1235 CALL ssbev( 'N', 'U', -1, 0, a, 1, x, z, 1, w, info )
1236 CALL chkxer( 'SSBEV ', infot, nout, lerr, ok )
1237 infot = 4
1238 CALL ssbev( 'N', 'U', 0, -1, a, 1, x, z, 1, w, info )
1239 CALL chkxer( 'SSBEV ', infot, nout, lerr, ok )
1240 infot = 6
1241 CALL ssbev( 'N', 'U', 2, 1, a, 1, x, z, 1, w, info )
1242 CALL chkxer( 'SSBEV ', infot, nout, lerr, ok )
1243 infot = 9
1244 CALL ssbev( 'V', 'U', 2, 0, a, 1, x, z, 1, w, info )
1245 CALL chkxer( 'SSBEV ', infot, nout, lerr, ok )
1246 nt = nt + 6
1247*
1248* SSBEV_2STAGE
1249*
1250 srnamt = 'SSBEV_2STAGE '
1251 infot = 1
1252 CALL ssbev_2stage( '/', 'U', 0, 0, a, 1, x, z, 1, w, 0, info )
1253 CALL chkxer( 'SSBEV_2STAGE ', infot, nout, lerr, ok )
1254 infot = 1
1255 CALL ssbev_2stage( 'V', 'U', 0, 0, a, 1, x, z, 1, w, 0, info )
1256 CALL chkxer( 'SSBEV_2STAGE ', infot, nout, lerr, ok )
1257 infot = 2
1258 CALL ssbev_2stage( 'N', '/', 0, 0, a, 1, x, z, 1, w, 0, info )
1259 CALL chkxer( 'SSBEV_2STAGE ', infot, nout, lerr, ok )
1260 infot = 3
1261 CALL ssbev_2stage( 'N', 'U', -1, 0, a, 1, x, z, 1, w, 0, info )
1262 CALL chkxer( 'SSBEV_2STAGE ', infot, nout, lerr, ok )
1263 infot = 4
1264 CALL ssbev_2stage( 'N', 'U', 0, -1, a, 1, x, z, 1, w, 0, info )
1265 CALL chkxer( 'SSBEV_2STAGE ', infot, nout, lerr, ok )
1266 infot = 6
1267 CALL ssbev_2stage( 'N', 'U', 2, 1, a, 1, x, z, 1, w, 0, info )
1268 CALL chkxer( 'SSBEV_2STAGE ', infot, nout, lerr, ok )
1269 infot = 9
1270 CALL ssbev_2stage( 'N', 'U', 2, 0, a, 1, x, z, 0, w, 0, info )
1271 CALL chkxer( 'SSBEV_2STAGE ', infot, nout, lerr, ok )
1272 infot = 11
1273 CALL ssbev_2stage( 'N', 'U', 0, 0, a, 1, x, z, 1, w, 0, info )
1274 CALL chkxer( 'SSBEV_2STAGE ', infot, nout, lerr, ok )
1275 nt = nt + 8
1276*
1277* SSBEVX
1278*
1279 srnamt = 'SSBEVX'
1280 infot = 1
1281 CALL ssbevx( '/', 'A', 'U', 0, 0, a, 1, q, 1, 0.0, 0.0, 0, 0,
1282 $ 0.0, m, x, z, 1, w, iw, i3, info )
1283 CALL chkxer( 'SSBEVX', infot, nout, lerr, ok )
1284 infot = 2
1285 CALL ssbevx( 'N', '/', 'U', 0, 0, a, 1, q, 1, 0.0, 0.0, 0, 0,
1286 $ 0.0, m, x, z, 1, w, iw, i3, info )
1287 CALL chkxer( 'SSBEVX', infot, nout, lerr, ok )
1288 infot = 3
1289 CALL ssbevx( 'N', 'A', '/', 0, 0, a, 1, q, 1, 0.0, 0.0, 0, 0,
1290 $ 0.0, m, x, z, 1, w, iw, i3, info )
1291 CALL chkxer( 'SSBEVX', infot, nout, lerr, ok )
1292 infot = 4
1293 CALL ssbevx( 'N', 'A', 'U', -1, 0, a, 1, q, 1, 0.0, 0.0, 0, 0,
1294 $ 0.0, m, x, z, 1, w, iw, i3, info )
1295 CALL chkxer( 'SSBEVX', infot, nout, lerr, ok )
1296 infot = 5
1297 CALL ssbevx( 'N', 'A', 'U', 0, -1, a, 1, q, 1, 0.0, 0.0, 0, 0,
1298 $ 0.0, m, x, z, 1, w, iw, i3, info )
1299 CALL chkxer( 'SSBEVX', infot, nout, lerr, ok )
1300 infot = 7
1301 CALL ssbevx( 'N', 'A', 'U', 2, 1, a, 1, q, 1, 0.0, 0.0, 0, 0,
1302 $ 0.0, m, x, z, 1, w, iw, i3, info )
1303 CALL chkxer( 'SSBEVX', infot, nout, lerr, ok )
1304 infot = 9
1305 CALL ssbevx( 'V', 'A', 'U', 2, 0, a, 1, q, 1, 0.0, 0.0, 0, 0,
1306 $ 0.0, m, x, z, 2, w, iw, i3, info )
1307 CALL chkxer( 'SSBEVX', infot, nout, lerr, ok )
1308 infot = 11
1309 CALL ssbevx( 'N', 'V', 'U', 1, 0, a, 1, q, 1, 0.0, 0.0, 0, 0,
1310 $ 0.0, m, x, z, 1, w, iw, i3, info )
1311 CALL chkxer( 'SSBEVX', infot, nout, lerr, ok )
1312 infot = 12
1313 CALL ssbevx( 'N', 'I', 'U', 1, 0, a, 1, q, 1, 0.0, 0.0, 0, 0,
1314 $ 0.0, m, x, z, 1, w, iw, i3, info )
1315 CALL chkxer( 'SSBEVX', infot, nout, lerr, ok )
1316 infot = 12
1317 CALL ssbevx( 'N', 'I', 'U', 1, 0, a, 1, q, 1, 0.0, 0.0, 2, 1,
1318 $ 0.0, m, x, z, 1, w, iw, i3, info )
1319 CALL chkxer( 'SSBEVX', infot, nout, lerr, ok )
1320 infot = 13
1321 CALL ssbevx( 'N', 'I', 'U', 2, 0, a, 1, q, 1, 0.0, 0.0, 2, 1,
1322 $ 0.0, m, x, z, 1, w, iw, i3, info )
1323 CALL chkxer( 'SSBEVX', infot, nout, lerr, ok )
1324 infot = 13
1325 CALL ssbevx( 'N', 'I', 'U', 1, 0, a, 1, q, 1, 0.0, 0.0, 1, 2,
1326 $ 0.0, m, x, z, 1, w, iw, i3, info )
1327 CALL chkxer( 'SSBEVX', infot, nout, lerr, ok )
1328 infot = 18
1329 CALL ssbevx( 'V', 'A', 'U', 2, 0, a, 1, q, 2, 0.0, 0.0, 0, 0,
1330 $ 0.0, m, x, z, 1, w, iw, i3, info )
1331 CALL chkxer( 'SSBEVX', infot, nout, lerr, ok )
1332 nt = nt + 13
1333*
1334* SSBEVX_2STAGE
1335*
1336 srnamt = 'SSBEVX_2STAGE'
1337 infot = 1
1338 CALL ssbevx_2stage( '/', 'A', 'U', 0, 0, a, 1, q, 1, 0.0e0,
1339 $ 0.0e0, 0, 0, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1340 CALL chkxer( 'SSBEVX_2STAGE', infot, nout, lerr, ok )
1341 infot = 1
1342 CALL ssbevx_2stage( 'V', 'A', 'U', 0, 0, a, 1, q, 1, 0.0e0,
1343 $ 0.0e0, 0, 0, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1344 CALL chkxer( 'SSBEVX_2STAGE', infot, nout, lerr, ok )
1345 infot = 2
1346 CALL ssbevx_2stage( 'N', '/', 'U', 0, 0, a, 1, q, 1, 0.0e0,
1347 $ 0.0e0, 0, 0, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1348 CALL chkxer( 'SSBEVX_2STAGE', infot, nout, lerr, ok )
1349 infot = 3
1350 CALL ssbevx_2stage( 'N', 'A', '/', 0, 0, a, 1, q, 1, 0.0e0,
1351 $ 0.0e0, 0, 0, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1352 CALL chkxer( 'SSBEVX_2STAGE', infot, nout, lerr, ok )
1353 infot = 4
1354 CALL ssbevx_2stage( 'N', 'A', 'U', -1, 0, a, 1, q, 1, 0.0e0,
1355 $ 0.0e0, 0, 0, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1356 CALL chkxer( 'SSBEVX_2STAGE', infot, nout, lerr, ok )
1357 infot = 5
1358 CALL ssbevx_2stage( 'N', 'A', 'U', 0, -1, a, 1, q, 1, 0.0e0,
1359 $ 0.0e0, 0, 0, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1360 CALL chkxer( 'SSBEVX_2STAGE', infot, nout, lerr, ok )
1361 infot = 7
1362 CALL ssbevx_2stage( 'N', 'A', 'U', 2, 1, a, 1, q, 1, 0.0e0,
1363 $ 0.0e0, 0, 0, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1364 CALL chkxer( 'SSBEVX_2STAGE', infot, nout, lerr, ok )
1365* INFOT = 9
1366* CALL SSBEVX_2STAGE( 'V', 'A', 'U', 2, 0, A, 1, Q, 1, 0.0E0,
1367* $ 0.0E0, 0, 0, 0.0E0, M, X, Z, 2, W, 0, IW, I3, INFO )
1368* CALL CHKXER( 'SSBEVX_2STAGE', INFOT, NOUT, LERR, OK )
1369 infot = 11
1370 CALL ssbevx_2stage( 'N', 'V', 'U', 1, 0, a, 1, q, 1, 0.0e0,
1371 $ 0.0e0, 0, 0, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1372 CALL chkxer( 'SSBEVX_2STAGE', infot, nout, lerr, ok )
1373 infot = 12
1374 CALL ssbevx_2stage( 'N', 'I', 'U', 1, 0, a, 1, q, 1, 0.0e0,
1375 $ 0.0e0, 0, 0, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1376 CALL chkxer( 'SSBEVX_2STAGE', infot, nout, lerr, ok )
1377 infot = 12
1378 CALL ssbevx_2stage( 'N', 'I', 'U', 1, 0, a, 1, q, 1, 0.0e0,
1379 $ 0.0e0, 2, 1, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1380 CALL chkxer( 'SSBEVX_2STAGE', infot, nout, lerr, ok )
1381 infot = 13
1382 CALL ssbevx_2stage( 'N', 'I', 'U', 2, 0, a, 1, q, 1, 0.0e0,
1383 $ 0.0e0, 2, 1, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1384 CALL chkxer( 'SSBEVX_2STAGE', infot, nout, lerr, ok )
1385 infot = 13
1386 CALL ssbevx_2stage( 'N', 'I', 'U', 1, 0, a, 1, q, 1, 0.0e0,
1387 $ 0.0e0, 1, 2, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1388 CALL chkxer( 'SSBEVX_2STAGE', infot, nout, lerr, ok )
1389* INFOT = 18
1390* CALL SSBEVX_2STAGE( 'V', 'A', 'U', 2, 0, A, 1, Q, 2, 0.0E0,
1391* $ 0.0E0, 0, 0, 0.0E0, M, X, Z, 1, W, 0, IW, I3, INFO )
1392* CALL CHKXER( 'SSBEVX_2STAGE', INFOT, NOUT, LERR, OK )
1393 infot = 20
1394 CALL ssbevx_2stage( 'N', 'A', 'U', 0, 0, a, 1, q, 1, 0.0e0,
1395 $ 0.0e0, 0, 0, 0.0e0, m, x, z, 1, w, 0, iw, i3, info )
1396 CALL chkxer( 'SSBEVX_2STAGE', infot, nout, lerr, ok )
1397* NT = NT + 15
1398 nt = nt + 13
1399 END IF
1400*
1401* Print a summary line.
1402*
1403 IF( ok ) THEN
1404 WRITE( nout, fmt = 9999 )path, nt
1405 ELSE
1406 WRITE( nout, fmt = 9998 )path
1407 END IF
1408*
1409 9999 FORMAT( 1x, a3, ' routines passed the tests of the error exits',
1410 $ ' (', i3, ' tests done)' )
1411 9998 FORMAT( ' *** ', a3, ' routines failed the tests of the error ',
1412 $ 'exits ***' )
1413*
1414 RETURN
1415*
1416* End of SERRST
1417*
1418 END
chkxer
subroutine chkxer(srnamt, infot, nout, lerr, ok)
Definition cblat2.f:3196
ssteqr
subroutine ssteqr(compz, n, d, e, z, ldz, work, info)
SSTEQR
Definition ssteqr.f:131
sstedc
subroutine sstedc(compz, n, d, e, z, ldz, work, lwork, iwork, liwork, info)
SSTEDC
Definition sstedc.f:188
sstebz
subroutine sstebz(range, order, n, vl, vu, il, iu, abstol, d, e, m, nsplit, w, iblock, isplit, work, iwork, info)
SSTEBZ
Definition sstebz.f:273
ssterf
subroutine ssterf(n, d, e, info)
SSTERF
Definition ssterf.f:86
sopgtr
subroutine sopgtr(uplo, n, ap, tau, q, ldq, work, info)
SOPGTR
Definition sopgtr.f:114
sorgtr
subroutine sorgtr(uplo, n, a, lda, tau, work, lwork, info)
SORGTR
Definition sorgtr.f:123
sormtr
subroutine sormtr(side, uplo, trans, m, n, a, lda, tau, c, ldc, work, lwork, info)
SORMTR
Definition sormtr.f:172
ssptrd
subroutine ssptrd(uplo, n, ap, d, e, tau, info)
SSPTRD
Definition ssptrd.f:150
sstein
subroutine sstein(n, d, e, m, w, iblock, isplit, z, ldz, work, iwork, ifail, info)
SSTEIN
Definition sstein.f:174
sopmtr
subroutine sopmtr(side, uplo, trans, m, n, ap, tau, c, ldc, work, info)
SOPMTR
Definition sopmtr.f:150
ssbtrd
subroutine ssbtrd(vect, uplo, n, kd, ab, ldab, d, e, q, ldq, work, info)
SSBTRD
Definition ssbtrd.f:163
sstevx
subroutine sstevx(jobz, range, n, d, e, vl, vu, il, iu, abstol, m, w, z, ldz, work, iwork, ifail, info)
SSTEVX computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
Definition sstevx.f:227
sspevd
subroutine sspevd(jobz, uplo, n, ap, w, z, ldz, work, lwork, iwork, liwork, info)
SSPEVD computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
Definition sspevd.f:178
ssbevx
subroutine ssbevx(jobz, range, uplo, n, kd, ab, ldab, q, ldq, vl, vu, il, iu, abstol, m, w, z, ldz, work, iwork, ifail, info)
SSBEVX computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
Definition ssbevx.f:265
ssbev_2stage
subroutine ssbev_2stage(jobz, uplo, n, kd, ab, ldab, w, z, ldz, work, lwork, info)
SSBEV_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER m...
Definition ssbev_2stage.f:204
ssbevd
subroutine ssbevd(jobz, uplo, n, kd, ab, ldab, w, z, ldz, work, lwork, iwork, liwork, info)
SSBEVD computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
Definition ssbevd.f:193
ssbev
subroutine ssbev(jobz, uplo, n, kd, ab, ldab, w, z, ldz, work, info)
SSBEV computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrices
Definition ssbev.f:146
sspevx
subroutine sspevx(jobz, range, uplo, n, ap, vl, vu, il, iu, abstol, m, w, z, ldz, work, iwork, ifail, info)
SSPEVX computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
Definition sspevx.f:234
ssbevx_2stage
subroutine ssbevx_2stage(jobz, range, uplo, n, kd, ab, ldab, q, ldq, vl, vu, il, iu, abstol, m, w, z, ldz, work, lwork, iwork, ifail, info)
SSBEVX_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER ...
Definition ssbevx_2stage.f:322
sspev
subroutine sspev(jobz, uplo, n, ap, w, z, ldz, work, info)
SSPEV computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrices
Definition sspev.f:130
sstev
subroutine sstev(jobz, n, d, e, z, ldz, work, info)
SSTEV computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrices
Definition sstev.f:116
ssbevd_2stage
subroutine ssbevd_2stage(jobz, uplo, n, kd, ab, ldab, w, z, ldz, work, lwork, iwork, liwork, info)
SSBEVD_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER ...
Definition ssbevd_2stage.f:234
sstevr
subroutine sstevr(jobz, range, n, d, e, vl, vu, il, iu, abstol, m, w, z, ldz, isuppz, work, lwork, iwork, liwork, info)
SSTEVR computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
Definition sstevr.f:306
sstevd
subroutine sstevd(jobz, n, d, e, z, ldz, work, lwork, iwork, liwork, info)
SSTEVD computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
Definition sstevd.f:163
spteqr
subroutine spteqr(compz, n, d, e, z, ldz, work, info)
SPTEQR
Definition spteqr.f:145
ssytrd_sy2sb
subroutine ssytrd_sy2sb(uplo, n, kd, a, lda, ab, ldab, tau, work, lwork, info)
SSYTRD_SY2SB
Definition ssytrd_sy2sb.f:243
ssytrd_2stage
subroutine ssytrd_2stage(vect, uplo, n, a, lda, d, e, tau, hous2, lhous2, work, lwork, info)
SSYTRD_2STAGE
Definition ssytrd_2stage.f:224
ssytrd
subroutine ssytrd(uplo, n, a, lda, d, e, tau, work, lwork, info)
SSYTRD
Definition ssytrd.f:192
ssyevd
subroutine ssyevd(jobz, uplo, n, a, lda, w, work, lwork, iwork, liwork, info)
SSYEVD computes the eigenvalues and, optionally, the left and/or right eigenvectors for SY matrices
Definition ssyevd.f:183
ssyevx
subroutine ssyevx(jobz, range, uplo, n, a, lda, vl, vu, il, iu, abstol, m, w, z, ldz, work, lwork, iwork, ifail, info)
SSYEVX computes the eigenvalues and, optionally, the left and/or right eigenvectors for SY matrices
Definition ssyevx.f:253
ssyevd_2stage
subroutine ssyevd_2stage(jobz, uplo, n, a, lda, w, work, lwork, iwork, liwork, info)
SSYEVD_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for SY mat...
Definition ssyevd_2stage.f:227
ssyevr_2stage
subroutine ssyevr_2stage(jobz, range, uplo, n, a, lda, vl, vu, il, iu, abstol, m, w, z, ldz, isuppz, work, lwork, iwork, liwork, info)
SSYEVR_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for SY mat...
Definition ssyevr_2stage.f:381
ssyevr
subroutine ssyevr(jobz, range, uplo, n, a, lda, vl, vu, il, iu, abstol, m, w, z, ldz, isuppz, work, lwork, iwork, liwork, info)
SSYEVR computes the eigenvalues and, optionally, the left and/or right eigenvectors for SY matrices
Definition ssyevr.f:336
ssyev_2stage
subroutine ssyev_2stage(jobz, uplo, n, a, lda, w, work, lwork, info)
SSYEV_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for SY matr...
Definition ssyev_2stage.f:183
ssyev
subroutine ssyev(jobz, uplo, n, a, lda, w, work, lwork, info)
SSYEV computes the eigenvalues and, optionally, the left and/or right eigenvectors for SY matrices
Definition ssyev.f:132
ssyevx_2stage
subroutine ssyevx_2stage(jobz, range, uplo, n, a, lda, vl, vu, il, iu, abstol, m, w, z, ldz, work, lwork, iwork, ifail, info)
SSYEVX_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for SY mat...
Definition ssyevx_2stage.f:300
serrst
subroutine serrst(path, nunit)
SERRST
Definition serrst.f:61
ssytrd_sb2st
subroutine ssytrd_sb2st(stage1, vect, uplo, n, kd, ab, ldab, d, e, hous, lhous, work, lwork, info)
SSYTRD_SB2ST reduces a real symmetric band matrix A to real symmetric tridiagonal form T
Definition ssytrd_sb2st.F:230