56 SUBROUTINE derrvx( PATH, NUNIT )
72 parameter ( nmax = 4 )
74 parameter ( one = 1.0d+0 )
79 INTEGER i, info, j, n_err_bnds, nparams
80 DOUBLE PRECISION rcond, rpvgrw, berr
83 INTEGER ip( nmax ), iw( nmax )
84 DOUBLE PRECISION a( nmax, nmax ), af( nmax, nmax ), b( nmax ),
85 $ c( nmax ), e( nmax ), r( nmax ), r1( nmax ),
86 $ r2( nmax ), w( 2*nmax ), x( nmax ),
87 $ err_bnds_n( nmax, 3 ), err_bnds_c( nmax, 3 ),
107 COMMON / infoc / infot, nout, ok, lerr
108 COMMON / srnamc / srnamt
116 WRITE( nout, fmt = * )
123 a( i, j ) = 1.d0 / dble( i+j )
124 af( i, j ) = 1.d0 / dble( i+j )
139 IF(
lsamen( 2, c2,
'GE' ) )
THEN
145 CALL dgesv( -1, 0, a, 1, ip, b, 1, info )
146 CALL chkxer(
'DGESV ', infot, nout, lerr, ok )
148 CALL dgesv( 0, -1, a, 1, ip, b, 1, info )
149 CALL chkxer(
'DGESV ', infot, nout, lerr, ok )
151 CALL dgesv( 2, 1, a, 1, ip, b, 2, info )
152 CALL chkxer(
'DGESV ', infot, nout, lerr, ok )
154 CALL dgesv( 2, 1, a, 2, ip, b, 1, info )
155 CALL chkxer(
'DGESV ', infot, nout, lerr, ok )
161 CALL dgesvx(
'/',
'N', 0, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
162 $ x, 1, rcond, r1, r2, w, iw, info )
163 CALL chkxer(
'DGESVX', infot, nout, lerr, ok )
165 CALL dgesvx(
'N',
'/', 0, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
166 $ x, 1, rcond, r1, r2, w, iw, info )
167 CALL chkxer(
'DGESVX', infot, nout, lerr, ok )
169 CALL dgesvx(
'N',
'N', -1, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
170 $ x, 1, rcond, r1, r2, w, iw, info )
171 CALL chkxer(
'DGESVX', infot, nout, lerr, ok )
173 CALL dgesvx(
'N',
'N', 0, -1, a, 1, af, 1, ip, eq, r, c, b, 1,
174 $ x, 1, rcond, r1, r2, w, iw, info )
175 CALL chkxer(
'DGESVX', infot, nout, lerr, ok )
177 CALL dgesvx(
'N',
'N', 2, 1, a, 1, af, 2, ip, eq, r, c, b, 2,
178 $ x, 2, rcond, r1, r2, w, iw, info )
179 CALL chkxer(
'DGESVX', infot, nout, lerr, ok )
181 CALL dgesvx(
'N',
'N', 2, 1, a, 2, af, 1, ip, eq, r, c, b, 2,
182 $ x, 2, rcond, r1, r2, w, iw, info )
183 CALL chkxer(
'DGESVX', infot, nout, lerr, ok )
186 CALL dgesvx(
'F',
'N', 0, 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 )
196 CALL dgesvx(
'F',
'N', 1, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
197 $ x, 1, rcond, r1, r2, w, iw, info )
198 CALL chkxer(
'DGESVX', infot, nout, lerr, ok )
200 CALL dgesvx(
'N',
'N', 2, 1, a, 2, af, 2, ip, eq, r, c, b, 1,
201 $ x, 2, rcond, r1, r2, w, iw, info )
202 CALL chkxer(
'DGESVX', infot, nout, lerr, ok )
204 CALL dgesvx(
'N',
'N', 2, 1, a, 2, af, 2, ip, eq, r, c, b, 2,
205 $ x, 1, rcond, r1, r2, w, iw, info )
206 CALL chkxer(
'DGESVX', 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',
'/', 0, 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', -1, 0, 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', 0, -1, a, 1, af, 1, ip, eq, r, c, b, 1,
230 $ x, 1, 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, 1, af, 2, 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 )
239 CALL dgesvxx(
'N',
'N', 2, 1, a, 2, af, 1, ip, eq, r, c, b, 2,
240 $ x, 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
241 $ err_bnds_c, nparams, params, w, iw, info )
242 CALL chkxer(
'DGESVXX', infot, nout, lerr, ok )
245 CALL dgesvxx(
'F',
'N', 0, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
246 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
247 $ err_bnds_c, nparams, params, w, iw, info )
248 CALL chkxer(
'DGESVXX', infot, nout, lerr, ok )
251 CALL dgesvxx(
'F',
'N', 1, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
252 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
253 $ err_bnds_c, nparams, params, w, iw, info )
254 CALL chkxer(
'DGESVXX', infot, nout, lerr, ok )
257 CALL dgesvxx(
'F',
'N', 1, 0, a, 1, af, 1, ip, eq, r, c, b, 1,
258 $ x, 1, 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, 1,
263 $ x, 2, 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 CALL dgesvxx(
'N',
'N', 2, 1, a, 2, af, 2, ip, eq, r, c, b, 2,
268 $ x, 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
269 $ err_bnds_c, nparams, params, w, iw, info )
270 CALL chkxer(
'DGESVXX', infot, nout, lerr, ok )
272 ELSE IF(
lsamen( 2, c2,
'GB' ) )
THEN
278 CALL dgbsv( -1, 0, 0, 0, a, 1, ip, b, 1, info )
279 CALL chkxer(
'DGBSV ', infot, nout, lerr, ok )
281 CALL dgbsv( 1, -1, 0, 0, a, 1, ip, b, 1, info )
282 CALL chkxer(
'DGBSV ', infot, nout, lerr, ok )
284 CALL dgbsv( 1, 0, -1, 0, a, 1, ip, b, 1, info )
285 CALL chkxer(
'DGBSV ', infot, nout, lerr, ok )
287 CALL dgbsv( 0, 0, 0, -1, a, 1, ip, b, 1, info )
288 CALL chkxer(
'DGBSV ', infot, nout, lerr, ok )
290 CALL dgbsv( 1, 1, 1, 0, a, 3, ip, b, 1, info )
291 CALL chkxer(
'DGBSV ', infot, nout, lerr, ok )
293 CALL dgbsv( 2, 0, 0, 0, a, 1, ip, b, 1, info )
294 CALL chkxer(
'DGBSV ', infot, nout, lerr, ok )
300 CALL dgbsvx(
'/',
'N', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
301 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
302 CALL chkxer(
'DGBSVX', infot, nout, lerr, ok )
304 CALL dgbsvx(
'N',
'/', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
305 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
306 CALL chkxer(
'DGBSVX', infot, nout, lerr, ok )
308 CALL dgbsvx(
'N',
'N', -1, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
309 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
310 CALL chkxer(
'DGBSVX', infot, nout, lerr, ok )
312 CALL dgbsvx(
'N',
'N', 1, -1, 0, 0, a, 1, af, 1, ip, eq, r, c,
313 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
314 CALL chkxer(
'DGBSVX', infot, nout, lerr, ok )
316 CALL dgbsvx(
'N',
'N', 1, 0, -1, 0, a, 1, af, 1, ip, eq, r, c,
317 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
318 CALL chkxer(
'DGBSVX', infot, nout, lerr, ok )
320 CALL dgbsvx(
'N',
'N', 0, 0, 0, -1, a, 1, af, 1, ip, eq, r, c,
321 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
322 CALL chkxer(
'DGBSVX', infot, nout, lerr, ok )
324 CALL dgbsvx(
'N',
'N', 1, 1, 1, 0, a, 2, af, 4, ip, eq, r, c,
325 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
326 CALL chkxer(
'DGBSVX', infot, nout, lerr, ok )
328 CALL dgbsvx(
'N',
'N', 1, 1, 1, 0, a, 3, af, 3, 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', 0, 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 )
343 CALL dgbsvx(
'F',
'N', 1, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
344 $ b, 1, x, 1, rcond, r1, r2, w, iw, info )
345 CALL chkxer(
'DGBSVX', infot, nout, lerr, ok )
347 CALL dgbsvx(
'N',
'N', 2, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
348 $ b, 1, x, 2, rcond, r1, r2, w, iw, info )
349 CALL chkxer(
'DGBSVX', infot, nout, lerr, ok )
351 CALL dgbsvx(
'N',
'N', 2, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
352 $ b, 2, x, 1, rcond, r1, r2, w, iw, info )
353 CALL chkxer(
'DGBSVX', infot, nout, lerr, ok )
361 CALL dgbsvxx(
'/',
'N', 0, 0, 0, 0, a, 1, af, 1, ip, eq, r, c,
362 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
363 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
365 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
367 CALL dgbsvxx(
'N',
'/', 0, 1, 1, 0, a, 1, af, 1, ip, eq, r, c,
368 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
369 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
371 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
373 CALL dgbsvxx(
'N',
'N', -1, 1, 1, 0, a, 1, af, 1, ip, eq, r, c,
374 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
375 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
377 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
379 CALL dgbsvxx(
'N',
'N', 2, -1, 1, 0, a, 1, af, 1, ip, eq,
380 $ r, c, b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
381 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
383 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
385 CALL dgbsvxx(
'N',
'N', 2, 1, -1, 0, a, 1, af, 1, ip, eq,
386 $ r, c, b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
387 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
389 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
391 CALL dgbsvxx(
'N',
'N', 0, 1, 1, -1, a, 1, af, 1, ip, eq, r, c,
392 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
393 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
395 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
397 CALL dgbsvxx(
'N',
'N', 2, 1, 1, 1, a, 2, af, 2, ip, eq, r, c,
398 $ b, 2, x, 2, rcond, rpvgrw, berr, n_err_bnds,
399 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
401 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
403 CALL dgbsvxx(
'N',
'N', 2, 1, 1, 1, a, 3, af, 3, ip, eq, r, c,
404 $ b, 2, x, 2, rcond, rpvgrw, berr, n_err_bnds,
405 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
407 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
410 CALL dgbsvxx(
'F',
'N', 0, 1, 1, 0, a, 3, af, 4, ip, eq, r, c,
411 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
412 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
414 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
417 CALL dgbsvxx(
'F',
'N', 1, 1, 1, 0, a, 3, af, 4, ip, eq, r, c,
418 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
419 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
421 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
424 CALL dgbsvxx(
'F',
'N', 1, 1, 1, 0, a, 3, af, 4, ip, eq, r, c,
425 $ b, 1, x, 1, rcond, rpvgrw, berr, n_err_bnds,
426 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
428 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
430 CALL dgbsvxx(
'N',
'N', 2, 1, 1, 1, a, 3, af, 4, ip, eq, r, c,
431 $ b, 1, x, 2, rcond, rpvgrw, berr, n_err_bnds,
432 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
434 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
436 CALL dgbsvxx(
'N',
'N', 2, 1, 1, 1, a, 3, af, 4, ip, eq, r, c,
437 $ b, 2, x, 1, rcond, rpvgrw, berr, n_err_bnds,
438 $ err_bnds_n, err_bnds_c, nparams, params, w, iw,
440 CALL chkxer(
'DGBSVXX', infot, nout, lerr, ok )
442 ELSE IF(
lsamen( 2, c2,
'GT' ) )
THEN
448 CALL dgtsv( -1, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ), b, 1,
450 CALL chkxer(
'DGTSV ', infot, nout, lerr, ok )
452 CALL dgtsv( 0, -1, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ), b, 1,
454 CALL chkxer(
'DGTSV ', infot, nout, lerr, ok )
456 CALL dgtsv( 2, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ), b, 1, info )
457 CALL chkxer(
'DGTSV ', 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',
'/', 0, 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', -1, 0, 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', 0, -1, 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, 1, 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, 1, x, 2, rcond, r1, r2, w, iw, info )
486 CALL chkxer(
'DGTSVX', infot, nout, lerr, ok )
488 CALL dgtsvx(
'N',
'N', 2, 0, a( 1, 1 ), a( 1, 2 ), a( 1, 3 ),
489 $ af( 1, 1 ), af( 1, 2 ), af( 1, 3 ), af( 1, 4 ),
490 $ ip, b, 2, x, 1, rcond, r1, r2, w, iw, info )
491 CALL chkxer(
'DGTSVX', infot, nout, lerr, ok )
493 ELSE IF(
lsamen( 2, c2,
'PO' ) )
THEN
499 CALL dposv(
'/', 0, 0, a, 1, b, 1, info )
500 CALL chkxer(
'DPOSV ', infot, nout, lerr, ok )
502 CALL dposv(
'U', -1, 0, a, 1, b, 1, info )
503 CALL chkxer(
'DPOSV ', infot, nout, lerr, ok )
505 CALL dposv(
'U', 0, -1, a, 1, b, 1, info )
506 CALL chkxer(
'DPOSV ', infot, nout, lerr, ok )
508 CALL dposv(
'U', 2, 0, a, 1, b, 2, info )
509 CALL chkxer(
'DPOSV ', infot, nout, lerr, ok )
511 CALL dposv(
'U', 2, 0, a, 2, b, 1, info )
512 CALL chkxer(
'DPOSV ', infot, nout, lerr, ok )
518 CALL dposvx(
'/',
'U', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
519 $ rcond, r1, r2, w, iw, info )
520 CALL chkxer(
'DPOSVX', infot, nout, lerr, ok )
522 CALL dposvx(
'N',
'/', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
523 $ rcond, r1, r2, w, iw, info )
524 CALL chkxer(
'DPOSVX', infot, nout, lerr, ok )
526 CALL dposvx(
'N',
'U', -1, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
527 $ rcond, r1, r2, w, iw, info )
528 CALL chkxer(
'DPOSVX', infot, nout, lerr, ok )
530 CALL dposvx(
'N',
'U', 0, -1, a, 1, af, 1, eq, c, b, 1, x, 1,
531 $ rcond, r1, r2, w, iw, info )
532 CALL chkxer(
'DPOSVX', infot, nout, lerr, ok )
534 CALL dposvx(
'N',
'U', 2, 0, a, 1, af, 2, eq, c, b, 2, x, 2,
535 $ rcond, r1, r2, w, iw, info )
536 CALL chkxer(
'DPOSVX', infot, nout, lerr, ok )
538 CALL dposvx(
'N',
'U', 2, 0, a, 2, af, 1, eq, c, b, 2, x, 2,
539 $ rcond, r1, r2, w, iw, info )
540 CALL chkxer(
'DPOSVX', infot, nout, lerr, ok )
543 CALL dposvx(
'F',
'U', 0, 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 )
548 CALL dposvx(
'F',
'U', 1, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
549 $ rcond, r1, r2, w, iw, info )
550 CALL chkxer(
'DPOSVX', infot, nout, lerr, ok )
552 CALL dposvx(
'N',
'U', 2, 0, a, 2, af, 2, eq, c, b, 1, x, 2,
553 $ rcond, r1, r2, w, iw, info )
554 CALL chkxer(
'DPOSVX', infot, nout, lerr, ok )
556 CALL dposvx(
'N',
'U', 2, 0, a, 2, af, 2, eq, c, b, 2, x, 1,
557 $ rcond, r1, r2, w, iw, info )
558 CALL chkxer(
'DPOSVX', infot, nout, lerr, ok )
566 CALL dposvxx(
'/',
'U', 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',
'/', 0, 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', -1, 0, 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', 0, -1, a, 1, af, 1, eq, c, b, 1, x, 1,
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, 1, af, 2, 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 )
591 CALL dposvxx(
'N',
'U', 2, 0, a, 2, af, 1, eq, c, b, 2, x, 2,
592 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
593 $ err_bnds_c, nparams, params, w, iw, info )
594 CALL chkxer(
'DPOSVXX', infot, nout, lerr, ok )
597 CALL dposvxx(
'F',
'U', 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
598 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
599 $ err_bnds_c, nparams, params, w, iw, info )
600 CALL chkxer(
'DPOSVXX', infot, nout, lerr, ok )
603 CALL dposvxx(
'F',
'U', 1, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
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, 1, x, 2,
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 CALL dposvxx(
'N',
'U', 2, 0, a, 2, af, 2, eq, c, b, 2, x, 1,
614 $ rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
615 $ err_bnds_c, nparams, params, w, iw, info )
616 CALL chkxer(
'DPOSVXX', infot, nout, lerr, ok )
618 ELSE IF(
lsamen( 2, c2,
'PP' ) )
THEN
624 CALL dppsv(
'/', 0, 0, a, b, 1, info )
625 CALL chkxer(
'DPPSV ', infot, nout, lerr, ok )
627 CALL dppsv(
'U', -1, 0, a, b, 1, info )
628 CALL chkxer(
'DPPSV ', infot, nout, lerr, ok )
630 CALL dppsv(
'U', 0, -1, a, b, 1, info )
631 CALL chkxer(
'DPPSV ', infot, nout, lerr, ok )
633 CALL dppsv(
'U', 2, 0, a, b, 1, info )
634 CALL chkxer(
'DPPSV ', infot, nout, lerr, ok )
640 CALL dppsvx(
'/',
'U', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
641 $ r1, r2, w, iw, info )
642 CALL chkxer(
'DPPSVX', infot, nout, lerr, ok )
644 CALL dppsvx(
'N',
'/', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
645 $ r1, r2, w, iw, info )
646 CALL chkxer(
'DPPSVX', infot, nout, lerr, ok )
648 CALL dppsvx(
'N',
'U', -1, 0, a, af, eq, c, b, 1, x, 1, rcond,
649 $ r1, r2, w, iw, info )
650 CALL chkxer(
'DPPSVX', infot, nout, lerr, ok )
652 CALL dppsvx(
'N',
'U', 0, -1, 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', 0, 0, a, af, eq, c, b, 1, x, 1, rcond,
658 $ r1, r2, w, iw, info )
659 CALL chkxer(
'DPPSVX', infot, nout, lerr, ok )
662 CALL dppsvx(
'F',
'U', 1, 0, a, af, eq, c, b, 1, x, 1, rcond,
663 $ r1, r2, w, iw, info )
664 CALL chkxer(
'DPPSVX', infot, nout, lerr, ok )
666 CALL dppsvx(
'N',
'U', 2, 0, a, af, eq, c, b, 1, x, 2, rcond,
667 $ r1, r2, w, iw, info )
668 CALL chkxer(
'DPPSVX', infot, nout, lerr, ok )
670 CALL dppsvx(
'N',
'U', 2, 0, a, af, eq, c, b, 2, x, 1, rcond,
671 $ r1, r2, w, iw, info )
672 CALL chkxer(
'DPPSVX', infot, nout, lerr, ok )
674 ELSE IF(
lsamen( 2, c2,
'PB' ) )
THEN
680 CALL dpbsv(
'/', 0, 0, 0, a, 1, b, 1, info )
681 CALL chkxer(
'DPBSV ', infot, nout, lerr, ok )
683 CALL dpbsv(
'U', -1, 0, 0, a, 1, b, 1, info )
684 CALL chkxer(
'DPBSV ', infot, nout, lerr, ok )
686 CALL dpbsv(
'U', 1, -1, 0, a, 1, b, 1, info )
687 CALL chkxer(
'DPBSV ', infot, nout, lerr, ok )
689 CALL dpbsv(
'U', 0, 0, -1, a, 1, b, 1, info )
690 CALL chkxer(
'DPBSV ', infot, nout, lerr, ok )
692 CALL dpbsv(
'U', 1, 1, 0, a, 1, b, 2, info )
693 CALL chkxer(
'DPBSV ', infot, nout, lerr, ok )
695 CALL dpbsv(
'U', 2, 0, 0, a, 1, b, 1, info )
696 CALL chkxer(
'DPBSV ', infot, nout, lerr, ok )
702 CALL dpbsvx(
'/',
'U', 0, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
703 $ rcond, r1, r2, w, iw, info )
704 CALL chkxer(
'DPBSVX', infot, nout, lerr, ok )
706 CALL dpbsvx(
'N',
'/', 0, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
707 $ rcond, r1, r2, w, iw, info )
708 CALL chkxer(
'DPBSVX', infot, nout, lerr, ok )
710 CALL dpbsvx(
'N',
'U', -1, 0, 0, a, 1, af, 1, eq, c, b, 1, x,
711 $ 1, rcond, r1, r2, w, iw, info )
712 CALL chkxer(
'DPBSVX', infot, nout, lerr, ok )
714 CALL dpbsvx(
'N',
'U', 1, -1, 0, a, 1, af, 1, eq, c, b, 1, x,
715 $ 1, rcond, r1, r2, w, iw, info )
716 CALL chkxer(
'DPBSVX', infot, nout, lerr, ok )
718 CALL dpbsvx(
'N',
'U', 0, 0, -1, a, 1, af, 1, eq, c, b, 1, x,
719 $ 1, rcond, r1, r2, w, iw, info )
720 CALL chkxer(
'DPBSVX', infot, nout, lerr, ok )
722 CALL dpbsvx(
'N',
'U', 1, 1, 0, a, 1, af, 2, eq, c, b, 2, x, 2,
723 $ rcond, r1, r2, w, iw, info )
724 CALL chkxer(
'DPBSVX', infot, nout, lerr, ok )
726 CALL dpbsvx(
'N',
'U', 1, 1, 0, a, 2, af, 1, eq, c, b, 2, x, 2,
727 $ rcond, r1, r2, w, iw, info )
728 CALL chkxer(
'DPBSVX', infot, nout, lerr, ok )
731 CALL dpbsvx(
'F',
'U', 0, 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 )
736 CALL dpbsvx(
'F',
'U', 1, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 1,
737 $ rcond, r1, r2, w, iw, info )
738 CALL chkxer(
'DPBSVX', infot, nout, lerr, ok )
740 CALL dpbsvx(
'N',
'U', 2, 0, 0, a, 1, af, 1, eq, c, b, 1, x, 2,
741 $ rcond, r1, r2, w, iw, info )
742 CALL chkxer(
'DPBSVX', infot, nout, lerr, ok )
744 CALL dpbsvx(
'N',
'U', 2, 0, 0, a, 1, af, 1, eq, c, b, 2, x, 1,
745 $ rcond, r1, r2, w, iw, info )
746 CALL chkxer(
'DPBSVX', infot, nout, lerr, ok )
748 ELSE IF(
lsamen( 2, c2,
'PT' ) )
THEN
754 CALL dptsv( -1, 0, a( 1, 1 ), a( 1, 2 ), b, 1, info )
755 CALL chkxer(
'DPTSV ', infot, nout, lerr, ok )
757 CALL dptsv( 0, -1, a( 1, 1 ), a( 1, 2 ), b, 1, info )
758 CALL chkxer(
'DPTSV ', infot, nout, lerr, ok )
760 CALL dptsv( 2, 0, a( 1, 1 ), a( 1, 2 ), b, 1, info )
761 CALL chkxer(
'DPTSV ', infot, nout, lerr, ok )
767 CALL dptsvx(
'/', 0, 0, a( 1, 1 ), a( 1, 2 ), af( 1, 1 ),
768 $ af( 1, 2 ), b, 1, x, 1, rcond, r1, r2, w, info )
769 CALL chkxer(
'DPTSVX', infot, nout, lerr, ok )
771 CALL dptsvx(
'N', -1, 0, a( 1, 1 ), a( 1, 2 ), af( 1, 1 ),
772 $ af( 1, 2 ), b, 1, x, 1, rcond, r1, r2, w, info )
773 CALL chkxer(
'DPTSVX', infot, nout, lerr, ok )
775 CALL dptsvx(
'N', 0, -1, a( 1, 1 ), a( 1, 2 ), af( 1, 1 ),
776 $ af( 1, 2 ), b, 1, x, 1, rcond, r1, r2, w, info )
777 CALL chkxer(
'DPTSVX', infot, nout, lerr, ok )
779 CALL dptsvx(
'N', 2, 0, a( 1, 1 ), a( 1, 2 ), af( 1, 1 ),
780 $ af( 1, 2 ), b, 1, x, 2, rcond, r1, r2, w, info )
781 CALL chkxer(
'DPTSVX', infot, nout, lerr, ok )
783 CALL dptsvx(
'N', 2, 0, a( 1, 1 ), a( 1, 2 ), af( 1, 1 ),
784 $ af( 1, 2 ), b, 2, x, 1, rcond, r1, r2, w, info )
785 CALL chkxer(
'DPTSVX', infot, nout, lerr, ok )
787 ELSE IF(
lsamen( 2, c2,
'SY' ) )
THEN
793 CALL dsysv(
'/', 0, 0, a, 1, ip, b, 1, w, 1, info )
794 CALL chkxer(
'DSYSV ', infot, nout, lerr, ok )
796 CALL dsysv(
'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
797 CALL chkxer(
'DSYSV ', infot, nout, lerr, ok )
799 CALL dsysv(
'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
800 CALL chkxer(
'DSYSV ', infot, nout, lerr, ok )
802 CALL dsysv(
'U', 2, 0, a, 1, ip, b, 2, w, 1, info )
803 CALL chkxer(
'DSYSV_ROOK', infot, nout, lerr, ok )
805 CALL dsysv(
'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
806 CALL chkxer(
'DSYSV ', infot, nout, lerr, ok )
808 CALL dsysv(
'U', 0, 0, a, 1, ip, b, 1, w, 0, info )
809 CALL chkxer(
'DSYSV ', infot, nout, lerr, ok )
811 CALL dsysv(
'U', 0, 0, a, 1, ip, b, 1, w, -2, info )
812 CALL chkxer(
'DSYSV ', infot, nout, lerr, ok )
818 CALL dsysvx(
'/',
'U', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
819 $ rcond, r1, r2, w, 1, iw, info )
820 CALL chkxer(
'DSYSVX', infot, nout, lerr, ok )
822 CALL dsysvx(
'N',
'/', 0, 0, a, 1, af, 1, ip, b, 1, x, 1,
823 $ rcond, r1, r2, w, 1, iw, info )
824 CALL chkxer(
'DSYSVX', infot, nout, lerr, ok )
826 CALL dsysvx(
'N',
'U', -1, 0, a, 1, af, 1, ip, b, 1, x, 1,
827 $ rcond, r1, r2, w, 1, iw, info )
828 CALL chkxer(
'DSYSVX', infot, nout, lerr, ok )
830 CALL dsysvx(
'N',
'U', 0, -1, a, 1, af, 1, ip, b, 1, x, 1,
831 $ rcond, r1, r2, w, 1, iw, info )
832 CALL chkxer(
'DSYSVX', infot, nout, lerr, ok )
834 CALL dsysvx(
'N',
'U', 2, 0, a, 1, af, 2, ip, b, 2, x, 2,
835 $ rcond, r1, r2, w, 4, iw, info )
836 CALL chkxer(
'DSYSVX', infot, nout, lerr, ok )
838 CALL dsysvx(
'N',
'U', 2, 0, a, 2, af, 1, ip, b, 2, x, 2,
839 $ rcond, r1, r2, w, 4, iw, info )
840 CALL chkxer(
'DSYSVX', infot, nout, lerr, ok )
842 CALL dsysvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 1, x, 2,
843 $ rcond, r1, r2, w, 4, iw, info )
844 CALL chkxer(
'DSYSVX', infot, nout, lerr, ok )
846 CALL dsysvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 1,
847 $ rcond, r1, r2, w, 4, iw, info )
848 CALL chkxer(
'DSYSVX', infot, nout, lerr, ok )
850 CALL dsysvx(
'N',
'U', 2, 0, a, 2, af, 2, ip, b, 2, x, 2,
851 $ rcond, r1, r2, w, 3, iw, info )
852 CALL chkxer(
'DSYSVX', infot, nout, lerr, ok )
861 CALL dsysvxx(
'/',
'U', 0, 0, a, 1, af, 1, ip, eq, r, b, 1, x,
862 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
863 $ err_bnds_c, nparams, params, w, iw, info )
864 CALL chkxer(
'DSYSVXX', infot, nout, lerr, ok )
866 CALL dsysvxx(
'N',
'/', 0, 0, a, 1, af, 1, ip, eq, r, b, 1, x,
867 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
868 $ err_bnds_c, nparams, params, w, iw, info )
869 CALL chkxer(
'DSYSVXX', infot, nout, lerr, ok )
871 CALL dsysvxx(
'N',
'U', -1, 0, a, 1, af, 1, ip, eq, r, b, 1, x,
872 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
873 $ err_bnds_c, nparams, params, w, iw, info )
874 CALL chkxer(
'DSYSVXX', infot, nout, lerr, ok )
877 CALL dsysvxx(
'N',
'U', 0, -1, a, 1, af, 1, ip, eq, r, b, 1, x,
878 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
879 $ err_bnds_c, nparams, params, w, iw, info )
880 CALL chkxer(
'DSYSVXX', infot, nout, lerr, ok )
883 CALL dsysvxx(
'N',
'U', 2, 0, a, 1, af, 2, ip, eq, r, b, 2, x,
884 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
885 $ err_bnds_c, nparams, params, w, iw, info )
886 CALL chkxer(
'DSYSVXX', infot, nout, lerr, ok )
888 CALL dsysvxx(
'N',
'U', 2, 0, a, 2, af, 1, ip, eq, r, b, 2, x,
889 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
890 $ err_bnds_c, nparams, params, w, iw, info )
891 CALL chkxer(
'DSYSVXX', infot, nout, lerr, ok )
893 CALL dsysvxx(
'F',
'U', 2, 0, a, 2, af, 2, ip,
'A', r, b, 2, x,
894 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
895 $ err_bnds_c, nparams, params, w, iw, info )
896 CALL chkxer(
'DSYSVXX', infot, nout, lerr, ok )
899 CALL dsysvxx(
'F',
'U', 2, 0, a, 2, af, 2, ip, eq, r, b, 2, x,
900 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
901 $ err_bnds_c, nparams, params, w, iw, info )
902 CALL chkxer(
'DSYSVXX', infot, nout, lerr, ok )
906 CALL dsysvxx(
'F',
'U', 2, 0, a, 2, af, 2, ip, eq, r, b, 2, x,
907 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
908 $ err_bnds_c, nparams, params, w, iw, info )
909 CALL chkxer(
'DSYSVXX', infot, nout, lerr, ok )
912 CALL dsysvxx(
'N',
'U', 2, 0, a, 2, af, 2, ip, eq, r, b, 1, x,
913 $ 2, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
914 $ err_bnds_c, nparams, params, w, iw, info )
915 CALL chkxer(
'DSYSVXX', infot, nout, lerr, ok )
917 CALL dsysvxx(
'N',
'U', 2, 0, a, 2, af, 2, ip, eq, r, b, 2, x,
918 $ 1, rcond, rpvgrw, berr, n_err_bnds, err_bnds_n,
919 $ err_bnds_c, nparams, params, w, iw, info )
920 CALL chkxer(
'DSYSVXX', infot, nout, lerr, ok )
922 ELSE IF(
lsamen( 2, c2,
'SR' ) )
THEN
926 srnamt =
'DSYSV_ROOK'
928 CALL dsysv_rook(
'/', 0, 0, a, 1, ip, b, 1, w, 1, info )
929 CALL chkxer(
'DSYSV_ROOK', infot, nout, lerr, ok )
931 CALL dsysv_rook(
'U', -1, 0, a, 1, ip, b, 1, w, 1, info )
932 CALL chkxer(
'DSYSV_ROOK', infot, nout, lerr, ok )
934 CALL dsysv_rook(
'U', 0, -1, a, 1, ip, b, 1, w, 1, info )
935 CALL chkxer(
'DSYSV_ROOK', infot, nout, lerr, ok )
937 CALL dsysv_rook(
'U', 2, 0, a, 1, ip, b, 2, w, 1, info )
938 CALL chkxer(
'DSYSV_ROOK', infot, nout, lerr, ok )
940 CALL dsysv_rook(
'U', 2, 0, a, 2, ip, b, 1, w, 1, info )
941 CALL chkxer(
'DSYSV_ROOK', infot, nout, lerr, ok )
943 CALL dsysv_rook(
'U', 0, 0, a, 1, ip, b, 1, w, 0, info )
944 CALL chkxer(
'DSYSV_ROOK', infot, nout, lerr, ok )
946 CALL dsysv_rook(
'U', 0, 0, a, 1, ip, b, 1, w, -2, info )
947 CALL chkxer(
'DSYSV_ROOK', infot, nout, lerr, ok )
949 ELSE IF(
lsamen( 2, c2,
'SK' ) )
THEN
963 CALL dsysv_rk(
'/', 0, 0, a, 1, e, ip, b, 1, w, 1, info )
964 CALL chkxer(
'DSYSV_RK', infot, nout, lerr, ok )
966 CALL dsysv_rk(
'U', -1, 0, a, 1, e, ip, b, 1, w, 1, info )
967 CALL chkxer(
'DSYSV_RK', infot, nout, lerr, ok )
969 CALL dsysv_rk(
'U', 0, -1, a, 1, e, ip, b, 1, w, 1, info )
970 CALL chkxer(
'DSYSV_RK', infot, nout, lerr, ok )
972 CALL dsysv_rk(
'U', 2, 0, a, 1, e, ip, b, 2, w, 1, info )
973 CALL chkxer(
'DSYSV_RK', infot, nout, lerr, ok )
975 CALL dsysv_rk(
'U', 2, 0, a, 2, e, ip, b, 1, w, 1, info )
976 CALL chkxer(
'DSYSV_RK', infot, nout, lerr, ok )
978 CALL dsysv_rk(
'U', 0, 0, a, 1, e, ip, b, 1, w, 0, info )
979 CALL chkxer(
'DSYSV_RK', infot, nout, lerr, ok )
981 CALL dsysv_rk(
'U', 0, 0, a, 1, e, ip, b, 1, w, -2, info )
982 CALL chkxer(
'DSYSV_RK', infot, nout, lerr, ok )
984 ELSE IF(
lsamen( 2, c2,
'SP' ) )
THEN
990 CALL dspsv(
'/', 0, 0, a, ip, b, 1, info )
991 CALL chkxer(
'DSPSV ', infot, nout, lerr, ok )
993 CALL dspsv(
'U', -1, 0, a, ip, b, 1, info )
994 CALL chkxer(
'DSPSV ', infot, nout, lerr, ok )
996 CALL dspsv(
'U', 0, -1, a, ip, b, 1, info )
997 CALL chkxer(
'DSPSV ', infot, nout, lerr, ok )
999 CALL dspsv(
'U', 2, 0, a, ip, b, 1, info )
1000 CALL chkxer(
'DSPSV ', infot, nout, lerr, ok )
1006 CALL dspsvx(
'/',
'U', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
1008 CALL chkxer(
'DSPSVX', infot, nout, lerr, ok )
1010 CALL dspsvx(
'N',
'/', 0, 0, a, af, ip, b, 1, x, 1, rcond, r1,
1012 CALL chkxer(
'DSPSVX', infot, nout, lerr, ok )
1014 CALL dspsvx(
'N',
'U', -1, 0, a, af, ip, b, 1, x, 1, rcond, r1,
1016 CALL chkxer(
'DSPSVX', infot, nout, lerr, ok )
1018 CALL dspsvx(
'N',
'U', 0, -1, a, af, ip, b, 1, x, 1, rcond, r1,
1020 CALL chkxer(
'DSPSVX', infot, nout, lerr, ok )
1022 CALL dspsvx(
'N',
'U', 2, 0, a, af, ip, b, 1, x, 2, rcond, r1,
1024 CALL chkxer(
'DSPSVX', infot, nout, lerr, ok )
1026 CALL dspsvx(
'N',
'U', 2, 0, a, af, ip, b, 2, x, 1, rcond, r1,
1028 CALL chkxer(
'DSPSVX', infot, nout, lerr, ok )
1034 WRITE( nout, fmt = 9999 )path
1036 WRITE( nout, fmt = 9998 )path
1039 9999
FORMAT( 1x, a3,
' drivers passed the tests of the error exits' )
1040 9998
FORMAT(
' *** ', a3,
' drivers failed the tests of the error ',
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 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 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 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 chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
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 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 ...
logical function lsamen(N, CA, CB)
LSAMEN
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 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 derrvx(PATH, NUNIT)
DERRVX
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 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 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 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 dsysv_rook(UPLO, N, NRHS, A, LDA, IPIV, B, LDB, WORK, LWORK, INFO)
DSYSV_ROOK 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 dppsv(UPLO, N, NRHS, AP, B, LDB, INFO)
DPPSV computes the solution to system of linear equations A * X = B for OTHER 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 dsysv_rk(UPLO, N, NRHS, A, LDA, E, IPIV, B, LDB, WORK, LWORK, INFO)
DSYSV_RK computes the solution to system of linear equations A * X = B for SY matrices ...
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 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 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 dptsv(N, NRHS, D, E, B, LDB, INFO)
DPTSV computes the solution to system of linear equations A * X = B for PT 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 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 ...