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 derrvx(PATH, NUNIT)
DERRVX
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 ...
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 ...