OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
serrvx.f
Go to the documentation of this file.
1*> \brief \b SERRVX
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 SERRVX( 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*> SERRVX tests the error exits for the REAL 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 single_lin
52*
53* =====================================================================
54 SUBROUTINE serrvx( 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 ), IW( NMAX )
79 REAL A( NMAX, NMAX ), AF( NMAX, NMAX ), B( NMAX ),
80 $ C( NMAX ), E( NMAX ), R( NMAX ), R1( NMAX ),
81 $ R2( NMAX ), W( 2*NMAX ), X( NMAX )
82* ..
83* .. External Functions ..
84 LOGICAL LSAMEN
85 EXTERNAL lsamen
86* ..
87* .. External Subroutines ..
88 EXTERNAL chkxer, sgbsv, sgbsvx, sgesv, sgesvx, sgtsv,
89 $ sgtsvx, spbsv, spbsvx, sposv, sposvx, sppsv,
90 $ sppsvx, sptsv, sptsvx, sspsv, sspsvx, ssysv,
91 $ ssysv_aa, ssysv_rk, ssysv_rook, ssysvx,
92 $ ssysv_aa_2stage
93* ..
94* .. Scalars in Common ..
95 LOGICAL LERR, OK
96 CHARACTER*32 SRNAMT
97 INTEGER INFOT, NOUT
98* ..
99* .. Common blocks ..
100 COMMON / infoc / infot, nout, ok, lerr
101 COMMON / srnamc / srnamt
102* ..
103* .. Intrinsic Functions ..
104 INTRINSIC real
105* ..
106* .. Executable Statements ..
107*
108 nout = nunit
109 WRITE( nout, fmt = * )
110 c2 = path( 2: 3 )
111*
112* Set the variables to innocuous values.
113*
114 DO 20 j = 1, nmax
115 DO 10 i = 1, nmax
116 a( i, j ) = 1. / real( i+j )
117 af( i, j ) = 1. / real( i+j )
118 10 CONTINUE
119 b( j ) = 0.e+0
120 e( j ) = 0.e+0
121 r1( j ) = 0.e+0
122 r2( j ) = 0.e+0
123 w( j ) = 0.e+0
124 x( j ) = 0.e+0
125 c( j ) = 0.e+0
126 r( j ) = 0.e+0
127 ip( j ) = j
128 20 CONTINUE
129 eq = ' '
130 ok = .true.
131*
132 IF( lsamen( 2, c2, 'GE' ) ) THEN
133*
134* SGESV
135*
136 srnamt = 'SGESV '
137 infot = 1
138 CALL sgesv( -1, 0, a, 1, ip, b, 1, info )
139 CALL chkxer( 'SGESV ', infot, nout, lerr, ok )
140 infot = 2
141 CALL sgesv( 0, -1, a, 1, ip, b, 1, info )
142 CALL chkxer( 'SGESV ', infot, nout, lerr, ok )
143 infot = 4
144 CALL sgesv( 2, 1, a, 1, ip, b, 2, info )
145 CALL chkxer( 'SGESV ', infot, nout, lerr, ok )
146 infot = 7
147 CALL sgesv( 2, 1, a, 2, ip, b, 1, info )
148 CALL chkxer( 'SGESV ', infot, nout, lerr, ok )
149*
150* SGESVX
151*
152 srnamt = 'SGESVX'
153 infot = 1
154 CALL sgesvx( '/', 'N', 0, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
155 $ x, 1, rcond, r1, r2, w, iw, info )
156 CALL chkxer( 'SGESVX', infot, nout, lerr, ok )
157 infot = 2
158 CALL sgesvx( 'N', '/', 0, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
159 $ x, 1, rcond, r1, r2, w, iw, info )
160 CALL chkxer( 'SGESVX', infot, nout, lerr, ok )
161 infot = 3
162 CALL sgesvx( 'N', 'N', -1, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
163 $ x, 1, rcond, r1, r2, w, iw, info )
164 CALL chkxer( 'SGESVX', infot, nout, lerr, ok )
165 infot = 4
166 CALL sgesvx( 'N', 'N', 0, -1, a, 1, af, 1, ip, eq, r, c, b, 1,
167 $ x, 1, rcond, r1, r2, w, iw, info )
168 CALL chkxer( 'SGESVX', infot, nout, lerr, ok )
169 infot = 6
170 CALL sgesvx( 'N', 'N', 2, 1, a, 1, af, 2, ip, eq, r, c, b, 2,
171 $ x, 2, rcond, r1, r2, w, iw, info )
172 CALL chkxer( 'SGESVX', infot, nout, lerr, ok )
173 infot = 8
174 CALL sgesvx( 'N', 'n', 2, 1, A, 2, AF, 1, IP, EQ, R, C, B, 2,
175 $ X, 2, RCOND, R1, R2, W, IW, INFO )
176 CALL CHKXER( 'sgesvx', INFOT, NOUT, LERR, OK )
177 INFOT = 10
178 EQ = '/'
179 CALL SGESVX( 'f', 'n', 0, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
180 $ X, 1, RCOND, R1, R2, W, IW, INFO )
181 CALL CHKXER( 'sgesvx', INFOT, NOUT, LERR, OK )
182 INFOT = 11
183 EQ = 'r'
184 CALL SGESVX( 'f', 'n', 1, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
185 $ X, 1, RCOND, R1, R2, W, IW, INFO )
186 CALL CHKXER( 'sgesvx', INFOT, NOUT, LERR, OK )
187 INFOT = 12
188 EQ = 'c'
189 CALL SGESVX( 'f', 'n', 1, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
190 $ X, 1, RCOND, R1, R2, W, IW, INFO )
191 CALL CHKXER( 'sgesvx', INFOT, NOUT, LERR, OK )
192 INFOT = 14
193 CALL SGESVX( 'n', 'n', 2, 1, A, 2, AF, 2, IP, EQ, R, C, B, 1,
194 $ X, 2, RCOND, R1, R2, W, IW, INFO )
195 CALL CHKXER( 'sgesvx', INFOT, NOUT, LERR, OK )
196 INFOT = 16
197 CALL SGESVX( 'n', 'n', 2, 1, A, 2, AF, 2, IP, EQ, R, C, B, 2,
198 $ X, 1, RCOND, R1, R2, W, IW, INFO )
199 CALL CHKXER( 'sgesvx', INFOT, NOUT, LERR, OK )
200*
201 ELSE IF( LSAMEN( 2, C2, 'gb' ) ) THEN
202*
203* SGBSV
204*
205 SRNAMT = 'sgbsv '
206 INFOT = 1
207 CALL SGBSV( -1, 0, 0, 0, A, 1, IP, B, 1, INFO )
208 CALL CHKXER( 'sgbsv ', INFOT, NOUT, LERR, OK )
209 INFOT = 2
210 CALL SGBSV( 1, -1, 0, 0, A, 1, IP, B, 1, INFO )
211 CALL CHKXER( 'sgbsv ', INFOT, NOUT, LERR, OK )
212 INFOT = 3
213 CALL SGBSV( 1, 0, -1, 0, A, 1, IP, B, 1, INFO )
214 CALL CHKXER( 'sgbsv ', INFOT, NOUT, LERR, OK )
215 INFOT = 4
216 CALL SGBSV( 0, 0, 0, -1, A, 1, IP, B, 1, INFO )
217 CALL CHKXER( 'sgbsv ', INFOT, NOUT, LERR, OK )
218 INFOT = 6
219 CALL SGBSV( 1, 1, 1, 0, A, 3, IP, B, 1, INFO )
220 CALL CHKXER( 'sgbsv ', INFOT, NOUT, LERR, OK )
221 INFOT = 9
222 CALL SGBSV( 2, 0, 0, 0, A, 1, IP, B, 1, INFO )
223 CALL CHKXER( 'sgbsv ', INFOT, NOUT, LERR, OK )
224*
225* SGBSVX
226*
227 SRNAMT = 'sgbsvx'
228 INFOT = 1
229 CALL SGBSVX( '/', 'n', 0, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
230 $ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
231 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
232 INFOT = 2
233 CALL SGBSVX( 'n', '/', 0, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
234 $ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
235 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
236 INFOT = 3
237 CALL SGBSVX( 'n', 'n', -1, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
238 $ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
239 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
240 INFOT = 4
241 CALL SGBSVX( 'n', 'n', 1, -1, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
242 $ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
243 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
244 INFOT = 5
245 CALL SGBSVX( 'n', 'n', 1, 0, -1, 0, A, 1, AF, 1, IP, EQ, R, C,
246 $ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
247 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
248 INFOT = 6
249 CALL SGBSVX( 'n', 'n', 0, 0, 0, -1, A, 1, AF, 1, IP, EQ, R, C,
250 $ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
251 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
252 INFOT = 8
253 CALL SGBSVX( 'n', 'n', 1, 1, 1, 0, A, 2, AF, 4, IP, EQ, R, C,
254 $ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
255 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
256 INFOT = 10
257 CALL SGBSVX( 'n', 'n', 1, 1, 1, 0, A, 3, AF, 3, IP, EQ, R, C,
258 $ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
259 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
260 INFOT = 12
261 EQ = '/'
262 CALL SGBSVX( 'f', 'n', 0, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
263 $ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
264 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
265 INFOT = 13
266 EQ = 'r'
267 CALL SGBSVX( 'f', 'n', 1, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
268 $ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
269 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
270 INFOT = 14
271 EQ = 'c'
272 CALL SGBSVX( 'f', 'n', 1, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
273 $ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
274 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
275 INFOT = 16
276 CALL SGBSVX( 'n', 'n', 2, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
277 $ B, 1, X, 2, RCOND, R1, R2, W, IW, INFO )
278 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
279 INFOT = 18
280 CALL SGBSVX( 'n', 'n', 2, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
281 $ B, 2, X, 1, RCOND, R1, R2, W, IW, INFO )
282 CALL CHKXER( 'sgbsvx', INFOT, NOUT, LERR, OK )
283*
284 ELSE IF( LSAMEN( 2, C2, 'gt' ) ) THEN
285*
286* SGTSV
287*
288 SRNAMT = 'sgtsv '
289 INFOT = 1
290 CALL SGTSV( -1, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), B, 1,
291 $ INFO )
292 CALL CHKXER( 'sgtsv ', INFOT, NOUT, LERR, OK )
293 INFOT = 2
294 CALL SGTSV( 0, -1, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), B, 1,
295 $ INFO )
296 CALL CHKXER( 'sgtsv ', INFOT, NOUT, LERR, OK )
297 INFOT = 7
298 CALL SGTSV( 2, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), B, 1, INFO )
299 CALL CHKXER( 'sgtsv ', INFOT, NOUT, LERR, OK )
300*
301* SGTSVX
302*
303 SRNAMT = 'sgtsvx'
304 INFOT = 1
305 CALL SGTSVX( '/', 'n', 0, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
306 $ AF( 1, 1 ), AF( 1, 2 ), AF( 1, 3 ), AF( 1, 4 ),
307 $ IP, B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
308 CALL CHKXER( 'sgtsvx', INFOT, NOUT, LERR, OK )
309 INFOT = 2
310 CALL SGTSVX( 'n', '/', 0, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
311 $ AF( 1, 1 ), AF( 1, 2 ), AF( 1, 3 ), AF( 1, 4 ),
312 $ IP, B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
313 CALL CHKXER( 'sgtsvx', INFOT, NOUT, LERR, OK )
314 INFOT = 3
315 CALL SGTSVX( 'n', 'n', -1, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
316 $ AF( 1, 1 ), AF( 1, 2 ), AF( 1, 3 ), AF( 1, 4 ),
317 $ IP, B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
318 CALL CHKXER( 'sgtsvx', INFOT, NOUT, LERR, OK )
319 INFOT = 4
320 CALL SGTSVX( 'n', 'n', 0, -1, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
321 $ AF( 1, 1 ), AF( 1, 2 ), AF( 1, 3 ), AF( 1, 4 ),
322 $ IP, B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
323 CALL CHKXER( 'sgtsvx', INFOT, NOUT, LERR, OK )
324 INFOT = 14
325 CALL SGTSVX( 'n', 'N', 2, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
326 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
327 $ ip, b, 1, x, 2, rcond, r1, r2, w, iw, info )
328 CALL chkxer( 'SGTSVX', infot, nout, lerr, ok )
329 infot = 16
330 CALL sgtsvx( 'N', 'N', 2, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
331 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
332 $ ip, b, 2, x, 1, rcond, r1, r2, w, iw, info )
333 CALL chkxer( 'SGTSVX', infot, nout, lerr, ok )
334*
335 ELSE IF( lsamen( 2, c2, 'PO' ) ) THEN
336*
337* SPOSV
338*
339 srnamt = 'SPOSV '
340 infot = 1
341 CALL sposv( '/', 0, 0, a, 1, b, 1, info )
342 CALL chkxer( 'SPOSV ', infot, nout, lerr, ok )
343 infot = 2
344 CALL sposv( 'U', -1, 0, a, 1, b, 1, info )
345 CALL chkxer( 'SPOSV ', infot, nout, lerr, ok )
346 infot = 3
347 CALL sposv( 'U', 0, -1, a, 1, b, 1, info )
348 CALL chkxer( 'SPOSV ', infot, nout, lerr, ok )
349 infot = 5
350 CALL sposv( 'U', 2, 0, a, 1, b, 2, info )
351 CALL chkxer( 'SPOSV ', infot, nout, lerr, ok )
352 infot = 7
353 CALL sposv( 'U', 2, 0, a, 2, b, 1, info )
354 CALL chkxer( 'SPOSV ', infot, nout, lerr, ok )
355*
356* SPOSVX
357*
358 srnamt = 'SPOSVX'
359 infot = 1
360 CALL sposvx( '/', 'U', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
361 $ rcond, r1, r2, w, iw, info )
362 CALL chkxer( 'sposvx', INFOT, NOUT, LERR, OK )
363 INFOT = 2
364 CALL SPOSVX( 'n', '/', 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
365 $ RCOND, R1, R2, W, IW, INFO )
366 CALL CHKXER( 'sposvx', INFOT, NOUT, LERR, OK )
367 INFOT = 3
368 CALL SPOSVX( 'n', 'u', -1, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
369 $ RCOND, R1, R2, W, IW, INFO )
370 CALL CHKXER( 'sposvx', INFOT, NOUT, LERR, OK )
371 INFOT = 4
372 CALL SPOSVX( 'n', 'u', 0, -1, A, 1, AF, 1, EQ, C, B, 1, X, 1,
373 $ RCOND, R1, R2, W, IW, INFO )
374 CALL CHKXER( 'sposvx', INFOT, NOUT, LERR, OK )
375 INFOT = 6
376 CALL SPOSVX( 'n', 'u', 2, 0, A, 1, AF, 2, EQ, C, B, 2, X, 2,
377 $ RCOND, R1, R2, W, IW, INFO )
378 CALL CHKXER( 'sposvx', INFOT, NOUT, LERR, OK )
379 INFOT = 8
380 CALL SPOSVX( 'n', 'u', 2, 0, A, 2, AF, 1, EQ, C, B, 2, X, 2,
381 $ RCOND, R1, R2, W, IW, INFO )
382 CALL CHKXER( 'sposvx', INFOT, NOUT, LERR, OK )
383 INFOT = 9
384 EQ = '/'
385 CALL SPOSVX( 'f', 'u', 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
386 $ RCOND, R1, R2, W, IW, INFO )
387 CALL CHKXER( 'sposvx', INFOT, NOUT, LERR, OK )
388 INFOT = 10
389 EQ = 'y'
390 CALL SPOSVX( 'f', 'u', 1, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
391 $ RCOND, R1, R2, W, IW, INFO )
392 CALL CHKXER( 'sposvx', INFOT, NOUT, LERR, OK )
393 INFOT = 12
394 CALL SPOSVX( 'n', 'u', 2, 0, A, 2, AF, 2, EQ, C, B, 1, X, 2,
395 $ RCOND, R1, R2, W, IW, INFO )
396 CALL CHKXER( 'sposvx', INFOT, NOUT, LERR, OK )
397 INFOT = 14
398 CALL SPOSVX( 'n', 'u', 2, 0, A, 2, AF, 2, EQ, C, B, 2, X, 1,
399 $ RCOND, R1, R2, W, IW, INFO )
400 CALL CHKXER( 'sposvx', INFOT, NOUT, LERR, OK )
401*
402 ELSE IF( LSAMEN( 2, C2, 'pp' ) ) THEN
403*
404* SPPSV
405*
406 SRNAMT = 'sppsv '
407 infot = 1
408 CALL sppsv( '/', 0, 0, a, b, 1, info )
409 CALL chkxer( 'SPPSV ', infot, nout, lerr, ok )
410 infot = 2
411 CALL sppsv( 'U', -1, 0, a, b, 1, info )
412 CALL chkxer( 'SPPSV ', infot, nout, lerr, ok )
413 infot = 3
414 CALL sppsv( 'U', 0, -1, a, b, 1, info )
415 CALL chkxer( 'SPPSV ', infot, nout, lerr, ok )
416 infot = 6
417 CALL sppsv( 'U', 2, 0, a, b, 1, info )
418 CALL chkxer( 'SPPSV ', infot, nout, lerr, ok )
419*
420* SPPSVX
421*
422 srnamt = 'SPPSVX'
423 infot = 1
424 CALL sppsvx( '/', 'U', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
425 $ r1, r2, w, iw, info )
426 CALL chkxer( 'SPPSVX', infot, nout, lerr, ok )
427 infot = 2
428 CALL sppsvx( 'N', '/', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
429 $ r1, r2, w, iw, info )
430 CALL chkxer( 'SPPSVX', infot, nout, lerr, ok )
431 infot = 3
432 CALL sppsvx( 'N', 'U', -1, 0, a, af, eq, c, b, 1, x, 1, rcond,
433 $ r1, r2, w, iw, info )
434 CALL chkxer( 'SPPSVX', infot, nout, lerr, ok )
435 infot = 4
436 CALL sppsvx( 'N', 'U', 0, -1, a, af, eq, c, b, 1, x, 1, rcond,
437 $ r1, r2, w, iw, info )
438 CALL chkxer( 'SPPSVX', infot, nout, lerr, ok )
439 infot = 7
440 eq = '/'
441 CALL sppsvx( 'F', 'U', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
442 $ r1, r2, w, iw, info )
443 CALL chkxer( 'SPPSVX', infot, nout, lerr, ok )
444 infot = 8
445 eq = 'Y'
446 CALL sppsvx( 'F', 'U', 1, 0, a, af, eq, c, b, 1, x, 1, rcond,
447 $ r1, r2, w, iw, info )
448 CALL chkxer( 'SPPSVX', infot, nout, lerr, ok )
449 infot = 10
450 CALL sppsvx( 'N', 'U', 2, 0, a, af, eq, c, b, 1, x, 2, rcond,
451 $ r1, r2, w, iw, info )
452 CALL chkxer( 'SPPSVX', infot, nout, lerr, ok )
453 infot = 12
454 CALL sppsvx( 'N', 'U', 2, 0, a, af, eq, c, b, 2, x, 1, rcond,
455 $ r1, r2, w, iw, info )
456 CALL chkxer( 'SPPSVX', infot, nout, lerr, ok )
457*
458 ELSE IF( lsamen( 2, c2, 'PB' ) ) THEN
459*
460* SPBSV
461*
462 srnamt = 'SPBSV '
463 infot = 1
464 CALL spbsv( '/', 0, 0, 0, a, 1, b, 1, info )
465 CALL chkxer( 'SPBSV ', infot, nout, lerr, ok )
466 infot = 2
467 CALL spbsv( 'U', -1, 0, 0, a, 1, b, 1, info )
468 CALL chkxer( 'SPBSV ', infot, nout, lerr, ok )
469 infot = 3
470 CALL spbsv( 'U', 1, -1, 0, a, 1, b, 1, info )
471 CALL chkxer( 'SPBSV ', infot, nout, lerr, ok )
472 infot = 4
473 CALL spbsv( 'U', 0, 0, -1, a, 1, b, 1, info )
474 CALL chkxer( 'SPBSV ', infot, nout, lerr, ok )
475 infot = 6
476 CALL spbsv( 'U', 1, 1, 0, a, 1, b, 2, info )
477 CALL chkxer( 'SPBSV ', infot, nout, lerr, ok )
478 infot = 8
479 CALL spbsv( 'U', 2, 0, 0, a, 1, b, 1, info )
480 CALL chkxer( 'SPBSV ', infot, nout, lerr, ok )
481*
482* SPBSVX
483*
484 srnamt = 'SPBSVX'
485 infot = 1
486 CALL spbsvx( '/', 'U', 0, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
487 $ rcond, r1, r2, w, iw, info )
488 CALL chkxer( 'SPBSVX', infot, nout, lerr, ok )
489 infot = 2
490 CALL spbsvx( 'N', '/', 0, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
491 $ rcond, r1, r2, w, iw, info )
492 CALL chkxer( 'SPBSVX', infot, nout, lerr, ok )
493 infot = 3
494 CALL spbsvx( 'N', 'U', -1, 0, 0, a, 1, af, 1, eq, c, b, 1, x,
495 $ 1, rcond, r1, r2, w, iw, info )
496 CALL chkxer( 'SPBSVX', infot, nout, lerr, ok )
497 infot = 4
498 CALL spbsvx( 'N', 'U', 1, -1, 0, a, 1, af, 1, eq, c, b, 1, x,
499 $ 1, rcond, r1, r2, w, iw, info )
500 CALL chkxer( 'SPBSVX', infot, nout, lerr, ok )
501 infot = 5
502 CALL spbsvx( 'N', 'U', 0, 0, -1, a, 1, af, 1, eq, c, b, 1, x,
503 $ 1, rcond, r1, r2, w, iw, info )
504 CALL chkxer( 'SPBSVX', infot, nout, lerr, ok )
505 infot = 7
506 CALL spbsvx( 'N', 'U', 1, 1, 0, a, 1, af, 2, eq, c, b, 2, x, 2,
507 $ rcond, r1, r2, w, iw, info )
508 CALL chkxer( 'SPBSVX', infot, nout, lerr, ok )
509 infot = 9
510 CALL spbsvx( 'N', 'U', 1, 1, 0, a, 2, af, 1, eq, c, b, 2, x, 2,
511 $ rcond, r1, r2, w, iw, info )
512 CALL chkxer( 'SPBSVX', infot, nout, lerr, ok )
513 infot = 10
514 eq = '/'
515 CALL spbsvx( 'F', 'U', 0, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
516 $ rcond, r1, r2, w, iw, info )
517 CALL chkxer( 'SPBSVX', infot, nout, lerr, ok )
518 infot = 11
519 eq = 'Y'
520 CALL spbsvx( 'F', 'U', 1, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
521 $ rcond, r1, r2, w, iw, info )
522 CALL chkxer( 'SPBSVX', infot, nout, lerr, ok )
523 infot = 13
524 CALL spbsvx( 'N', 'U', 2, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 2,
525 $ rcond, r1, r2, w, iw, info )
526 CALL chkxer( 'SPBSVX', infot, nout, lerr, ok )
527 infot = 15
528 CALL spbsvx( 'N', 'U', 2, 0, 0, a, 1, af, 1, eq, c, b, 2, x, 1,
529 $ rcond, r1, r2, w, iw, info )
530 CALL chkxer( 'SPBSVX', infot, nout, lerr, ok )
531*
532 ELSE IF( lsamen( 2, c2, 'PT' ) ) THEN
533*
534* SPTSV
535*
536 srnamt = 'SPTSV '
537 infot = 1
538 CALL sptsv( -1, 0, a( 1, 1 ), a( 1, 2 ), b, 1, info )
539 CALL chkxer( 'SPTSV ', infot, nout, lerr, ok )
540 infot = 2
541 CALL sptsv( 0, -1, a( 1, 1 ), a( 1, 2 ), b, 1, info )
542 CALL chkxer( 'SPTSV ', infot, nout, lerr, ok )
543 infot = 6
544 CALL sptsv( 2, 0, a( 1, 1 ), a( 1, 2 ), b, 1, info )
545 CALL chkxer( 'SPTSV ', infot, nout, lerr, ok )
546*
547* SPTSVX
548*
549 srnamt = 'SPTSVX'
550 infot = 1
551 CALL sptsvx( '/', 0, 0, a( 1, 1 ), a( 1, 2 ), af( 1, 1 ),
552 $ af( 1, 2 ), b, 1, x, 1, rcond, r1, r2, w, info )
553 CALL chkxer( 'SPTSVX', infot, nout, lerr, ok )
554 infot = 2
555 CALL sptsvx( 'N', -1, 0, a( 1, 1 ), a( 1, 2 ), af( 1, 1 ),
556 $ af( 1, 2 ), b, 1, x, 1, rcond, r1, r2, w, info )
557 CALL chkxer( 'SPTSVX', infot, nout, lerr, ok )
558 infot = 3
559 CALL sptsvx( 'N', 0, -1, a( 1, 1 ), a( 1, 2 ), af( 1, 1 ),
560 $ af( 1, 2 ), b, 1, x, 1, rcond, r1, r2, w, info )
561 CALL chkxer( 'SPTSVX', infot, nout, lerr, ok )
562 infot = 9
563 CALL sptsvx( 'N', 2, 0, a( 1, 1 ), a( 1, 2 ), af( 1, 1 ),
564 $ af( 1, 2 ), b, 1, x, 2, rcond, r1, r2, w, info )
565 CALL chkxer( 'SPTSVX', infot, nout, lerr, ok )
566 infot = 11
567 CALL sptsvx( 'N', 2, 0, a( 1, 1 ), a( 1, 2 ), af( 1, 1 ),
568 $ af( 1, 2 ), b, 2, x, 1, rcond, r1, r2, w, info )
569 CALL chkxer( 'SPTSVX', infot, nout, lerr, ok )
570*
571 ELSE IF( lsamen( 2, c2, 'SY' ) ) THEN
572*
573* SSYSV
574*
575 srnamt = 'SSYSV '
576 infot = 1
577 CALL ssysv( '/', 0, 0, a, 1, ip, b, 1, w, 1, info )
578 CALL chkxer( 'SSYSV ', infot, nout, lerr, ok )
579 infot = 2
580 CALL ssysv( 'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
581 CALL chkxer( 'SSYSV ', infot, nout, lerr, ok )
582 infot = 3
583 CALL ssysv( 'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
584 CALL chkxer( 'SSYSV ', infot, nout, lerr, ok )
585 infot = 8
586 CALL ssysv( 'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
587 CALL chkxer( 'SSYSV ', infot, nout, lerr, ok )
588 infot = 10
589 CALL ssysv( 'U', 0, 0, a, 1, ip, b, 1, w, 0, info )
590 CALL chkxer( 'SSYSV ', infot, nout, lerr, ok )
591 infot = 10
592 CALL ssysv( 'U', 0, 0, a, 1, ip, b, 1, w, -2, info )
593 CALL chkxer( 'SSYSV ', infot, nout, lerr, ok )
594*
595* SSYSVX
596*
597 srnamt = 'SSYSVX'
598 infot = 1
599 CALL ssysvx( '/', 'U', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
600 $ rcond, r1, r2, w, 1, iw, info )
601 CALL chkxer( 'SSYSVX', infot, nout, lerr, ok )
602 infot = 2
603 CALL ssysvx( 'N', '/', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
604 $ rcond, r1, r2, w, 1, iw, info )
605 CALL chkxer( 'SSYSVX', infot, nout, lerr, ok )
606 infot = 3
607 CALL ssysvx( 'N', 'U', -1, 0, a, 1, af, 1, ip, b, 1, x, 1,
608 $ rcond, r1, r2, w, 1, iw, info )
609 CALL chkxer( 'SSYSVX', infot, nout, lerr, ok )
610 infot = 4
611 CALL ssysvx( 'N', 'U', 0, -1, a, 1, af, 1, ip, b, 1, x, 1,
612 $ rcond, r1, r2, w, 1, iw, info )
613 CALL chkxer( 'SSYSVX', infot, nout, lerr, ok )
614 infot = 6
615 CALL ssysvx( 'N', 'U', 2, 0, a, 1, af, 2, ip, b, 2, x, 2,
616 $ rcond, r1, r2, w, 4, iw, info )
617 CALL chkxer( 'SSYSVX', infot, nout, lerr, ok )
618 infot = 8
619 CALL ssysvx( 'N', 'U', 2, 0, a, 2, af, 1, ip, b, 2, x, 2,
620 $ rcond, r1, r2, w, 4, iw, info )
621 CALL chkxer( 'SSYSVX', infot, nout, lerr, ok )
622 infot = 11
623 CALL ssysvx( 'N', 'U', 2, 0, a, 2, af, 2, ip, b, 1, x, 2,
624 $ rcond, r1, r2, w, 4, iw, info )
625 CALL chkxer( 'SSYSVX', infot, nout, lerr, ok )
626 infot = 13
627 CALL ssysvx( 'N', 'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 1,
628 $ rcond, r1, r2, w, 4, iw, info )
629 CALL chkxer( 'SSYSVX', infot, nout, lerr, ok )
630 infot = 18
631 CALL ssysvx( 'N', 'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 2,
632 $ rcond, r1, r2, w, 3, iw, info )
633 CALL chkxer( 'SSYSVX', infot, nout, lerr, ok )
634*
635*
636 ELSE IF( lsamen( 2, c2, 'SR' ) ) THEN
637*
638* SSYSV_ROOK
639*
640 srnamt = 'SSYSV_ROOK'
641 infot = 1
642 CALL ssysv_rook( '/', 0, 0, a, 1, ip, b, 1, w, 1, info )
643 CALL chkxer( 'SSYSV_ROOK', infot, nout, lerr, ok )
644 infot = 2
645 CALL ssysv_rook( 'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
646 CALL chkxer( 'SSYSV_ROOK', infot, nout, lerr, ok )
647 infot = 3
648 CALL ssysv_rook( 'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
649 CALL chkxer( 'SSYSV_ROOK', infot, nout, lerr, ok )
650 infot = 8
651 CALL ssysv_rook( 'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
652 CALL chkxer( 'SSYSV_ROOK', infot, nout, lerr, ok )
653 infot = 10
654 CALL ssysv_rook( 'U', 0, 0, a, 1, ip, b, 1, w, 0, info )
655 CALL chkxer( 'SSYSV_ROOK', infot, nout, lerr, ok )
656 infot = 10
657 CALL ssysv_rook( 'U', 0, 0, a, 1, ip, b, 1, w, -2, info )
658 CALL chkxer( 'SSYSV_ROOK', infot, nout, lerr, ok )
659*
660 ELSE IF( lsamen( 2, c2, 'SK' ) ) THEN
661*
662* SSYSV_RK
663*
664* Test error exits of the driver that uses factorization
665* of a symmetric indefinite matrix with rook
666* (bounded Bunch-Kaufman) pivoting with the new storage
667* format for factors L ( or U) and D.
668*
669* L (or U) is stored in A, diagonal of D is stored on the
670* diagonal of A, subdiagonal of D is stored in a separate array E.
671*
672 srnamt = 'SSYSV_RK'
673 infot = 1
674 CALL ssysv_rk( '/', 0, 0, a, 1, e, ip, b, 1, w, 1, info )
675 CALL chkxer( 'SSYSV_RK', infot, nout, lerr, ok )
676 infot = 2
677 CALL ssysv_rk( 'U', -1, 0, a, 1, e, ip, b, 1, w, 1, info )
678 CALL chkxer( 'SSYSV_RK', infot, nout, lerr, ok )
679 infot = 3
680 CALL ssysv_rk( 'U', 0, -1, a, 1, e, ip, b, 1, w, 1, info )
681 CALL chkxer( 'SSYSV_RK', infot, nout, lerr, ok )
682 infot = 5
683 CALL ssysv_rk( 'U', 2, 0, a, 1, e, ip, b, 2, w, 1, info )
684 CALL chkxer( 'SSYSV_RK', infot, nout, lerr, ok )
685 infot = 9
686 CALL ssysv_rk( 'U', 2, 0, a, 2, e, ip, b, 1, w, 1, info )
687 CALL chkxer( 'SSYSV_RK', infot, nout, lerr, ok )
688 infot = 11
689 CALL ssysv_rk( 'U', 0, 0, a, 1, e, ip, b, 1, w, 0, info )
690 CALL chkxer( 'SSYSV_RK', infot, nout, lerr, ok )
691 infot = 11
692 CALL ssysv_rk( 'U', 0, 0, a, 1, e, ip, b, 1, w, -2, info )
693 CALL chkxer( 'SSYSV_RK', infot, nout, lerr, ok )
694*
695 ELSE IF( lsamen( 2, c2, 'SA' ) ) THEN
696*
697* SSYSV_AA
698*
699 srnamt = 'SSYSV_AA'
700 infot = 1
701 CALL ssysv_aa( '/', 0, 0, a, 1, ip, b, 1, w, 1, info )
702 CALL chkxer( 'ssysv_aa', INFOT, NOUT, LERR, OK )
703 INFOT = 2
704 CALL SSYSV_AA( 'u', -1, 0, A, 1, IP, B, 1, W, 1, INFO )
705 CALL CHKXER( 'ssysv_aa', INFOT, NOUT, LERR, OK )
706 INFOT = 3
707 CALL SSYSV_AA( 'u', 0, -1, A, 1, IP, B, 1, W, 1, INFO )
708 CALL CHKXER( 'ssysv_aa', INFOT, NOUT, LERR, OK )
709 INFOT = 8
710 CALL SSYSV_AA( 'u', 2, 0, A, 2, IP, B, 1, W, 1, INFO )
711 CALL CHKXER( 'ssysv_aa', INFOT, NOUT, LERR, OK )
712*
713 ELSE IF( LSAMEN( 2, C2, 's2' ) ) THEN
714*
715* DSYSV_AASEN_2STAGE
716*
717 SRNAMT = 'ssysv_aa_2stage'
718 INFOT = 1
719 CALL SSYSV_AA_2STAGE( '/', 0, 0, A, 1, A, 1, IP, IP, B, 1,
720 $ W, 1, INFO )
721 CALL CHKXER( 'ssysv_aa_2stage', INFOT, NOUT, LERR, OK )
722 INFOT = 2
723 CALL SSYSV_AA_2STAGE( 'u', -1, 0, A, 1, A, 1, IP, IP, B, 1,
724 $ W, 1, INFO )
725 CALL CHKXER( 'ssysv_aa_2stage', INFOT, NOUT, LERR, OK )
726 INFOT = 3
727 CALL SSYSV_AA_2STAGE( 'u', 0, -1, A, 1, A, 1, IP, IP, B, 1,
728 $ W, 1, INFO )
729 CALL CHKXER( 'ssysv_aa_2stage', INFOT, NOUT, LERR, OK )
730 INFOT = 5
731 CALL SSYSV_AA_2STAGE( 'u', 2, 1, A, 1, A, 1, IP, IP, B, 1,
732 $ W, 1, INFO )
733 CALL CHKXER( 'ssysv_aa_2stage', INFOT, NOUT, LERR, OK )
734 INFOT = 11
735 CALL SSYSV_AA_2STAGE( 'u', 2, 1, A, 2, A, 8, IP, IP, B, 1,
736 $ W, 1, INFO )
737 CALL CHKXER( 'ssysv_aa_2stage', INFOT, NOUT, LERR, OK )
738 INFOT = 7
739 CALL SSYSV_AA_2STAGE( 'u', 2, 1, A, 2, A, 1, IP, IP, B, 2,
740 $ W, 1, INFO )
741 CALL CHKXER( 'ssysv_aa_2stage', INFOT, NOUT, LERR, OK )
742*
743 ELSE IF( LSAMEN( 2, C2, 'sp' ) ) THEN
744*
745* SSPSV
746*
747 SRNAMT = 'sspsv '
748 INFOT = 1
749 CALL SSPSV( '/', 0, 0, A, IP, B, 1, INFO )
750 CALL CHKXER( 'sspsv ', INFOT, NOUT, LERR, OK )
751 INFOT = 2
752 CALL SSPSV( 'u', -1, 0, A, IP, B, 1, INFO )
753 CALL CHKXER( 'sspsv ', INFOT, NOUT, LERR, OK )
754 INFOT = 3
755 CALL SSPSV( 'u', 0, -1, a, ip, b, 1, info )
756 CALL chkxer( 'SSPSV ', infot, nout, lerr, ok )
757 infot = 7
758 CALL sspsv( 'U', 2, 0, a, ip, b, 1, info )
759 CALL chkxer( 'SSPSV ', infot, nout, lerr, ok )
760*
761* SSPSVX
762*
763 srnamt = 'SSPSVX'
764 infot = 1
765 CALL sspsvx( '/', 'U', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
766 $ r2, w, iw, info )
767 CALL chkxer( 'SSPSVX', infot, nout, lerr, ok )
768 infot = 2
769 CALL sspsvx( 'N', '/', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
770 $ r2, w, iw, info )
771 CALL chkxer( 'SSPSVX', infot, nout, lerr, ok )
772 infot = 3
773 CALL sspsvx( 'N', 'U', -1, 0, a, af, ip, b, 1, x, 1, rcond, r1,
774 $ r2, w, iw, info )
775 CALL chkxer( 'sspsvx', INFOT, NOUT, LERR, OK )
776 INFOT = 4
777 CALL SSPSVX( 'n', 'u', 0, -1, A, AF, IP, B, 1, X, 1, RCOND, R1,
778 $ R2, W, IW, INFO )
779 CALL CHKXER( 'sspsvx', INFOT, NOUT, LERR, OK )
780 INFOT = 9
781 CALL SSPSVX( 'n', 'u', 2, 0, A, AF, IP, B, 1, X, 2, RCOND, R1,
782 $ R2, W, IW, INFO )
783 CALL CHKXER( 'sspsvx', INFOT, NOUT, LERR, OK )
784 INFOT = 11
785 CALL SSPSVX( 'n', 'u', 2, 0, A, AF, IP, B, 2, X, 1, RCOND, R1,
786 $ R2, W, IW, INFO )
787 CALL CHKXER( 'sspsvx', INFOT, NOUT, LERR, OK )
788 END IF
789*
790* Print a summary line.
791*
792 IF( OK ) THEN
793 WRITE( NOUT, FMT = 9999 )PATH
794 ELSE
795 WRITE( NOUT, FMT = 9998 )PATH
796 END IF
797*
798 9999 FORMAT( 1X, A3, ' drivers passed the tests of the error exits' )
799 9998 FORMAT( ' *** ', a3, ' drivers failed the tests of the error ',
800 $ 'exits ***' )
801*
802 RETURN
803*
804* End of SERRVX
805*
806 END
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
sgbsvx
subroutine sgbsvx(fact, trans, n, kl, ku, nrhs, ab, ldab, afb, ldafb, ipiv, equed, r, c, b, ldb, x, ldx, rcond, ferr, berr, work, iwork, info)
SGBSVX computes the solution to system of linear equations A * X = B for GB matrices
Definition sgbsvx.f:368
sgbsv
subroutine sgbsv(n, kl, ku, nrhs, ab, ldab, ipiv, b, ldb, info)
SGBSV computes the solution to system of linear equations A * X = B for GB matrices (simple driver)
Definition sgbsv.f:162
sgesv
subroutine sgesv(n, nrhs, a, lda, ipiv, b, ldb, info)
SGESV computes the solution to system of linear equations A * X = B for GE matrices (simple driver)
Definition sgesv.f:122
sgesvx
subroutine sgesvx(fact, trans, n, nrhs, a, lda, af, ldaf, ipiv, equed, r, c, b, ldb, x, ldx, rcond, ferr, berr, work, iwork, info)
SGESVX computes the solution to system of linear equations A * X = B for GE matrices
Definition sgesvx.f:349
sgtsv
subroutine sgtsv(n, nrhs, dl, d, du, b, ldb, info)
SGTSV computes the solution to system of linear equations A * X = B for GT matrices
Definition sgtsv.f:127
sgtsvx
subroutine sgtsvx(fact, trans, n, nrhs, dl, d, du, dlf, df, duf, du2, ipiv, b, ldb, x, ldx, rcond, ferr, berr, work, iwork, info)
SGTSVX computes the solution to system of linear equations A * X = B for GT matrices
Definition sgtsvx.f:293
sppsvx
subroutine sppsvx(fact, uplo, n, nrhs, ap, afp, equed, s, b, ldb, x, ldx, rcond, ferr, berr, work, iwork, info)
SPPSVX computes the solution to system of linear equations A * X = B for OTHER matrices
Definition sppsvx.f:311
sspsvx
subroutine sspsvx(fact, uplo, n, nrhs, ap, afp, ipiv, b, ldb, x, ldx, rcond, ferr, berr, work, iwork, info)
SSPSVX computes the solution to system of linear equations A * X = B for OTHER matrices
Definition sspsvx.f:276
spbsv
subroutine spbsv(uplo, n, kd, nrhs, ab, ldab, b, ldb, info)
SPBSV computes the solution to system of linear equations A * X = B for OTHER matrices
Definition spbsv.f:164
sppsv
subroutine sppsv(uplo, n, nrhs, ap, b, ldb, info)
SPPSV computes the solution to system of linear equations A * X = B for OTHER matrices
Definition sppsv.f:144
spbsvx
subroutine spbsvx(fact, uplo, n, kd, nrhs, ab, ldab, afb, ldafb, equed, s, b, ldb, x, ldx, rcond, ferr, berr, work, iwork, info)
SPBSVX computes the solution to system of linear equations A * X = B for OTHER matrices
Definition spbsvx.f:343
sspsv
subroutine sspsv(uplo, n, nrhs, ap, ipiv, b, ldb, info)
SSPSV computes the solution to system of linear equations A * X = B for OTHER matrices
Definition sspsv.f:162
sposv
subroutine sposv(uplo, n, nrhs, a, lda, b, ldb, info)
SPOSV computes the solution to system of linear equations A * X = B for PO matrices
Definition sposv.f:130
sposvx
subroutine sposvx(fact, uplo, n, nrhs, a, lda, af, ldaf, equed, s, b, ldb, x, ldx, rcond, ferr, berr, work, iwork, info)
SPOSVX computes the solution to system of linear equations A * X = B for PO matrices
Definition sposvx.f:307
sptsv
subroutine sptsv(n, nrhs, d, e, b, ldb, info)
SPTSV computes the solution to system of linear equations A * X = B for PT matrices
Definition sptsv.f:114
sptsvx
subroutine sptsvx(fact, n, nrhs, d, e, df, ef, b, ldb, x, ldx, rcond, ferr, berr, work, info)
SPTSVX computes the solution to system of linear equations A * X = B for PT matrices
Definition sptsvx.f:228
ssysv_aa_2stage
subroutine ssysv_aa_2stage(uplo, n, nrhs, a, lda, tb, ltb, ipiv, ipiv2, b, ldb, work, lwork, info)
SSYSV_AA_2STAGE computes the solution to system of linear equations A * X = B for SY matrices
Definition ssysv_aa_2stage.f:187
ssysv
subroutine ssysv(uplo, n, nrhs, a, lda, ipiv, b, ldb, work, lwork, info)
SSYSV computes the solution to system of linear equations A * X = B for SY matrices
Definition ssysv.f:171
ssysvx
subroutine ssysvx(fact, uplo, n, nrhs, a, lda, af, ldaf, ipiv, b, ldb, x, ldx, rcond, ferr, berr, work, lwork, iwork, info)
SSYSVX computes the solution to system of linear equations A * X = B for SY matrices
Definition ssysvx.f:284
ssysv_rook
subroutine ssysv_rook(uplo, n, nrhs, a, lda, ipiv, b, ldb, work, lwork, info)
SSYSV_ROOK computes the solution to system of linear equations A * X = B for SY matrices
Definition ssysv_rook.f:204
ssysv_aa
subroutine ssysv_aa(uplo, n, nrhs, a, lda, ipiv, b, ldb, work, lwork, info)
SSYSV_AA computes the solution to system of linear equations A * X = B for SY matrices
Definition ssysv_aa.f:162
serrvx
subroutine serrvx(path, nunit)
SERRVX
Definition serrvx.f:55
ssysv_rk
subroutine ssysv_rk(uplo, n, nrhs, a, lda, e, ipiv, b, ldb, work, lwork, info)
SSYSV_RK computes the solution to system of linear equations A * X = B for SY matrices
Definition ssysv_rk.f:228