LAPACK 3.12.0
LAPACK: Linear Algebra PACKage
Loading...
Searching...
No Matches

◆ serrvx()

subroutine serrvx ( character*3  path,
integer  nunit 
)

SERRVX

Purpose:
 SERRVX tests the error exits for the REAL driver routines
 for solving linear systems of equations.
Parameters
[in]PATH
          PATH is CHARACTER*3
          The LAPACK path name for the routines to be tested.
[in]NUNIT
          NUNIT is INTEGER
          The unit number for output.
Author
Univ. of Tennessee
Univ. of California Berkeley
Univ. of Colorado Denver
NAG Ltd.

Definition at line 54 of file serrvx.f.

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,
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 = 5
586 CALL ssysv( 'U', 2, 0, a, 1, ip, b, 2, w, 1, info )
587 CALL chkxer( 'SSYSV ', infot, nout, lerr, ok )
588 infot = 8
589 CALL ssysv( 'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
590 CALL chkxer( 'SSYSV ', infot, nout, lerr, ok )
591 infot = 10
592 CALL ssysv( 'U', 0, 0, a, 1, ip, b, 1, w, 0, info )
593 CALL chkxer( 'SSYSV ', infot, nout, lerr, ok )
594 infot = 10
595 CALL ssysv( 'U', 0, 0, a, 1, ip, b, 1, w, -2, info )
596 CALL chkxer( 'SSYSV ', infot, nout, lerr, ok )
597*
598* SSYSVX
599*
600 srnamt = 'SSYSVX'
601 infot = 1
602 CALL ssysvx( '/', 'U', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
603 $ rcond, r1, r2, w, 1, iw, info )
604 CALL chkxer( 'SSYSVX', infot, nout, lerr, ok )
605 infot = 2
606 CALL ssysvx( 'N', '/', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
607 $ rcond, r1, r2, w, 1, iw, info )
608 CALL chkxer( 'SSYSVX', infot, nout, lerr, ok )
609 infot = 3
610 CALL ssysvx( 'N', 'U', -1, 0, a, 1, af, 1, ip, b, 1, x, 1,
611 $ rcond, r1, r2, w, 1, iw, info )
612 CALL chkxer( 'SSYSVX', infot, nout, lerr, ok )
613 infot = 4
614 CALL ssysvx( 'N', 'U', 0, -1, a, 1, af, 1, ip, b, 1, x, 1,
615 $ rcond, r1, r2, w, 1, iw, info )
616 CALL chkxer( 'SSYSVX', infot, nout, lerr, ok )
617 infot = 6
618 CALL ssysvx( 'N', 'U', 2, 0, a, 1, af, 2, ip, b, 2, x, 2,
619 $ rcond, r1, r2, w, 4, iw, info )
620 CALL chkxer( 'SSYSVX', infot, nout, lerr, ok )
621 infot = 8
622 CALL ssysvx( 'N', 'U', 2, 0, a, 2, af, 1, ip, b, 2, x, 2,
623 $ rcond, r1, r2, w, 4, iw, info )
624 CALL chkxer( 'SSYSVX', infot, nout, lerr, ok )
625 infot = 11
626 CALL ssysvx( 'N', 'U', 2, 0, a, 2, af, 2, ip, b, 1, x, 2,
627 $ rcond, r1, r2, w, 4, iw, info )
628 CALL chkxer( 'SSYSVX', infot, nout, lerr, ok )
629 infot = 13
630 CALL ssysvx( 'N', 'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 1,
631 $ rcond, r1, r2, w, 4, iw, info )
632 CALL chkxer( 'SSYSVX', infot, nout, lerr, ok )
633 infot = 18
634 CALL ssysvx( 'N', 'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 2,
635 $ rcond, r1, r2, w, 3, iw, info )
636 CALL chkxer( 'SSYSVX', infot, nout, lerr, ok )
637*
638*
639 ELSE IF( lsamen( 2, c2, 'SR' ) ) THEN
640*
641* SSYSV_ROOK
642*
643 srnamt = 'SSYSV_ROOK'
644 infot = 1
645 CALL ssysv_rook( '/', 0, 0, a, 1, ip, b, 1, w, 1, info )
646 CALL chkxer( 'SSYSV_ROOK', infot, nout, lerr, ok )
647 infot = 2
648 CALL ssysv_rook( 'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
649 CALL chkxer( 'SSYSV_ROOK', infot, nout, lerr, ok )
650 infot = 3
651 CALL ssysv_rook( 'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
652 CALL chkxer( 'SSYSV_ROOK', infot, nout, lerr, ok )
653 infot = 5
654 CALL ssysv_rook( 'U', 2, 0, a, 1, ip, b, 2, w, 1, info )
655 CALL chkxer( 'SSYSV_ROOK', infot, nout, lerr, ok )
656 infot = 8
657 CALL ssysv_rook( 'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
658 CALL chkxer( 'SSYSV_ROOK', infot, nout, lerr, ok )
659 infot = 10
660 CALL ssysv_rook( 'U', 0, 0, a, 1, ip, b, 1, w, 0, info )
661 CALL chkxer( 'SSYSV_ROOK', infot, nout, lerr, ok )
662 infot = 10
663 CALL ssysv_rook( 'U', 0, 0, a, 1, ip, b, 1, w, -2, info )
664 CALL chkxer( 'SSYSV_ROOK', infot, nout, lerr, ok )
665*
666 ELSE IF( lsamen( 2, c2, 'SK' ) ) THEN
667*
668* SSYSV_RK
669*
670* Test error exits of the driver that uses factorization
671* of a symmetric indefinite matrix with rook
672* (bounded Bunch-Kaufman) pivoting with the new storage
673* format for factors L ( or U) and D.
674*
675* L (or U) is stored in A, diagonal of D is stored on the
676* diagonal of A, subdiagonal of D is stored in a separate array E.
677*
678 srnamt = 'SSYSV_RK'
679 infot = 1
680 CALL ssysv_rk( '/', 0, 0, a, 1, e, ip, b, 1, w, 1, info )
681 CALL chkxer( 'SSYSV_RK', infot, nout, lerr, ok )
682 infot = 2
683 CALL ssysv_rk( 'U', -1, 0, a, 1, e, ip, b, 1, w, 1, info )
684 CALL chkxer( 'SSYSV_RK', infot, nout, lerr, ok )
685 infot = 3
686 CALL ssysv_rk( 'U', 0, -1, a, 1, e, ip, b, 1, w, 1, info )
687 CALL chkxer( 'SSYSV_RK', infot, nout, lerr, ok )
688 infot = 5
689 CALL ssysv_rk( 'U', 2, 0, a, 1, e, ip, b, 2, w, 1, info )
690 CALL chkxer( 'SSYSV_RK', infot, nout, lerr, ok )
691 infot = 9
692 CALL ssysv_rk( 'U', 2, 0, a, 2, e, ip, b, 1, w, 1, info )
693 CALL chkxer( 'SSYSV_RK', infot, nout, lerr, ok )
694 infot = 11
695 CALL ssysv_rk( 'U', 0, 0, a, 1, e, ip, b, 1, w, 0, info )
696 CALL chkxer( 'SSYSV_RK', infot, nout, lerr, ok )
697 infot = 11
698 CALL ssysv_rk( 'U', 0, 0, a, 1, e, ip, b, 1, w, -2, info )
699 CALL chkxer( 'SSYSV_RK', infot, nout, lerr, ok )
700*
701 ELSE IF( lsamen( 2, c2, 'SA' ) ) THEN
702*
703* SSYSV_AASEN
704*
705 srnamt = 'SSYSV_AA'
706 infot = 1
707 CALL ssysv_aa( '/', 0, 0, a, 1, ip, b, 1, w, 1, info )
708 CALL chkxer( 'SSYSV_AA', infot, nout, lerr, ok )
709 infot = 2
710 CALL ssysv_aa( 'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
711 CALL chkxer( 'SSYSV_AA', infot, nout, lerr, ok )
712 infot = 3
713 CALL ssysv_aa( 'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
714 CALL chkxer( 'SSYSV_AA', infot, nout, lerr, ok )
715 infot = 5
716 CALL ssysv_aa( 'U', 2, 0, a, 1, ip, b, 2, w, 1, info )
717 CALL chkxer( 'SSYSV_AA', infot, nout, lerr, ok )
718 infot = 8
719 CALL ssysv_aa( 'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
720 CALL chkxer( 'SSYSV_AA', infot, nout, lerr, ok )
721 infot = 10
722 CALL ssysv_aa( 'U', 3, 1, a, 3, ip, b, 3, w, 6, info )
723 CALL chkxer( 'SSYSV_AA', infot, nout, lerr, ok )
724*
725 ELSE IF( lsamen( 2, c2, 'S2' ) ) THEN
726*
727* SSYSV_AASEN_2STAGE
728*
729 srnamt = 'SSYSV_AA_2STAGE'
730 infot = 1
731 CALL ssysv_aa_2stage( '/', 0, 0, a, 1, a, 1, ip, ip, b, 1,
732 $ w, 1, info )
733 CALL chkxer( 'SSYSV_AA_2STAGE', infot, nout, lerr, ok )
734 infot = 2
735 CALL ssysv_aa_2stage( 'U', -1, 0, a, 1, a, 1, ip, ip, b, 1,
736 $ w, 1, info )
737 CALL chkxer( 'SSYSV_AA_2STAGE', infot, nout, lerr, ok )
738 infot = 3
739 CALL ssysv_aa_2stage( 'U', 0, -1, a, 1, a, 1, ip, ip, b, 1,
740 $ w, 1, info )
741 CALL chkxer( 'SSYSV_AA_2STAGE', infot, nout, lerr, ok )
742 infot = 5
743 CALL ssysv_aa_2stage( 'U', 2, 1, a, 1, a, 1, ip, ip, b, 1,
744 $ w, 1, info )
745 CALL chkxer( 'SSYSV_AA_2STAGE', infot, nout, lerr, ok )
746 infot = 7
747 CALL ssysv_aa_2stage( 'U', 2, 1, a, 2, a, 1, ip, ip, b, 2,
748 $ w, 1, info )
749 CALL chkxer( 'SSYSV_AA_2STAGE', infot, nout, lerr, ok )
750 infot = 11
751 CALL ssysv_aa_2stage( 'U', 2, 1, a, 2, a, 8, ip, ip, b, 1,
752 $ w, 1, info )
753 CALL chkxer( 'SSYSV_AA_2STAGE', infot, nout, lerr, ok )
754 infot = 13
755 CALL ssysv_aa_2stage( 'U', 2, 1, a, 2, a, 8, ip, ip, b, 2,
756 $ w, 1, info )
757 CALL chkxer( 'SSYSV_AA_2STAGE', infot, nout, lerr, ok )
758*
759 ELSE IF( lsamen( 2, c2, 'SP' ) ) THEN
760*
761* SSPSV
762*
763 srnamt = 'SSPSV '
764 infot = 1
765 CALL sspsv( '/', 0, 0, a, ip, b, 1, info )
766 CALL chkxer( 'SSPSV ', infot, nout, lerr, ok )
767 infot = 2
768 CALL sspsv( 'U', -1, 0, a, ip, b, 1, info )
769 CALL chkxer( 'SSPSV ', infot, nout, lerr, ok )
770 infot = 3
771 CALL sspsv( 'U', 0, -1, a, ip, b, 1, info )
772 CALL chkxer( 'SSPSV ', infot, nout, lerr, ok )
773 infot = 7
774 CALL sspsv( 'U', 2, 0, a, ip, b, 1, info )
775 CALL chkxer( 'SSPSV ', infot, nout, lerr, ok )
776*
777* SSPSVX
778*
779 srnamt = 'SSPSVX'
780 infot = 1
781 CALL sspsvx( '/', 'U', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
782 $ r2, w, iw, info )
783 CALL chkxer( 'SSPSVX', infot, nout, lerr, ok )
784 infot = 2
785 CALL sspsvx( 'N', '/', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
786 $ r2, w, iw, info )
787 CALL chkxer( 'SSPSVX', infot, nout, lerr, ok )
788 infot = 3
789 CALL sspsvx( 'N', 'U', -1, 0, a, af, ip, b, 1, x, 1, rcond, r1,
790 $ r2, w, iw, info )
791 CALL chkxer( 'SSPSVX', infot, nout, lerr, ok )
792 infot = 4
793 CALL sspsvx( 'N', 'U', 0, -1, a, af, ip, b, 1, x, 1, rcond, r1,
794 $ r2, w, iw, info )
795 CALL chkxer( 'SSPSVX', infot, nout, lerr, ok )
796 infot = 9
797 CALL sspsvx( 'N', 'U', 2, 0, a, af, ip, b, 1, x, 2, rcond, r1,
798 $ r2, w, iw, info )
799 CALL chkxer( 'SSPSVX', infot, nout, lerr, ok )
800 infot = 11
801 CALL sspsvx( 'N', 'U', 2, 0, a, af, ip, b, 2, x, 1, rcond, r1,
802 $ r2, w, iw, info )
803 CALL chkxer( 'SSPSVX', infot, nout, lerr, ok )
804 END IF
805*
806* Print a summary line.
807*
808 IF( ok ) THEN
809 WRITE( nout, fmt = 9999 )path
810 ELSE
811 WRITE( nout, fmt = 9998 )path
812 END IF
813*
814 9999 FORMAT( 1x, a3, ' drivers passed the tests of the error exits' )
815 9998 FORMAT( ' *** ', a3, ' drivers failed the tests of the error ',
816 $ 'exits ***' )
817*
818 RETURN
819*
820* End of SERRVX
821*
subroutine chkxer(srnamt, infot, nout, lerr, ok)
Definition cblat2.f:3224
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
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
subroutine sgesv(n, nrhs, a, lda, ipiv, b, ldb, info)
Download SGESV + dependencies <a href="http://www.netlib.org/cgi-bin/netlibfiles....
Definition sgesv.f:124
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
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
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
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
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
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
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
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
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
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
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
logical function lsamen(n, ca, cb)
LSAMEN
Definition lsamen.f:74
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
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
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
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
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
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
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
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
Here is the call graph for this function:
Here is the caller graph for this function: