61
62
63
64
65
66
67 CHARACTER*3 PATH
68 INTEGER NUNIT
69
70
71
72
73
74
75 INTEGER NMAX, LIW, LW
76 parameter( nmax = 3, liw = 12*nmax, lw = 20*nmax )
77
78
79 CHARACTER*2 C2
80 INTEGER I, INFO, J, M, N, NSPLIT, NT
81
82
83 INTEGER I1( NMAX ), I2( NMAX ), I3( NMAX ), IW( LIW )
84 DOUBLE PRECISION 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
90 LOGICAL LSAMEN
92
93
103
104
105 LOGICAL LERR, OK
106 CHARACTER*32 SRNAMT
107 INTEGER INFOT, NOUT
108
109
110 COMMON / infoc / infot, nout, ok, lerr
111 COMMON / srnamc / srnamt
112
113
114 INTRINSIC dble
115
116
117
118 nout = nunit
119 WRITE( nout, fmt = * )
120 c2 = path( 2: 3 )
121
122
123
124 DO 20 j = 1, nmax
125 DO 10 i = 1, nmax
126 a( i, j ) = 1.d0 / dble( i+j )
127 10 CONTINUE
128 20 CONTINUE
129 DO 30 j = 1, nmax
130 d( j ) = dble( j )
131 e( j ) = 0.0d0
132 i1( j ) = j
133 i2( j ) = j
134 tau( j ) = 1.d0
135 30 CONTINUE
136 ok = .true.
137 nt = 0
138
139
140
141 IF(
lsamen( 2, c2,
'ST' ) )
THEN
142
143
144
145 srnamt = 'DSYTRD'
146 infot = 1
147 CALL dsytrd(
'/', 0, a, 1, d, e, tau, w, 1, info )
148 CALL chkxer(
'DSYTRD', infot, nout, lerr, ok )
149 infot = 2
150 CALL dsytrd(
'U', -1, a, 1, d, e, tau, w, 1, info )
151 CALL chkxer(
'DSYTRD', infot, nout, lerr, ok )
152 infot = 4
153 CALL dsytrd(
'U', 2, a, 1, d, e, tau, w, 1, info )
154 CALL chkxer(
'DSYTRD', infot, nout, lerr, ok )
155 infot = 9
156 CALL dsytrd(
'U', 0, a, 1, d, e, tau, w, 0, info )
157 CALL chkxer(
'DSYTRD', infot, nout, lerr, ok )
158 nt = nt + 4
159
160
161
162 srnamt = 'DSYTRD_2STAGE'
163 infot = 1
165 $ c, 1, w, 1, info )
166 CALL chkxer(
'DSYTRD_2STAGE', infot, nout, lerr, ok )
167 infot = 1
169 $ c, 1, w, 1, info )
170 CALL chkxer(
'DSYTRD_2STAGE', infot, nout, lerr, ok )
171 infot = 2
173 $ c, 1, w, 1, info )
174 CALL chkxer(
'DSYTRD_2STAGE', infot, nout, lerr, ok )
175 infot = 3
177 $ c, 1, w, 1, info )
178 CALL chkxer(
'DSYTRD_2STAGE', infot, nout, lerr, ok )
179 infot = 5
181 $ c, 1, w, 1, info )
182 CALL chkxer(
'DSYTRD_2STAGE', infot, nout, lerr, ok )
183 infot = 10
185 $ c, 0, w, 1, info )
186 CALL chkxer(
'DSYTRD_2STAGE', infot, nout, lerr, ok )
187 infot = 12
189 $ c, 1, w, 0, info )
190 CALL chkxer(
'DSYTRD_2STAGE', infot, nout, lerr, ok )
191 nt = nt + 7
192
193
194
195 srnamt = 'DSYTRD_SY2SB'
196 infot = 1
197 CALL dsytrd_sy2sb(
'/', 0, 0, a, 1, c, 1, tau, w, 1, info )
198 CALL chkxer(
'DSYTRD_SY2SB', infot, nout, lerr, ok )
199 infot = 2
200 CALL dsytrd_sy2sb(
'U', -1, 0, a, 1, c, 1, tau, w, 1, info )
201 CALL chkxer(
'DSYTRD_SY2SB', infot, nout, lerr, ok )
202 infot = 3
203 CALL dsytrd_sy2sb(
'U', 0, -1, a, 1, c, 1, tau, w, 1, info )
204 CALL chkxer(
'DSYTRD_SY2SB', infot, nout, lerr, ok )
205 infot = 5
206 CALL dsytrd_sy2sb(
'U', 2, 0, a, 1, c, 1, tau, w, 1, info )
207 CALL chkxer(
'DSYTRD_SY2SB', infot, nout, lerr, ok )
208 infot = 7
209 CALL dsytrd_sy2sb(
'U', 0, 2, a, 1, c, 1, tau, w, 1, info )
210 CALL chkxer(
'DSYTRD_SY2SB', infot, nout, lerr, ok )
211 infot = 10
212 CALL dsytrd_sy2sb(
'U', 0, 0, a, 1, c, 1, tau, w, 0, info )
213 CALL chkxer(
'DSYTRD_SY2SB', infot, nout, lerr, ok )
214 nt = nt + 6
215
216
217
218 srnamt = 'DSYTRD_SB2ST'
219 infot = 1
221 $ c, 1, w, 1, info )
222 CALL chkxer(
'DSYTRD_SB2ST', infot, nout, lerr, ok )
223 infot = 2
225 $ c, 1, w, 1, info )
226 CALL chkxer(
'DSYTRD_SB2ST', infot, nout, lerr, ok )
227 infot = 2
229 $ c, 1, w, 1, info )
230 CALL chkxer(
'DSYTRD_SB2ST', infot, nout, lerr, ok )
231 infot = 3
233 $ c, 1, w, 1, info )
234 CALL chkxer(
'DSYTRD_SB2ST', infot, nout, lerr, ok )
235 infot = 4
237 $ c, 1, w, 1, info )
238 CALL chkxer(
'DSYTRD_SB2ST', infot, nout, lerr, ok )
239 infot = 5
241 $ c, 1, w, 1, info )
242 CALL chkxer(
'DSYTRD_SB2ST', infot, nout, lerr, ok )
243 infot = 7
245 $ c, 1, w, 1, info )
246 CALL chkxer(
'DSYTRD_SB2ST', infot, nout, lerr, ok )
247 infot = 11
249 $ c, 0, w, 1, info )
250 CALL chkxer(
'DSYTRD_SB2ST', infot, nout, lerr, ok )
251 infot = 13
253 $ c, 1, w, 0, info )
254 CALL chkxer(
'DSYTRD_SB2ST', infot, nout, lerr, ok )
255 nt = nt + 9
256
257
258
259 srnamt = 'DORGTR'
260 infot = 1
261 CALL dorgtr(
'/', 0, a, 1, tau, w, 1, info )
262 CALL chkxer(
'DORGTR', infot, nout, lerr, ok )
263 infot = 2
264 CALL dorgtr(
'U', -1, a, 1, tau, w, 1, info )
265 CALL chkxer(
'DORGTR', infot, nout, lerr, ok )
266 infot = 4
267 CALL dorgtr(
'U', 2, a, 1, tau, w, 1, info )
268 CALL chkxer(
'DORGTR', infot, nout, lerr, ok )
269 infot = 7
270 CALL dorgtr(
'U', 3, a, 3, tau, w, 1, info )
271 CALL chkxer(
'DORGTR', infot, nout, lerr, ok )
272 nt = nt + 4
273
274
275
276 srnamt = 'DORMTR'
277 infot = 1
278 CALL dormtr(
'/',
'U',
'N', 0, 0, a, 1, tau, c, 1, w, 1, info )
279 CALL chkxer(
'DORMTR', infot, nout, lerr, ok )
280 infot = 2
281 CALL dormtr(
'L',
'/',
'N', 0, 0, a, 1, tau, c, 1, w, 1, info )
282 CALL chkxer(
'DORMTR', infot, nout, lerr, ok )
283 infot = 3
284 CALL dormtr(
'L',
'U',
'/', 0, 0, a, 1, tau, c, 1, w, 1, info )
285 CALL chkxer(
'DORMTR', infot, nout, lerr, ok )
286 infot = 4
287 CALL dormtr(
'L',
'U',
'N', -1, 0, a, 1, tau, c, 1, w, 1,
288 $ info )
289 CALL chkxer(
'DORMTR', infot, nout, lerr, ok )
290 infot = 5
291 CALL dormtr(
'L',
'U',
'N', 0, -1, a, 1, tau, c, 1, w, 1,
292 $ info )
293 CALL chkxer(
'DORMTR', infot, nout, lerr, ok )
294 infot = 7
295 CALL dormtr(
'L',
'U',
'N', 2, 0, a, 1, tau, c, 2, w, 1, info )
296 CALL chkxer(
'DORMTR', infot, nout, lerr, ok )
297 infot = 7
298 CALL dormtr(
'R',
'U',
'N', 0, 2, a, 1, tau, c, 1, w, 1, info )
299 CALL chkxer(
'DORMTR', infot, nout, lerr, ok )
300 infot = 10
301 CALL dormtr(
'L',
'U',
'N', 2, 0, a, 2, tau, c, 1, w, 1, info )
302 CALL chkxer(
'DORMTR', infot, nout, lerr, ok )
303 infot = 12
304 CALL dormtr(
'L',
'U',
'N', 0, 2, a, 1, tau, c, 1, w, 1, info )
305 CALL chkxer(
'DORMTR', infot, nout, lerr, ok )
306 infot = 12
307 CALL dormtr(
'R',
'U',
'N', 2, 0, a, 1, tau, c, 2, w, 1, info )
308 CALL chkxer(
'DORMTR', infot, nout, lerr, ok )
309 nt = nt + 10
310
311
312
313 srnamt = 'DSPTRD'
314 infot = 1
315 CALL dsptrd(
'/', 0, a, d, e, tau, info )
316 CALL chkxer(
'DSPTRD', infot, nout, lerr, ok )
317 infot = 2
318 CALL dsptrd(
'U', -1, a, d, e, tau, info )
319 CALL chkxer(
'DSPTRD', infot, nout, lerr, ok )
320 nt = nt + 2
321
322
323
324 srnamt = 'DOPGTR'
325 infot = 1
326 CALL dopgtr(
'/', 0, a, tau, z, 1, w, info )
327 CALL chkxer(
'DOPGTR', infot, nout, lerr, ok )
328 infot = 2
329 CALL dopgtr(
'U', -1, a, tau, z, 1, w, info )
330 CALL chkxer(
'DOPGTR', infot, nout, lerr, ok )
331 infot = 6
332 CALL dopgtr(
'U', 2, a, tau, z, 1, w, info )
333 CALL chkxer(
'DOPGTR', infot, nout, lerr, ok )
334 nt = nt + 3
335
336
337
338 srnamt = 'DOPMTR'
339 infot = 1
340 CALL dopmtr(
'/',
'U',
'N', 0, 0, a, tau, c, 1, w, info )
341 CALL chkxer(
'DOPMTR', infot, nout, lerr, ok )
342 infot = 2
343 CALL dopmtr(
'L',
'/',
'N', 0, 0, a, tau, c, 1, w, info )
344 CALL chkxer(
'DOPMTR', infot, nout, lerr, ok )
345 infot = 3
346 CALL dopmtr(
'L',
'U',
'/', 0, 0, a, tau, c, 1, w, info )
347 CALL chkxer(
'DOPMTR', infot, nout, lerr, ok )
348 infot = 4
349 CALL dopmtr(
'L',
'U',
'N', -1, 0, a, tau, c, 1, w, info )
350 CALL chkxer(
'DOPMTR', infot, nout, lerr, ok )
351 infot = 5
352 CALL dopmtr(
'L',
'U',
'N', 0, -1, a, tau, c, 1, w, info )
353 CALL chkxer(
'DOPMTR', infot, nout, lerr, ok )
354 infot = 9
355 CALL dopmtr(
'L',
'U',
'N', 2, 0, a, tau, c, 1, w, info )
356 CALL chkxer(
'DOPMTR', infot, nout, lerr, ok )
357 nt = nt + 6
358
359
360
361 srnamt = 'DPTEQR'
362 infot = 1
363 CALL dpteqr(
'/', 0, d, e, z, 1, w, info )
364 CALL chkxer(
'DPTEQR', infot, nout, lerr, ok )
365 infot = 2
366 CALL dpteqr(
'N', -1, d, e, z, 1, w, info )
367 CALL chkxer(
'DPTEQR', infot, nout, lerr, ok )
368 infot = 6
369 CALL dpteqr(
'V', 2, d, e, z, 1, w, info )
370 CALL chkxer(
'DPTEQR', infot, nout, lerr, ok )
371 nt = nt + 3
372
373
374
375 srnamt = 'DSTEBZ'
376 infot = 1
377 CALL dstebz(
'/',
'E', 0, 0.0d0, 1.0d0, 1, 0, 0.0d0, d, e, m,
378 $ nsplit, x, i1, i2, w, iw, info )
379 CALL chkxer(
'DSTEBZ', infot, nout, lerr, ok )
380 infot = 2
381 CALL dstebz(
'A',
'/', 0, 0.0d0, 0.0d0, 0, 0, 0.0d0, d, e, m,
382 $ nsplit, x, i1, i2, w, iw, info )
383 CALL chkxer(
'DSTEBZ', infot, nout, lerr, ok )
384 infot = 3
385 CALL dstebz(
'A',
'E', -1, 0.0d0, 0.0d0, 0, 0, 0.0d0, d, e, m,
386 $ nsplit, x, i1, i2, w, iw, info )
387 CALL chkxer(
'DSTEBZ', infot, nout, lerr, ok )
388 infot = 5
389 CALL dstebz(
'V',
'E', 0, 0.0d0, 0.0d0, 0, 0, 0.0d0, d, e, m,
390 $ nsplit, x, i1, i2, w, iw, info )
391 CALL chkxer(
'DSTEBZ', infot, nout, lerr, ok )
392 infot = 6
393 CALL dstebz(
'I',
'E', 0, 0.0d0, 0.0d0, 0, 0, 0.0d0, d, e, m,
394 $ nsplit, x, i1, i2, w, iw, info )
395 CALL chkxer(
'DSTEBZ', infot, nout, lerr, ok )
396 infot = 6
397 CALL dstebz(
'I',
'E', 1, 0.0d0, 0.0d0, 2, 1, 0.0d0, d, e, m,
398 $ nsplit, x, i1, i2, w, iw, info )
399 CALL chkxer(
'DSTEBZ', infot, nout, lerr, ok )
400 infot = 7
401 CALL dstebz(
'I',
'E', 1, 0.0d0, 0.0d0, 1, 0, 0.0d0, d, e, m,
402 $ nsplit, x, i1, i2, w, iw, info )
403 CALL chkxer(
'DSTEBZ', infot, nout, lerr, ok )
404 infot = 7
405 CALL dstebz(
'I',
'E', 1, 0.0d0, 0.0d0, 1, 2, 0.0d0, d, e, m,
406 $ nsplit, x, i1, i2, w, iw, info )
407 CALL chkxer(
'DSTEBZ', infot, nout, lerr, ok )
408 nt = nt + 8
409
410
411
412 srnamt = 'DSTEIN'
413 infot = 1
414 CALL dstein( -1, d, e, 0, x, i1, i2, z, 1, w, iw, i3, info )
415 CALL chkxer(
'DSTEIN', infot, nout, lerr, ok )
416 infot = 4
417 CALL dstein( 0, d, e, -1, x, i1, i2, z, 1, w, iw, i3, info )
418 CALL chkxer(
'DSTEIN', infot, nout, lerr, ok )
419 infot = 4
420 CALL dstein( 0, d, e, 1, x, i1, i2, z, 1, w, iw, i3, info )
421 CALL chkxer(
'DSTEIN', infot, nout, lerr, ok )
422 infot = 9
423 CALL dstein( 2, d, e, 0, x, i1, i2, z, 1, w, iw, i3, info )
424 CALL chkxer(
'DSTEIN', infot, nout, lerr, ok )
425 nt = nt + 4
426
427
428
429 srnamt = 'DSTEQR'
430 infot = 1
431 CALL dsteqr(
'/', 0, d, e, z, 1, w, info )
432 CALL chkxer(
'DSTEQR', infot, nout, lerr, ok )
433 infot = 2
434 CALL dsteqr(
'N', -1, d, e, z, 1, w, info )
435 CALL chkxer(
'DSTEQR', infot, nout, lerr, ok )
436 infot = 6
437 CALL dsteqr(
'V', 2, d, e, z, 1, w, info )
438 CALL chkxer(
'DSTEQR', infot, nout, lerr, ok )
439 nt = nt + 3
440
441
442
443 srnamt = 'DSTERF'
444 infot = 1
445 CALL dsterf( -1, d, e, info )
446 CALL chkxer(
'DSTERF', infot, nout, lerr, ok )
447 nt = nt + 1
448
449
450
451 srnamt = 'DSTEDC'
452 infot = 1
453 CALL dstedc(
'/', 0, d, e, z, 1, w, 1, iw, 1, info )
454 CALL chkxer(
'DSTEDC', infot, nout, lerr, ok )
455 infot = 2
456 CALL dstedc(
'N', -1, d, e, z, 1, w, 1, iw, 1, info )
457 CALL chkxer(
'DSTEDC', infot, nout, lerr, ok )
458 infot = 6
459 CALL dstedc(
'V', 2, d, e, z, 1, w, 23, iw, 28, info )
460 CALL chkxer(
'DSTEDC', infot, nout, lerr, ok )
461 infot = 8
462 CALL dstedc(
'N', 1, d, e, z, 1, w, 0, iw, 1, info )
463 CALL chkxer(
'DSTEDC', infot, nout, lerr, ok )
464 infot = 8
465 CALL dstedc(
'I', 2, d, e, z, 2, w, 0, iw, 12, info )
466 CALL chkxer(
'DSTEDC', infot, nout, lerr, ok )
467 infot = 8
468 CALL dstedc(
'V', 2, d, e, z, 2, w, 0, iw, 28, info )
469 CALL chkxer(
'DSTEDC', infot, nout, lerr, ok )
470 infot = 10
471 CALL dstedc(
'N', 1, d, e, z, 1, w, 1, iw, 0, info )
472 CALL chkxer(
'DSTEDC', infot, nout, lerr, ok )
473 infot = 10
474 CALL dstedc(
'I', 2, d, e, z, 2, w, 19, iw, 0, info )
475 CALL chkxer(
'DSTEDC', infot, nout, lerr, ok )
476 infot = 10
477 CALL dstedc(
'V', 2, d, e, z, 2, w, 23, iw, 0, info )
478 CALL chkxer(
'DSTEDC', infot, nout, lerr, ok )
479 nt = nt + 9
480
481
482
483 srnamt = 'DSTEVD'
484 infot = 1
485 CALL dstevd(
'/', 0, d, e, z, 1, w, 1, iw, 1, info )
486 CALL chkxer(
'DSTEVD', infot, nout, lerr, ok )
487 infot = 2
488 CALL dstevd(
'N', -1, d, e, z, 1, w, 1, iw, 1, info )
489 CALL chkxer(
'DSTEVD', infot, nout, lerr, ok )
490 infot = 6
491 CALL dstevd(
'V', 2, d, e, z, 1, w, 19, iw, 12, info )
492 CALL chkxer(
'DSTEVD', infot, nout, lerr, ok )
493 infot = 8
494 CALL dstevd(
'N', 1, d, e, z, 1, w, 0, iw, 1, info )
495 CALL chkxer(
'DSTEVD', infot, nout, lerr, ok )
496 infot = 8
497 CALL dstevd(
'V', 2, d, e, z, 2, w, 12, iw, 12, info )
498 CALL chkxer(
'DSTEVD', infot, nout, lerr, ok )
499 infot = 10
500 CALL dstevd(
'N', 0, d, e, z, 1, w, 1, iw, 0, info )
501 CALL chkxer(
'DSTEVD', infot, nout, lerr, ok )
502 infot = 10
503 CALL dstevd(
'V', 2, d, e, z, 2, w, 19, iw, 11, info )
504 CALL chkxer(
'DSTEVD', infot, nout, lerr, ok )
505 nt = nt + 7
506
507
508
509 srnamt = 'DSTEV '
510 infot = 1
511 CALL dstev(
'/', 0, d, e, z, 1, w, info )
512 CALL chkxer(
'DSTEV ', infot, nout, lerr, ok )
513 infot = 2
514 CALL dstev(
'N', -1, d, e, z, 1, w, info )
515 CALL chkxer(
'DSTEV ', infot, nout, lerr, ok )
516 infot = 6
517 CALL dstev(
'V', 2, d, e, z, 1, w, info )
518 CALL chkxer(
'DSTEV ', infot, nout, lerr, ok )
519 nt = nt + 3
520
521
522
523 srnamt = 'DSTEVX'
524 infot = 1
525 CALL dstevx(
'/',
'A', 0, d, e, 0.0d0, 0.0d0, 0, 0, 0.0d0, m,
526 $ x, z, 1, w, iw, i3, info )
527 CALL chkxer(
'DSTEVX', infot, nout, lerr, ok )
528 infot = 2
529 CALL dstevx(
'N',
'/', 0, d, e, 0.0d0, 1.0d0, 1, 0, 0.0d0, m,
530 $ x, z, 1, w, iw, i3, info )
531 CALL chkxer(
'DSTEVX', infot, nout, lerr, ok )
532 infot = 3
533 CALL dstevx(
'N',
'A', -1, d, e, 0.0d0, 0.0d0, 0, 0, 0.0d0, m,
534 $ x, z, 1, w, iw, i3, info )
535 CALL chkxer(
'DSTEVX', infot, nout, lerr, ok )
536 infot = 7
537 CALL dstevx(
'N',
'V', 1, d, e, 0.0d0, 0.0d0, 0, 0, 0.0d0, m,
538 $ x, z, 1, w, iw, i3, info )
539 CALL chkxer(
'DSTEVX', infot, nout, lerr, ok )
540 infot = 8
541 CALL dstevx(
'N',
'I', 1, d, e, 0.0d0, 0.0d0, 0, 0, 0.0d0, m,
542 $ x, z, 1, w, iw, i3, info )
543 CALL chkxer(
'DSTEVX', infot, nout, lerr, ok )
544 infot = 8
545 CALL dstevx(
'N',
'I', 1, d, e, 0.0d0, 0.0d0, 2, 1, 0.0d0, m,
546 $ x, z, 1, w, iw, i3, info )
547 CALL chkxer(
'DSTEVX', infot, nout, lerr, ok )
548 infot = 9
549 CALL dstevx(
'N',
'I', 2, d, e, 0.0d0, 0.0d0, 2, 1, 0.0d0, m,
550 $ x, z, 1, w, iw, i3, info )
551 CALL chkxer(
'DSTEVX', infot, nout, lerr, ok )
552 infot = 9
553 CALL dstevx(
'N',
'I', 1, d, e, 0.0d0, 0.0d0, 1, 2, 0.0d0, m,
554 $ x, z, 1, w, iw, i3, info )
555 CALL chkxer(
'DSTEVX', infot, nout, lerr, ok )
556 infot = 14
557 CALL dstevx(
'V',
'A', 2, d, e, 0.0d0, 0.0d0, 0, 0, 0.0d0, m,
558 $ x, z, 1, w, iw, i3, info )
559 CALL chkxer(
'DSTEVX', infot, nout, lerr, ok )
560 nt = nt + 9
561
562
563
564 n = 1
565 srnamt = 'DSTEVR'
566 infot = 1
567 CALL dstevr(
'/',
'A', 0, d, e, 0.0d0, 0.0d0, 1, 1, 0.0d0, m,
568 $ r, z, 1, iw, x, 20*n, iw( 2*n+1 ), 10*n, info )
569 CALL chkxer(
'DSTEVR', infot, nout, lerr, ok )
570 infot = 2
571 CALL dstevr(
'V',
'/', 0, d, e, 0.0d0, 0.0d0, 1, 1, 0.0d0, m,
572 $ r, z, 1, iw, x, 20*n, iw( 2*n+1 ), 10*n, info )
573 CALL chkxer(
'DSTEVR', infot, nout, lerr, ok )
574 infot = 3
575 CALL dstevr(
'V',
'A', -1, d, e, 0.0d0, 0.0d0, 1, 1, 0.0d0, m,
576 $ r, z, 1, iw, x, 20*n, iw( 2*n+1 ), 10*n, info )
577 CALL chkxer(
'DSTEVR', infot, nout, lerr, ok )
578 infot = 7
579 CALL dstevr(
'V',
'V', 1, d, e, 0.0d0, 0.0d0, 1, 1, 0.0d0, m,
580 $ r, z, 1, iw, x, 20*n, iw( 2*n+1 ), 10*n, info )
581 CALL chkxer(
'DSTEVR', infot, nout, lerr, ok )
582 infot = 8
583 CALL dstevr(
'V',
'I', 1, d, e, 0.0d0, 0.0d0, 0, 1, 0.0d0, m,
584 $ w, z, 1, iw, x, 20*n, iw( 2*n+1 ), 10*n, info )
585 CALL chkxer(
'DSTEVR', infot, nout, lerr, ok )
586 infot = 9
587 n = 2
588 CALL dstevr(
'V',
'I', 2, d, e, 0.0d0, 0.0d0, 2, 1, 0.0d0, m,
589 $ w, z, 1, iw, x, 20*n, iw( 2*n+1 ), 10*n, info )
590 CALL chkxer(
'DSTEVR', infot, nout, lerr, ok )
591 infot = 14
592 n = 1
593 CALL dstevr(
'V',
'I', 1, d, e, 0.0d0, 0.0d0, 1, 1, 0.0d0, m,
594 $ w, z, 0, iw, x, 20*n, iw( 2*n+1 ), 10*n, info )
595 CALL chkxer(
'DSTEVR', infot, nout, lerr, ok )
596 infot = 17
597 CALL dstevr(
'V',
'I', 1, d, e, 0.0d0, 0.0d0, 1, 1, 0.0d0, m,
598 $ w, z, 1, iw, x, 20*n-1, iw( 2*n+1 ), 10*n, info )
599 CALL chkxer(
'DSTEVR', infot, nout, lerr, ok )
600 infot = 19
601 CALL dstevr(
'V',
'I', 1, d, e, 0.0d0, 0.0d0, 1, 1, 0.0d0, m,
602 $ w, z, 1, iw, x, 20*n, iw( 2*n+1 ), 10*n-1, info )
603 CALL chkxer(
'DSTEVR', infot, nout, lerr, ok )
604 nt = nt + 9
605
606
607
608 srnamt = 'DSYEVD'
609 infot = 1
610 CALL dsyevd(
'/',
'U', 0, a, 1, x, w, 1, iw, 1, info )
611 CALL chkxer(
'DSYEVD', infot, nout, lerr, ok )
612 infot = 2
613 CALL dsyevd(
'N',
'/', 0, a, 1, x, w, 1, iw, 1, info )
614 CALL chkxer(
'DSYEVD', infot, nout, lerr, ok )
615 infot = 3
616 CALL dsyevd(
'N',
'U', -1, a, 1, x, w, 1, iw, 1, info )
617 CALL chkxer(
'DSYEVD', infot, nout, lerr, ok )
618 infot = 5
619 CALL dsyevd(
'N',
'U', 2, a, 1, x, w, 3, iw, 1, info )
620 CALL chkxer(
'DSYEVD', infot, nout, lerr, ok )
621 infot = 8
622 CALL dsyevd(
'N',
'U', 1, a, 1, x, w, 0, iw, 1, info )
623 CALL chkxer(
'DSYEVD', infot, nout, lerr, ok )
624 infot = 8
625 CALL dsyevd(
'N',
'U', 2, a, 2, x, w, 4, iw, 1, info )
626 CALL chkxer(
'DSYEVD', infot, nout, lerr, ok )
627 infot = 8
628 CALL dsyevd(
'V',
'U', 2, a, 2, x, w, 20, iw, 12, info )
629 CALL chkxer(
'DSYEVD', infot, nout, lerr, ok )
630 infot = 10
631 CALL dsyevd(
'N',
'U', 1, a, 1, x, w, 1, iw, 0, info )
632 CALL chkxer(
'DSYEVD', infot, nout, lerr, ok )
633 infot = 10
634 CALL dsyevd(
'N',
'U', 2, a, 2, x, w, 5, iw, 0, info )
635 CALL chkxer(
'DSYEVD', infot, nout, lerr, ok )
636 infot = 10
637 CALL dsyevd(
'V',
'U', 2, a, 2, x, w, 27, iw, 11, info )
638 CALL chkxer(
'DSYEVD', infot, nout, lerr, ok )
639 nt = nt + 10
640
641
642
643 srnamt = 'DSYEVD_2STAGE'
644 infot = 1
645 CALL dsyevd_2stage(
'/',
'U', 0, a, 1, x, w, 1, iw, 1, info )
646 CALL chkxer(
'DSYEVD_2STAGE', infot, nout, lerr, ok )
647 infot = 1
648 CALL dsyevd_2stage(
'V',
'U', 0, a, 1, x, w, 1, iw, 1, info )
649 CALL chkxer(
'DSYEVD_2STAGE', infot, nout, lerr, ok )
650 infot = 2
651 CALL dsyevd_2stage(
'N',
'/', 0, a, 1, x, w, 1, iw, 1, info )
652 CALL chkxer(
'DSYEVD_2STAGE', infot, nout, lerr, ok )
653 infot = 3
654 CALL dsyevd_2stage(
'N',
'U', -1, a, 1, x, w, 1, iw, 1, info )
655 CALL chkxer(
'DSYEVD_2STAGE', infot, nout, lerr, ok )
656 infot = 5
657 CALL dsyevd_2stage(
'N',
'U', 2, a, 1, x, w, 3, iw, 1, info )
658 CALL chkxer(
'DSYEVD_2STAGE', infot, nout, lerr, ok )
659 infot = 8
660 CALL dsyevd_2stage(
'N',
'U', 1, a, 1, x, w, 0, iw, 1, info )
661 CALL chkxer(
'DSYEVD_2STAGE', infot, nout, lerr, ok )
662 infot = 8
663 CALL dsyevd_2stage(
'N',
'U', 2, a, 2, x, w, 4, iw, 1, info )
664 CALL chkxer(
'DSYEVD_2STAGE', infot, nout, lerr, ok )
665
666
667
668 infot = 10
669 CALL dsyevd_2stage(
'N',
'U', 1, a, 1, x, w, 1, iw, 0, info )
670 CALL chkxer(
'DSYEVD_2STAGE', infot, nout, lerr, ok )
671 infot = 10
672 CALL dsyevd_2stage(
'N',
'U', 2, a, 2, x, w, 25, iw, 0, info )
673 CALL chkxer(
'DSYEVD_2STAGE', infot, nout, lerr, ok )
674
675
676
677 nt = nt + 9
678
679
680
681 srnamt = 'DSYEVR'
682 n = 1
683 infot = 1
684 CALL dsyevr(
'/',
'A',
'U', 0, a, 1, 0.0d0, 0.0d0, 1, 1, 0.0d0,
685 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
686 CALL chkxer(
'DSYEVR', infot, nout, lerr, ok )
687 infot = 2
688 CALL dsyevr(
'V',
'/',
'U', 0, a, 1, 0.0d0, 0.0d0, 1, 1, 0.0d0,
689 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
690 CALL chkxer(
'DSYEVR', infot, nout, lerr, ok )
691 infot = 3
692 CALL dsyevr(
'V',
'A',
'/', -1, a, 1, 0.0d0, 0.0d0, 1, 1,
693 $ 0.0d0, m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n,
694 $ info )
695 CALL chkxer(
'DSYEVR', infot, nout, lerr, ok )
696 infot = 4
697 CALL dsyevr(
'V',
'A',
'U', -1, a, 1, 0.0d0, 0.0d0, 1, 1,
698 $ 0.0d0, m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n,
699 $ info )
700 CALL chkxer(
'DSYEVR', infot, nout, lerr, ok )
701 infot = 6
702 CALL dsyevr(
'V',
'A',
'U', 2, a, 1, 0.0d0, 0.0d0, 1, 1, 0.0d0,
703 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
704 CALL chkxer(
'DSYEVR', infot, nout, lerr, ok )
705 infot = 8
706 CALL dsyevr(
'V',
'V',
'U', 1, a, 1, 0.0d0, 0.0d0, 1, 1, 0.0d0,
707 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
708 CALL chkxer(
'DSYEVR', infot, nout, lerr, ok )
709 infot = 9
710 CALL dsyevr(
'V',
'I',
'U', 1, a, 1, 0.0d0, 0.0d0, 0, 1, 0.0d0,
711 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
712 CALL chkxer(
'DSYEVR', infot, nout, lerr, ok )
713 infot = 10
714
715 CALL dsyevr(
'V',
'I',
'U', 2, a, 2, 0.0d0, 0.0d0, 2, 1, 0.0d0,
716 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
717 CALL chkxer(
'DSYEVR', infot, nout, lerr, ok )
718 infot = 15
719 CALL dsyevr(
'V',
'I',
'U', 1, a, 1, 0.0d0, 0.0d0, 1, 1, 0.0d0,
720 $ m, r, z, 0, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
721 CALL chkxer(
'DSYEVR', infot, nout, lerr, ok )
722 infot = 18
723 CALL dsyevr(
'V',
'I',
'U', 1, a, 1, 0.0d0, 0.0d0, 1, 1, 0.0d0,
724 $ m, r, z, 1, iw, q, 26*n-1, iw( 2*n+1 ), 10*n,
725 $ info )
726 CALL chkxer(
'DSYEVR', infot, nout, lerr, ok )
727 infot = 20
728 CALL dsyevr(
'V',
'I',
'U', 1, a, 1, 0.0d0, 0.0d0, 1, 1, 0.0d0,
729 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n-1,
730 $ info )
731 CALL chkxer(
'DSYEVR', infot, nout, lerr, ok )
732 nt = nt + 11
733
734
735
736 srnamt = 'DSYEVR_2STAGE'
737 n = 1
738 infot = 1
740 $ 0.0d0, 0.0d0, 1, 1, 0.0d0,
741 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
742 CALL chkxer(
'DSYEVR_2STAGE', infot, nout, lerr, ok )
743 infot = 1
745 $ 0.0d0, 0.0d0, 1, 1, 0.0d0,
746 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
747 CALL chkxer(
'DSYEVR_2STAGE', infot, nout, lerr, ok )
748 infot = 2
750 $ 0.0d0, 0.0d0, 1, 1, 0.0d0,
751 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
752 CALL chkxer(
'DSYEVR_2STAGE', infot, nout, lerr, ok )
753 infot = 3
755 $ 0.0d0, 0.0d0, 1, 1, 0.0d0,
756 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
757 CALL chkxer(
'DSYEVR_2STAGE', infot, nout, lerr, ok )
758 infot = 4
760 $ 0.0d0, 0.0d0, 1, 1, 0.0d0,
761 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
762 CALL chkxer(
'DSYEVR_2STAGE', infot, nout, lerr, ok )
763 infot = 6
765 $ 0.0d0, 0.0d0, 1, 1, 0.0d0,
766 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
767 CALL chkxer(
'DSYEVR_2STAGE', infot, nout, lerr, ok )
768 infot = 8
770 $ 0.0d0, 0.0d0, 1, 1, 0.0d0,
771 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
772 CALL chkxer(
'DSYEVR_2STAGE', infot, nout, lerr, ok )
773 infot = 9
775 $ 0.0d0, 0.0d0, 0, 1, 0.0d0,
776 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
777 CALL chkxer(
'DSYEVR_2STAGE', infot, nout, lerr, ok )
778 infot = 10
780 $ 0.0d0, 0.0d0, 2, 1, 0.0d0,
781 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
782 CALL chkxer(
'DSYEVR_2STAGE', infot, nout, lerr, ok )
783 infot = 15
785 $ 0.0d0, 0.0d0, 1, 1, 0.0d0,
786 $ m, r, z, 0, iw, q, 26*n, iw( 2*n+1 ), 10*n, info )
787 CALL chkxer(
'DSYEVR_2STAGE', infot, nout, lerr, ok )
788 infot = 18
790 $ 0.0d0, 0.0d0, 1, 1, 0.0d0,
791 $ m, r, z, 1, iw, q, 0, iw( 2*n+1 ), 10*n,
792 $ info )
793 CALL chkxer(
'DSYEVR_2STAGE', infot, nout, lerr, ok )
794 infot = 20
796 $ 0.0d0, 0.0d0, 1, 1, 0.0d0,
797 $ m, r, z, 1, iw, q, 26*n, iw( 2*n+1 ), 0,
798 $ info )
799 CALL chkxer(
'DSYEVR_2STAGE', infot, nout, lerr, ok )
800 nt = nt + 12
801
802
803
804 srnamt = 'DSYEV '
805 infot = 1
806 CALL dsyev(
'/',
'U', 0, a, 1, x, w, 1, info )
807 CALL chkxer(
'DSYEV ', infot, nout, lerr, ok )
808 infot = 2
809 CALL dsyev(
'N',
'/', 0, a, 1, x, w, 1, info )
810 CALL chkxer(
'DSYEV ', infot, nout, lerr, ok )
811 infot = 3
812 CALL dsyev(
'N',
'U', -1, a, 1, x, w, 1, info )
813 CALL chkxer(
'DSYEV ', infot, nout, lerr, ok )
814 infot = 5
815 CALL dsyev(
'N',
'U', 2, a, 1, x, w, 3, info )
816 CALL chkxer(
'DSYEV ', infot, nout, lerr, ok )
817 infot = 8
818 CALL dsyev(
'N',
'U', 1, a, 1, x, w, 1, info )
819 CALL chkxer(
'DSYEV ', infot, nout, lerr, ok )
820 nt = nt + 5
821
822
823
824 srnamt = 'DSYEV_2STAGE '
825 infot = 1
827 CALL chkxer(
'DSYEV_2STAGE ', infot, nout, lerr, ok )
828 infot = 1
830 CALL chkxer(
'DSYEV_2STAGE ', infot, nout, lerr, ok )
831 infot = 2
833 CALL chkxer(
'DSYEV_2STAGE ', infot, nout, lerr, ok )
834 infot = 3
836 CALL chkxer(
'DSYEV_2STAGE ', infot, nout, lerr, ok )
837 infot = 5
839 CALL chkxer(
'DSYEV_2STAGE ', infot, nout, lerr, ok )
840 infot = 8
842 CALL chkxer(
'DSYEV_2STAGE ', infot, nout, lerr, ok )
843 nt = nt + 6
844
845
846
847 srnamt = 'DSYEVX'
848 infot = 1
849 CALL dsyevx(
'/',
'A',
'U', 0, a, 1, 0.0d0, 0.0d0, 0, 0, 0.0d0,
850 $ m, x, z, 1, w, 1, iw, i3, info )
851 CALL chkxer(
'DSYEVX', infot, nout, lerr, ok )
852 infot = 2
853 CALL dsyevx(
'N',
'/',
'U', 0, a, 1, 0.0d0, 1.0d0, 1, 0, 0.0d0,
854 $ m, x, z, 1, w, 1, iw, i3, info )
855 CALL chkxer(
'DSYEVX', infot, nout, lerr, ok )
856 infot = 3
857 CALL dsyevx(
'N',
'A',
'/', 0, a, 1, 0.0d0, 0.0d0, 0, 0, 0.0d0,
858 $ m, x, z, 1, w, 1, iw, i3, info )
859 infot = 4
860 CALL dsyevx(
'N',
'A',
'U', -1, a, 1, 0.0d0, 0.0d0, 0, 0,
861 $ 0.0d0, m, x, z, 1, w, 1, iw, i3, info )
862 CALL chkxer(
'DSYEVX', infot, nout, lerr, ok )
863 infot = 6
864 CALL dsyevx(
'N',
'A',
'U', 2, a, 1, 0.0d0, 0.0d0, 0, 0, 0.0d0,
865 $ m, x, z, 1, w, 16, iw, i3, info )
866 CALL chkxer(
'DSYEVX', infot, nout, lerr, ok )
867 infot = 8
868 CALL dsyevx(
'N',
'V',
'U', 1, a, 1, 0.0d0, 0.0d0, 0, 0, 0.0d0,
869 $ m, x, z, 1, w, 8, iw, i3, info )
870 CALL chkxer(
'DSYEVX', infot, nout, lerr, ok )
871 infot = 9
872 CALL dsyevx(
'N',
'I',
'U', 1, a, 1, 0.0d0, 0.0d0, 0, 0, 0.0d0,
873 $ m, x, z, 1, w, 8, iw, i3, info )
874 CALL chkxer(
'DSYEVX', infot, nout, lerr, ok )
875 infot = 9
876 CALL dsyevx(
'N',
'I',
'U', 1, a, 1, 0.0d0, 0.0d0, 2, 1, 0.0d0,
877 $ m, x, z, 1, w, 8, iw, i3, info )
878 CALL chkxer(
'DSYEVX', infot, nout, lerr, ok )
879 infot = 10
880 CALL dsyevx(
'N',
'I',
'U', 2, a, 2, 0.0d0, 0.0d0, 2, 1, 0.0d0,
881 $ m, x, z, 1, w, 16, iw, i3, info )
882 CALL chkxer(
'DSYEVX', infot, nout, lerr, ok )
883 infot = 10
884 CALL dsyevx(
'N',
'I',
'U', 1, a, 1, 0.0d0, 0.0d0, 1, 2, 0.0d0,
885 $ m, x, z, 1, w, 8, iw, i3, info )
886 CALL chkxer(
'DSYEVX', infot, nout, lerr, ok )
887 infot = 15
888 CALL dsyevx(
'V',
'A',
'U', 2, a, 2, 0.0d0, 0.0d0, 0, 0, 0.0d0,
889 $ m, x, z, 1, w, 16, iw, i3, info )
890 CALL chkxer(
'DSYEVX', infot, nout, lerr, ok )
891 infot = 17
892 CALL dsyevx(
'V',
'A',
'U', 1, a, 1, 0.0d0, 0.0d0, 0, 0, 0.0d0,
893 $ m, x, z, 1, w, 0, iw, i3, info )
894 CALL chkxer(
'DSYEVX', infot, nout, lerr, ok )
895 nt = nt + 12
896
897
898
899 srnamt = 'DSYEVX_2STAGE'
900 infot = 1
902 $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
903 $ m, x, z, 1, w, 1, iw, i3, info )
904 CALL chkxer(
'DSYEVX_2STAGE', infot, nout, lerr, ok )
905 infot = 1
907 $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
908 $ m, x, z, 1, w, 1, iw, i3, info )
909 CALL chkxer(
'DSYEVX_2STAGE', infot, nout, lerr, ok )
910 infot = 2
912 $ 0.0d0, 1.0d0, 1, 0, 0.0d0,
913 $ m, x, z, 1, w, 1, iw, i3, info )
914 CALL chkxer(
'DSYEVX_2STAGE', infot, nout, lerr, ok )
915 infot = 3
917 $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
918 $ m, x, z, 1, w, 1, iw, i3, info )
919 infot = 4
921 $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
922 $ m, x, z, 1, w, 1, iw, i3, info )
923 CALL chkxer(
'DSYEVX_2STAGE', infot, nout, lerr, ok )
924 infot = 6
926 $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
927 $ m, x, z, 1, w, 16, iw, i3, info )
928 CALL chkxer(
'DSYEVX_2STAGE', infot, nout, lerr, ok )
929 infot = 8
931 $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
932 $ m, x, z, 1, w, 8, iw, i3, info )
933 CALL chkxer(
'DSYEVX_2STAGE', infot, nout, lerr, ok )
934 infot = 9
936 $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
937 $ m, x, z, 1, w, 8, iw, i3, info )
938 CALL chkxer(
'DSYEVX_2STAGE', infot, nout, lerr, ok )
939 infot = 9
941 $ 0.0d0, 0.0d0, 2, 1, 0.0d0,
942 $ m, x, z, 1, w, 8, iw, i3, info )
943 CALL chkxer(
'DSYEVX_2STAGE', infot, nout, lerr, ok )
944 infot = 10
946 $ 0.0d0, 0.0d0, 2, 1, 0.0d0,
947 $ m, x, z, 1, w, 16, iw, i3, info )
948 CALL chkxer(
'DSYEVX_2STAGE', infot, nout, lerr, ok )
949 infot = 10
951 $ 0.0d0, 0.0d0, 1, 2, 0.0d0,
952 $ m, x, z, 1, w, 8, iw, i3, info )
953 CALL chkxer(
'DSYEVX_2STAGE', infot, nout, lerr, ok )
954 infot = 15
956 $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
957 $ m, x, z, 0, w, 16, iw, i3, info )
958 CALL chkxer(
'DSYEVX_2STAGE', infot, nout, lerr, ok )
959 infot = 17
961 $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
962 $ m, x, z, 1, w, 0, iw, i3, info )
963 CALL chkxer(
'DSYEVX_2STAGE', infot, nout, lerr, ok )
964 nt = nt + 13
965
966
967
968 srnamt = 'DSPEVD'
969 infot = 1
970 CALL dspevd(
'/',
'U', 0, a, x, z, 1, w, 1, iw, 1, info )
971 CALL chkxer(
'DSPEVD', infot, nout, lerr, ok )
972 infot = 2
973 CALL dspevd(
'N',
'/', 0, a, x, z, 1, w, 1, iw, 1, info )
974 CALL chkxer(
'DSPEVD', infot, nout, lerr, ok )
975 infot = 3
976 CALL dspevd(
'N',
'U', -1, a, x, z, 1, w, 1, iw, 1, info )
977 CALL chkxer(
'DSPEVD', infot, nout, lerr, ok )
978 infot = 7
979 CALL dspevd(
'V',
'U', 2, a, x, z, 1, w, 23, iw, 12, info )
980 CALL chkxer(
'DSPEVD', infot, nout, lerr, ok )
981 infot = 9
982 CALL dspevd(
'N',
'U', 1, a, x, z, 1, w, 0, iw, 1, info )
983 CALL chkxer(
'DSPEVD', infot, nout, lerr, ok )
984 infot = 9
985 CALL dspevd(
'N',
'U', 2, a, x, z, 1, w, 3, iw, 1, info )
986 CALL chkxer(
'DSPEVD', infot, nout, lerr, ok )
987 infot = 9
988 CALL dspevd(
'V',
'U', 2, a, x, z, 2, w, 16, iw, 12, info )
989 CALL chkxer(
'DSPEVD', infot, nout, lerr, ok )
990 infot = 11
991 CALL dspevd(
'N',
'U', 1, a, x, z, 1, w, 1, iw, 0, info )
992 CALL chkxer(
'DSPEVD', infot, nout, lerr, ok )
993 infot = 11
994 CALL dspevd(
'N',
'U', 2, a, x, z, 1, w, 4, iw, 0, info )
995 CALL chkxer(
'DSPEVD', infot, nout, lerr, ok )
996 infot = 11
997 CALL dspevd(
'V',
'U', 2, a, x, z, 2, w, 23, iw, 11, info )
998 CALL chkxer(
'DSPEVD', infot, nout, lerr, ok )
999 nt = nt + 10
1000
1001
1002
1003 srnamt = 'DSPEV '
1004 infot = 1
1005 CALL dspev(
'/',
'U', 0, a, w, z, 1, x, info )
1006 CALL chkxer(
'DSPEV ', infot, nout, lerr, ok )
1007 infot = 2
1008 CALL dspev(
'N',
'/', 0, a, w, z, 1, x, info )
1009 CALL chkxer(
'DSPEV ', infot, nout, lerr, ok )
1010 infot = 3
1011 CALL dspev(
'N',
'U', -1, a, w, z, 1, x, info )
1012 CALL chkxer(
'DSPEV ', infot, nout, lerr, ok )
1013 infot = 7
1014 CALL dspev(
'V',
'U', 2, a, w, z, 1, x, info )
1015 CALL chkxer(
'DSPEV ', infot, nout, lerr, ok )
1016 nt = nt + 4
1017
1018
1019
1020 srnamt = 'DSPEVX'
1021 infot = 1
1022 CALL dspevx(
'/',
'A',
'U', 0, a, 0.0d0, 0.0d0, 0, 0, 0.0d0, m,
1023 $ x, z, 1, w, iw, i3, info )
1024 CALL chkxer(
'DSPEVX', infot, nout, lerr, ok )
1025 infot = 2
1026 CALL dspevx(
'N',
'/',
'U', 0, a, 0.0d0, 0.0d0, 0, 0, 0.0d0, m,
1027 $ x, z, 1, w, iw, i3, info )
1028 CALL chkxer(
'DSPEVX', infot, nout, lerr, ok )
1029 infot = 3
1030 CALL dspevx(
'N',
'A',
'/', 0, a, 0.0d0, 0.0d0, 0, 0, 0.0d0, m,
1031 $ x, z, 1, w, iw, i3, info )
1032 infot = 4
1033 CALL dspevx(
'N',
'A',
'U', -1, a, 0.0d0, 0.0d0, 0, 0, 0.0d0,
1034 $ m, x, z, 1, w, iw, i3, info )
1035 CALL chkxer(
'DSPEVX', infot, nout, lerr, ok )
1036 infot = 7
1037 CALL dspevx(
'N',
'V',
'U', 1, a, 0.0d0, 0.0d0, 0, 0, 0.0d0, m,
1038 $ x, z, 1, w, iw, i3, info )
1039 CALL chkxer(
'DSPEVX', infot, nout, lerr, ok )
1040 infot = 8
1041 CALL dspevx(
'N',
'I',
'U', 1, a, 0.0d0, 0.0d0, 0, 0, 0.0d0, m,
1042 $ x, z, 1, w, iw, i3, info )
1043 CALL chkxer(
'DSPEVX', infot, nout, lerr, ok )
1044 infot = 8
1045 CALL dspevx(
'N',
'I',
'U', 1, a, 0.0d0, 0.0d0, 2, 1, 0.0d0, m,
1046 $ x, z, 1, w, iw, i3, info )
1047 CALL chkxer(
'DSPEVX', infot, nout, lerr, ok )
1048 infot = 9
1049 CALL dspevx(
'N',
'I',
'U', 2, a, 0.0d0, 0.0d0, 2, 1, 0.0d0, m,
1050 $ x, z, 1, w, iw, i3, info )
1051 CALL chkxer(
'DSPEVX', infot, nout, lerr, ok )
1052 infot = 9
1053 CALL dspevx(
'N',
'I',
'U', 1, a, 0.0d0, 0.0d0, 1, 2, 0.0d0, m,
1054 $ x, z, 1, w, iw, i3, info )
1055 CALL chkxer(
'DSPEVX', infot, nout, lerr, ok )
1056 infot = 14
1057 CALL dspevx(
'V',
'A',
'U', 2, a, 0.0d0, 0.0d0, 0, 0, 0.0d0, m,
1058 $ x, z, 1, w, iw, i3, info )
1059 CALL chkxer(
'DSPEVX', infot, nout, lerr, ok )
1060 nt = nt + 10
1061
1062
1063
1064 ELSE IF(
lsamen( 2, c2,
'SB' ) )
THEN
1065
1066
1067
1068 srnamt = 'DSBTRD'
1069 infot = 1
1070 CALL dsbtrd(
'/',
'U', 0, 0, a, 1, d, e, z, 1, w, info )
1071 CALL chkxer(
'DSBTRD', infot, nout, lerr, ok )
1072 infot = 2
1073 CALL dsbtrd(
'N',
'/', 0, 0, a, 1, d, e, z, 1, w, info )
1074 CALL chkxer(
'DSBTRD', infot, nout, lerr, ok )
1075 infot = 3
1076 CALL dsbtrd(
'N',
'U', -1, 0, a, 1, d, e, z, 1, w, info )
1077 CALL chkxer(
'DSBTRD', infot, nout, lerr, ok )
1078 infot = 4
1079 CALL dsbtrd(
'N',
'U', 0, -1, a, 1, d, e, z, 1, w, info )
1080 CALL chkxer(
'DSBTRD', infot, nout, lerr, ok )
1081 infot = 6
1082 CALL dsbtrd(
'N',
'U', 1, 1, a, 1, d, e, z, 1, w, info )
1083 CALL chkxer(
'DSBTRD', infot, nout, lerr, ok )
1084 infot = 10
1085 CALL dsbtrd(
'V',
'U', 2, 0, a, 1, d, e, z, 1, w, info )
1086 CALL chkxer(
'DSBTRD', infot, nout, lerr, ok )
1087 nt = nt + 6
1088
1089
1090
1091 srnamt = 'DSYTRD_SB2ST'
1092 infot = 1
1094 $ c, 1, w, 1, info )
1095 CALL chkxer(
'DSYTRD_SB2ST', infot, nout, lerr, ok )
1096 infot = 2
1098 $ c, 1, w, 1, info )
1099 CALL chkxer(
'DSYTRD_SB2ST', infot, nout, lerr, ok )
1100 infot = 2
1102 $ c, 1, w, 1, info )
1103 CALL chkxer(
'DSYTRD_SB2ST', infot, nout, lerr, ok )
1104 infot = 3
1106 $ c, 1, w, 1, info )
1107 CALL chkxer(
'DSYTRD_SB2ST', infot, nout, lerr, ok )
1108 infot = 4
1110 $ c, 1, w, 1, info )
1111 CALL chkxer(
'DSYTRD_SB2ST', infot, nout, lerr, ok )
1112 infot = 5
1114 $ c, 1, w, 1, info )
1115 CALL chkxer(
'DSYTRD_SB2ST', infot, nout, lerr, ok )
1116 infot = 7
1118 $ c, 1, w, 1, info )
1119 CALL chkxer(
'DSYTRD_SB2ST', infot, nout, lerr, ok )
1120 infot = 11
1122 $ c, 0, w, 1, info )
1123 CALL chkxer(
'DSYTRD_SB2ST', infot, nout, lerr, ok )
1124 infot = 13
1126 $ c, 1, w, 0, info )
1127 CALL chkxer(
'DSYTRD_SB2ST', infot, nout, lerr, ok )
1128 nt = nt + 9
1129
1130
1131
1132 srnamt = 'DSBEVD'
1133 infot = 1
1134 CALL dsbevd(
'/',
'U', 0, 0, a, 1, x, z, 1, w, 1, iw, 1, info )
1135 CALL chkxer(
'DSBEVD', infot, nout, lerr, ok )
1136 infot = 2
1137 CALL dsbevd(
'N',
'/', 0, 0, a, 1, x, z, 1, w, 1, iw, 1, info )
1138 CALL chkxer(
'DSBEVD', infot, nout, lerr, ok )
1139 infot = 3
1140 CALL dsbevd(
'N',
'U', -1, 0, a, 1, x, z, 1, w, 1, iw, 1,
1141 $ info )
1142 CALL chkxer(
'DSBEVD', infot, nout, lerr, ok )
1143 infot = 4
1144 CALL dsbevd(
'N',
'U', 0, -1, a, 1, x, z, 1, w, 1, iw, 1,
1145 $ info )
1146 CALL chkxer(
'DSBEVD', infot, nout, lerr, ok )
1147 infot = 6
1148 CALL dsbevd(
'N',
'U', 2, 1, a, 1, x, z, 1, w, 4, iw, 1, info )
1149 CALL chkxer(
'DSBEVD', infot, nout, lerr, ok )
1150 infot = 9
1151 CALL dsbevd(
'V',
'U', 2, 1, a, 2, x, z, 1, w, 25, iw, 12,
1152 $ info )
1153 CALL chkxer(
'DSBEVD', infot, nout, lerr, ok )
1154 infot = 11
1155 CALL dsbevd(
'N',
'U', 1, 0, a, 1, x, z, 1, w, 0, iw, 1, info )
1156 CALL chkxer(
'DSBEVD', infot, nout, lerr, ok )
1157 infot = 11
1158 CALL dsbevd(
'N',
'U', 2, 0, a, 1, x, z, 1, w, 3, iw, 1, info )
1159 CALL chkxer(
'DSBEVD', infot, nout, lerr, ok )
1160 infot = 11
1161 CALL dsbevd(
'V',
'U', 2, 0, a, 1, x, z, 2, w, 18, iw, 12,
1162 $ info )
1163 CALL chkxer(
'DSBEVD', infot, nout, lerr, ok )
1164 infot = 13
1165 CALL dsbevd(
'N',
'U', 1, 0, a, 1, x, z, 1, w, 1, iw, 0, info )
1166 CALL chkxer(
'DSBEVD', infot, nout, lerr, ok )
1167 infot = 13
1168 CALL dsbevd(
'V',
'U', 2, 0, a, 1, x, z, 2, w, 25, iw, 11,
1169 $ info )
1170 CALL chkxer(
'DSBEVD', infot, nout, lerr, ok )
1171 nt = nt + 11
1172
1173
1174
1175 srnamt = 'DSBEVD_2STAGE'
1176 infot = 1
1177 CALL dsbevd_2stage(
'/',
'U', 0, 0, a, 1, x, z, 1, w,
1178 $ 1, iw, 1, info )
1179 CALL chkxer(
'DSBEVD_2STAGE', infot, nout, lerr, ok )
1180 infot = 1
1181 CALL dsbevd_2stage(
'V',
'U', 0, 0, a, 1, x, z, 1, w,
1182 $ 1, iw, 1, info )
1183 CALL chkxer(
'DSBEVD_2STAGE', infot, nout, lerr, ok )
1184 infot = 2
1185 CALL dsbevd_2stage(
'N',
'/', 0, 0, a, 1, x, z, 1, w,
1186 $ 1, iw, 1, info )
1187 CALL chkxer(
'DSBEVD_2STAGE', infot, nout, lerr, ok )
1188 infot = 3
1189 CALL dsbevd_2stage(
'N',
'U', -1, 0, a, 1, x, z, 1, w,
1190 $ 1, iw, 1, info )
1191 CALL chkxer(
'DSBEVD_2STAGE', infot, nout, lerr, ok )
1192 infot = 4
1193 CALL dsbevd_2stage(
'N',
'U', 0, -1, a, 1, x, z, 1, w,
1194 $ 1, iw, 1, info )
1195 CALL chkxer(
'DSBEVD_2STAGE', infot, nout, lerr, ok )
1196 infot = 6
1197 CALL dsbevd_2stage(
'N',
'U', 2, 1, a, 1, x, z, 1, w,
1198 $ 4, iw, 1, info )
1199 CALL chkxer(
'DSBEVD_2STAGE', infot, nout, lerr, ok )
1200
1201
1202
1203
1204 infot = 11
1205 CALL dsbevd_2stage(
'N',
'U', 1, 0, a, 1, x, z, 1, w,
1206 $ 0, iw, 1, info )
1207 CALL chkxer(
'DSBEVD_2STAGE', infot, nout, lerr, ok )
1208 infot = 11
1209 CALL dsbevd_2stage(
'N',
'U', 2, 0, a, 1, x, z, 1, w,
1210 $ 3, iw, 1, info )
1211 CALL chkxer(
'DSBEVD_2STAGE', infot, nout, lerr, ok )
1212
1213
1214
1215
1216 infot = 13
1217 CALL dsbevd_2stage(
'N',
'U', 1, 0, a, 1, x, z, 1, w,
1218 $ 1, iw, 0, info )
1219 CALL chkxer(
'DSBEVD_2STAGE', infot, nout, lerr, ok )
1220
1221
1222
1223
1224
1225 nt = nt + 9
1226
1227
1228
1229 srnamt = 'DSBEV '
1230 infot = 1
1231 CALL dsbev(
'/',
'U', 0, 0, a, 1, x, z, 1, w, info )
1232 CALL chkxer(
'DSBEV ', infot, nout, lerr, ok )
1233 infot = 2
1234 CALL dsbev(
'N',
'/', 0, 0, a, 1, x, z, 1, w, info )
1235 CALL chkxer(
'DSBEV ', infot, nout, lerr, ok )
1236 infot = 3
1237 CALL dsbev(
'N',
'U', -1, 0, a, 1, x, z, 1, w, info )
1238 CALL chkxer(
'DSBEV ', infot, nout, lerr, ok )
1239 infot = 4
1240 CALL dsbev(
'N',
'U', 0, -1, a, 1, x, z, 1, w, info )
1241 CALL chkxer(
'DSBEV ', infot, nout, lerr, ok )
1242 infot = 6
1243 CALL dsbev(
'N',
'U', 2, 1, a, 1, x, z, 1, w, info )
1244 CALL chkxer(
'DSBEV ', infot, nout, lerr, ok )
1245 infot = 9
1246 CALL dsbev(
'V',
'U', 2, 0, a, 1, x, z, 1, w, info )
1247 CALL chkxer(
'DSBEV ', infot, nout, lerr, ok )
1248 nt = nt + 6
1249
1250
1251
1252 srnamt = 'DSBEV_2STAGE '
1253 infot = 1
1254 CALL dsbev_2stage(
'/',
'U', 0, 0, a, 1, x, z, 1, w, 0, info )
1255 CALL chkxer(
'DSBEV_2STAGE ', infot, nout, lerr, ok )
1256 infot = 1
1257 CALL dsbev_2stage(
'V',
'U', 0, 0, a, 1, x, z, 1, w, 0, info )
1258 CALL chkxer(
'DSBEV_2STAGE ', infot, nout, lerr, ok )
1259 infot = 2
1260 CALL dsbev_2stage(
'N',
'/', 0, 0, a, 1, x, z, 1, w, 0, info )
1261 CALL chkxer(
'DSBEV_2STAGE ', infot, nout, lerr, ok )
1262 infot = 3
1263 CALL dsbev_2stage(
'N',
'U', -1, 0, a, 1, x, z, 1, w, 0, info )
1264 CALL chkxer(
'DSBEV_2STAGE ', infot, nout, lerr, ok )
1265 infot = 4
1266 CALL dsbev_2stage(
'N',
'U', 0, -1, a, 1, x, z, 1, w, 0, info )
1267 CALL chkxer(
'DSBEV_2STAGE ', infot, nout, lerr, ok )
1268 infot = 6
1269 CALL dsbev_2stage(
'N',
'U', 2, 1, a, 1, x, z, 1, w, 0, info )
1270 CALL chkxer(
'DSBEV_2STAGE ', infot, nout, lerr, ok )
1271 infot = 9
1272 CALL dsbev_2stage(
'N',
'U', 2, 0, a, 1, x, z, 0, w, 0, info )
1273 CALL chkxer(
'DSBEV_2STAGE ', infot, nout, lerr, ok )
1274 infot = 11
1275 CALL dsbev_2stage(
'N',
'U', 0, 0, a, 1, x, z, 1, w, 0, info )
1276 CALL chkxer(
'DSBEV_2STAGE ', infot, nout, lerr, ok )
1277 nt = nt + 8
1278
1279
1280
1281 srnamt = 'DSBEVX'
1282 infot = 1
1283 CALL dsbevx(
'/',
'A',
'U', 0, 0, a, 1, q, 1, 0.0d0, 0.0d0, 0,
1284 $ 0, 0.0d0, m, x, z, 1, w, iw, i3, info )
1285 CALL chkxer(
'DSBEVX', infot, nout, lerr, ok )
1286 infot = 2
1287 CALL dsbevx(
'N',
'/',
'U', 0, 0, a, 1, q, 1, 0.0d0, 0.0d0, 0,
1288 $ 0, 0.0d0, m, x, z, 1, w, iw, i3, info )
1289 CALL chkxer(
'DSBEVX', infot, nout, lerr, ok )
1290 infot = 3
1291 CALL dsbevx(
'N',
'A',
'/', 0, 0, a, 1, q, 1, 0.0d0, 0.0d0, 0,
1292 $ 0, 0.0d0, m, x, z, 1, w, iw, i3, info )
1293 CALL chkxer(
'DSBEVX', infot, nout, lerr, ok )
1294 infot = 4
1295 CALL dsbevx(
'N',
'A',
'U', -1, 0, a, 1, q, 1, 0.0d0, 0.0d0, 0,
1296 $ 0, 0.0d0, m, x, z, 1, w, iw, i3, info )
1297 CALL chkxer(
'DSBEVX', infot, nout, lerr, ok )
1298 infot = 5
1299 CALL dsbevx(
'N',
'A',
'U', 0, -1, a, 1, q, 1, 0.0d0, 0.0d0, 0,
1300 $ 0, 0.0d0, m, x, z, 1, w, iw, i3, info )
1301 CALL chkxer(
'DSBEVX', infot, nout, lerr, ok )
1302 infot = 7
1303 CALL dsbevx(
'N',
'A',
'U', 2, 1, a, 1, q, 1, 0.0d0, 0.0d0, 0,
1304 $ 0, 0.0d0, m, x, z, 1, w, iw, i3, info )
1305 CALL chkxer(
'DSBEVX', infot, nout, lerr, ok )
1306 infot = 9
1307 CALL dsbevx(
'V',
'A',
'U', 2, 0, a, 1, q, 1, 0.0d0, 0.0d0, 0,
1308 $ 0, 0.0d0, m, x, z, 2, w, iw, i3, info )
1309 CALL chkxer(
'DSBEVX', infot, nout, lerr, ok )
1310 infot = 11
1311 CALL dsbevx(
'N',
'V',
'U', 1, 0, a, 1, q, 1, 0.0d0, 0.0d0, 0,
1312 $ 0, 0.0d0, m, x, z, 1, w, iw, i3, info )
1313 CALL chkxer(
'DSBEVX', infot, nout, lerr, ok )
1314 infot = 12
1315 CALL dsbevx(
'N',
'I',
'U', 1, 0, a, 1, q, 1, 0.0d0, 0.0d0, 0,
1316 $ 0, 0.0d0, m, x, z, 1, w, iw, i3, info )
1317 CALL chkxer(
'DSBEVX', infot, nout, lerr, ok )
1318 infot = 12
1319 CALL dsbevx(
'N',
'I',
'U', 1, 0, a, 1, q, 1, 0.0d0, 0.0d0, 2,
1320 $ 1, 0.0d0, m, x, z, 1, w, iw, i3, info )
1321 CALL chkxer(
'DSBEVX', infot, nout, lerr, ok )
1322 infot = 13
1323 CALL dsbevx(
'N',
'I',
'U', 2, 0, a, 1, q, 1, 0.0d0, 0.0d0, 2,
1324 $ 1, 0.0d0, m, x, z, 1, w, iw, i3, info )
1325 CALL chkxer(
'DSBEVX', infot, nout, lerr, ok )
1326 infot = 13
1327 CALL dsbevx(
'N',
'I',
'U', 1, 0, a, 1, q, 1, 0.0d0, 0.0d0, 1,
1328 $ 2, 0.0d0, m, x, z, 1, w, iw, i3, info )
1329 CALL chkxer(
'DSBEVX', infot, nout, lerr, ok )
1330 infot = 18
1331 CALL dsbevx(
'V',
'A',
'U', 2, 0, a, 1, q, 2, 0.0d0, 0.0d0, 0,
1332 $ 0, 0.0d0, m, x, z, 1, w, iw, i3, info )
1333 CALL chkxer(
'DSBEVX', infot, nout, lerr, ok )
1334 nt = nt + 13
1335
1336
1337
1338 srnamt = 'DSBEVX_2STAGE'
1339 infot = 1
1340 CALL dsbevx_2stage(
'/',
'A',
'U', 0, 0, a, 1, q, 1, 0.0d0,
1341 $ 0.0d0, 0, 0, 0.0d0, m, x, z, 1, w, 0, iw, i3, info )
1342 CALL chkxer(
'DSBEVX_2STAGE', infot, nout, lerr, ok )
1343 infot = 1
1344 CALL dsbevx_2stage(
'V',
'A',
'U', 0, 0, a, 1, q, 1, 0.0d0,
1345 $ 0.0d0, 0, 0, 0.0d0, m, x, z, 1, w, 0, iw, i3, info )
1346 CALL chkxer(
'DSBEVX_2STAGE', infot, nout, lerr, ok )
1347 infot = 2
1348 CALL dsbevx_2stage(
'N',
'/',
'U', 0, 0, a, 1, q, 1, 0.0d0,
1349 $ 0.0d0, 0, 0, 0.0d0, m, x, z, 1, w, 0, iw, i3, info )
1350 CALL chkxer(
'DSBEVX_2STAGE', infot, nout, lerr, ok )
1351 infot = 3
1352 CALL dsbevx_2stage(
'N',
'A',
'/', 0, 0, a, 1, q, 1, 0.0d0,
1353 $ 0.0d0, 0, 0, 0.0d0, m, x, z, 1, w, 0, iw, i3, info )
1354 CALL chkxer(
'DSBEVX_2STAGE', infot, nout, lerr, ok )
1355 infot = 4
1356 CALL dsbevx_2stage(
'N',
'A',
'U', -1, 0, a, 1, q, 1, 0.0d0,
1357 $ 0.0d0, 0, 0, 0.0d0, m, x, z, 1, w, 0, iw, i3, info )
1358 CALL chkxer(
'DSBEVX_2STAGE', infot, nout, lerr, ok )
1359 infot = 5
1360 CALL dsbevx_2stage(
'N',
'A',
'U', 0, -1, a, 1, q, 1, 0.0d0,
1361 $ 0.0d0, 0, 0, 0.0d0, m, x, z, 1, w, 0, iw, i3, info )
1362 CALL chkxer(
'DSBEVX_2STAGE', infot, nout, lerr, ok )
1363 infot = 7
1364 CALL dsbevx_2stage(
'N',
'A',
'U', 2, 1, a, 1, q, 1, 0.0d0,
1365 $ 0.0d0, 0, 0, 0.0d0, m, x, z, 1, w, 0, iw, i3, info )
1366 CALL chkxer(
'DSBEVX_2STAGE', infot, nout, lerr, ok )
1367
1368
1369
1370
1371 infot = 11
1372 CALL dsbevx_2stage(
'N',
'V',
'U', 1, 0, a, 1, q, 1, 0.0d0,
1373 $ 0.0d0, 0, 0, 0.0d0, m, x, z, 1, w, 0, iw, i3, info )
1374 CALL chkxer(
'DSBEVX_2STAGE', infot, nout, lerr, ok )
1375 infot = 12
1376 CALL dsbevx_2stage(
'N',
'I',
'U', 1, 0, a, 1, q, 1, 0.0d0,
1377 $ 0.0d0, 0, 0, 0.0d0, m, x, z, 1, w, 0, iw, i3, info )
1378 CALL chkxer(
'DSBEVX_2STAGE', infot, nout, lerr, ok )
1379 infot = 12
1380 CALL dsbevx_2stage(
'N',
'I',
'U', 1, 0, a, 1, q, 1, 0.0d0,
1381 $ 0.0d0, 2, 1, 0.0d0, m, x, z, 1, w, 0, iw, i3, info )
1382 CALL chkxer(
'DSBEVX_2STAGE', infot, nout, lerr, ok )
1383 infot = 13
1384 CALL dsbevx_2stage(
'N',
'I',
'U', 2, 0, a, 1, q, 1, 0.0d0,
1385 $ 0.0d0, 2, 1, 0.0d0, m, x, z, 1, w, 0, iw, i3, info )
1386 CALL chkxer(
'DSBEVX_2STAGE', infot, nout, lerr, ok )
1387 infot = 13
1388 CALL dsbevx_2stage(
'N',
'I',
'U', 1, 0, a, 1, q, 1, 0.0d0,
1389 $ 0.0d0, 1, 2, 0.0d0, m, x, z, 1, w, 0, iw, i3, info )
1390 CALL chkxer(
'DSBEVX_2STAGE', infot, nout, lerr, ok )
1391
1392
1393
1394
1395 infot = 20
1396 CALL dsbevx_2stage(
'N',
'A',
'U', 0, 0, a, 1, q, 1, 0.0d0,
1397 $ 0.0d0, 0, 0, 0.0d0, m, x, z, 1, w, 0, iw, i3, info )
1398 CALL chkxer(
'DSBEVX_2STAGE', infot, nout, lerr, ok )
1399
1400 nt = nt + 13
1401 END IF
1402
1403
1404
1405 IF( ok ) THEN
1406 WRITE( nout, fmt = 9999 )path, nt
1407 ELSE
1408 WRITE( nout, fmt = 9998 )path
1409 END IF
1410
1411 9999 FORMAT( 1x, a3, ' routines passed the tests of the error exits',
1412 $ ' (', i3, ' tests done)' )
1413 9998 FORMAT( ' *** ', a3, ' routines failed the tests of the error ',
1414 $ 'exits ***' )
1415
1416 RETURN
1417
1418
1419
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine dsytrd_sb2st(STAGE1, VECT, UPLO, N, KD, AB, LDAB, D, E, HOUS, LHOUS, WORK, LWORK, INFO)
DSYTRD_SB2ST reduces a real symmetric band matrix A to real symmetric tridiagonal form T
logical function lsamen(N, CA, CB)
LSAMEN
subroutine dstebz(RANGE, ORDER, N, VL, VU, IL, IU, ABSTOL, D, E, M, NSPLIT, W, IBLOCK, ISPLIT, WORK, IWORK, INFO)
DSTEBZ
subroutine dsteqr(COMPZ, N, D, E, Z, LDZ, WORK, INFO)
DSTEQR
subroutine dstedc(COMPZ, N, D, E, Z, LDZ, WORK, LWORK, IWORK, LIWORK, INFO)
DSTEDC
subroutine dsterf(N, D, E, INFO)
DSTERF
subroutine dsptrd(UPLO, N, AP, D, E, TAU, INFO)
DSPTRD
subroutine dstein(N, D, E, M, W, IBLOCK, ISPLIT, Z, LDZ, WORK, IWORK, IFAIL, INFO)
DSTEIN
subroutine dopgtr(UPLO, N, AP, TAU, Q, LDQ, WORK, INFO)
DOPGTR
subroutine dsbtrd(VECT, UPLO, N, KD, AB, LDAB, D, E, Q, LDQ, WORK, INFO)
DSBTRD
subroutine dormtr(SIDE, UPLO, TRANS, M, N, A, LDA, TAU, C, LDC, WORK, LWORK, INFO)
DORMTR
subroutine dorgtr(UPLO, N, A, LDA, TAU, WORK, LWORK, INFO)
DORGTR
subroutine dopmtr(SIDE, UPLO, TRANS, M, N, AP, TAU, C, LDC, WORK, INFO)
DOPMTR
subroutine dspevx(JOBZ, RANGE, UPLO, N, AP, VL, VU, IL, IU, ABSTOL, M, W, Z, LDZ, WORK, IWORK, IFAIL, INFO)
DSPEVX computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
subroutine dstevd(JOBZ, N, D, E, Z, LDZ, WORK, LWORK, IWORK, LIWORK, INFO)
DSTEVD computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
subroutine dstevr(JOBZ, RANGE, N, D, E, VL, VU, IL, IU, ABSTOL, M, W, Z, LDZ, ISUPPZ, WORK, LWORK, IWORK, LIWORK, INFO)
DSTEVR computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
subroutine dsbevx(JOBZ, RANGE, UPLO, N, KD, AB, LDAB, Q, LDQ, VL, VU, IL, IU, ABSTOL, M, W, Z, LDZ, WORK, IWORK, IFAIL, INFO)
DSBEVX computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
subroutine dspevd(JOBZ, UPLO, N, AP, W, Z, LDZ, WORK, LWORK, IWORK, LIWORK, INFO)
DSPEVD computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
subroutine dstevx(JOBZ, RANGE, N, D, E, VL, VU, IL, IU, ABSTOL, M, W, Z, LDZ, WORK, IWORK, IFAIL, INFO)
DSTEVX computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
subroutine dsbevd_2stage(JOBZ, UPLO, N, KD, AB, LDAB, W, Z, LDZ, WORK, LWORK, IWORK, LIWORK, INFO)
DSBEVD_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER ...
subroutine dspev(JOBZ, UPLO, N, AP, W, Z, LDZ, WORK, INFO)
DSPEV computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrices
subroutine dsbev_2stage(JOBZ, UPLO, N, KD, AB, LDAB, W, Z, LDZ, WORK, LWORK, INFO)
DSBEV_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER m...
subroutine dstev(JOBZ, N, D, E, Z, LDZ, WORK, INFO)
DSTEV computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrices
subroutine dsbevx_2stage(JOBZ, RANGE, UPLO, N, KD, AB, LDAB, Q, LDQ, VL, VU, IL, IU, ABSTOL, M, W, Z, LDZ, WORK, LWORK, IWORK, IFAIL, INFO)
DSBEVX_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER ...
subroutine dsbev(JOBZ, UPLO, N, KD, AB, LDAB, W, Z, LDZ, WORK, INFO)
DSBEV computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrices
subroutine dsbevd(JOBZ, UPLO, N, KD, AB, LDAB, W, Z, LDZ, WORK, LWORK, IWORK, LIWORK, INFO)
DSBEVD computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
subroutine dpteqr(COMPZ, N, D, E, Z, LDZ, WORK, INFO)
DPTEQR
subroutine dsytrd_sy2sb(UPLO, N, KD, A, LDA, AB, LDAB, TAU, WORK, LWORK, INFO)
DSYTRD_SY2SB
subroutine dsytrd_2stage(VECT, UPLO, N, A, LDA, D, E, TAU, HOUS2, LHOUS2, WORK, LWORK, INFO)
DSYTRD_2STAGE
subroutine dsytrd(UPLO, N, A, LDA, D, E, TAU, WORK, LWORK, INFO)
DSYTRD
subroutine dsyev(JOBZ, UPLO, N, A, LDA, W, WORK, LWORK, INFO)
DSYEV computes the eigenvalues and, optionally, the left and/or right eigenvectors for SY matrices
subroutine dsyevx(JOBZ, RANGE, UPLO, N, A, LDA, VL, VU, IL, IU, ABSTOL, M, W, Z, LDZ, WORK, LWORK, IWORK, IFAIL, INFO)
DSYEVX computes the eigenvalues and, optionally, the left and/or right eigenvectors for SY matrices
subroutine dsyevd(JOBZ, UPLO, N, A, LDA, W, WORK, LWORK, IWORK, LIWORK, INFO)
DSYEVD computes the eigenvalues and, optionally, the left and/or right eigenvectors for SY matrices
subroutine dsyevr_2stage(JOBZ, RANGE, UPLO, N, A, LDA, VL, VU, IL, IU, ABSTOL, M, W, Z, LDZ, ISUPPZ, WORK, LWORK, IWORK, LIWORK, INFO)
DSYEVR_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for SY mat...
subroutine dsyevx_2stage(JOBZ, RANGE, UPLO, N, A, LDA, VL, VU, IL, IU, ABSTOL, M, W, Z, LDZ, WORK, LWORK, IWORK, IFAIL, INFO)
DSYEVX_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for SY mat...
subroutine dsyev_2stage(JOBZ, UPLO, N, A, LDA, W, WORK, LWORK, INFO)
DSYEV_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for SY matr...
subroutine dsyevr(JOBZ, RANGE, UPLO, N, A, LDA, VL, VU, IL, IU, ABSTOL, M, W, Z, LDZ, ISUPPZ, WORK, LWORK, IWORK, LIWORK, INFO)
DSYEVR computes the eigenvalues and, optionally, the left and/or right eigenvectors for SY matrices
subroutine dsyevd_2stage(JOBZ, UPLO, N, A, LDA, W, WORK, LWORK, IWORK, LIWORK, INFO)
DSYEVD_2STAGE computes the eigenvalues and, optionally, the left and/or right eigenvectors for SY mat...