OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
derrst.f
Go to the documentation of this file.
1*> \brief \b DERRST
2*
3* =========== DOCUMENTATION ===========
4*
5* Online html documentation available at
6* http://www.netlib.org/lapack/explore-html/
7*
8* Definition:
9* ===========
10*
11* SUBROUTINE DERRST( PATH, NUNIT )
12*
13* .. Scalar Arguments ..
14* CHARACTER*3 PATH
15* INTEGER NUNIT
16* ..
17*
18*
19*> \par Purpose:
20* =============
21*>
22*> \verbatim
23*>
24*> DERRST tests the error exits for DSYTRD, DORGTR, DORMTR, DSPTRD,
25*> DOPGTR, DOPMTR, DSTEQR, SSTERF, SSTEBZ, SSTEIN, DPTEQR, DSBTRD,
26*> DSYEV, SSYEVX, SSYEVD, DSBEV, SSBEVX, SSBEVD,
27*> DSPEV, SSPEVX, SSPEVD, DSTEV, SSTEVX, SSTEVD, and SSTEDC.
28*> DSYEVD_2STAGE, DSYEVR_2STAGE, DSYEVX_2STAGE,
29*> DSYEV_2STAGE, DSBEV_2STAGE, DSBEVD_2STAGE,
30*> DSBEVX_2STAGE, DSYTRD_2STAGE, DSYTRD_SY2SB,
31*> DSYTRD_SB2ST
32*> \endverbatim
33*
34* Arguments:
35* ==========
36*
37*> \param[in] PATH
38*> \verbatim
39*> PATH is CHARACTER*3
40*> The LAPACK path name for the routines to be tested.
41*> \endverbatim
42*>
43*> \param[in] NUNIT
44*> \verbatim
45*> NUNIT is INTEGER
46*> The unit number for output.
47*> \endverbatim
48*
49* Authors:
50* ========
51*
52*> \author Univ. of Tennessee
53*> \author Univ. of California Berkeley
54*> \author Univ. of Colorado Denver
55*> \author NAG Ltd.
56*
57*> \ingroup double_eig
58*
59* =====================================================================
60 SUBROUTINE derrst( PATH, NUNIT )
61*
62* -- LAPACK test routine --
63* -- LAPACK is a software package provided by Univ. of Tennessee, --
64* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
65*
66* .. Scalar Arguments ..
67 CHARACTER*3 PATH
68 INTEGER NUNIT
69* ..
70*
71* =====================================================================
72*
73* NMAX has to be at least 3 or LIW may be too small
74* .. Parameters ..
75 INTEGER NMAX, LIW, LW
76 parameter( nmax = 3, liw = 12*nmax, lw = 20*nmax )
77* ..
78* .. Local Scalars ..
79 CHARACTER*2 C2
80 INTEGER I, INFO, J, M, N, NSPLIT, NT
81* ..
82* .. Local Arrays ..
83 INTEGER I1( NMAX ), I2( NMAX ), I3( NMAX ), IW( LIW )
84 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*
1420 END
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
dsteqr
subroutine dsteqr(compz, n, d, e, z, ldz, work, info)
DSTEQR
Definition dsteqr.f:131
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
dsterf
subroutine dsterf(n, d, e, info)
DSTERF
Definition dsterf.f:86
dstedc
subroutine dstedc(compz, n, d, e, z, ldz, work, lwork, iwork, liwork, info)
DSTEDC
Definition dstedc.f:188
dormtr
subroutine dormtr(side, uplo, trans, m, n, a, lda, tau, c, ldc, work, lwork, info)
DORMTR
Definition dormtr.f:171
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
dopmtr
subroutine dopmtr(side, uplo, trans, m, n, ap, tau, c, ldc, work, info)
DOPMTR
Definition dopmtr.f:150
dsbtrd
subroutine dsbtrd(vect, uplo, n, kd, ab, ldab, d, e, q, ldq, work, info)
DSBTRD
Definition dsbtrd.f:163
dopgtr
subroutine dopgtr(uplo, n, ap, tau, q, ldq, work, info)
DOPGTR
Definition dopgtr.f:114
dorgtr
subroutine dorgtr(uplo, n, a, lda, tau, work, lwork, info)
DORGTR
Definition dorgtr.f:123
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
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
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
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
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
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
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
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
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_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
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
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
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
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
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
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
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
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
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
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
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
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
derrst
subroutine derrst(path, nunit)
DERRST
Definition derrst.f:61