OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
cerrvx.f
Go to the documentation of this file.
1*> \brief \b CERRVX
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 CERRVX( 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*> CERRVX tests the error exits for the COMPLEX driver routines
25*> for solving linear systems of equations.
26*> \endverbatim
27*
28* Arguments:
29* ==========
30*
31*> \param[in] PATH
32*> \verbatim
33*> PATH is CHARACTER*3
34*> The LAPACK path name for the routines to be tested.
35*> \endverbatim
36*>
37*> \param[in] NUNIT
38*> \verbatim
39*> NUNIT is INTEGER
40*> The unit number for output.
41*> \endverbatim
42*
43* Authors:
44* ========
45*
46*> \author Univ. of Tennessee
47*> \author Univ. of California Berkeley
48*> \author Univ. of Colorado Denver
49*> \author NAG Ltd.
50*
51*> \ingroup complex_lin
52*
53* =====================================================================
54 SUBROUTINE cerrvx( PATH, NUNIT )
55*
56* -- LAPACK test routine --
57* -- LAPACK is a software package provided by Univ. of Tennessee, --
58* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
59*
60* .. Scalar Arguments ..
61 CHARACTER*3 PATH
62 INTEGER NUNIT
63* ..
64*
65* =====================================================================
66*
67* .. Parameters ..
68 INTEGER NMAX
69 parameter( nmax = 4 )
70* ..
71* .. Local Scalars ..
72 CHARACTER EQ
73 CHARACTER*2 C2
74 INTEGER I, INFO, J
75 REAL RCOND
76* ..
77* .. Local Arrays ..
78 INTEGER IP( NMAX )
79 REAL C( NMAX ), R( NMAX ), R1( NMAX ), R2( NMAX ),
80 $ RF( NMAX ), RW( NMAX )
81 COMPLEX A( NMAX, NMAX ), AF( NMAX, NMAX ), B( NMAX ),
82 $ E( NMAX ), W( 2*NMAX ), X( NMAX )
83* ..
84* .. External Functions ..
85 LOGICAL LSAMEN
86 EXTERNAL lsamen
87* ..
88* .. External Subroutines ..
89 EXTERNAL cgbsv, cgbsvx, cgesv, cgesvx, cgtsv, cgtsvx,
90 $ chesv, chesv_rk ,chesv_rook, chesvx, chkxer,
91 $ chpsv, chpsvx, cpbsv, cpbsvx, cposv, cposvx,
92 $ cppsv, cppsvx, cptsv, cptsvx, cspsv, cspsvx,
93 $ csysv, csysv_aa, csysv_rk, csysv_rook,
94 $ csysvx, csysv_aa_2stage
95* ..
96* .. Scalars in Common ..
97 LOGICAL LERR, OK
98 CHARACTER*32 SRNAMT
99 INTEGER INFOT, NOUT
100* ..
101* .. Common blocks ..
102 COMMON / infoc / infot, nout, ok, lerr
103 COMMON / srnamc / srnamt
104* ..
105* .. Intrinsic Functions ..
106 INTRINSIC cmplx, real
107* ..
108* .. Executable Statements ..
109*
110 nout = nunit
111 WRITE( nout, fmt = * )
112 c2 = path( 2: 3 )
113*
114* Set the variables to innocuous values.
115*
116 DO 20 j = 1, nmax
117 DO 10 i = 1, nmax
118 a( i, j ) = cmplx( 1. / real( i+j ), -1. / real( i+j ) )
119 af( i, j ) = cmplx( 1. / real( i+j ), -1. / real( i+j ) )
120 10 CONTINUE
121 b( j ) = 0.e+0
122 e( j ) = 0.e+0
123 r1( j ) = 0.e+0
124 r2( j ) = 0.e+0
125 w( j ) = 0.e+0
126 x( j ) = 0.e+0
127 c( j ) = 0.e+0
128 r( j ) = 0.e+0
129 ip( j ) = j
130 20 CONTINUE
131 eq = ' '
132 ok = .true.
133*
134 IF( lsamen( 2, c2, 'GE' ) ) THEN
135*
136* CGESV
137*
138 srnamt = 'CGESV '
139 infot = 1
140 CALL cgesv( -1, 0, a, 1, ip, b, 1, info )
141 CALL chkxer( 'CGESV ', infot, nout, lerr, ok )
142 infot = 2
143 CALL cgesv( 0, -1, a, 1, ip, b, 1, info )
144 CALL chkxer( 'CGESV ', infot, nout, lerr, ok )
145 infot = 4
146 CALL cgesv( 2, 1, a, 1, ip, b, 2, info )
147 CALL chkxer( 'CGESV ', infot, nout, lerr, ok )
148 infot = 7
149 CALL cgesv( 2, 1, a, 2, ip, b, 1, info )
150 CALL chkxer( 'CGESV ', infot, nout, lerr, ok )
151*
152* CGESVX
153*
154 srnamt = 'CGESVX'
155 infot = 1
156 CALL cgesvx( '/', 'N', 0, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
157 $ x, 1, rcond, r1, r2, w, rw, info )
158 CALL chkxer( 'CGESVX', infot, nout, lerr, ok )
159 infot = 2
160 CALL cgesvx( 'n', '/', 0, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
161 $ X, 1, RCOND, R1, R2, W, RW, INFO )
162 CALL CHKXER( 'cgesvx', INFOT, NOUT, LERR, OK )
163 INFOT = 3
164 CALL CGESVX( 'n', 'n', -1, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
165 $ X, 1, RCOND, R1, R2, W, RW, INFO )
166 CALL CHKXER( 'cgesvx', INFOT, NOUT, LERR, OK )
167 INFOT = 4
168 CALL CGESVX( 'n', 'n', 0, -1, A, 1, AF, 1, IP, EQ, R, C, B, 1,
169 $ X, 1, RCOND, R1, R2, W, RW, INFO )
170 CALL CHKXER( 'cgesvx', INFOT, NOUT, LERR, OK )
171 INFOT = 6
172 CALL CGESVX( 'n', 'n', 2, 1, A, 1, AF, 2, IP, EQ, R, C, B, 2,
173 $ X, 2, RCOND, R1, R2, W, RW, INFO )
174 CALL CHKXER( 'cgesvx', INFOT, NOUT, LERR, OK )
175 INFOT = 8
176 CALL CGESVX( 'n', 'n', 2, 1, A, 2, AF, 1, IP, EQ, R, C, B, 2,
177 $ X, 2, RCOND, R1, R2, W, RW, INFO )
178 CALL CHKXER( 'cgesvx', INFOT, NOUT, LERR, OK )
179 INFOT = 10
180 EQ = '/'
181 CALL CGESVX( 'f', 'n', 0, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
182 $ X, 1, RCOND, R1, R2, W, RW, INFO )
183 CALL CHKXER( 'cgesvx', INFOT, NOUT, LERR, OK )
184 INFOT = 11
185 EQ = 'r'
186 CALL CGESVX( 'f', 'n', 1, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
187 $ X, 1, RCOND, R1, R2, W, RW, INFO )
188 CALL CHKXER( 'cgesvx', infot, nout, lerr, ok )
189 infot = 12
190 eq = 'C'
191 CALL cgesvx( 'F', 'N', 1, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
192 $ x, 1, rcond, r1, r2, w, rw, info )
193 CALL chkxer( 'CGESVX', infot, nout, lerr, ok )
194 infot = 14
195 CALL cgesvx( 'N', 'N', 2, 1, a, 2, af, 2, ip, eq, r, c, b, 1,
196 $ x, 2, rcond, r1, r2, w, rw, info )
197 CALL chkxer( 'CGESVX', infot, nout, lerr, ok )
198 infot = 16
199 CALL cgesvx( 'N', 'N', 2, 1, a, 2, af, 2, ip, eq, r, c, b, 2,
200 $ x, 1, rcond, r1, r2, w, rw, info )
201 CALL chkxer( 'CGESVX', infot, nout, lerr, ok )
202*
203 ELSE IF( lsamen( 2, c2, 'GB' ) ) THEN
204*
205* CGBSV
206*
207 srnamt = 'CGBSV '
208 infot = 1
209 CALL cgbsv( -1, 0, 0, 0, a, 1, ip, b, 1, info )
210 CALL chkxer( 'CGBSV ', infot, nout, lerr, ok )
211 infot = 2
212 CALL cgbsv( 1, -1, 0, 0, a, 1, ip, b, 1, info )
213 CALL chkxer( 'CGBSV ', infot, nout, lerr, ok )
214 infot = 3
215 CALL cgbsv( 1, 0, -1, 0, a, 1, ip, b, 1, info )
216 CALL chkxer( 'CGBSV ', infot, nout, lerr, ok )
217 infot = 4
218 CALL cgbsv( 0, 0, 0, -1, a, 1, ip, b, 1, info )
219 CALL chkxer( 'CGBSV ', infot, nout, lerr, ok )
220 infot = 6
221 CALL cgbsv( 1, 1, 1, 0, a, 3, ip, b, 1, info )
222 CALL chkxer( 'CGBSV ', infot, nout, lerr, ok )
223 infot = 9
224 CALL cgbsv( 2, 0, 0, 0, a, 1, ip, b, 1, info )
225 CALL chkxer( 'CGBSV ', infot, nout, lerr, ok )
226*
227* CGBSVX
228*
229 srnamt = 'CGBSVX'
230 infot = 1
231 CALL cgbsvx( '/', 'N', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
232 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
233 CALL chkxer( 'CGBSVX', infot, nout, lerr, ok )
234 infot = 2
235 CALL cgbsvx( 'N', '/', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
236 $ b, 1, x, 1, rcond, r1, r2, w, rw, info )
237 CALL chkxer( 'CGBSVX', infot, nout, lerr, ok )
238 infot = 3
239 CALL cgbsvx( 'N', 'n', -1, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
240 $ B, 1, X, 1, RCOND, R1, R2, W, RW, INFO )
241 CALL CHKXER( 'cgbsvx', INFOT, NOUT, LERR, OK )
242 INFOT = 4
243 CALL CGBSVX( 'n', 'n', 1, -1, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
244 $ B, 1, X, 1, RCOND, R1, R2, W, RW, INFO )
245 CALL CHKXER( 'cgbsvx', INFOT, NOUT, LERR, OK )
246 INFOT = 5
247 CALL CGBSVX( 'n', 'n', 1, 0, -1, 0, A, 1, AF, 1, IP, EQ, R, C,
248 $ B, 1, X, 1, RCOND, R1, R2, W, RW, INFO )
249 CALL CHKXER( 'cgbsvx', INFOT, NOUT, LERR, OK )
250 INFOT = 6
251 CALL CGBSVX( 'n', 'n', 0, 0, 0, -1, A, 1, AF, 1, IP, EQ, R, C,
252 $ B, 1, X, 1, RCOND, R1, R2, W, RW, INFO )
253 CALL CHKXER( 'cgbsvx', INFOT, NOUT, LERR, OK )
254 INFOT = 8
255 CALL CGBSVX( 'n', 'n', 1, 1, 1, 0, A, 2, AF, 4, IP, EQ, R, C,
256 $ B, 1, X, 1, RCOND, R1, R2, W, RW, INFO )
257 CALL CHKXER( 'cgbsvx', INFOT, NOUT, LERR, OK )
258 INFOT = 10
259 CALL CGBSVX( 'n', 'n', 1, 1, 1, 0, A, 3, AF, 3, IP, EQ, R, C,
260 $ B, 1, X, 1, RCOND, R1, R2, W, RW, INFO )
261 CALL CHKXER( 'cgbsvx', INFOT, NOUT, LERR, OK )
262 INFOT = 12
263 EQ = '/'
264 CALL CGBSVX( 'f', 'n', 0, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
265 $ B, 1, X, 1, RCOND, R1, R2, W, RW, INFO )
266 CALL CHKXER( 'cgbsvx', INFOT, NOUT, LERR, OK )
267 INFOT = 13
268 EQ = 'r'
269 CALL CGBSVX( 'f', 'n', 1, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
270 $ B, 1, X, 1, RCOND, R1, R2, W, RW, INFO )
271 CALL CHKXER( 'cgbsvx', INFOT, NOUT, LERR, OK )
272 INFOT = 14
273 EQ = 'c'
274 CALL CGBSVX( 'f', 'n', 1, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
275 $ B, 1, X, 1, RCOND, R1, R2, W, RW, INFO )
276 CALL CHKXER( 'cgbsvx', INFOT, NOUT, LERR, OK )
277 INFOT = 16
278 CALL CGBSVX( 'n', 'n', 2, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
279 $ B, 1, X, 2, RCOND, R1, R2, W, RW, INFO )
280 CALL CHKXER( 'cgbsvx', INFOT, NOUT, LERR, OK )
281 INFOT = 18
282 CALL CGBSVX( 'n', 'n', 2, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
283 $ B, 2, X, 1, RCOND, R1, R2, W, RW, INFO )
284 CALL CHKXER( 'cgbsvx', INFOT, NOUT, LERR, OK )
285*
286 ELSE IF( LSAMEN( 2, C2, 'gt' ) ) THEN
287*
288* CGTSV
289*
290 SRNAMT = 'cgtsv '
291 INFOT = 1
292 CALL CGTSV( -1, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), B, 1,
293 $ INFO )
294 CALL CHKXER( 'cgtsv ', INFOT, NOUT, LERR, OK )
295 INFOT = 2
296 CALL CGTSV( 0, -1, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), B, 1,
297 $ INFO )
298 CALL CHKXER( 'cgtsv ', infot, nout, lerr, ok )
299 infot = 7
300 CALL cgtsv( 2, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ), b, 1, info )
301 CALL chkxer( 'CGTSV ', infot, nout, lerr, ok )
302*
303* CGTSVX
304*
305 srnamt = 'CGTSVX'
306 infot = 1
307 CALL cgtsvx( '/', 'N', 0, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
308 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
309 $ ip, b, 1, x, 1, rcond, r1, r2, w, rw, info )
310 CALL chkxer( 'CGTSVX', infot, nout, lerr, ok )
311 infot = 2
312 CALL cgtsvx( 'N', '/', 0, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
313 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
314 $ ip, b, 1, x, 1, rcond, r1, r2, w, rw, info )
315 CALL chkxer( 'CGTSVX', infot, nout, lerr, ok )
316 infot = 3
317 CALL cgtsvx( 'N', 'N', -1, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
318 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
319 $ ip, b, 1, x, 1, rcond, r1, r2, w, rw, info )
320 CALL chkxer( 'CGTSVX', infot, nout, lerr, ok )
321 infot = 4
322 CALL cgtsvx( 'N', 'N', 0, -1, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
323 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
324 $ ip, b, 1, x, 1, rcond, r1, r2, w, rw, info )
325 CALL chkxer( 'CGTSVX', infot, nout, lerr, ok )
326 infot = 14
327 CALL cgtsvx( 'N', 'N', 2, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
328 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
329 $ ip, b, 1, x, 2, rcond, r1, r2, w, rw, info )
330 CALL chkxer( 'cgtsvx', INFOT, NOUT, LERR, OK )
331 INFOT = 16
332 CALL CGTSVX( 'n', 'n', 2, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
333 $ AF( 1, 1 ), AF( 1, 2 ), AF( 1, 3 ), AF( 1, 4 ),
334 $ IP, B, 2, X, 1, RCOND, R1, R2, W, RW, INFO )
335 CALL CHKXER( 'cgtsvx', INFOT, NOUT, LERR, OK )
336*
337 ELSE IF( LSAMEN( 2, C2, 'po' ) ) THEN
338*
339* CPOSV
340*
341 SRNAMT = 'cposv '
342 INFOT = 1
343 CALL CPOSV( '/', 0, 0, A, 1, B, 1, INFO )
344 CALL CHKXER( 'cposv ', INFOT, NOUT, LERR, OK )
345 INFOT = 2
346 CALL CPOSV( 'u', -1, 0, A, 1, B, 1, INFO )
347 CALL CHKXER( 'cposv ', INFOT, NOUT, LERR, OK )
348 INFOT = 3
349 CALL CPOSV( 'u', 0, -1, A, 1, B, 1, INFO )
350 CALL CHKXER( 'cposv ', INFOT, NOUT, LERR, OK )
351 INFOT = 5
352 CALL CPOSV( 'u', 2, 0, A, 1, B, 2, INFO )
353 CALL CHKXER( 'cposv ', INFOT, NOUT, LERR, OK )
354 INFOT = 7
355 CALL CPOSV( 'u', 2, 0, A, 2, B, 1, INFO )
356 CALL CHKXER( 'cposv ', INFOT, NOUT, LERR, OK )
357*
358* CPOSVX
359*
360 SRNAMT = 'cposvx'
361 INFOT = 1
362 CALL CPOSVX( '/', 'u', 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
363 $ RCOND, R1, R2, W, RW, INFO )
364 CALL CHKXER( 'cposvx', INFOT, NOUT, LERR, OK )
365 INFOT = 2
366 CALL CPOSVX( 'n', '/', 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
367 $ RCOND, R1, R2, W, RW, INFO )
368 CALL CHKXER( 'cposvx', INFOT, NOUT, LERR, OK )
369 INFOT = 3
370 CALL CPOSVX( 'n', 'u', -1, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
371 $ RCOND, R1, R2, W, RW, INFO )
372 CALL CHKXER( 'cposvx', INFOT, NOUT, LERR, OK )
373 INFOT = 4
374 CALL CPOSVX( 'n', 'u', 0, -1, A, 1, AF, 1, EQ, C, B, 1, X, 1,
375 $ RCOND, R1, R2, W, RW, INFO )
376 CALL CHKXER( 'cposvx', INFOT, NOUT, LERR, OK )
377 INFOT = 6
378 CALL CPOSVX( 'n', 'u', 2, 0, A, 1, AF, 2, EQ, C, B, 2, X, 2,
379 $ RCOND, R1, R2, W, RW, INFO )
380 CALL CHKXER( 'cposvx', INFOT, NOUT, LERR, OK )
381 INFOT = 8
382 CALL CPOSVX( 'n', 'u', 2, 0, A, 2, AF, 1, EQ, C, B, 2, X, 2,
383 $ RCOND, R1, R2, W, RW, INFO )
384 CALL CHKXER( 'cposvx', INFOT, NOUT, LERR, OK )
385 INFOT = 9
386 EQ = '/'
387 CALL CPOSVX( 'f', 'u', 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
388 $ RCOND, R1, R2, W, RW, INFO )
389 CALL CHKXER( 'cposvx', INFOT, NOUT, LERR, OK )
390 INFOT = 10
391 EQ = 'y'
392 CALL CPOSVX( 'f', 'u', 1, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
393 $ RCOND, R1, R2, W, RW, INFO )
394 CALL CHKXER( 'cposvx', INFOT, NOUT, LERR, OK )
395 INFOT = 12
396 CALL CPOSVX( 'n', 'u', 2, 0, A, 2, AF, 2, EQ, C, B, 1, X, 2,
397 $ RCOND, R1, R2, W, RW, INFO )
398 CALL CHKXER( 'cposvx', INFOT, NOUT, LERR, OK )
399 INFOT = 14
400 CALL CPOSVX( 'n', 'u', 2, 0, A, 2, AF, 2, EQ, C, B, 2, X, 1,
401 $ RCOND, R1, R2, W, RW, INFO )
402 CALL CHKXER( 'cposvx', INFOT, NOUT, LERR, OK )
403*
404 ELSE IF( LSAMEN( 2, C2, 'pp' ) ) THEN
405*
406* CPPSV
407*
408 SRNAMT = 'cppsv '
409 INFOT = 1
410 CALL CPPSV( '/', 0, 0, A, B, 1, INFO )
411 CALL CHKXER( 'cppsv ', INFOT, NOUT, LERR, OK )
412 INFOT = 2
413 CALL CPPSV( 'u', -1, 0, A, B, 1, INFO )
414 CALL CHKXER( 'cppsv ', INFOT, NOUT, LERR, OK )
415 INFOT = 3
416 CALL CPPSV( 'u', 0, -1, A, B, 1, INFO )
417 CALL CHKXER( 'cppsv ', INFOT, NOUT, LERR, OK )
418 INFOT = 6
419 CALL CPPSV( 'u', 2, 0, a, b, 1, info )
420 CALL chkxer( 'CPPSV ', infot, nout, lerr, ok )
421*
422* CPPSVX
423*
424 srnamt = 'CPPSVX'
425 infot = 1
426 CALL cppsvx( '/', 'U', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
427 $ r1, r2, w, rw, info )
428 CALL chkxer( 'CPPSVX', infot, nout, lerr, ok )
429 infot = 2
430 CALL cppsvx( 'N', '/', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
431 $ r1, r2, w, rw, info )
432 CALL chkxer( 'CPPSVX', infot, nout, lerr, ok )
433 infot = 3
434 CALL cppsvx( 'N', 'U', -1, 0, a, af, eq, c, b, 1, x, 1, rcond,
435 $ r1, r2, w, rw, info )
436 CALL chkxer( 'cppsvx', INFOT, NOUT, LERR, OK )
437 INFOT = 4
438 CALL CPPSVX( 'n', 'u', 0, -1, A, AF, EQ, C, B, 1, X, 1, RCOND,
439 $ R1, R2, W, RW, INFO )
440 CALL CHKXER( 'cppsvx', INFOT, NOUT, LERR, OK )
441 INFOT = 7
442 EQ = '/'
443 CALL CPPSVX( 'f', 'u', 0, 0, A, AF, EQ, C, B, 1, X, 1, RCOND,
444 $ R1, R2, W, RW, INFO )
445 CALL CHKXER( 'cppsvx', INFOT, NOUT, LERR, OK )
446 INFOT = 8
447 EQ = 'y'
448 CALL CPPSVX( 'f', 'u', 1, 0, A, AF, EQ, C, B, 1, X, 1, RCOND,
449 $ R1, R2, W, RW, INFO )
450 CALL CHKXER( 'cppsvx', INFOT, NOUT, LERR, OK )
451 INFOT = 10
452 CALL CPPSVX( 'n', 'u', 2, 0, A, AF, EQ, C, B, 1, X, 2, RCOND,
453 $ R1, R2, W, RW, INFO )
454 CALL CHKXER( 'cppsvx', INFOT, NOUT, LERR, OK )
455 INFOT = 12
456 CALL CPPSVX( 'n', 'u', 2, 0, A, AF, EQ, C, B, 2, X, 1, RCOND,
457 $ R1, R2, W, RW, INFO )
458 CALL CHKXER( 'cppsvx', INFOT, NOUT, LERR, OK )
459*
460 ELSE IF( LSAMEN( 2, C2, 'pb' ) ) THEN
461*
462* CPBSV
463*
464 SRNAMT = 'cpbsv '
465 INFOT = 1
466 CALL CPBSV( '/', 0, 0, 0, A, 1, B, 1, INFO )
467 CALL CHKXER( 'cpbsv ', INFOT, NOUT, LERR, OK )
468 INFOT = 2
469 CALL CPBSV( 'u', -1, 0, 0, A, 1, B, 1, INFO )
470 CALL CHKXER( 'cpbsv ', INFOT, NOUT, LERR, OK )
471 INFOT = 3
472 CALL CPBSV( 'u', 1, -1, 0, A, 1, B, 1, INFO )
473 CALL CHKXER( 'cpbsv ', INFOT, NOUT, LERR, OK )
474 INFOT = 4
475 CALL CPBSV( 'u', 0, 0, -1, A, 1, B, 1, INFO )
476 CALL CHKXER( 'cpbsv ', INFOT, NOUT, LERR, OK )
477 INFOT = 6
478 CALL CPBSV( 'u', 1, 1, 0, A, 1, B, 2, INFO )
479 CALL CHKXER( 'cpbsv ', INFOT, NOUT, LERR, OK )
480 INFOT = 8
481 CALL CPBSV( 'u', 2, 0, 0, A, 1, B, 1, INFO )
482 CALL CHKXER( 'cpbsv ', INFOT, NOUT, LERR, OK )
483*
484* CPBSVX
485*
486 SRNAMT = 'cpbsvx'
487 INFOT = 1
488 CALL CPBSVX( '/', 'u', 0, 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
489 $ RCOND, R1, R2, W, RW, INFO )
490 CALL CHKXER( 'cpbsvx', INFOT, NOUT, LERR, OK )
491 INFOT = 2
492 CALL CPBSVX( 'n', '/', 0, 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
493 $ RCOND, R1, R2, W, RW, INFO )
494 CALL CHKXER( 'cpbsvx', INFOT, NOUT, LERR, OK )
495 INFOT = 3
496 CALL CPBSVX( 'n', 'u', -1, 0, 0, A, 1, AF, 1, EQ, C, B, 1, X,
497 $ 1, RCOND, R1, R2, W, RW, INFO )
498 CALL CHKXER( 'cpbsvx', INFOT, NOUT, LERR, OK )
499 INFOT = 4
500 CALL CPBSVX( 'n', 'u', 1, -1, 0, A, 1, AF, 1, EQ, C, B, 1, X,
501 $ 1, RCOND, R1, R2, W, RW, INFO )
502 CALL CHKXER( 'cpbsvx', INFOT, NOUT, LERR, OK )
503 INFOT = 5
504 CALL CPBSVX( 'n', 'u', 0, 0, -1, A, 1, AF, 1, EQ, C, B, 1, X,
505 $ 1, RCOND, R1, R2, W, RW, INFO )
506 CALL CHKXER( 'cpbsvx', INFOT, NOUT, LERR, OK )
507 INFOT = 7
508 CALL CPBSVX( 'n', 'u', 1, 1, 0, A, 1, AF, 2, EQ, C, B, 2, X, 2,
509 $ RCOND, R1, R2, W, RW, INFO )
510 CALL CHKXER( 'cpbsvx', INFOT, NOUT, LERR, OK )
511 INFOT = 9
512 CALL CPBSVX( 'n', 'u', 1, 1, 0, A, 2, AF, 1, EQ, C, B, 2, X, 2,
513 $ RCOND, R1, R2, W, RW, INFO )
514 CALL CHKXER( 'cpbsvx', INFOT, NOUT, LERR, OK )
515 INFOT = 10
516 EQ = '/'
517 CALL CPBSVX( 'f', 'u', 0, 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
518 $ RCOND, R1, R2, W, RW, INFO )
519 CALL CHKXER( 'cpbsvx', INFOT, NOUT, LERR, OK )
520 INFOT = 11
521 EQ = 'y'
522 CALL CPBSVX( 'f', 'u', 1, 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
523 $ RCOND, R1, R2, W, RW, INFO )
524 CALL CHKXER( 'cpbsvx', INFOT, NOUT, LERR, OK )
525 INFOT = 13
526 CALL CPBSVX( 'n', 'u', 2, 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 2,
527 $ RCOND, R1, R2, W, RW, INFO )
528 CALL CHKXER( 'cpbsvx', infot, nout, lerr, ok )
529 infot = 15
530 CALL cpbsvx( 'N', 'U', 2, 0, 0, a, 1, af, 1, eq, c, b, 2, x, 1,
531 $ rcond, r1, r2, w, rw, info )
532 CALL chkxer( 'CPBSVX', infot, nout, lerr, ok )
533*
534 ELSE IF( lsamen( 2, c2, 'PT' ) ) THEN
535*
536* CPTSV
537*
538 srnamt = 'CPTSV '
539 infot = 1
540 CALL cptsv( -1, 0, r, a( 1, 1 ), b, 1, info )
541 CALL chkxer( 'cptsv ', INFOT, NOUT, LERR, OK )
542 INFOT = 2
543 CALL CPTSV( 0, -1, R, A( 1, 1 ), B, 1, INFO )
544 CALL CHKXER( 'cptsv ', INFOT, NOUT, LERR, OK )
545 INFOT = 6
546 CALL CPTSV( 2, 0, R, A( 1, 1 ), B, 1, INFO )
547 CALL CHKXER( 'cptsv ', INFOT, NOUT, LERR, OK )
548*
549* CPTSVX
550*
551 SRNAMT = 'cptsvx'
552 INFOT = 1
553 CALL CPTSVX( '/', 0, 0, R, A( 1, 1 ), RF, AF( 1, 1 ), B, 1, X,
554 $ 1, RCOND, R1, R2, W, RW, INFO )
555 CALL CHKXER( 'cptsvx', INFOT, NOUT, LERR, OK )
556 INFOT = 2
557 CALL CPTSVX( 'n', -1, 0, R, A( 1, 1 ), RF, AF( 1, 1 ), B, 1, X,
558 $ 1, RCOND, R1, R2, W, RW, INFO )
559 CALL CHKXER( 'cptsvx', INFOT, NOUT, LERR, OK )
560 INFOT = 3
561 CALL CPTSVX( 'n', 0, -1, R, A( 1, 1 ), RF, AF( 1, 1 ), B, 1, X,
562 $ 1, RCOND, R1, R2, W, RW, INFO )
563 CALL CHKXER( 'cptsvx', infot, nout, lerr, ok )
564 infot = 9
565 CALL cptsvx( 'N', 2, 0, r, a( 1, 1 ), rf, af( 1, 1 ), b, 1, x,
566 $ 2, rcond, r1, r2, w, rw, info )
567 CALL chkxer( 'CPTSVX', infot, nout, lerr, ok )
568 infot = 11
569 CALL cptsvx( 'N', 2, 0, r, a( 1, 1 ), rf, af( 1, 1 ), b, 2, x,
570 $ 1, rcond, r1, r2, w, rw, info )
571 CALL chkxer( 'CPTSVX', infot, nout, lerr, ok )
572*
573 ELSE IF( lsamen( 2, c2, 'HE' ) ) THEN
574*
575* CHESV
576*
577 srnamt = 'CHESV '
578 infot = 1
579 CALL chesv( '/', 0, 0, a, 1, ip, b, 1, w, 1, info )
580 CALL chkxer( 'CHESV ', infot, nout, lerr, ok )
581 infot = 2
582 CALL chesv( 'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
583 CALL chkxer( 'CHESV ', infot, nout, lerr, ok )
584 infot = 3
585 CALL chesv( 'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
586 CALL chkxer( 'CHESV ', infot, nout, lerr, ok )
587 infot = 5
588 CALL chesv( 'U', 2, 0, a, 1, ip, b, 2, w, 1, info )
589 CALL chkxer( 'chesv ', INFOT, NOUT, LERR, OK )
590 INFOT = 8
591 CALL CHESV( 'u', 2, 0, A, 2, IP, B, 1, W, 1, INFO )
592 CALL CHKXER( 'chesv ', INFOT, NOUT, LERR, OK )
593 INFOT = 10
594 CALL CHESV( 'u', 0, 0, A, 1, IP, B, 1, W, 0, INFO )
595 CALL CHKXER( 'chesv ', INFOT, NOUT, LERR, OK )
596 INFOT = 10
597 CALL CHESV( 'u', 0, 0, A, 1, IP, B, 1, W, -2, INFO )
598 CALL CHKXER( 'chesv ', INFOT, NOUT, LERR, OK )
599*
600* CHESVX
601*
602 SRNAMT = 'chesvx'
603 INFOT = 1
604 CALL CHESVX( '/', 'u', 0, 0, A, 1, AF, 1, IP, B, 1, X, 1,
605 $ RCOND, R1, R2, W, 1, RW, INFO )
606 CALL CHKXER( 'chesvx', INFOT, NOUT, LERR, OK )
607 INFOT = 2
608 CALL CHESVX( 'n', '/', 0, 0, A, 1, AF, 1, IP, B, 1, X, 1,
609 $ RCOND, R1, R2, W, 1, RW, INFO )
610 CALL CHKXER( 'chesvx', INFOT, NOUT, LERR, OK )
611 INFOT = 3
612 CALL CHESVX( 'n', 'u', -1, 0, A, 1, AF, 1, IP, B, 1, X, 1,
613 $ RCOND, R1, R2, W, 1, RW, INFO )
614 CALL CHKXER( 'chesvx', INFOT, NOUT, LERR, OK )
615 INFOT = 4
616 CALL CHESVX( 'n', 'u', 0, -1, A, 1, AF, 1, IP, B, 1, X, 1,
617 $ RCOND, R1, R2, W, 1, RW, INFO )
618 CALL CHKXER( 'chesvx', INFOT, NOUT, LERR, OK )
619 INFOT = 6
620 CALL CHESVX( 'n', 'u', 2, 0, A, 1, AF, 2, IP, B, 2, X, 2,
621 $ RCOND, R1, R2, W, 4, RW, INFO )
622 CALL CHKXER( 'chesvx', INFOT, NOUT, LERR, OK )
623 INFOT = 8
624 CALL CHESVX( 'n', 'u', 2, 0, A, 2, AF, 1, IP, B, 2, X, 2,
625 $ RCOND, R1, R2, W, 4, RW, INFO )
626 CALL CHKXER( 'chesvx', INFOT, NOUT, LERR, OK )
627 INFOT = 11
628 CALL CHESVX( 'n', 'u', 2, 0, A, 2, AF, 2, IP, B, 1, X, 2,
629 $ RCOND, R1, R2, W, 4, RW, INFO )
630 CALL CHKXER( 'chesvx', INFOT, NOUT, LERR, OK )
631 INFOT = 13
632 CALL CHESVX( 'n', 'u', 2, 0, A, 2, AF, 2, IP, B, 2, X, 1,
633 $ RCOND, R1, R2, W, 4, RW, INFO )
634 CALL CHKXER( 'chesvx', INFOT, NOUT, LERR, OK )
635 INFOT = 18
636 CALL CHESVX( 'n', 'u', 2, 0, A, 2, AF, 2, IP, B, 2, X, 2,
637 $ RCOND, R1, R2, W, 3, RW, INFO )
638 CALL CHKXER( 'chesvx', INFOT, NOUT, LERR, OK )
639*
640 ELSE IF( LSAMEN( 2, C2, 'hr' ) ) THEN
641*
642* CHESV_ROOK
643*
644 SRNAMT = 'chesv_rook'
645 INFOT = 1
646 CALL CHESV_ROOK( '/', 0, 0, A, 1, IP, B, 1, W, 1, INFO )
647 CALL CHKXER( 'chesv_rook', INFOT, NOUT, LERR, OK )
648 INFOT = 2
649 CALL CHESV_ROOK( 'u', -1, 0, A, 1, IP, B, 1, W, 1, INFO )
650 CALL CHKXER( 'chesv_rook', INFOT, NOUT, LERR, OK )
651 INFOT = 3
652 CALL CHESV_ROOK( 'u', 0, -1, A, 1, IP, B, 1, W, 1, INFO )
653 CALL CHKXER( 'chesv_rook', INFOT, NOUT, LERR, OK )
654 INFOT = 8
655 CALL CHESV_ROOK( 'u', 2, 0, A, 2, IP, B, 1, W, 1, INFO )
656 CALL CHKXER( 'chesv_rook', INFOT, NOUT, LERR, OK )
657 INFOT = 10
658 CALL CHESV_ROOK( 'u', 0, 0, A, 1, IP, B, 1, W, 0, INFO )
659 CALL CHKXER( 'chesv_rook', INFOT, NOUT, LERR, OK )
660 INFOT = 10
661 CALL CHESV_ROOK( 'u', 0, 0, A, 1, IP, B, 1, W, -2, INFO )
662 CALL CHKXER( 'chesv_rook', INFOT, NOUT, LERR, OK )
663*
664 ELSE IF( LSAMEN( 2, C2, 'hk' ) ) THEN
665*
666* CHESV_RK
667*
668* Test error exits of the driver that uses factorization
669* of a symmetric indefinite matrix with rook
670* (bounded Bunch-Kaufman) pivoting with the new storage
671* format for factors L ( or U) and D.
672*
673* L (or U) is stored in A, diagonal of D is stored on the
674* diagonal of A, subdiagonal of D is stored in a separate array E.
675*
676 SRNAMT = 'chesv_rk'
677 INFOT = 1
678 CALL CHESV_RK( '/', 0, 0, A, 1, E, IP, B, 1, W, 1, INFO )
679 CALL CHKXER( 'chesv_rk', INFOT, NOUT, LERR, OK )
680 INFOT = 2
681 CALL CHESV_RK( 'u', -1, 0, A, 1, E, IP, B, 1, W, 1, INFO )
682 CALL CHKXER( 'chesv_rk', INFOT, NOUT, LERR, OK )
683 INFOT = 3
684 CALL CHESV_RK( 'u', 0, -1, A, 1, E, IP, B, 1, W, 1, INFO )
685 CALL CHKXER( 'chesv_rk', INFOT, NOUT, LERR, OK )
686 INFOT = 5
687 CALL CHESV_RK( 'u', 2, 0, A, 1, E, IP, B, 2, W, 1, INFO )
688 CALL CHKXER( 'chesv_rk', INFOT, NOUT, LERR, OK )
689 INFOT = 9
690 CALL CHESV_RK( 'u', 2, 0, A, 2, E, IP, B, 1, W, 1, INFO )
691 CALL CHKXER( 'chesv_rk', INFOT, NOUT, LERR, OK )
692 INFOT = 11
693 CALL CHESV_RK( 'u', 0, 0, A, 1, E, IP, B, 1, W, 0, INFO )
694 CALL CHKXER( 'chesv_rk', INFOT, NOUT, LERR, OK )
695 INFOT = 11
696 CALL CHESV_RK( 'u', 0, 0, A, 1, E, IP, B, 1, W, -2, INFO )
697 CALL CHKXER( 'chesv_rk', INFOT, NOUT, LERR, OK )
698*
699 ELSE IF( LSAMEN( 2, C2, 'ha' ) ) THEN
700*
701* CHESV_AASEN
702*
703 SRNAMT = 'chesv_aa'
704 INFOT = 1
705 CALL CHESV_AA( '/', 0, 0, A, 1, IP, B, 1, W, 1, INFO )
706 CALL CHKXER( 'chesv_aa', INFOT, NOUT, LERR, OK )
707 INFOT = 2
708 CALL CHESV_AA( 'u', -1, 0, A, 1, IP, B, 1, W, 1, INFO )
709 CALL CHKXER( 'chesv_aa', INFOT, NOUT, LERR, OK )
710 INFOT = 3
711 CALL CHESV_AA( 'u', 0, -1, A, 1, IP, B, 1, W, 1, INFO )
712 CALL CHKXER( 'chesv_aa', INFOT, NOUT, LERR, OK )
713 INFOT = 8
714 CALL CHESV_AA( 'u', 2, 0, A, 2, IP, B, 1, W, 1, INFO )
715 CALL CHKXER( 'chesv_aa', INFOT, NOUT, LERR, OK )
716*
717 ELSE IF( LSAMEN( 2, C2, 'h2' ) ) THEN
718*
719* CHESV_AASEN_2STAGE
720*
721 SRNAMT = 'chesv_aa_2stage'
722 INFOT = 1
723 CALL CHESV_AA_2STAGE( '/', 0, 0, A, 1, A, 1, IP, IP, B, 1,
724 $ W, 1, INFO )
725 CALL CHKXER( 'chesv_aa_2stage', INFOT, NOUT, LERR, OK )
726 INFOT = 2
727 CALL CHESV_AA_2STAGE( 'u', -1, 0, A, 1, A, 1, IP, IP, B, 1,
728 $ W, 1, INFO )
729 CALL CHKXER( 'chesv_aa_2stage', INFOT, NOUT, LERR, OK )
730 INFOT = 3
731 CALL CHESV_AA_2STAGE( 'u', 0, -1, A, 1, A, 1, IP, IP, B, 1,
732 $ W, 1, INFO )
733 CALL CHKXER( 'chesv_aa_2stage', INFOT, NOUT, LERR, OK )
734 INFOT = 5
735 CALL CHESV_AA_2STAGE( 'u', 2, 1, A, 1, A, 1, IP, IP, B, 1,
736 $ W, 1, INFO )
737 CALL CHKXER( 'chesv_aa_2stage', INFOT, NOUT, LERR, OK )
738 INFOT = 11
739 CALL CHESV_AA_2STAGE( 'u', 2, 1, A, 2, A, 8, IP, IP, B, 1,
740 $ W, 1, INFO )
741 CALL CHKXER( 'chesv_aa_2stage', INFOT, NOUT, LERR, OK )
742 INFOT = 7
743 CALL CHESV_AA_2STAGE( 'u', 2, 1, a, 2, a, 1, ip, ip, b, 2,
744 $ w, 1, info )
745 CALL chkxer( 'CHESV_AA_2STAGE', infot, nout, lerr, ok )
746*
747 ELSE IF( lsamen( 2, c2, 'S2' ) ) THEN
748*
749* CSYSV_AASEN_2STAGE
750*
751 srnamt = 'CSYSV_AA_2STAGE'
752 infot = 1
753 CALL csysv_aa_2stage( '/', 0, 0, a, 1, a, 1, ip, ip, b, 1,
754 $ w, 1, info )
755 CALL chkxer( 'CSYSV_AA_2STAGE', infot, nout, lerr, ok )
756 infot = 2
757 CALL csysv_aa_2stage( 'U', -1, 0, a, 1, a, 1, ip, ip, b, 1,
758 $ w, 1, info )
759 CALL chkxer( 'CSYSV_AA_2STAGE', infot, nout, lerr, ok )
760 infot = 3
761 CALL csysv_aa_2stage( 'U', 0, -1, a, 1, a, 1, ip, ip, b, 1,
762 $ w, 1, info )
763 CALL chkxer( 'CSYSV_AA_2STAGE', infot, nout, lerr, ok )
764 infot = 5
765 CALL csysv_aa_2stage( 'U', 2, 1, a, 1, a, 1, ip, ip, b, 1,
766 $ w, 1, info )
767 CALL chkxer( 'CSYSV_AA_2STAGE', infot, nout, lerr, ok )
768 infot = 11
769 CALL csysv_aa_2stage( 'U', 2, 1, a, 2, a, 8, ip, ip, b, 1,
770 $ w, 1, info )
771 CALL chkxer( 'CSYSV_AA_2STAGE', infot, nout, lerr, ok )
772 infot = 7
773 CALL csysv_aa_2stage( 'U', 2, 1, a, 2, a, 1, ip, ip, b, 2,
774 $ w, 1, info )
775 CALL chkxer( 'CSYSV_AA_2STAGE', infot, nout, lerr, ok )
776*
777 ELSE IF( lsamen( 2, c2, 'HP' ) ) THEN
778*
779* CHPSV
780*
781 srnamt = 'CHPSV '
782 infot = 1
783 CALL chpsv( '/', 0, 0, a, ip, b, 1, info )
784 CALL chkxer( 'CHPSV ', infot, nout, lerr, ok )
785 infot = 2
786 CALL chpsv( 'U', -1, 0, a, ip, b, 1, info )
787 CALL chkxer( 'CHPSV ', infot, nout, lerr, ok )
788 infot = 3
789 CALL chpsv( 'U', 0, -1, a, ip, b, 1, info )
790 CALL chkxer( 'CHPSV ', infot, nout, lerr, ok )
791 infot = 7
792 CALL chpsv( 'U', 2, 0, a, ip, b, 1, info )
793 CALL chkxer( 'CHPSV ', infot, nout, lerr, ok )
794*
795* CHPSVX
796*
797 srnamt = 'CHPSVX'
798 infot = 1
799 CALL chpsvx( '/', 'U', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
800 $ r2, w, rw, info )
801 CALL chkxer( 'CHPSVX', infot, nout, lerr, ok )
802 infot = 2
803 CALL chpsvx( 'N', '/', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
804 $ r2, w, rw, info )
805 CALL chkxer( 'CHPSVX', infot, nout, lerr, ok )
806 infot = 3
807 CALL chpsvx( 'N', 'U', -1, 0, a, af, ip, b, 1, x, 1, rcond, r1,
808 $ r2, w, rw, info )
809 CALL chkxer( 'CHPSVX', infot, nout, lerr, ok )
810 infot = 4
811 CALL chpsvx( 'N', 'U', 0, -1, a, af, ip, b, 1, x, 1, rcond, r1,
812 $ r2, w, rw, info )
813 CALL chkxer( 'CHPSVX', infot, nout, lerr, ok )
814 infot = 9
815 CALL chpsvx( 'N', 'U', 2, 0, a, af, ip, b, 1, x, 2, rcond, r1,
816 $ r2, w, rw, info )
817 CALL chkxer( 'CHPSVX', infot, nout, lerr, ok )
818 infot = 11
819 CALL chpsvx( 'N', 'U', 2, 0, a, af, ip, b, 2, x, 1, rcond, r1,
820 $ r2, w, rw, info )
821 CALL chkxer( 'CHPSVX', infot, nout, lerr, ok )
822*
823 ELSE IF( lsamen( 2, c2, 'SY' ) ) THEN
824*
825* CSYSV
826*
827 srnamt = 'CSYSV '
828 infot = 1
829 CALL csysv( '/', 0, 0, a, 1, ip, b, 1, w, 1, info )
830 CALL chkxer( 'CSYSV ', infot, nout, lerr, ok )
831 infot = 2
832 CALL csysv( 'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
833 CALL chkxer( 'CSYSV ', infot, nout, lerr, ok )
834 infot = 3
835 CALL csysv( 'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
836 CALL chkxer( 'CSYSV ', infot, nout, lerr, ok )
837 infot = 8
838 CALL csysv( 'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
839 CALL chkxer( 'csysv ', INFOT, NOUT, LERR, OK )
840 INFOT = 10
841 CALL CSYSV( 'u', 0, 0, A, 1, IP, B, 1, W, 0, INFO )
842 CALL CHKXER( 'csysv ', INFOT, NOUT, LERR, OK )
843 INFOT = 10
844 CALL CSYSV( 'u', 0, 0, A, 1, IP, B, 1, W, -2, INFO )
845 CALL CHKXER( 'csysv ', INFOT, NOUT, LERR, OK )
846*
847* CSYSVX
848*
849 SRNAMT = 'csysvx'
850 INFOT = 1
851 CALL CSYSVX( '/', 'u', 0, 0, A, 1, AF, 1, IP, B, 1, X, 1,
852 $ RCOND, R1, R2, W, 1, RW, INFO )
853 CALL CHKXER( 'csysvx', INFOT, NOUT, LERR, OK )
854 INFOT = 2
855 CALL CSYSVX( 'n', '/', 0, 0, A, 1, AF, 1, IP, B, 1, X, 1,
856 $ RCOND, R1, R2, W, 1, RW, INFO )
857 CALL CHKXER( 'csysvx', INFOT, NOUT, LERR, OK )
858 INFOT = 3
859 CALL CSYSVX( 'n', 'u', -1, 0, a, 1, af, 1, ip, b, 1, x, 1,
860 $ rcond, r1, r2, w, 1, rw, info )
861 CALL chkxer( 'CSYSVX', infot, nout, lerr, ok )
862 infot = 4
863 CALL csysvx( 'N', 'U', 0, -1, a, 1, af, 1, ip, b, 1, x, 1,
864 $ rcond, r1, r2, w, 1, rw, info )
865 CALL chkxer( 'CSYSVX', infot, nout, lerr, ok )
866 infot = 6
867 CALL csysvx( 'N', 'U', 2, 0, a, 1, af, 2, ip, b, 2, x, 2,
868 $ rcond, r1, r2, w, 4, rw, info )
869 CALL chkxer( 'CSYSVX', infot, nout, lerr, ok )
870 infot = 8
871 CALL csysvx( 'N', 'U', 2, 0, a, 2, af, 1, ip, b, 2, x, 2,
872 $ rcond, r1, r2, w, 4, rw, info )
873 CALL chkxer( 'CSYSVX', infot, nout, lerr, ok )
874 infot = 11
875 CALL csysvx( 'N', 'U', 2, 0, a, 2, af, 2, ip, b, 1, x, 2,
876 $ rcond, r1, r2, w, 4, rw, info )
877 CALL chkxer( 'CSYSVX', infot, nout, lerr, ok )
878 infot = 13
879 CALL csysvx( 'N', 'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 1,
880 $ rcond, r1, r2, w, 4, rw, info )
881 CALL chkxer( 'CSYSVX', infot, nout, lerr, ok )
882 infot = 18
883 CALL csysvx( 'N', 'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 2,
884 $ rcond, r1, r2, w, 3, rw, info )
885 CALL chkxer( 'CSYSVX', infot, nout, lerr, ok )
886*
887 ELSE IF( lsamen( 2, c2, 'SR' ) ) THEN
888*
889* CSYSV_ROOK
890*
891 srnamt = 'CSYSV_ROOK'
892 infot = 1
893 CALL csysv_rook( '/', 0, 0, a, 1, ip, b, 1, w, 1, info )
894 CALL chkxer( 'CSYSV_ROOK', infot, nout, lerr, ok )
895 infot = 2
896 CALL csysv_rook( 'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
897 CALL chkxer( 'CSYSV_ROOK', infot, nout, lerr, ok )
898 infot = 3
899 CALL csysv_rook( 'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
900 CALL chkxer( 'CSYSV_ROOK', infot, nout, lerr, ok )
901 infot = 8
902 CALL csysv_rook( 'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
903 CALL chkxer( 'CSYSV_ROOK', infot, nout, lerr, ok )
904 infot = 10
905 CALL csysv_rook( 'U', 0, 0, a, 1, ip, b, 1, w, 0, info )
906 CALL chkxer( 'CSYSV_ROOK', infot, nout, lerr, ok )
907 infot = 10
908 CALL csysv_rook( 'U', 0, 0, a, 1, ip, b, 1, w, -2, info )
909 CALL chkxer( 'csysv_rook', INFOT, NOUT, LERR, OK )
910*
911 ELSE IF( LSAMEN( 2, C2, 'sk' ) ) THEN
912*
913* CSYSV_RK
914*
915* Test error exits of the driver that uses factorization
916* of a symmetric indefinite matrix with rook
917* (bounded Bunch-Kaufman) pivoting with the new storage
918* format for factors L ( or U) and D.
919*
920* L (or U) is stored in A, diagonal of D is stored on the
921* diagonal of A, subdiagonal of D is stored in a separate array E.
922*
923 SRNAMT = 'csysv_rk'
924 INFOT = 1
925 CALL CSYSV_RK( '/', 0, 0, A, 1, E, IP, B, 1, W, 1, INFO )
926 CALL CHKXER( 'csysv_rk', INFOT, NOUT, LERR, OK )
927 INFOT = 2
928 CALL CSYSV_RK( 'u', -1, 0, A, 1, E, IP, B, 1, W, 1, INFO )
929 CALL CHKXER( 'csysv_rk', INFOT, NOUT, LERR, OK )
930 INFOT = 3
931 CALL CSYSV_RK( 'u', 0, -1, A, 1, E, IP, B, 1, W, 1, INFO )
932 CALL CHKXER( 'csysv_rk', INFOT, NOUT, LERR, OK )
933 INFOT = 5
934 CALL CSYSV_RK( 'u', 2, 0, A, 1, E, IP, B, 2, W, 1, INFO )
935 CALL CHKXER( 'csysv_rk', INFOT, NOUT, LERR, OK )
936 INFOT = 9
937 CALL CSYSV_RK( 'u', 2, 0, A, 2, E, IP, B, 1, W, 1, INFO )
938 CALL CHKXER( 'csysv_rk', INFOT, NOUT, LERR, OK )
939 INFOT = 11
940 CALL CSYSV_RK( 'u', 0, 0, A, 1, E, IP, B, 1, W, 0, INFO )
941 CALL CHKXER( 'csysv_rk', INFOT, NOUT, LERR, OK )
942 INFOT = 11
943 CALL CSYSV_RK( 'u', 0, 0, A, 1, E, IP, B, 1, W, -2, INFO )
944 CALL CHKXER( 'csysv_rk', INFOT, NOUT, LERR, OK )
945*
946 ELSE IF( LSAMEN( 2, C2, 'sp' ) ) THEN
947*
948* CSPSV
949*
950 SRNAMT = 'cspsv '
951 INFOT = 1
952 CALL CSPSV( '/', 0, 0, A, IP, B, 1, INFO )
953 CALL CHKXER( 'cspsv ', INFOT, NOUT, LERR, OK )
954 INFOT = 2
955 CALL CSPSV( 'u', -1, 0, A, IP, B, 1, INFO )
956 CALL CHKXER( 'cspsv ', INFOT, NOUT, LERR, OK )
957 INFOT = 3
958 CALL CSPSV( 'u', 0, -1, A, IP, B, 1, INFO )
959 CALL CHKXER( 'cspsv ', INFOT, NOUT, LERR, OK )
960 INFOT = 7
961 CALL CSPSV( 'u', 2, 0, A, IP, B, 1, INFO )
962 CALL CHKXER( 'cspsv ', INFOT, NOUT, LERR, OK )
963*
964* CSPSVX
965*
966 SRNAMT = 'cspsvx'
967 INFOT = 1
968 CALL CSPSVX( '/', 'u', 0, 0, A, AF, IP, B, 1, X, 1, RCOND, R1,
969 $ R2, W, RW, INFO )
970 CALL CHKXER( 'cspsvx', INFOT, NOUT, LERR, OK )
971 INFOT = 2
972 CALL CSPSVX( 'n', '/', 0, 0, A, AF, IP, B, 1, X, 1, RCOND, R1,
973 $ R2, W, RW, INFO )
974 CALL CHKXER( 'cspsvx', INFOT, NOUT, LERR, OK )
975 INFOT = 3
976 CALL CSPSVX( 'n', 'u', -1, 0, A, AF, IP, B, 1, X, 1, RCOND, R1,
977 $ R2, W, RW, INFO )
978 CALL CHKXER( 'cspsvx', INFOT, NOUT, LERR, OK )
979 INFOT = 4
980 CALL CSPSVX( 'n', 'u', 0, -1, A, AF, IP, B, 1, X, 1, RCOND, R1,
981 $ R2, W, RW, INFO )
982 CALL CHKXER( 'cspsvx', INFOT, NOUT, LERR, OK )
983 INFOT = 9
984 CALL CSPSVX( 'n', 'u', 2, 0, A, AF, IP, B, 1, X, 2, RCOND, R1,
985 $ R2, W, RW, INFO )
986 CALL CHKXER( 'cspsvx', INFOT, NOUT, LERR, OK )
987 INFOT = 11
988 CALL CSPSVX( 'n', 'u', 2, 0, A, AF, IP, B, 2, X, 1, RCOND, R1,
989 $ R2, W, RW, INFO )
990 CALL CHKXER( 'cspsvx', INFOT, NOUT, LERR, OK )
991 END IF
992*
993* Print a summary line.
994*
995 IF( OK ) THEN
996 WRITE( NOUT, FMT = 9999 )PATH
997 ELSE
998 WRITE( NOUT, FMT = 9998 )PATH
999 END IF
1000*
1001 9999 FORMAT( 1X, A3, ' drivers passed the tests of the error exits' )
1002 9998 FORMAT( ' *** ', A3, ' drivers failed the tests of the error ',
1003 $ 'exits ***' )
1004*
1005 RETURN
1006*
1007* End of CERRVX
1008*
1009 END
cmplx
float cmplx[2]
Definition pblas.h:136
chkxer
subroutine chkxer(srnamt, infot, nout, lerr, ok)
Definition cblat2.f:3196
the
end diagonal values have been computed in the(sparse) matrix id.SOL
cgbsvx
subroutine cgbsvx(fact, trans, n, kl, ku, nrhs, ab, ldab, afb, ldafb, ipiv, equed, r, c, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
CGBSVX computes the solution to system of linear equations A * X = B for GB matrices
Definition cgbsvx.f:370
cgbsv
subroutine cgbsv(n, kl, ku, nrhs, ab, ldab, ipiv, b, ldb, info)
CGBSV computes the solution to system of linear equations A * X = B for GB matrices (simple driver)
Definition cgbsv.f:162
cgesv
subroutine cgesv(n, nrhs, a, lda, ipiv, b, ldb, info)
CGESV computes the solution to system of linear equations A * X = B for GE matrices (simple driver)
Definition cgesv.f:122
cgesvx
subroutine cgesvx(fact, trans, n, nrhs, a, lda, af, ldaf, ipiv, equed, r, c, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
CGESVX computes the solution to system of linear equations A * X = B for GE matrices
Definition cgesvx.f:350
cgtsv
subroutine cgtsv(n, nrhs, dl, d, du, b, ldb, info)
CGTSV computes the solution to system of linear equations A * X = B for GT matrices
Definition cgtsv.f:124
cgtsvx
subroutine cgtsvx(fact, trans, n, nrhs, dl, d, du, dlf, df, duf, du2, ipiv, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
CGTSVX computes the solution to system of linear equations A * X = B for GT matrices
Definition cgtsvx.f:294
chesv_rk
subroutine chesv_rk(uplo, n, nrhs, a, lda, e, ipiv, b, ldb, work, lwork, info)
CHESV_RK computes the solution to system of linear equations A * X = B for SY matrices
Definition chesv_rk.f:228
chesv_rook
subroutine chesv_rook(uplo, n, nrhs, a, lda, ipiv, b, ldb, work, lwork, info)
CHESV_ROOK computes the solution to a system of linear equations A * X = B for HE matrices using the ...
Definition chesv_rook.f:205
chesvx
subroutine chesvx(fact, uplo, n, nrhs, a, lda, af, ldaf, ipiv, b, ldb, x, ldx, rcond, ferr, berr, work, lwork, rwork, info)
CHESVX computes the solution to system of linear equations A * X = B for HE matrices
Definition chesvx.f:285
chesv_aa
subroutine chesv_aa(uplo, n, nrhs, a, lda, ipiv, b, ldb, work, lwork, info)
CHESV_AA computes the solution to system of linear equations A * X = B for HE matrices
Definition chesv_aa.f:162
chesv
subroutine chesv(uplo, n, nrhs, a, lda, ipiv, b, ldb, work, lwork, info)
CHESV computes the solution to system of linear equations A * X = B for HE matrices
Definition chesv.f:171
chpsv
subroutine chpsv(uplo, n, nrhs, ap, ipiv, b, ldb, info)
CHPSV computes the solution to system of linear equations A * X = B for OTHER matrices
Definition chpsv.f:162
cspsvx
subroutine cspsvx(fact, uplo, n, nrhs, ap, afp, ipiv, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
CSPSVX computes the solution to system of linear equations A * X = B for OTHER matrices
Definition cspsvx.f:277
cpbsv
subroutine cpbsv(uplo, n, kd, nrhs, ab, ldab, b, ldb, info)
CPBSV computes the solution to system of linear equations A * X = B for OTHER matrices
Definition cpbsv.f:164
cppsvx
subroutine cppsvx(fact, uplo, n, nrhs, ap, afp, equed, s, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
CPPSVX computes the solution to system of linear equations A * X = B for OTHER matrices
Definition cppsvx.f:311
cpbsvx
subroutine cpbsvx(fact, uplo, n, kd, nrhs, ab, ldab, afb, ldafb, equed, s, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
CPBSVX computes the solution to system of linear equations A * X = B for OTHER matrices
Definition cpbsvx.f:342
cspsv
subroutine cspsv(uplo, n, nrhs, ap, ipiv, b, ldb, info)
CSPSV computes the solution to system of linear equations A * X = B for OTHER matrices
Definition cspsv.f:162
cppsv
subroutine cppsv(uplo, n, nrhs, ap, b, ldb, info)
CPPSV computes the solution to system of linear equations A * X = B for OTHER matrices
Definition cppsv.f:144
chpsvx
subroutine chpsvx(fact, uplo, n, nrhs, ap, afp, ipiv, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
CHPSVX computes the solution to system of linear equations A * X = B for OTHER matrices
Definition chpsvx.f:277
cposv
subroutine cposv(uplo, n, nrhs, a, lda, b, ldb, info)
CPOSV computes the solution to system of linear equations A * X = B for PO matrices
Definition cposv.f:130
cposvx
subroutine cposvx(fact, uplo, n, nrhs, a, lda, af, ldaf, equed, s, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
CPOSVX computes the solution to system of linear equations A * X = B for PO matrices
Definition cposvx.f:306
cptsvx
subroutine cptsvx(fact, n, nrhs, d, e, df, ef, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
CPTSVX computes the solution to system of linear equations A * X = B for PT matrices
Definition cptsvx.f:234
cptsv
subroutine cptsv(n, nrhs, d, e, b, ldb, info)
CPTSV computes the solution to system of linear equations A * X = B for PT matrices
Definition cptsv.f:115
chesv_aa_2stage
subroutine chesv_aa_2stage(uplo, n, nrhs, a, lda, tb, ltb, ipiv, ipiv2, b, ldb, work, lwork, info)
CHESV_AA_2STAGE computes the solution to system of linear equations A * X = B for HE matrices
Definition chesv_aa_2stage.f:186
csysv_aa_2stage
subroutine csysv_aa_2stage(uplo, n, nrhs, a, lda, tb, ltb, ipiv, ipiv2, b, ldb, work, lwork, info)
CSYSV_AA_2STAGE computes the solution to system of linear equations A * X = B for SY matrices
Definition csysv_aa_2stage.f:186
csysv_rk
subroutine csysv_rk(uplo, n, nrhs, a, lda, e, ipiv, b, ldb, work, lwork, info)
CSYSV_RK computes the solution to system of linear equations A * X = B for SY matrices
Definition csysv_rk.f:228
csysv
subroutine csysv(uplo, n, nrhs, a, lda, ipiv, b, ldb, work, lwork, info)
CSYSV computes the solution to system of linear equations A * X = B for SY matrices
Definition csysv.f:171
csysv_aa
subroutine csysv_aa(uplo, n, nrhs, a, lda, ipiv, b, ldb, work, lwork, info)
CSYSV_AA computes the solution to system of linear equations A * X = B for SY matrices
Definition csysv_aa.f:162
csysv_rook
subroutine csysv_rook(uplo, n, nrhs, a, lda, ipiv, b, ldb, work, lwork, info)
CSYSV_ROOK computes the solution to system of linear equations A * X = B for SY matrices
Definition csysv_rook.f:204
csysvx
subroutine csysvx(fact, uplo, n, nrhs, a, lda, af, ldaf, ipiv, b, ldb, x, ldx, rcond, ferr, berr, work, lwork, rwork, info)
CSYSVX computes the solution to system of linear equations A * X = B for SY matrices
Definition csysvx.f:285
cerrvx
subroutine cerrvx(path, nunit)
CERRVX
Definition cerrvx.f:55