swiftnav: derrpo.c Source File

Go to the documentation of this file.
00001 /* derrpo.f -- translated by f2c (version 20061008).
00002    You must link the resulting object file with libf2c:
00003         on Microsoft Windows system, link with libf2c.lib;
00004         on Linux or Unix systems, link with .../path/to/libf2c.a -lm
00005         or, if you install libf2c.a in a standard place, with -lf2c -lm
00006         -- in that order, at the end of the command line, as in
00007                 cc *.o -lf2c -lm
00008         Source for libf2c is in /netlib/f2c/libf2c.zip, e.g.,
00009 
00010                 http://www.netlib.org/f2c/libf2c.zip
00011 */
00012 
00013 #include "f2c.h"
00014 #include "blaswrap.h"
00015 
00016 /* Common Block Declarations */
00017 
00018 struct {
00019     integer infot, nout;
00020     logical ok, lerr;
00021 } infoc_;
00022 
00023 #define infoc_1 infoc_
00024 
00025 struct {
00026     char srnamt[32];
00027 } srnamc_;
00028 
00029 #define srnamc_1 srnamc_
00030 
00031 /* Table of constant values */
00032 
00033 static integer c__2 = 2;
00034 static integer c__0 = 0;
00035 static integer c__1 = 1;
00036 static integer c_n1 = -1;
00037 
00038 /* Subroutine */ int derrpo_(char *path, integer *nunit)
00039 {
00040     /* Builtin functions */
00041     integer s_wsle(cilist *), e_wsle(void);
00042     /* Subroutine */ int s_copy(char *, char *, ftnlen, ftnlen);
00043 
00044     /* Local variables */
00045     doublereal a[16]    /* was [4][4] */, b[4];
00046     integer i__, j;
00047     doublereal w[12], x[4];
00048     char c2[2];
00049     doublereal r1[4], r2[4], af[16]     /* was [4][4] */;
00050     integer iw[4], info;
00051     doublereal anrm, rcond;
00052     extern /* Subroutine */ int dpbtf2_(char *, integer *, integer *, 
00053             doublereal *, integer *, integer *), dpotf2_(char *, 
00054             integer *, doublereal *, integer *, integer *), alaesm_(
00055             char *, logical *, integer *), dpbcon_(char *, integer *, 
00056             integer *, doublereal *, integer *, doublereal *, doublereal *, 
00057             doublereal *, integer *, integer *);
00058     extern logical lsamen_(integer *, char *, char *);
00059     extern /* Subroutine */ int dpbequ_(char *, integer *, integer *, 
00060             doublereal *, integer *, doublereal *, doublereal *, doublereal *, 
00061              integer *), dpbrfs_(char *, integer *, integer *, 
00062             integer *, doublereal *, integer *, doublereal *, integer *, 
00063             doublereal *, integer *, doublereal *, integer *, doublereal *, 
00064             doublereal *, doublereal *, integer *, integer *), 
00065             dpbtrf_(char *, integer *, integer *, doublereal *, integer *, 
00066             integer *), dpocon_(char *, integer *, doublereal *, 
00067             integer *, doublereal *, doublereal *, doublereal *, integer *, 
00068             integer *), chkxer_(char *, integer *, integer *, logical 
00069             *, logical *), dppcon_(char *, integer *, doublereal *, 
00070             doublereal *, doublereal *, doublereal *, integer *, integer *), dpoequ_(integer *, doublereal *, integer *, doublereal *, 
00071              doublereal *, doublereal *, integer *), dpbtrs_(char *, integer *
00072 , integer *, integer *, doublereal *, integer *, doublereal *, 
00073             integer *, integer *), dporfs_(char *, integer *, integer 
00074             *, doublereal *, integer *, doublereal *, integer *, doublereal *, 
00075              integer *, doublereal *, integer *, doublereal *, doublereal *, 
00076             doublereal *, integer *, integer *), dpotrf_(char *, 
00077             integer *, doublereal *, integer *, integer *), dpotri_(
00078             char *, integer *, doublereal *, integer *, integer *), 
00079             dppequ_(char *, integer *, doublereal *, doublereal *, doublereal 
00080             *, doublereal *, integer *), dpprfs_(char *, integer *, 
00081             integer *, doublereal *, doublereal *, doublereal *, integer *, 
00082             doublereal *, integer *, doublereal *, doublereal *, doublereal *, 
00083              integer *, integer *), dpptrf_(char *, integer *, 
00084             doublereal *, integer *), dpptri_(char *, integer *, 
00085             doublereal *, integer *), dpotrs_(char *, integer *, 
00086             integer *, doublereal *, integer *, doublereal *, integer *, 
00087             integer *), dpptrs_(char *, integer *, integer *, 
00088             doublereal *, doublereal *, integer *, integer *);
00089 
00090     /* Fortran I/O blocks */
00091     static cilist io___1 = { 0, 0, 0, 0, 0 };
00092 
00093 
00094 
00095 /*  -- LAPACK test routine (version 3.1) -- */
00096 /*     Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. */
00097 /*     November 2006 */
00098 
00099 /*     .. Scalar Arguments .. */
00100 /*     .. */
00101 
00102 /*  Purpose */
00103 /*  ======= */
00104 
00105 /*  DERRPO tests the error exits for the DOUBLE PRECISION routines */
00106 /*  for symmetric positive definite matrices. */
00107 
00108 /*  Arguments */
00109 /*  ========= */
00110 
00111 /*  PATH    (input) CHARACTER*3 */
00112 /*          The LAPACK path name for the routines to be tested. */
00113 
00114 /*  NUNIT   (input) INTEGER */
00115 /*          The unit number for output. */
00116 
00117 /*  ===================================================================== */
00118 
00119 /*     .. Parameters .. */
00120 /*     .. */
00121 /*     .. Local Scalars .. */
00122 /*     .. */
00123 /*     .. Local Arrays .. */
00124 /*     .. */
00125 /*     .. External Functions .. */
00126 /*     .. */
00127 /*     .. External Subroutines .. */
00128 /*     .. */
00129 /*     .. Scalars in Common .. */
00130 /*     .. */
00131 /*     .. Common blocks .. */
00132 /*     .. */
00133 /*     .. Intrinsic Functions .. */
00134 /*     .. */
00135 /*     .. Executable Statements .. */
00136 
00137     infoc_1.nout = *nunit;
00138     io___1.ciunit = infoc_1.nout;
00139     s_wsle(&io___1);
00140     e_wsle();
00141     s_copy(c2, path + 1, (ftnlen)2, (ftnlen)2);
00142 
00143 /*     Set the variables to innocuous values. */
00144 
00145     for (j = 1; j <= 4; ++j) {
00146         for (i__ = 1; i__ <= 4; ++i__) {
00147             a[i__ + (j << 2) - 5] = 1. / (doublereal) (i__ + j);
00148             af[i__ + (j << 2) - 5] = 1. / (doublereal) (i__ + j);
00149 /* L10: */
00150         }
00151         b[j - 1] = 0.;
00152         r1[j - 1] = 0.;
00153         r2[j - 1] = 0.;
00154         w[j - 1] = 0.;
00155         x[j - 1] = 0.;
00156         iw[j - 1] = j;
00157 /* L20: */
00158     }
00159     infoc_1.ok = TRUE_;
00160 
00161     if (lsamen_(&c__2, c2, "PO")) {
00162 
00163 /*        Test error exits of the routines that use the Cholesky */
00164 /*        decomposition of a symmetric positive definite matrix. */
00165 
00166 /*        DPOTRF */
00167 
00168         s_copy(srnamc_1.srnamt, "DPOTRF", (ftnlen)32, (ftnlen)6);
00169         infoc_1.infot = 1;
00170         dpotrf_("/", &c__0, a, &c__1, &info);
00171         chkxer_("DPOTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00172                 infoc_1.ok);
00173         infoc_1.infot = 2;
00174         dpotrf_("U", &c_n1, a, &c__1, &info);
00175         chkxer_("DPOTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00176                 infoc_1.ok);
00177         infoc_1.infot = 4;
00178         dpotrf_("U", &c__2, a, &c__1, &info);
00179         chkxer_("DPOTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00180                 infoc_1.ok);
00181 
00182 /*        DPOTF2 */
00183 
00184         s_copy(srnamc_1.srnamt, "DPOTF2", (ftnlen)32, (ftnlen)6);
00185         infoc_1.infot = 1;
00186         dpotf2_("/", &c__0, a, &c__1, &info);
00187         chkxer_("DPOTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00188                 infoc_1.ok);
00189         infoc_1.infot = 2;
00190         dpotf2_("U", &c_n1, a, &c__1, &info);
00191         chkxer_("DPOTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00192                 infoc_1.ok);
00193         infoc_1.infot = 4;
00194         dpotf2_("U", &c__2, a, &c__1, &info);
00195         chkxer_("DPOTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00196                 infoc_1.ok);
00197 
00198 /*        DPOTRI */
00199 
00200         s_copy(srnamc_1.srnamt, "DPOTRI", (ftnlen)32, (ftnlen)6);
00201         infoc_1.infot = 1;
00202         dpotri_("/", &c__0, a, &c__1, &info);
00203         chkxer_("DPOTRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00204                 infoc_1.ok);
00205         infoc_1.infot = 2;
00206         dpotri_("U", &c_n1, a, &c__1, &info);
00207         chkxer_("DPOTRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00208                 infoc_1.ok);
00209         infoc_1.infot = 4;
00210         dpotri_("U", &c__2, a, &c__1, &info);
00211         chkxer_("DPOTRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00212                 infoc_1.ok);
00213 
00214 /*        DPOTRS */
00215 
00216         s_copy(srnamc_1.srnamt, "DPOTRS", (ftnlen)32, (ftnlen)6);
00217         infoc_1.infot = 1;
00218         dpotrs_("/", &c__0, &c__0, a, &c__1, b, &c__1, &info);
00219         chkxer_("DPOTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00220                 infoc_1.ok);
00221         infoc_1.infot = 2;
00222         dpotrs_("U", &c_n1, &c__0, a, &c__1, b, &c__1, &info);
00223         chkxer_("DPOTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00224                 infoc_1.ok);
00225         infoc_1.infot = 3;
00226         dpotrs_("U", &c__0, &c_n1, a, &c__1, b, &c__1, &info);
00227         chkxer_("DPOTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00228                 infoc_1.ok);
00229         infoc_1.infot = 5;
00230         dpotrs_("U", &c__2, &c__1, a, &c__1, b, &c__2, &info);
00231         chkxer_("DPOTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00232                 infoc_1.ok);
00233         infoc_1.infot = 7;
00234         dpotrs_("U", &c__2, &c__1, a, &c__2, b, &c__1, &info);
00235         chkxer_("DPOTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00236                 infoc_1.ok);
00237 
00238 /*        DPORFS */
00239 
00240         s_copy(srnamc_1.srnamt, "DPORFS", (ftnlen)32, (ftnlen)6);
00241         infoc_1.infot = 1;
00242         dporfs_("/", &c__0, &c__0, a, &c__1, af, &c__1, b, &c__1, x, &c__1, 
00243                 r1, r2, w, iw, &info);
00244         chkxer_("DPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00245                 infoc_1.ok);
00246         infoc_1.infot = 2;
00247         dporfs_("U", &c_n1, &c__0, a, &c__1, af, &c__1, b, &c__1, x, &c__1, 
00248                 r1, r2, w, iw, &info);
00249         chkxer_("DPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00250                 infoc_1.ok);
00251         infoc_1.infot = 3;
00252         dporfs_("U", &c__0, &c_n1, a, &c__1, af, &c__1, b, &c__1, x, &c__1, 
00253                 r1, r2, w, iw, &info);
00254         chkxer_("DPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00255                 infoc_1.ok);
00256         infoc_1.infot = 5;
00257         dporfs_("U", &c__2, &c__1, a, &c__1, af, &c__2, b, &c__2, x, &c__2, 
00258                 r1, r2, w, iw, &info);
00259         chkxer_("DPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00260                 infoc_1.ok);
00261         infoc_1.infot = 7;
00262         dporfs_("U", &c__2, &c__1, a, &c__2, af, &c__1, b, &c__2, x, &c__2, 
00263                 r1, r2, w, iw, &info);
00264         chkxer_("DPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00265                 infoc_1.ok);
00266         infoc_1.infot = 9;
00267         dporfs_("U", &c__2, &c__1, a, &c__2, af, &c__2, b, &c__1, x, &c__2, 
00268                 r1, r2, w, iw, &info);
00269         chkxer_("DPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00270                 infoc_1.ok);
00271         infoc_1.infot = 11;
00272         dporfs_("U", &c__2, &c__1, a, &c__2, af, &c__2, b, &c__2, x, &c__1, 
00273                 r1, r2, w, iw, &info);
00274         chkxer_("DPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00275                 infoc_1.ok);
00276 
00277 /*        DPOCON */
00278 
00279         s_copy(srnamc_1.srnamt, "DPOCON", (ftnlen)32, (ftnlen)6);
00280         infoc_1.infot = 1;
00281         dpocon_("/", &c__0, a, &c__1, &anrm, &rcond, w, iw, &info);
00282         chkxer_("DPOCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00283                 infoc_1.ok);
00284         infoc_1.infot = 2;
00285         dpocon_("U", &c_n1, a, &c__1, &anrm, &rcond, w, iw, &info);
00286         chkxer_("DPOCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00287                 infoc_1.ok);
00288         infoc_1.infot = 4;
00289         dpocon_("U", &c__2, a, &c__1, &anrm, &rcond, w, iw, &info);
00290         chkxer_("DPOCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00291                 infoc_1.ok);
00292 
00293 /*        DPOEQU */
00294 
00295         s_copy(srnamc_1.srnamt, "DPOEQU", (ftnlen)32, (ftnlen)6);
00296         infoc_1.infot = 1;
00297         dpoequ_(&c_n1, a, &c__1, r1, &rcond, &anrm, &info);
00298         chkxer_("DPOEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00299                 infoc_1.ok);
00300         infoc_1.infot = 3;
00301         dpoequ_(&c__2, a, &c__1, r1, &rcond, &anrm, &info);
00302         chkxer_("DPOEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00303                 infoc_1.ok);
00304 
00305     } else if (lsamen_(&c__2, c2, "PP")) {
00306 
00307 /*        Test error exits of the routines that use the Cholesky */
00308 /*        decomposition of a symmetric positive definite packed matrix. */
00309 
00310 /*        DPPTRF */
00311 
00312         s_copy(srnamc_1.srnamt, "DPPTRF", (ftnlen)32, (ftnlen)6);
00313         infoc_1.infot = 1;
00314         dpptrf_("/", &c__0, a, &info);
00315         chkxer_("DPPTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00316                 infoc_1.ok);
00317         infoc_1.infot = 2;
00318         dpptrf_("U", &c_n1, a, &info);
00319         chkxer_("DPPTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00320                 infoc_1.ok);
00321 
00322 /*        DPPTRI */
00323 
00324         s_copy(srnamc_1.srnamt, "DPPTRI", (ftnlen)32, (ftnlen)6);
00325         infoc_1.infot = 1;
00326         dpptri_("/", &c__0, a, &info);
00327         chkxer_("DPPTRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00328                 infoc_1.ok);
00329         infoc_1.infot = 2;
00330         dpptri_("U", &c_n1, a, &info);
00331         chkxer_("DPPTRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00332                 infoc_1.ok);
00333 
00334 /*        DPPTRS */
00335 
00336         s_copy(srnamc_1.srnamt, "DPPTRS", (ftnlen)32, (ftnlen)6);
00337         infoc_1.infot = 1;
00338         dpptrs_("/", &c__0, &c__0, a, b, &c__1, &info);
00339         chkxer_("DPPTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00340                 infoc_1.ok);
00341         infoc_1.infot = 2;
00342         dpptrs_("U", &c_n1, &c__0, a, b, &c__1, &info);
00343         chkxer_("DPPTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00344                 infoc_1.ok);
00345         infoc_1.infot = 3;
00346         dpptrs_("U", &c__0, &c_n1, a, b, &c__1, &info);
00347         chkxer_("DPPTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00348                 infoc_1.ok);
00349         infoc_1.infot = 6;
00350         dpptrs_("U", &c__2, &c__1, a, b, &c__1, &info);
00351         chkxer_("DPPTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00352                 infoc_1.ok);
00353 
00354 /*        DPPRFS */
00355 
00356         s_copy(srnamc_1.srnamt, "DPPRFS", (ftnlen)32, (ftnlen)6);
00357         infoc_1.infot = 1;
00358         dpprfs_("/", &c__0, &c__0, a, af, b, &c__1, x, &c__1, r1, r2, w, iw, &
00359                 info);
00360         chkxer_("DPPRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00361                 infoc_1.ok);
00362         infoc_1.infot = 2;
00363         dpprfs_("U", &c_n1, &c__0, a, af, b, &c__1, x, &c__1, r1, r2, w, iw, &
00364                 info);
00365         chkxer_("DPPRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00366                 infoc_1.ok);
00367         infoc_1.infot = 3;
00368         dpprfs_("U", &c__0, &c_n1, a, af, b, &c__1, x, &c__1, r1, r2, w, iw, &
00369                 info);
00370         chkxer_("DPPRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00371                 infoc_1.ok);
00372         infoc_1.infot = 7;
00373         dpprfs_("U", &c__2, &c__1, a, af, b, &c__1, x, &c__2, r1, r2, w, iw, &
00374                 info);
00375         chkxer_("DPPRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00376                 infoc_1.ok);
00377         infoc_1.infot = 9;
00378         dpprfs_("U", &c__2, &c__1, a, af, b, &c__2, x, &c__1, r1, r2, w, iw, &
00379                 info);
00380         chkxer_("DPPRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00381                 infoc_1.ok);
00382 
00383 /*        DPPCON */
00384 
00385         s_copy(srnamc_1.srnamt, "DPPCON", (ftnlen)32, (ftnlen)6);
00386         infoc_1.infot = 1;
00387         dppcon_("/", &c__0, a, &anrm, &rcond, w, iw, &info);
00388         chkxer_("DPPCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00389                 infoc_1.ok);
00390         infoc_1.infot = 2;
00391         dppcon_("U", &c_n1, a, &anrm, &rcond, w, iw, &info);
00392         chkxer_("DPPCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00393                 infoc_1.ok);
00394 
00395 /*        DPPEQU */
00396 
00397         s_copy(srnamc_1.srnamt, "DPPEQU", (ftnlen)32, (ftnlen)6);
00398         infoc_1.infot = 1;
00399         dppequ_("/", &c__0, a, r1, &rcond, &anrm, &info);
00400         chkxer_("DPPEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00401                 infoc_1.ok);
00402         infoc_1.infot = 2;
00403         dppequ_("U", &c_n1, a, r1, &rcond, &anrm, &info);
00404         chkxer_("DPPEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00405                 infoc_1.ok);
00406 
00407     } else if (lsamen_(&c__2, c2, "PB")) {
00408 
00409 /*        Test error exits of the routines that use the Cholesky */
00410 /*        decomposition of a symmetric positive definite band matrix. */
00411 
00412 /*        DPBTRF */
00413 
00414         s_copy(srnamc_1.srnamt, "DPBTRF", (ftnlen)32, (ftnlen)6);
00415         infoc_1.infot = 1;
00416         dpbtrf_("/", &c__0, &c__0, a, &c__1, &info);
00417         chkxer_("DPBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00418                 infoc_1.ok);
00419         infoc_1.infot = 2;
00420         dpbtrf_("U", &c_n1, &c__0, a, &c__1, &info);
00421         chkxer_("DPBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00422                 infoc_1.ok);
00423         infoc_1.infot = 3;
00424         dpbtrf_("U", &c__1, &c_n1, a, &c__1, &info);
00425         chkxer_("DPBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00426                 infoc_1.ok);
00427         infoc_1.infot = 5;
00428         dpbtrf_("U", &c__2, &c__1, a, &c__1, &info);
00429         chkxer_("DPBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00430                 infoc_1.ok);
00431 
00432 /*        DPBTF2 */
00433 
00434         s_copy(srnamc_1.srnamt, "DPBTF2", (ftnlen)32, (ftnlen)6);
00435         infoc_1.infot = 1;
00436         dpbtf2_("/", &c__0, &c__0, a, &c__1, &info);
00437         chkxer_("DPBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00438                 infoc_1.ok);
00439         infoc_1.infot = 2;
00440         dpbtf2_("U", &c_n1, &c__0, a, &c__1, &info);
00441         chkxer_("DPBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00442                 infoc_1.ok);
00443         infoc_1.infot = 3;
00444         dpbtf2_("U", &c__1, &c_n1, a, &c__1, &info);
00445         chkxer_("DPBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00446                 infoc_1.ok);
00447         infoc_1.infot = 5;
00448         dpbtf2_("U", &c__2, &c__1, a, &c__1, &info);
00449         chkxer_("DPBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00450                 infoc_1.ok);
00451 
00452 /*        DPBTRS */
00453 
00454         s_copy(srnamc_1.srnamt, "DPBTRS", (ftnlen)32, (ftnlen)6);
00455         infoc_1.infot = 1;
00456         dpbtrs_("/", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, &info);
00457         chkxer_("DPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00458                 infoc_1.ok);
00459         infoc_1.infot = 2;
00460         dpbtrs_("U", &c_n1, &c__0, &c__0, a, &c__1, b, &c__1, &info);
00461         chkxer_("DPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00462                 infoc_1.ok);
00463         infoc_1.infot = 3;
00464         dpbtrs_("U", &c__1, &c_n1, &c__0, a, &c__1, b, &c__1, &info);
00465         chkxer_("DPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00466                 infoc_1.ok);
00467         infoc_1.infot = 4;
00468         dpbtrs_("U", &c__0, &c__0, &c_n1, a, &c__1, b, &c__1, &info);
00469         chkxer_("DPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00470                 infoc_1.ok);
00471         infoc_1.infot = 6;
00472         dpbtrs_("U", &c__2, &c__1, &c__1, a, &c__1, b, &c__1, &info);
00473         chkxer_("DPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00474                 infoc_1.ok);
00475         infoc_1.infot = 8;
00476         dpbtrs_("U", &c__2, &c__0, &c__1, a, &c__1, b, &c__1, &info);
00477         chkxer_("DPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00478                 infoc_1.ok);
00479 
00480 /*        DPBRFS */
00481 
00482         s_copy(srnamc_1.srnamt, "DPBRFS", (ftnlen)32, (ftnlen)6);
00483         infoc_1.infot = 1;
00484         dpbrfs_("/", &c__0, &c__0, &c__0, a, &c__1, af, &c__1, b, &c__1, x, &
00485                 c__1, r1, r2, w, iw, &info);
00486         chkxer_("DPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00487                 infoc_1.ok);
00488         infoc_1.infot = 2;
00489         dpbrfs_("U", &c_n1, &c__0, &c__0, a, &c__1, af, &c__1, b, &c__1, x, &
00490                 c__1, r1, r2, w, iw, &info);
00491         chkxer_("DPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00492                 infoc_1.ok);
00493         infoc_1.infot = 3;
00494         dpbrfs_("U", &c__1, &c_n1, &c__0, a, &c__1, af, &c__1, b, &c__1, x, &
00495                 c__1, r1, r2, w, iw, &info);
00496         chkxer_("DPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00497                 infoc_1.ok);
00498         infoc_1.infot = 4;
00499         dpbrfs_("U", &c__0, &c__0, &c_n1, a, &c__1, af, &c__1, b, &c__1, x, &
00500                 c__1, r1, r2, w, iw, &info);
00501         chkxer_("DPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00502                 infoc_1.ok);
00503         infoc_1.infot = 6;
00504         dpbrfs_("U", &c__2, &c__1, &c__1, a, &c__1, af, &c__2, b, &c__2, x, &
00505                 c__2, r1, r2, w, iw, &info);
00506         chkxer_("DPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00507                 infoc_1.ok);
00508         infoc_1.infot = 8;
00509         dpbrfs_("U", &c__2, &c__1, &c__1, a, &c__2, af, &c__1, b, &c__2, x, &
00510                 c__2, r1, r2, w, iw, &info);
00511         chkxer_("DPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00512                 infoc_1.ok);
00513         infoc_1.infot = 10;
00514         dpbrfs_("U", &c__2, &c__0, &c__1, a, &c__1, af, &c__1, b, &c__1, x, &
00515                 c__2, r1, r2, w, iw, &info);
00516         chkxer_("DPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00517                 infoc_1.ok);
00518         infoc_1.infot = 12;
00519         dpbrfs_("U", &c__2, &c__0, &c__1, a, &c__1, af, &c__1, b, &c__2, x, &
00520                 c__1, r1, r2, w, iw, &info);
00521         chkxer_("DPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00522                 infoc_1.ok);
00523 
00524 /*        DPBCON */
00525 
00526         s_copy(srnamc_1.srnamt, "DPBCON", (ftnlen)32, (ftnlen)6);
00527         infoc_1.infot = 1;
00528         dpbcon_("/", &c__0, &c__0, a, &c__1, &anrm, &rcond, w, iw, &info);
00529         chkxer_("DPBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00530                 infoc_1.ok);
00531         infoc_1.infot = 2;
00532         dpbcon_("U", &c_n1, &c__0, a, &c__1, &anrm, &rcond, w, iw, &info);
00533         chkxer_("DPBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00534                 infoc_1.ok);
00535         infoc_1.infot = 3;
00536         dpbcon_("U", &c__1, &c_n1, a, &c__1, &anrm, &rcond, w, iw, &info);
00537         chkxer_("DPBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00538                 infoc_1.ok);
00539         infoc_1.infot = 5;
00540         dpbcon_("U", &c__2, &c__1, a, &c__1, &anrm, &rcond, w, iw, &info);
00541         chkxer_("DPBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00542                 infoc_1.ok);
00543 
00544 /*        DPBEQU */
00545 
00546         s_copy(srnamc_1.srnamt, "DPBEQU", (ftnlen)32, (ftnlen)6);
00547         infoc_1.infot = 1;
00548         dpbequ_("/", &c__0, &c__0, a, &c__1, r1, &rcond, &anrm, &info);
00549         chkxer_("DPBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00550                 infoc_1.ok);
00551         infoc_1.infot = 2;
00552         dpbequ_("U", &c_n1, &c__0, a, &c__1, r1, &rcond, &anrm, &info);
00553         chkxer_("DPBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00554                 infoc_1.ok);
00555         infoc_1.infot = 3;
00556         dpbequ_("U", &c__1, &c_n1, a, &c__1, r1, &rcond, &anrm, &info);
00557         chkxer_("DPBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00558                 infoc_1.ok);
00559         infoc_1.infot = 5;
00560         dpbequ_("U", &c__2, &c__1, a, &c__1, r1, &rcond, &anrm, &info);
00561         chkxer_("DPBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00562                 infoc_1.ok);
00563     }
00564 
00565 /*     Print a summary line. */
00566 
00567     alaesm_(path, &infoc_1.ok, &infoc_1.nout);
00568 
00569     return 0;
00570 
00571 /*     End of DERRPO */
00572 
00573 } /* derrpo_ */