77 INTEGER i, info, j, n_err_bnds, nparams
78 DOUBLE PRECISION rcond, rpvgrw, berr
81 INTEGER ip( nmax ), iw( nmax )
82 DOUBLE PRECISION a( nmax, nmax ), af( nmax, nmax ), b( nmax ),
83 $ c( nmax ), r( nmax ), r1( nmax ), r2( nmax ),
84 $ w( 2*nmax ), x( nmax ), err_bnds_n( nmax, 3 ),
85 $ err_bnds_c( nmax, 3 ), params( 1 )
103 COMMON / infoc / infot, nout, ok, lerr
104 COMMON / srnamc / srnamt
112 WRITE( nout, fmt = * )
119 a( i, j ) = 1.d0 / dble( i+j )
120 af( i, j ) = 1.d0 / dble( i+j )
134 IF(
lsamen( 2, c2,
'GE' ) )
THEN
140 CALL dgesv( -1, 0, a, 1, ip, b, 1, info )
141 CALL chkxer(
'DGESV ', infot, nout, lerr, ok )
143 CALL dgesv( 0, -1, a, 1, ip, b, 1, info )
144 CALL chkxer(
'DGESV ', infot, nout, lerr, ok )
146 CALL dgesv( 2, 1, a, 1, ip, b, 2, info )
147 CALL chkxer(
'DGESV ', infot, nout, lerr, ok )
149 CALL dgesv( 2, 1, a, 2, ip, b, 1, info )
150 CALL chkxer(
'DGESV ', infot, nout, lerr, ok )
156 CALL dgesvx(
'/',
'N', 0, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
157 $ x, 1, rcond, r1, r2, w, iw, info )
158 CALL chkxer(
'DGESVX', infot, nout, lerr, ok )
160 CALL dgesvx(
'N',
'/', 0, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
161 $ x, 1, rcond, r1, r2, w, iw, info )
162 CALL chkxer(
'DGESVX', infot, nout, lerr, ok )
164 CALL dgesvx(
'N',
'N', -1, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
165 $ x, 1, rcond, r1, r2, w, iw, info )
166 CALL chkxer(
'DGESVX', infot, nout, lerr, ok )
168 CALL dgesvx(
'N',
'N', 0, -1, a, 1, af, 1, ip, eq, r, c, b, 1,
169 $ x, 1, rcond, r1, r2, w, iw, info )
170 CALL chkxer(
'DGESVX', infot, nout, lerr, ok )
172 CALL dgesvx(
'N',
'N', 2, 1, a, 1, af, 2, ip, eq, r, c, b, 2,
173 $ x, 2, rcond, r1, r2, w, iw, info )
174 CALL chkxer(
'DGESVX', infot, nout, lerr, ok )
176 CALL dgesvx(
'N',
'N', 2, 1, a, 2, af, 1, ip, eq, r, c, b, 2,
177 $ x, 2, rcond, r1, r2, w, iw, info )
178 CALL chkxer(
'DGESVX', infot, nout, lerr, ok )
181 CALL dgesvx(
'F',
'N', 0, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
182 $ x, 1, rcond, r1, r2, w, iw, info )
183 CALL chkxer(
'DGESVX', infot, nout, lerr, ok )
186 CALL dgesvx(
'F',
'N', 1, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
187 $ x, 1, rcond, r1, r2, w, iw, info )
188 CALL chkxer(
'DGESVX', infot, nout, lerr, ok )
191 CALL dgesvx(
'F',
'N', 1, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
192 $ x, 1, rcond, r1, r2, w, iw, info )
193 CALL chkxer(
'DGESVX', infot, nout, lerr, ok )
195 CALL dgesvx(
'N',
'N', 2, 1, a, 2, af, 2, ip, eq, r, c, b, 1,
196 $ x, 2, rcond, r1, r2, w, iw, info )
197 CALL chkxer(
'DGESVX', infot, nout, lerr, ok )
199 CALL dgesvx(
'N',
'N', 2, 1, a, 2, af, 2, ip, eq, r, c, b, 2,
200 $ x, 1, rcond, r1, r2, w, iw, info )
201 CALL chkxer(
'DGESVX', infot, nout, lerr, ok )
209 CALL dgesvxx(
'/',
'N', 0, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
210 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
211 $ err_bnds_c, nparams, params, w, iw, info )
212 CALL chkxer(
'DGESVXX', infot, nout, lerr, ok )
214 CALL dgesvxx(
'N',
'/', 0, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
215 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
216 $ err_bnds_c, nparams, params, w, iw, info )
217 CALL chkxer(
'DGESVXX', infot, nout, lerr, ok )
219 CALL dgesvxx(
'N',
'N', -1, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
220 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
221 $ err_bnds_c, nparams, params, w, iw, info )
222 CALL chkxer(
'DGESVXX', infot, nout, lerr, ok )
224 CALL dgesvxx(
'N',
'N', 0, -1, a, 1, af, 1, ip, eq, r, c, b, 1,
225 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
226 $ err_bnds_c, nparams, params, w, iw, info )
227 CALL chkxer(
'DGESVXX', infot, nout, lerr, ok )
229 CALL dgesvxx(
'N',
'N', 2, 1, a, 1, af, 2, ip, eq, r, c, b, 2,
230 $ x, 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
231 $ err_bnds_c, nparams, params, w, iw, info )
232 CALL chkxer(
'DGESVXX', infot, nout, lerr, ok )
234 CALL dgesvxx(
'N',
'N', 2, 1, a, 2, af, 1, ip, eq, r, c, b, 2,
235 $ x, 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
236 $ err_bnds_c, nparams, params, w, iw, info )
237 CALL chkxer(
'DGESVXX', infot, nout, lerr, ok )
240 CALL dgesvxx(
'F',
'N', 0, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
241 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
242 $ err_bnds_c, nparams, params, w, iw, info )
243 CALL chkxer(
'DGESVXX', infot, nout, lerr, ok )
246 CALL dgesvxx(
'F',
'N', 1, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
247 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
248 $ err_bnds_c, nparams, params, w, iw, info )
249 CALL chkxer(
'DGESVXX', infot, nout, lerr, ok )
252 CALL dgesvxx(
'F',
'N', 1, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
253 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
254 $ err_bnds_c, nparams, params, w, iw, info )
255 CALL chkxer(
'DGESVXX', infot, nout, lerr, ok )
257 CALL dgesvxx(
'N',
'N', 2, 1, a, 2, af, 2, ip, eq, r, c, b, 1,
258 $ x, 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
259 $ err_bnds_c, nparams, params, w, iw, info )
260 CALL chkxer(
'DGESVXX', infot, nout, lerr, ok )
262 CALL dgesvxx(
'N',
'N', 2, 1, a, 2, af, 2, ip, eq, r, c, b, 2,
263 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
264 $ err_bnds_c, nparams, params, w, iw, info )
265 CALL chkxer(
'DGESVXX', infot, nout, lerr, ok )
267 ELSE IF(
lsamen( 2, c2,
'GB' ) )
THEN
273 CALL dgbsv( -1, 0, 0, 0, a, 1, ip, b, 1, info )
274 CALL chkxer(
'DGBSV ', infot, nout, lerr, ok )
276 CALL dgbsv( 1, -1, 0, 0, a, 1, ip, b, 1, info )
277 CALL chkxer(
'DGBSV ', infot, nout, lerr, ok )
279 CALL dgbsv( 1, 0, -1, 0, a, 1, ip, b, 1, info )
280 CALL chkxer(
'DGBSV ', infot, nout, lerr, ok )
282 CALL dgbsv( 0, 0, 0, -1, a, 1, ip, b, 1, info )
283 CALL chkxer(
'DGBSV ', infot, nout, lerr, ok )
285 CALL dgbsv( 1, 1, 1, 0, a, 3, ip, b, 1, info )
286 CALL chkxer(
'DGBSV ', infot, nout, lerr, ok )
288 CALL dgbsv( 2, 0, 0, 0, a, 1, ip, b, 1, info )
289 CALL chkxer(
'DGBSV ', infot, nout, lerr, ok )
295 CALL dgbsvx(
'/',
'N', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
296 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
297 CALL chkxer(
'DGBSVX', infot, nout, lerr, ok )
299 CALL dgbsvx(
'N',
'/', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
300 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
301 CALL chkxer(
'DGBSVX', infot, nout, lerr, ok )
303 CALL dgbsvx(
'N',
'N', -1, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
304 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
305 CALL chkxer(
'DGBSVX', infot, nout, lerr, ok )
307 CALL dgbsvx(
'N',
'N', 1, -1, 0, 0, a, 1, af, 1, ip, eq, r, c,
308 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
309 CALL chkxer(
'DGBSVX', infot, nout, lerr, ok )
311 CALL dgbsvx(
'N',
'N', 1, 0, -1, 0, a, 1, af, 1, ip, eq, r, c,
312 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
313 CALL chkxer(
'DGBSVX', infot, nout, lerr, ok )
315 CALL dgbsvx(
'N',
'N', 0, 0, 0, -1, a, 1, af, 1, ip, eq, r, c,
316 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
317 CALL chkxer(
'DGBSVX', infot, nout, lerr, ok )
319 CALL dgbsvx(
'N',
'N', 1, 1, 1, 0, a, 2, af, 4, ip, eq, r, c,
320 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
321 CALL chkxer(
'DGBSVX', infot, nout, lerr, ok )
323 CALL dgbsvx(
'N',
'N', 1, 1, 1, 0, a, 3, af, 3, ip, eq, r, c,
324 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
325 CALL chkxer(
'DGBSVX', infot, nout, lerr, ok )
328 CALL dgbsvx(
'F',
'N', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
329 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
330 CALL chkxer(
'DGBSVX', infot, nout, lerr, ok )
333 CALL dgbsvx(
'F',
'N', 1, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
334 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
335 CALL chkxer(
'DGBSVX', infot, nout, lerr, ok )
338 CALL dgbsvx(
'F',
'N', 1, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
339 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
340 CALL chkxer(
'DGBSVX', infot, nout, lerr, ok )
342 CALL dgbsvx(
'N',
'N', 2, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
343 $ b, 1, x, 2, rcond, r1, r2, w, iw, info )
344 CALL chkxer(
'DGBSVX', infot, nout, lerr, ok )
346 CALL dgbsvx(
'N',
'N', 2, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
347 $ b, 2, x, 1, rcond, r1, r2, w, iw, info )
348 CALL chkxer(
'DGBSVX', infot, nout, lerr, ok )
356 CALL dgbsvxx(
'/',
'N', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
357 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
358 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
360 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
362 CALL dgbsvxx(
'N',
'/', 0, 1, 1, 0, a, 1, af, 1, ip, eq, r, c,
363 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
364 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
366 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
368 CALL dgbsvxx(
'N',
'N', -1, 1, 1, 0, a, 1, af, 1, ip, eq, r, c,
369 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
370 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
372 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
374 CALL dgbsvxx(
'N',
'N', 2, -1, 1, 0, a, 1, af, 1, ip, eq,
375 $ r, c, b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
376 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
378 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
380 CALL dgbsvxx(
'N',
'N', 2, 1, -1, 0, a, 1, af, 1, ip, eq,
381 $ r, c, b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
382 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
384 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
386 CALL dgbsvxx(
'N',
'N', 0, 1, 1, -1, a, 1, af, 1, ip, eq, r, c,
387 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
388 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
390 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
392 CALL dgbsvxx(
'N',
'N', 2, 1, 1, 1, a, 2, af, 2, ip, eq, r, c,
393 $ b, 2, x, 2, rcond, rpvgrw, berr, n_err_bnds,
394 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
396 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
398 CALL dgbsvxx(
'N',
'N', 2, 1, 1, 1, a, 3, af, 3, ip, eq, r, c,
399 $ b, 2, x, 2, rcond, rpvgrw, berr, n_err_bnds,
400 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
402 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
405 CALL dgbsvxx(
'F',
'N', 0, 1, 1, 0, a, 3, af, 4, ip, eq, r, c,
406 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
407 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
409 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
412 CALL dgbsvxx(
'F',
'N', 1, 1, 1, 0, a, 3, af, 4, ip, eq, r, c,
413 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
414 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
416 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
419 CALL dgbsvxx(
'F',
'N', 1, 1, 1, 0, a, 3, af, 4, ip, eq, r, c,
420 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
421 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
423 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
425 CALL dgbsvxx(
'N',
'N', 2, 1, 1, 1, a, 3, af, 4, ip, eq, r, c,
426 $ b, 1, x, 2, rcond, rpvgrw, berr, n_err_bnds,
427 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
429 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
431 CALL dgbsvxx(
'N',
'N', 2, 1, 1, 1, a, 3, af, 4, ip, eq, r, c,
432 $ b, 2, x, 1, rcond, rpvgrw, berr, n_err_bnds,
433 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
435 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
437 ELSE IF(
lsamen( 2, c2,
'GT' ) )
THEN
443 CALL dgtsv( -1, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ), b, 1,
445 CALL chkxer(
'DGTSV ', infot, nout, lerr, ok )
447 CALL dgtsv( 0, -1, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ), b, 1,
449 CALL chkxer(
'DGTSV ', infot, nout, lerr, ok )
451 CALL dgtsv( 2, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ), b, 1, info )
452 CALL chkxer(
'DGTSV ', infot, nout, lerr, ok )
458 CALL dgtsvx(
'/',
'N', 0, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
459 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
460 $ ip, b, 1, x, 1, rcond, r1, r2, w, iw, info )
461 CALL chkxer(
'DGTSVX', infot, nout, lerr, ok )
463 CALL dgtsvx(
'N',
'/', 0, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
464 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
465 $ ip, b, 1, x, 1, rcond, r1, r2, w, iw, info )
466 CALL chkxer(
'DGTSVX', infot, nout, lerr, ok )
468 CALL dgtsvx(
'N',
'N', -1, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
469 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
470 $ ip, b, 1, x, 1, rcond, r1, r2, w, iw, info )
471 CALL chkxer(
'DGTSVX', infot, nout, lerr, ok )
473 CALL dgtsvx(
'N',
'N', 0, -1, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
474 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
475 $ ip, b, 1, x, 1, rcond, r1, r2, w, iw, info )
476 CALL chkxer(
'DGTSVX', infot, nout, lerr, ok )
478 CALL dgtsvx(
'N',
'N', 2, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
479 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
480 $ ip, b, 1, x, 2, rcond, r1, r2, w, iw, info )
481 CALL chkxer(
'DGTSVX', infot, nout, lerr, ok )
483 CALL dgtsvx(
'N',
'N', 2, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
484 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
485 $ ip, b, 2, x, 1, rcond, r1, r2, w, iw, info )
486 CALL chkxer(
'DGTSVX', infot, nout, lerr, ok )
488 ELSE IF(
lsamen( 2, c2,
'PO' ) )
THEN
494 CALL dposv(
'/', 0, 0, a, 1, b, 1, info )
495 CALL chkxer(
'DPOSV ', infot, nout, lerr, ok )
497 CALL dposv(
'U', -1, 0, a, 1, b, 1, info )
498 CALL chkxer(
'DPOSV ', infot, nout, lerr, ok )
500 CALL dposv(
'U', 0, -1, a, 1, b, 1, info )
501 CALL chkxer(
'DPOSV ', infot, nout, lerr, ok )
503 CALL dposv(
'U', 2, 0, a, 1, b, 2, info )
504 CALL chkxer(
'DPOSV ', infot, nout, lerr, ok )
506 CALL dposv(
'U', 2, 0, a, 2, b, 1, info )
507 CALL chkxer(
'DPOSV ', infot, nout, lerr, ok )
513 CALL dposvx(
'/',
'U', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
514 $ rcond, r1, r2, w, iw, info )
515 CALL chkxer(
'DPOSVX', infot, nout, lerr, ok )
517 CALL dposvx(
'N',
'/', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
518 $ rcond, r1, r2, w, iw, info )
519 CALL chkxer(
'DPOSVX', infot, nout, lerr, ok )
521 CALL dposvx(
'N',
'U', -1, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
522 $ rcond, r1, r2, w, iw, info )
523 CALL chkxer(
'DPOSVX', infot, nout, lerr, ok )
525 CALL dposvx(
'N',
'U', 0, -1, a, 1, af, 1, eq, c, b, 1, x, 1,
526 $ rcond, r1, r2, w, iw, info )
527 CALL chkxer(
'DPOSVX', infot, nout, lerr, ok )
529 CALL dposvx(
'N',
'U', 2, 0, a, 1, af, 2, eq, c, b, 2, x, 2,
530 $ rcond, r1, r2, w, iw, info )
531 CALL chkxer(
'DPOSVX', infot, nout, lerr, ok )
533 CALL dposvx(
'N',
'U', 2, 0, a, 2, af, 1, eq, c, b, 2, x, 2,
534 $ rcond, r1, r2, w, iw, info )
535 CALL chkxer(
'DPOSVX', infot, nout, lerr, ok )
538 CALL dposvx(
'F',
'U', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
539 $ rcond, r1, r2, w, iw, info )
540 CALL chkxer(
'DPOSVX', infot, nout, lerr, ok )
543 CALL dposvx(
'F',
'U', 1, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
544 $ rcond, r1, r2, w, iw, info )
545 CALL chkxer(
'DPOSVX', infot, nout, lerr, ok )
547 CALL dposvx(
'N',
'U', 2, 0, a, 2, af, 2, eq, c, b, 1, x, 2,
548 $ rcond, r1, r2, w, iw, info )
549 CALL chkxer(
'DPOSVX', infot, nout, lerr, ok )
551 CALL dposvx(
'N',
'U', 2, 0, a, 2, af, 2, eq, c, b, 2, x, 1,
552 $ rcond, r1, r2, w, iw, info )
553 CALL chkxer(
'DPOSVX', infot, nout, lerr, ok )
561 CALL dposvxx(
'/',
'U', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
562 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
563 $ err_bnds_c, nparams, params, w, iw, info )
564 CALL chkxer(
'DPOSVXX', infot, nout, lerr, ok )
566 CALL dposvxx(
'N',
'/', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
567 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
568 $ err_bnds_c, nparams, params, w, iw, info )
569 CALL chkxer(
'DPOSVXX', infot, nout, lerr, ok )
571 CALL dposvxx(
'N',
'U', -1, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
572 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
573 $ err_bnds_c, nparams, params, w, iw, info )
574 CALL chkxer(
'DPOSVXX', infot, nout, lerr, ok )
576 CALL dposvxx(
'N',
'U', 0, -1, a, 1, af, 1, eq, c, b, 1, x, 1,
577 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
578 $ err_bnds_c, nparams, params, w, iw, info )
579 CALL chkxer(
'DPOSVXX', infot, nout, lerr, ok )
581 CALL dposvxx(
'N',
'U', 2, 0, a, 1, af, 2, eq, c, b, 2, x, 2,
582 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
583 $ err_bnds_c, nparams, params, w, iw, info )
584 CALL chkxer(
'DPOSVXX', infot, nout, lerr, ok )
586 CALL dposvxx(
'N',
'U', 2, 0, a, 2, af, 1, eq, c, b, 2, x, 2,
587 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
588 $ err_bnds_c, nparams, params, w, iw, info )
589 CALL chkxer(
'DPOSVXX', infot, nout, lerr, ok )
592 CALL dposvxx(
'F',
'U', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
593 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
594 $ err_bnds_c, nparams, params, w, iw, info )
595 CALL chkxer(
'DPOSVXX', infot, nout, lerr, ok )
598 CALL dposvxx(
'F',
'U', 1, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
599 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
600 $ err_bnds_c, nparams, params, w, iw, info )
601 CALL chkxer(
'DPOSVXX', infot, nout, lerr, ok )
603 CALL dposvxx(
'N',
'U', 2, 0, a, 2, af, 2, eq, c, b, 1, x, 2,
604 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
605 $ err_bnds_c, nparams, params, w, iw, info )
606 CALL chkxer(
'DPOSVXX', infot, nout, lerr, ok )
608 CALL dposvxx(
'N',
'U', 2, 0, a, 2, af, 2, eq, c, b, 2, x, 1,
609 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
610 $ err_bnds_c, nparams, params, w, iw, info )
611 CALL chkxer(
'DPOSVXX', infot, nout, lerr, ok )
613 ELSE IF(
lsamen( 2, c2,
'PP' ) )
THEN
619 CALL dppsv(
'/', 0, 0, a, b, 1, info )
620 CALL chkxer(
'DPPSV ', infot, nout, lerr, ok )
622 CALL dppsv(
'U', -1, 0, a, b, 1, info )
623 CALL chkxer(
'DPPSV ', infot, nout, lerr, ok )
625 CALL dppsv(
'U', 0, -1, a, b, 1, info )
626 CALL chkxer(
'DPPSV ', infot, nout, lerr, ok )
628 CALL dppsv(
'U', 2, 0, a, b, 1, info )
629 CALL chkxer(
'DPPSV ', infot, nout, lerr, ok )
635 CALL dppsvx(
'/',
'U', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
636 $ r1, r2, w, iw, info )
637 CALL chkxer(
'DPPSVX', infot, nout, lerr, ok )
639 CALL dppsvx(
'N',
'/', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
640 $ r1, r2, w, iw, info )
641 CALL chkxer(
'DPPSVX', infot, nout, lerr, ok )
643 CALL dppsvx(
'N',
'U', -1, 0, a, af, eq, c, b, 1, x, 1, rcond,
644 $ r1, r2, w, iw, info )
645 CALL chkxer(
'DPPSVX', infot, nout, lerr, ok )
647 CALL dppsvx(
'N',
'U', 0, -1, a, af, eq, c, b, 1, x, 1, rcond,
648 $ r1, r2, w, iw, info )
649 CALL chkxer(
'DPPSVX', infot, nout, lerr, ok )
652 CALL dppsvx(
'F',
'U', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
653 $ r1, r2, w, iw, info )
654 CALL chkxer(
'DPPSVX', infot, nout, lerr, ok )
657 CALL dppsvx(
'F',
'U', 1, 0, a, af, eq, c, b, 1, x, 1, rcond,
658 $ r1, r2, w, iw, info )
659 CALL chkxer(
'DPPSVX', infot, nout, lerr, ok )
661 CALL dppsvx(
'N',
'U', 2, 0, a, af, eq, c, b, 1, x, 2, rcond,
662 $ r1, r2, w, iw, info )
663 CALL chkxer(
'DPPSVX', infot, nout, lerr, ok )
665 CALL dppsvx(
'N',
'U', 2, 0, a, af, eq, c, b, 2, x, 1, rcond,
666 $ r1, r2, w, iw, info )
667 CALL chkxer(
'DPPSVX', infot, nout, lerr, ok )
669 ELSE IF(
lsamen( 2, c2,
'PB' ) )
THEN
675 CALL dpbsv(
'/', 0, 0, 0, a, 1, b, 1, info )
676 CALL chkxer(
'DPBSV ', infot, nout, lerr, ok )
678 CALL dpbsv(
'U', -1, 0, 0, a, 1, b, 1, info )
679 CALL chkxer(
'DPBSV ', infot, nout, lerr, ok )
681 CALL dpbsv(
'U', 1, -1, 0, a, 1, b, 1, info )
682 CALL chkxer(
'DPBSV ', infot, nout, lerr, ok )
684 CALL dpbsv(
'U', 0, 0, -1, a, 1, b, 1, info )
685 CALL chkxer(
'DPBSV ', infot, nout, lerr, ok )
687 CALL dpbsv(
'U', 1, 1, 0, a, 1, b, 2, info )
688 CALL chkxer(
'DPBSV ', infot, nout, lerr, ok )
690 CALL dpbsv(
'U', 2, 0, 0, a, 1, b, 1, info )
691 CALL chkxer(
'DPBSV ', infot, nout, lerr, ok )
697 CALL dpbsvx(
'/',
'U', 0, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
698 $ rcond, r1, r2, w, iw, info )
699 CALL chkxer(
'DPBSVX', infot, nout, lerr, ok )
701 CALL dpbsvx(
'N',
'/', 0, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
702 $ rcond, r1, r2, w, iw, info )
703 CALL chkxer(
'DPBSVX', infot, nout, lerr, ok )
705 CALL dpbsvx(
'N',
'U', -1, 0, 0, a, 1, af, 1, eq, c, b, 1, x,
706 $ 1, rcond, r1, r2, w, iw, info )
707 CALL chkxer(
'DPBSVX', infot, nout, lerr, ok )
709 CALL dpbsvx(
'N',
'U', 1, -1, 0, a, 1, af, 1, eq, c, b, 1, x,
710 $ 1, rcond, r1, r2, w, iw, info )
711 CALL chkxer(
'DPBSVX', infot, nout, lerr, ok )
713 CALL dpbsvx(
'N',
'U', 0, 0, -1, a, 1, af, 1, eq, c, b, 1, x,
714 $ 1, rcond, r1, r2, w, iw, info )
715 CALL chkxer(
'DPBSVX', infot, nout, lerr, ok )
717 CALL dpbsvx(
'N',
'U', 1, 1, 0, a, 1, af, 2, eq, c, b, 2, x, 2,
718 $ rcond, r1, r2, w, iw, info )
719 CALL chkxer(
'DPBSVX', infot, nout, lerr, ok )
721 CALL dpbsvx(
'N',
'U', 1, 1, 0, a, 2, af, 1, eq, c, b, 2, x, 2,
722 $ rcond, r1, r2, w, iw, info )
723 CALL chkxer(
'DPBSVX', infot, nout, lerr, ok )
726 CALL dpbsvx(
'F',
'U', 0, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
727 $ rcond, r1, r2, w, iw, info )
728 CALL chkxer(
'DPBSVX', infot, nout, lerr, ok )
731 CALL dpbsvx(
'F',
'U', 1, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
732 $ rcond, r1, r2, w, iw, info )
733 CALL chkxer(
'DPBSVX', infot, nout, lerr, ok )
735 CALL dpbsvx(
'N',
'U', 2, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 2,
736 $ rcond, r1, r2, w, iw, info )
737 CALL chkxer(
'DPBSVX', infot, nout, lerr, ok )
739 CALL dpbsvx(
'N',
'U', 2, 0, 0, a, 1, af, 1, eq, c, b, 2, x, 1,
740 $ rcond, r1, r2, w, iw, info )
741 CALL chkxer(
'DPBSVX', infot, nout, lerr, ok )
743 ELSE IF(
lsamen( 2, c2,
'PT' ) )
THEN
749 CALL dptsv( -1, 0, a( 1, 1 ), a( 1, 2 ), b, 1, info )
750 CALL chkxer(
'DPTSV ', infot, nout, lerr, ok )
752 CALL dptsv( 0, -1, a( 1, 1 ), a( 1, 2 ), b, 1, info )
753 CALL chkxer(
'DPTSV ', infot, nout, lerr, ok )
755 CALL dptsv( 2, 0, a( 1, 1 ), a( 1, 2 ), b, 1, info )
756 CALL chkxer(
'DPTSV ', infot, nout, lerr, ok )
762 CALL dptsvx(
'/', 0, 0, a( 1, 1 ), a( 1, 2 ), af( 1, 1 ),
763 $ af( 1, 2 ), b, 1, x, 1, rcond, r1, r2, w, info )
764 CALL chkxer(
'DPTSVX', infot, nout, lerr, ok )
766 CALL dptsvx(
'N', -1, 0, a( 1, 1 ), a( 1, 2 ), af( 1, 1 ),
767 $ af( 1, 2 ), b, 1, x, 1, rcond, r1, r2, w, info )
768 CALL chkxer(
'DPTSVX', infot, nout, lerr, ok )
770 CALL dptsvx(
'N', 0, -1, a( 1, 1 ), a( 1, 2 ), af( 1, 1 ),
771 $ af( 1, 2 ), b, 1, x, 1, rcond, r1, r2, w, info )
772 CALL chkxer(
'DPTSVX', infot, nout, lerr, ok )
774 CALL dptsvx(
'N', 2, 0, a( 1, 1 ), a( 1, 2 ), af( 1, 1 ),
775 $ af( 1, 2 ), b, 1, x, 2, rcond, r1, r2, w, info )
776 CALL chkxer(
'DPTSVX', infot, nout, lerr, ok )
778 CALL dptsvx(
'N', 2, 0, a( 1, 1 ), a( 1, 2 ), af( 1, 1 ),
779 $ af( 1, 2 ), b, 2, x, 1, rcond, r1, r2, w, info )
780 CALL chkxer(
'DPTSVX', infot, nout, lerr, ok )
782 ELSE IF(
lsamen( 2, c2,
'SY' ) )
THEN
788 CALL dsysv(
'/', 0, 0, a, 1, ip, b, 1, w, 1, info )
789 CALL chkxer(
'DSYSV ', infot, nout, lerr, ok )
791 CALL dsysv(
'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
792 CALL chkxer(
'DSYSV ', infot, nout, lerr, ok )
794 CALL dsysv(
'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
795 CALL chkxer(
'DSYSV ', infot, nout, lerr, ok )
797 CALL dsysv(
'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
798 CALL chkxer(
'DSYSV ', infot, nout, lerr, ok )
804 CALL dsysvx(
'/',
'U', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
805 $ rcond, r1, r2, w, 1, iw, info )
806 CALL chkxer(
'DSYSVX', infot, nout, lerr, ok )
808 CALL dsysvx(
'N',
'/', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
809 $ rcond, r1, r2, w, 1, iw, info )
810 CALL chkxer(
'DSYSVX', infot, nout, lerr, ok )
812 CALL dsysvx(
'N',
'U', -1, 0, a, 1, af, 1, ip, b, 1, x, 1,
813 $ rcond, r1, r2, w, 1, iw, info )
814 CALL chkxer(
'DSYSVX', infot, nout, lerr, ok )
816 CALL dsysvx(
'N',
'U', 0, -1, a, 1, af, 1, ip, b, 1, x, 1,
817 $ rcond, r1, r2, w, 1, iw, info )
818 CALL chkxer(
'DSYSVX', infot, nout, lerr, ok )
820 CALL dsysvx(
'N',
'U', 2, 0, a, 1, af, 2, ip, b, 2, x, 2,
821 $ rcond, r1, r2, w, 4, iw, info )
822 CALL chkxer(
'DSYSVX', infot, nout, lerr, ok )
824 CALL dsysvx(
'N',
'U', 2, 0, a, 2, af, 1, ip, b, 2, x, 2,
825 $ rcond, r1, r2, w, 4, iw, info )
826 CALL chkxer(
'DSYSVX', infot, nout, lerr, ok )
828 CALL dsysvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 1, x, 2,
829 $ rcond, r1, r2, w, 4, iw, info )
830 CALL chkxer(
'DSYSVX', infot, nout, lerr, ok )
832 CALL dsysvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 1,
833 $ rcond, r1, r2, w, 4, iw, info )
834 CALL chkxer(
'DSYSVX', infot, nout, lerr, ok )
836 CALL dsysvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 2,
837 $ rcond, r1, r2, w, 3, iw, info )
838 CALL chkxer(
'DSYSVX', infot, nout, lerr, ok )
847 CALL dsysvxx(
'/',
'U', 0, 0, a, 1, af, 1, ip, eq, r, b, 1, x,
848 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
849 $ err_bnds_c, nparams, params, w, iw, info )
850 CALL chkxer(
'DSYSVXX', infot, nout, lerr, ok )
852 CALL dsysvxx(
'N',
'/', 0, 0, a, 1, af, 1, ip, eq, r, b, 1, x,
853 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
854 $ err_bnds_c, nparams, params, w, iw, info )
855 CALL chkxer(
'DSYSVXX', infot, nout, lerr, ok )
857 CALL dsysvxx(
'N',
'U', -1, 0, a, 1, af, 1, ip, eq, r, b, 1, x,
858 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
859 $ err_bnds_c, nparams, params, w, iw, info )
860 CALL chkxer(
'DSYSVXX', infot, nout, lerr, ok )
863 CALL dsysvxx(
'N',
'U', 0, -1, a, 1, af, 1, ip, eq, r, b, 1, x,
864 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
865 $ err_bnds_c, nparams, params, w, iw, info )
866 CALL chkxer(
'DSYSVXX', infot, nout, lerr, ok )
869 CALL dsysvxx(
'N',
'U', 2, 0, a, 1, af, 2, ip, eq, r, b, 2, x,
870 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
871 $ err_bnds_c, nparams, params, w, iw, info )
872 CALL chkxer(
'DSYSVXX', infot, nout, lerr, ok )
874 CALL dsysvxx(
'N',
'U', 2, 0, a, 2, af, 1, ip, eq, r, b, 2, x,
875 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
876 $ err_bnds_c, nparams, params, w, iw, info )
877 CALL chkxer(
'DSYSVXX', infot, nout, lerr, ok )
880 CALL dsysvxx(
'N',
'U', 2, 0, a, 2, af, 2, ip, eq, r, b, 1, x,
881 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
882 $ err_bnds_c, nparams, params, w, iw, info )
883 CALL chkxer(
'DSYSVXX', infot, nout, lerr, ok )
885 CALL dsysvxx(
'N',
'U', 2, 0, a, 2, af, 2, ip, eq, r, b, 2, x,
886 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
887 $ err_bnds_c, nparams, params, w, iw, info )
888 CALL chkxer(
'DSYSVXX', infot, nout, lerr, ok )
890 ELSE IF(
lsamen( 2, c2,
'SP' ) )
THEN
896 CALL dspsv(
'/', 0, 0, a, ip, b, 1, info )
897 CALL chkxer(
'DSPSV ', infot, nout, lerr, ok )
899 CALL dspsv(
'U', -1, 0, a, ip, b, 1, info )
900 CALL chkxer(
'DSPSV ', infot, nout, lerr, ok )
902 CALL dspsv(
'U', 0, -1, a, ip, b, 1, info )
903 CALL chkxer(
'DSPSV ', infot, nout, lerr, ok )
905 CALL dspsv(
'U', 2, 0, a, ip, b, 1, info )
906 CALL chkxer(
'DSPSV ', infot, nout, lerr, ok )
912 CALL dspsvx(
'/',
'U', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
914 CALL chkxer(
'DSPSVX', infot, nout, lerr, ok )
916 CALL dspsvx(
'N',
'/', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
918 CALL chkxer(
'DSPSVX', infot, nout, lerr, ok )
920 CALL dspsvx(
'N',
'U', -1, 0, a, af, ip, b, 1, x, 1, rcond, r1,
922 CALL chkxer(
'DSPSVX', infot, nout, lerr, ok )
924 CALL dspsvx(
'N',
'U', 0, -1, a, af, ip, b, 1, x, 1, rcond, r1,
926 CALL chkxer(
'DSPSVX', infot, nout, lerr, ok )
928 CALL dspsvx(
'N',
'U', 2, 0, a, af, ip, b, 1, x, 2, rcond, r1,
930 CALL chkxer(
'DSPSVX', infot, nout, lerr, ok )
932 CALL dspsvx(
'N',
'U', 2, 0, a, af, ip, b, 2, x, 1, rcond, r1,
934 CALL chkxer(
'DSPSVX', infot, nout, lerr, ok )
940 WRITE( nout, fmt = 9999 )path
942 WRITE( nout, fmt = 9998 )path
945 9999
FORMAT( 1x, a3,
' drivers passed the tests of the error exits' )
946 9998
FORMAT(
' *** ', a3,
' drivers failed the tests of the error ',
subroutine dpbsvx(FACT, UPLO, N, KD, NRHS, AB, LDAB, AFB, LDAFB, EQUED, S, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, IWORK, INFO)
DPBSVX computes the solution to system of linear equations A * X = B for OTHER matrices ...
subroutine dsysvx(FACT, UPLO, N, NRHS, A, LDA, AF, LDAF, IPIV, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, LWORK, IWORK, INFO)
DSYSVX computes the solution to system of linear equations A * X = B for SY matrices ...
subroutine dgbsv(N, KL, KU, NRHS, AB, LDAB, IPIV, B, LDB, INFO)
DGBSV computes the solution to system of linear equations A * X = B for GB matrices (simple driver) ...
subroutine dppsvx(FACT, UPLO, N, NRHS, AP, AFP, EQUED, S, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, IWORK, INFO)
DPPSVX computes the solution to system of linear equations A * X = B for OTHER matrices ...
subroutine dgtsvx(FACT, TRANS, N, NRHS, DL, D, DU, DLF, DF, DUF, DU2, IPIV, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, IWORK, INFO)
DGTSVX computes the solution to system of linear equations A * X = B for GT matrices ...
subroutine dposvx(FACT, UPLO, N, NRHS, A, LDA, AF, LDAF, EQUED, S, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, IWORK, INFO)
DPOSVX computes the solution to system of linear equations A * X = B for PO matrices ...
subroutine dposv(UPLO, N, NRHS, A, LDA, B, LDB, INFO)
DPOSV computes the solution to system of linear equations A * X = B for PO matrices ...
subroutine dgesvxx(FACT, TRANS, N, NRHS, A, LDA, AF, LDAF, IPIV, EQUED, R, C, B, LDB, X, LDX, RCOND, RPVGRW, BERR, N_ERR_BNDS, ERR_BNDS_NORM, ERR_BNDS_COMP, NPARAMS, PARAMS, WORK, IWORK, INFO)
DGESVXX computes the solution to system of linear equations A * X = B for GE matrices ...
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine dgesv(N, NRHS, A, LDA, IPIV, B, LDB, INFO)
DGESV computes the solution to system of linear equations A * X = B for GE matrices ...
subroutine dptsvx(FACT, N, NRHS, D, E, DF, EF, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, INFO)
DPTSVX computes the solution to system of linear equations A * X = B for PT matrices ...
subroutine dgtsv(N, NRHS, DL, D, DU, B, LDB, INFO)
DGTSV computes the solution to system of linear equations A * X = B for GT matrices ...
subroutine dposvxx(FACT, UPLO, N, NRHS, A, LDA, AF, LDAF, EQUED, S, B, LDB, X, LDX, RCOND, RPVGRW, BERR, N_ERR_BNDS, ERR_BNDS_NORM, ERR_BNDS_COMP, NPARAMS, PARAMS, WORK, IWORK, INFO)
DPOSVXX computes the solution to system of linear equations A * X = B for PO matrices ...
subroutine dsysv(UPLO, N, NRHS, A, LDA, IPIV, B, LDB, WORK, LWORK, INFO)
DSYSV computes the solution to system of linear equations A * X = B for SY matrices ...
subroutine dgbsvx(FACT, TRANS, N, KL, KU, NRHS, AB, LDAB, AFB, LDAFB, IPIV, EQUED, R, C, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, IWORK, INFO)
DGBSVX computes the solution to system of linear equations A * X = B for GB matrices ...
logical function lsamen(N, CA, CB)
LSAMEN
subroutine dsysvxx(FACT, UPLO, N, NRHS, A, LDA, AF, LDAF, IPIV, EQUED, S, B, LDB, X, LDX, RCOND, RPVGRW, BERR, N_ERR_BNDS, ERR_BNDS_NORM, ERR_BNDS_COMP, NPARAMS, PARAMS, WORK, IWORK, INFO)
DSYSVXX
subroutine dspsvx(FACT, UPLO, N, NRHS, AP, AFP, IPIV, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, IWORK, INFO)
DSPSVX computes the solution to system of linear equations A * X = B for OTHER matrices ...
subroutine dpbsv(UPLO, N, KD, NRHS, AB, LDAB, B, LDB, INFO)
DPBSV computes the solution to system of linear equations A * X = B for OTHER matrices ...
subroutine dgesvx(FACT, TRANS, N, NRHS, A, LDA, AF, LDAF, IPIV, EQUED, R, C, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, IWORK, INFO)
DGESVX computes the solution to system of linear equations A * X = B for GE matrices ...
subroutine dgbsvxx(FACT, TRANS, N, KL, KU, NRHS, AB, LDAB, AFB, LDAFB, IPIV, EQUED, R, C, B, LDB, X, LDX, RCOND, RPVGRW, BERR, N_ERR_BNDS, ERR_BNDS_NORM, ERR_BNDS_COMP, NPARAMS, PARAMS, WORK, IWORK, INFO)
DGBSVXX computes the solution to system of linear equations A * X = B for GB matrices ...
subroutine dspsv(UPLO, N, NRHS, AP, IPIV, B, LDB, INFO)
DSPSV computes the solution to system of linear equations A * X = B for OTHER matrices ...
subroutine dppsv(UPLO, N, NRHS, AP, B, LDB, INFO)
DPPSV computes the solution to system of linear equations A * X = B for OTHER matrices ...
subroutine dptsv(N, NRHS, D, E, B, LDB, INFO)
DPTSV computes the solution to system of linear equations A * X = B for PT matrices ...