LAPACK  3.10.1
LAPACK: Linear Algebra PACKage

◆ derrst()

subroutine derrst ( character*3  PATH,
integer  NUNIT 
)

DERRST

Purpose:
 DERRST tests the error exits for DSYTRD, DORGTR, DORMTR, DSPTRD,
 DOPGTR, DOPMTR, DSTEQR, SSTERF, SSTEBZ, SSTEIN, DPTEQR, DSBTRD,
 DSYEV, SSYEVX, SSYEVD, DSBEV, SSBEVX, SSBEVD,
 DSPEV, SSPEVX, SSPEVD, DSTEV, SSTEVX, SSTEVD, and SSTEDC.
 DSYEVD_2STAGE, DSYEVR_2STAGE, DSYEVX_2STAGE,
 DSYEV_2STAGE, DSBEV_2STAGE, DSBEVD_2STAGE,
 DSBEVX_2STAGE, DSYTRD_2STAGE, DSYTRD_SY2SB,
 DSYTRD_SB2ST
Parameters
[in]PATH
          PATH is CHARACTER*3
          The LAPACK path name for the routines to be tested.
[in]NUNIT
          NUNIT is INTEGER
          The unit number for output.
Author
Univ. of Tennessee
Univ. of California Berkeley
Univ. of Colorado Denver
NAG Ltd.

Definition at line 60 of file derrst.f.

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  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 * .. External Functions ..
90  LOGICAL LSAMEN
91  EXTERNAL lsamen
92 * ..
93 * .. External Subroutines ..
94  EXTERNAL chkxer, dopgtr, dopmtr, dorgtr, dormtr, dpteqr,
95  $ dsbev, dsbevd, dsbevx, dsbtrd, dspev, dspevd,
96  $ dspevx, dsptrd, dstebz, dstedc, dstein, dsteqr,
97  $ dsterf, dstev, dstevd, dstevr, dstevx, dsyev,
98  $ dsyevd, dsyevr, dsyevx, dsytrd,
99  $ dsyevd_2stage, dsyevr_2stage, dsyevx_2stage,
100  $ dsyev_2stage, dsbev_2stage, dsbevd_2stage,
101  $ dsbevx_2stage, dsytrd_2stage, dsytrd_sy2sb,
102  $ dsytrd_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 dble
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.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 * Test error exits for the ST path.
140 *
141  IF( lsamen( 2, c2, 'ST' ) ) THEN
142 *
143 * DSYTRD
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 * DSYTRD_2STAGE
161 *
162  srnamt = 'DSYTRD_2STAGE'
163  infot = 1
164  CALL dsytrd_2stage( '/', 'U', 0, a, 1, d, e, tau,
165  $ c, 1, w, 1, info )
166  CALL chkxer( 'DSYTRD_2STAGE', infot, nout, lerr, ok )
167  infot = 1
168  CALL dsytrd_2stage( 'H', 'U', 0, a, 1, d, e, tau,
169  $ c, 1, w, 1, info )
170  CALL chkxer( 'DSYTRD_2STAGE', infot, nout, lerr, ok )
171  infot = 2
172  CALL dsytrd_2stage( 'N', '/', 0, a, 1, d, e, tau,
173  $ c, 1, w, 1, info )
174  CALL chkxer( 'DSYTRD_2STAGE', infot, nout, lerr, ok )
175  infot = 3
176  CALL dsytrd_2stage( 'N', 'U', -1, a, 1, d, e, tau,
177  $ c, 1, w, 1, info )
178  CALL chkxer( 'DSYTRD_2STAGE', infot, nout, lerr, ok )
179  infot = 5
180  CALL dsytrd_2stage( 'N', 'U', 2, a, 1, d, e, tau,
181  $ c, 1, w, 1, info )
182  CALL chkxer( 'DSYTRD_2STAGE', infot, nout, lerr, ok )
183  infot = 10
184  CALL dsytrd_2stage( 'N', 'U', 0, a, 1, d, e, tau,
185  $ c, 0, w, 1, info )
186  CALL chkxer( 'DSYTRD_2STAGE', infot, nout, lerr, ok )
187  infot = 12
188  CALL dsytrd_2stage( 'N', 'U', 0, a, 1, d, e, tau,
189  $ c, 1, w, 0, info )
190  CALL chkxer( 'DSYTRD_2STAGE', infot, nout, lerr, ok )
191  nt = nt + 7
192 *
193 * DSYTRD_SY2SB
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 * DSYTRD_SB2ST
217 *
218  srnamt = 'DSYTRD_SB2ST'
219  infot = 1
220  CALL dsytrd_sb2st( '/', 'N', 'U', 0, 0, a, 1, d, e,
221  $ c, 1, w, 1, info )
222  CALL chkxer( 'DSYTRD_SB2ST', infot, nout, lerr, ok )
223  infot = 2
224  CALL dsytrd_sb2st( 'Y', '/', 'U', 0, 0, a, 1, d, e,
225  $ c, 1, w, 1, info )
226  CALL chkxer( 'DSYTRD_SB2ST', infot, nout, lerr, ok )
227  infot = 2
228  CALL dsytrd_sb2st( 'Y', 'H', 'U', 0, 0, a, 1, d, e,
229  $ c, 1, w, 1, info )
230  CALL chkxer( 'DSYTRD_SB2ST', infot, nout, lerr, ok )
231  infot = 3
232  CALL dsytrd_sb2st( 'Y', 'N', '/', 0, 0, a, 1, d, e,
233  $ c, 1, w, 1, info )
234  CALL chkxer( 'DSYTRD_SB2ST', infot, nout, lerr, ok )
235  infot = 4
236  CALL dsytrd_sb2st( 'Y', 'N', 'U', -1, 0, a, 1, d, e,
237  $ c, 1, w, 1, info )
238  CALL chkxer( 'DSYTRD_SB2ST', infot, nout, lerr, ok )
239  infot = 5
240  CALL dsytrd_sb2st( 'Y', 'N', 'U', 0, -1, a, 1, d, e,
241  $ c, 1, w, 1, info )
242  CALL chkxer( 'DSYTRD_SB2ST', infot, nout, lerr, ok )
243  infot = 7
244  CALL dsytrd_sb2st( 'Y', 'N', 'U', 0, 1, a, 1, d, e,
245  $ c, 1, w, 1, info )
246  CALL chkxer( 'DSYTRD_SB2ST', infot, nout, lerr, ok )
247  infot = 11
248  CALL dsytrd_sb2st( 'Y', 'N', 'U', 0, 0, a, 1, d, e,
249  $ c, 0, w, 1, info )
250  CALL chkxer( 'DSYTRD_SB2ST', infot, nout, lerr, ok )
251  infot = 13
252  CALL dsytrd_sb2st( 'Y', 'N', 'U', 0, 0, a, 1, d, e,
253  $ c, 1, w, 0, info )
254  CALL chkxer( 'DSYTRD_SB2ST', infot, nout, lerr, ok )
255  nt = nt + 9
256 *
257 * DORGTR
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 * DORMTR
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 * DSPTRD
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 * DOPGTR
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 * DOPMTR
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 * DPTEQR
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 * DSTEBZ
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 * DSTEIN
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 * DSTEQR
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 * DSTERF
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 * DSTEDC
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 * DSTEVD
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 * DSTEV
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 * DSTEVX
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 * DSTEVR
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 * DSYEVD
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 * DSYEVD_2STAGE
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 * INFOT = 8
666 * CALL DSYEVD_2STAGE( 'V', 'U', 2, A, 2, X, W, 20, IW, 12, INFO )
667 * CALL CHKXER( 'DSYEVD_2STAGE', INFOT, NOUT, LERR, OK )
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 * INFOT = 10
675 * CALL DSYEVD_2STAGE( 'V', 'U', 2, A, 2, X, W, 27, IW, 11, INFO )
676 * CALL CHKXER( 'DSYEVD_2STAGE', INFOT, NOUT, LERR, OK )
677  nt = nt + 9
678 *
679 * DSYEVR
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 * DSYEVR_2STAGE
735 *
736  srnamt = 'DSYEVR_2STAGE'
737  n = 1
738  infot = 1
739  CALL dsyevr_2stage( '/', 'A', 'U', 0, a, 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
744  CALL dsyevr_2stage( 'V', 'A', 'U', 0, a, 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
749  CALL dsyevr_2stage( 'N', '/', 'U', 0, a, 1,
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
754  CALL dsyevr_2stage( 'N', 'A', '/', -1, a, 1,
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
759  CALL dsyevr_2stage( 'N', 'A', 'U', -1, a, 1,
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
764  CALL dsyevr_2stage( 'N', 'A', 'U', 2, a, 1,
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
769  CALL dsyevr_2stage( 'N', 'V', 'U', 1, a, 1,
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
774  CALL dsyevr_2stage( 'N', 'I', 'U', 1, a, 1,
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
779  CALL dsyevr_2stage( 'N', 'I', 'U', 2, a, 2,
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
784  CALL dsyevr_2stage( 'N', 'I', 'U', 1, a, 1,
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
789  CALL dsyevr_2stage( 'N', 'I', 'U', 1, a, 1,
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
795  CALL dsyevr_2stage( 'N', 'I', 'U', 1, a, 1,
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 * DSYEV
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 * DSYEV_2STAGE
823 *
824  srnamt = 'DSYEV_2STAGE '
825  infot = 1
826  CALL dsyev_2stage( '/', 'U', 0, a, 1, x, w, 1, info )
827  CALL chkxer( 'DSYEV_2STAGE ', infot, nout, lerr, ok )
828  infot = 1
829  CALL dsyev_2stage( 'V', 'U', 0, a, 1, x, w, 1, info )
830  CALL chkxer( 'DSYEV_2STAGE ', infot, nout, lerr, ok )
831  infot = 2
832  CALL dsyev_2stage( 'N', '/', 0, a, 1, x, w, 1, info )
833  CALL chkxer( 'DSYEV_2STAGE ', infot, nout, lerr, ok )
834  infot = 3
835  CALL dsyev_2stage( 'N', 'U', -1, a, 1, x, w, 1, info )
836  CALL chkxer( 'DSYEV_2STAGE ', infot, nout, lerr, ok )
837  infot = 5
838  CALL dsyev_2stage( 'N', 'U', 2, a, 1, x, w, 3, info )
839  CALL chkxer( 'DSYEV_2STAGE ', infot, nout, lerr, ok )
840  infot = 8
841  CALL dsyev_2stage( 'N', 'U', 1, a, 1, x, w, 1, info )
842  CALL chkxer( 'DSYEV_2STAGE ', infot, nout, lerr, ok )
843  nt = nt + 6
844 *
845 * DSYEVX
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 * DSYEVX_2STAGE
898 *
899  srnamt = 'DSYEVX_2STAGE'
900  infot = 1
901  CALL dsyevx_2stage( '/', 'A', 'U', 0, a, 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
906  CALL dsyevx_2stage( 'V', 'A', 'U', 0, a, 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
911  CALL dsyevx_2stage( 'N', '/', 'U', 0, a, 1,
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
916  CALL dsyevx_2stage( 'N', 'A', '/', 0, a, 1,
917  $ 0.0d0, 0.0d0, 0, 0, 0.0d0,
918  $ m, x, z, 1, w, 1, iw, i3, info )
919  infot = 4
920  CALL dsyevx_2stage( 'N', 'A', 'U', -1, a, 1,
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
925  CALL dsyevx_2stage( 'N', 'A', 'U', 2, a, 1,
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
930  CALL dsyevx_2stage( 'N', 'V', 'U', 1, a, 1,
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
935  CALL dsyevx_2stage( 'N', 'I', 'U', 1, a, 1,
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
940  CALL dsyevx_2stage( 'N', 'I', 'U', 1, a, 1,
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
945  CALL dsyevx_2stage( 'N', 'I', 'U', 2, a, 2,
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
950  CALL dsyevx_2stage( 'N', 'I', 'U', 1, a, 1,
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
955  CALL dsyevx_2stage( 'N', 'A', 'U', 2, a, 2,
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
960  CALL dsyevx_2stage( 'N', 'A', 'U', 1, a, 1,
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 * DSPEVD
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 * DSPEV
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 * DSPEVX
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 * Test error exits for the SB path.
1063 *
1064  ELSE IF( lsamen( 2, c2, 'SB' ) ) THEN
1065 *
1066 * DSBTRD
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 * DSYTRD_SB2ST
1090 *
1091  srnamt = 'DSYTRD_SB2ST'
1092  infot = 1
1093  CALL dsytrd_sb2st( '/', 'N', 'U', 0, 0, a, 1, d, e,
1094  $ c, 1, w, 1, info )
1095  CALL chkxer( 'DSYTRD_SB2ST', infot, nout, lerr, ok )
1096  infot = 2
1097  CALL dsytrd_sb2st( 'N', '/', 'U', 0, 0, a, 1, d, e,
1098  $ c, 1, w, 1, info )
1099  CALL chkxer( 'DSYTRD_SB2ST', infot, nout, lerr, ok )
1100  infot = 2
1101  CALL dsytrd_sb2st( 'N', 'H', 'U', 0, 0, a, 1, d, e,
1102  $ c, 1, w, 1, info )
1103  CALL chkxer( 'DSYTRD_SB2ST', infot, nout, lerr, ok )
1104  infot = 3
1105  CALL dsytrd_sb2st( 'N', 'N', '/', 0, 0, a, 1, d, e,
1106  $ c, 1, w, 1, info )
1107  CALL chkxer( 'DSYTRD_SB2ST', infot, nout, lerr, ok )
1108  infot = 4
1109  CALL dsytrd_sb2st( 'N', 'N', 'U', -1, 0, a, 1, d, e,
1110  $ c, 1, w, 1, info )
1111  CALL chkxer( 'DSYTRD_SB2ST', infot, nout, lerr, ok )
1112  infot = 5
1113  CALL dsytrd_sb2st( 'N', 'N', 'U', 0, -1, a, 1, d, e,
1114  $ c, 1, w, 1, info )
1115  CALL chkxer( 'DSYTRD_SB2ST', infot, nout, lerr, ok )
1116  infot = 7
1117  CALL dsytrd_sb2st( 'N', 'N', 'U', 0, 1, a, 1, d, e,
1118  $ c, 1, w, 1, info )
1119  CALL chkxer( 'DSYTRD_SB2ST', infot, nout, lerr, ok )
1120  infot = 11
1121  CALL dsytrd_sb2st( 'N', 'N', 'U', 0, 0, a, 1, d, e,
1122  $ c, 0, w, 1, info )
1123  CALL chkxer( 'DSYTRD_SB2ST', infot, nout, lerr, ok )
1124  infot = 13
1125  CALL dsytrd_sb2st( 'N', 'N', 'U', 0, 0, a, 1, d, e,
1126  $ c, 1, w, 0, info )
1127  CALL chkxer( 'DSYTRD_SB2ST', infot, nout, lerr, ok )
1128  nt = nt + 9
1129 *
1130 * DSBEVD
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 * DSBEVD_2STAGE
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 * INFOT = 9
1201 * CALL DSBEVD_2STAGE( 'V', 'U', 2, 1, A, 2, X, Z, 1, W,
1202 * $ 25, IW, 12, INFO )
1203 * CALL CHKXER( 'DSBEVD_2STAGE', INFOT, NOUT, LERR, OK )
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 * INFOT = 11
1213 * CALL DSBEVD_2STAGE( 'V', 'U', 2, 0, A, 1, X, Z, 2, W,
1214 * $ 18, IW, 12, INFO )
1215 * CALL CHKXER( 'DSBEVD_2STAGE', INFOT, NOUT, LERR, OK )
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 * INFOT = 13
1221 * CALL DSBEVD_2STAGE( 'V', 'U', 2, 0, A, 1, X, Z, 2, W,
1222 * $ 25, IW, 11, INFO )
1223 * CALL CHKXER( 'DSBEVD_2STAGE', INFOT, NOUT, LERR, OK )
1224 * NT = NT + 12
1225  nt = nt + 9
1226 *
1227 * DSBEV
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 * DSBEV_2STAGE
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 * DSBEVX
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 * DSBEVX_2STAGE
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 * INFOT = 9
1368 * CALL DSBEVX_2STAGE( 'V', 'A', 'U', 2, 0, A, 1, Q, 1, 0.0D0,
1369 * $ 0.0D0, 0, 0, 0.0D0, M, X, Z, 2, W, 0, IW, I3, INFO )
1370 * CALL CHKXER( 'DSBEVX_2STAGE', INFOT, NOUT, LERR, OK )
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 * INFOT = 18
1392 * CALL DSBEVX_2STAGE( 'V', 'A', 'U', 2, 0, A, 1, Q, 2, 0.0D0,
1393 * $ 0.0D0, 0, 0, 0.0D0, M, X, Z, 1, W, 0, IW, I3, INFO )
1394 * CALL CHKXER( 'DSBEVX_2STAGE', INFOT, NOUT, LERR, OK )
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 * NT = NT + 15
1400  nt = nt + 13
1401  END IF
1402 *
1403 * Print a summary line.
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 * End of DERRST
1419 *
chkxer
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
Definition: cblat2.f:3196
dsytrd_sb2st
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
Definition: dsytrd_sb2st.F:230
lsamen
logical function lsamen(N, CA, CB)
LSAMEN
Definition: lsamen.f:74
dstebz
subroutine dstebz(RANGE, ORDER, N, VL, VU, IL, IU, ABSTOL, D, E, M, NSPLIT, W, IBLOCK, ISPLIT, WORK, IWORK, INFO)
DSTEBZ
Definition: dstebz.f:273
dsteqr
subroutine dsteqr(COMPZ, N, D, E, Z, LDZ, WORK, INFO)
DSTEQR
Definition: dsteqr.f:131
dstedc
subroutine dstedc(COMPZ, N, D, E, Z, LDZ, WORK, LWORK, IWORK, LIWORK, INFO)
DSTEDC
Definition: dstedc.f:188
dsterf
subroutine dsterf(N, D, E, INFO)
DSTERF
Definition: dsterf.f:86
dsptrd
subroutine dsptrd(UPLO, N, AP, D, E, TAU, INFO)
DSPTRD
Definition: dsptrd.f:150
dstein
subroutine dstein(N, D, E, M, W, IBLOCK, ISPLIT, Z, LDZ, WORK, IWORK, IFAIL, INFO)
DSTEIN
Definition: dstein.f:174
dopgtr
subroutine dopgtr(UPLO, N, AP, TAU, Q, LDQ, WORK, INFO)
DOPGTR
Definition: dopgtr.f:114
dsbtrd
subroutine dsbtrd(VECT, UPLO, N, KD, AB, LDAB, D, E, Q, LDQ, WORK, INFO)
DSBTRD
Definition: dsbtrd.f:163
dormtr
subroutine dormtr(SIDE, UPLO, TRANS, M, N, A, LDA, TAU, C, LDC, WORK, LWORK, INFO)
DORMTR
Definition: dormtr.f:171
dorgtr
subroutine dorgtr(UPLO, N, A, LDA, TAU, WORK, LWORK, INFO)
DORGTR
Definition: dorgtr.f:123
dopmtr
subroutine dopmtr(SIDE, UPLO, TRANS, M, N, AP, TAU, C, LDC, WORK, INFO)
DOPMTR
Definition: dopmtr.f:150
dspevx
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...
Definition: dspevx.f:234
dstevd
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...
Definition: dstevd.f:163
dstevr
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...
Definition: dstevr.f:304
dsbevx
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...
Definition: dsbevx.f:265
dspevd
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...
Definition: dspevd.f:178
dstevx
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...
Definition: dstevx.f:227
dsbevd_2stage
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 ...
Definition: dsbevd_2stage.f:234
dspev
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
Definition: dspev.f:130
dsbev_2stage
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...
Definition: dsbev_2stage.f:204
dstev
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
Definition: dstev.f:116
dsbevx_2stage
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 ...
Definition: dsbevx_2stage.f:322
dsbev
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
Definition: dsbev.f:146
dsbevd
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...
Definition: dsbevd.f:193
dpteqr
subroutine dpteqr(COMPZ, N, D, E, Z, LDZ, WORK, INFO)
DPTEQR
Definition: dpteqr.f:145
dsytrd_sy2sb
subroutine dsytrd_sy2sb(UPLO, N, KD, A, LDA, AB, LDAB, TAU, WORK, LWORK, INFO)
DSYTRD_SY2SB
Definition: dsytrd_sy2sb.f:243
dsytrd_2stage
subroutine dsytrd_2stage(VECT, UPLO, N, A, LDA, D, E, TAU, HOUS2, LHOUS2, WORK, LWORK, INFO)
DSYTRD_2STAGE
Definition: dsytrd_2stage.f:224
dsytrd
subroutine dsytrd(UPLO, N, A, LDA, D, E, TAU, WORK, LWORK, INFO)
DSYTRD
Definition: dsytrd.f:192
dsyev
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
Definition: dsyev.f:132
dsyevx
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
Definition: dsyevx.f:253
dsyevd
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
Definition: dsyevd.f:185
dsyevr_2stage
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...
Definition: dsyevr_2stage.f:381
dsyevx_2stage
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...
Definition: dsyevx_2stage.f:300
dsyev_2stage
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...
Definition: dsyev_2stage.f:183
dsyevr
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
Definition: dsyevr.f:334
dsyevd_2stage
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...
Definition: dsyevd_2stage.f:227
Here is the call graph for this function:
Here is the caller graph for this function: