swiftnav: derrpox.c Source File

Go to the documentation of this file.
00001 /* derrpox.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 s[4], w[12], x[4];
00048     char c2[2];
00049     doublereal r1[4], r2[4], af[16]     /* was [4][4] */;
00050     char eq[1];
00051     integer iw[4];
00052     doublereal err_bnds_c__[12] /* was [4][3] */;
00053     integer n_err_bnds__;
00054     doublereal err_bnds_n__[12] /* was [4][3] */, berr;
00055     integer info;
00056     doublereal anrm, rcond;
00057     extern /* Subroutine */ int dpbtf2_(char *, integer *, integer *, 
00058             doublereal *, integer *, integer *), dpotf2_(char *, 
00059             integer *, doublereal *, integer *, integer *), alaesm_(
00060             char *, logical *, integer *), dpbcon_(char *, integer *, 
00061             integer *, doublereal *, integer *, doublereal *, doublereal *, 
00062             doublereal *, integer *, integer *);
00063     extern logical lsamen_(integer *, char *, char *);
00064     extern /* Subroutine */ int dpbequ_(char *, integer *, integer *, 
00065             doublereal *, integer *, doublereal *, doublereal *, doublereal *, 
00066              integer *), dpbrfs_(char *, integer *, integer *, 
00067             integer *, doublereal *, integer *, doublereal *, integer *, 
00068             doublereal *, integer *, doublereal *, integer *, doublereal *, 
00069             doublereal *, doublereal *, integer *, integer *), 
00070             dpbtrf_(char *, integer *, integer *, doublereal *, integer *, 
00071             integer *);
00072     doublereal params;
00073     extern /* Subroutine */ int chkxer_(char *, integer *, integer *, logical 
00074             *, logical *), dpocon_(char *, integer *, doublereal *, 
00075             integer *, doublereal *, doublereal *, doublereal *, integer *, 
00076             integer *), dppcon_(char *, integer *, doublereal *, 
00077             doublereal *, doublereal *, doublereal *, integer *, integer *), dpoequ_(integer *, doublereal *, integer *, doublereal *, 
00078              doublereal *, doublereal *, integer *), dpbtrs_(char *, integer *
00079 , integer *, integer *, doublereal *, integer *, doublereal *, 
00080             integer *, integer *), dporfs_(char *, integer *, integer 
00081             *, doublereal *, integer *, doublereal *, integer *, doublereal *, 
00082              integer *, doublereal *, integer *, doublereal *, doublereal *, 
00083             doublereal *, integer *, integer *), dpotrf_(char *, 
00084             integer *, doublereal *, integer *, integer *), dpotri_(
00085             char *, integer *, doublereal *, integer *, integer *), 
00086             dppequ_(char *, integer *, doublereal *, doublereal *, doublereal 
00087             *, doublereal *, integer *), dpprfs_(char *, integer *, 
00088             integer *, doublereal *, doublereal *, doublereal *, integer *, 
00089             doublereal *, integer *, doublereal *, doublereal *, doublereal *, 
00090              integer *, integer *), dpptrf_(char *, integer *, 
00091             doublereal *, integer *), dpptri_(char *, integer *, 
00092             doublereal *, integer *), dpotrs_(char *, integer *, 
00093             integer *, doublereal *, integer *, doublereal *, integer *, 
00094             integer *), dpptrs_(char *, integer *, integer *, 
00095             doublereal *, doublereal *, integer *, integer *), 
00096             dpoequb_(integer *, doublereal *, integer *, doublereal *, 
00097             doublereal *, doublereal *, integer *);
00098     integer nparams;
00099     extern /* Subroutine */ int dporfsx_(char *, char *, integer *, integer *, 
00100              doublereal *, integer *, doublereal *, integer *, doublereal *, 
00101             doublereal *, integer *, doublereal *, integer *, doublereal *, 
00102             doublereal *, integer *, doublereal *, doublereal *, integer *, 
00103             doublereal *, doublereal *, integer *, integer *);
00104 
00105     /* Fortran I/O blocks */
00106     static cilist io___1 = { 0, 0, 0, 0, 0 };
00107 
00108 
00109 
00110 /*  -- LAPACK test routine (version 3.1) -- */
00111 /*     Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. */
00112 /*     November 2006 */
00113 
00114 /*     .. Scalar Arguments .. */
00115 /*     .. */
00116 
00117 /*  Purpose */
00118 /*  ======= */
00119 
00120 /*  DERRPO tests the error exits for the DOUBLE PRECISION routines */
00121 /*  for symmetric positive definite matrices. */
00122 
00123 /*  Arguments */
00124 /*  ========= */
00125 
00126 /*  PATH    (input) CHARACTER*3 */
00127 /*          The LAPACK path name for the routines to be tested. */
00128 
00129 /*  NUNIT   (input) INTEGER */
00130 /*          The unit number for output. */
00131 
00132 /*  ===================================================================== */
00133 
00134 /*     .. Parameters .. */
00135 /*     .. */
00136 /*     .. Local Scalars .. */
00137 /*     .. */
00138 /*     .. Local Arrays .. */
00139 /*     .. */
00140 /*     .. External Functions .. */
00141 /*     .. */
00142 /*     .. External Subroutines .. */
00143 /*     .. */
00144 /*     .. Scalars in Common .. */
00145 /*     .. */
00146 /*     .. Common blocks .. */
00147 /*     .. */
00148 /*     .. Intrinsic Functions .. */
00149 /*     .. */
00150 /*     .. Executable Statements .. */
00151 
00152     infoc_1.nout = *nunit;
00153     io___1.ciunit = infoc_1.nout;
00154     s_wsle(&io___1);
00155     e_wsle();
00156     s_copy(c2, path + 1, (ftnlen)2, (ftnlen)2);
00157 
00158 /*     Set the variables to innocuous values. */
00159 
00160     for (j = 1; j <= 4; ++j) {
00161         for (i__ = 1; i__ <= 4; ++i__) {
00162             a[i__ + (j << 2) - 5] = 1. / (doublereal) (i__ + j);
00163             af[i__ + (j << 2) - 5] = 1. / (doublereal) (i__ + j);
00164 /* L10: */
00165         }
00166         b[j - 1] = 0.;
00167         r1[j - 1] = 0.;
00168         r2[j - 1] = 0.;
00169         w[j - 1] = 0.;
00170         x[j - 1] = 0.;
00171         s[j - 1] = 0.;
00172         iw[j - 1] = j;
00173 /* L20: */
00174     }
00175     infoc_1.ok = TRUE_;
00176 
00177     if (lsamen_(&c__2, c2, "PO")) {
00178 
00179 /*        Test error exits of the routines that use the Cholesky */
00180 /*        decomposition of a symmetric positive definite matrix. */
00181 
00182 /*        DPOTRF */
00183 
00184         s_copy(srnamc_1.srnamt, "DPOTRF", (ftnlen)32, (ftnlen)6);
00185         infoc_1.infot = 1;
00186         dpotrf_("/", &c__0, a, &c__1, &info);
00187         chkxer_("DPOTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00188                 infoc_1.ok);
00189         infoc_1.infot = 2;
00190         dpotrf_("U", &c_n1, a, &c__1, &info);
00191         chkxer_("DPOTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00192                 infoc_1.ok);
00193         infoc_1.infot = 4;
00194         dpotrf_("U", &c__2, a, &c__1, &info);
00195         chkxer_("DPOTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00196                 infoc_1.ok);
00197 
00198 /*        DPOTF2 */
00199 
00200         s_copy(srnamc_1.srnamt, "DPOTF2", (ftnlen)32, (ftnlen)6);
00201         infoc_1.infot = 1;
00202         dpotf2_("/", &c__0, a, &c__1, &info);
00203         chkxer_("DPOTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00204                 infoc_1.ok);
00205         infoc_1.infot = 2;
00206         dpotf2_("U", &c_n1, a, &c__1, &info);
00207         chkxer_("DPOTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00208                 infoc_1.ok);
00209         infoc_1.infot = 4;
00210         dpotf2_("U", &c__2, a, &c__1, &info);
00211         chkxer_("DPOTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00212                 infoc_1.ok);
00213 
00214 /*        DPOTRI */
00215 
00216         s_copy(srnamc_1.srnamt, "DPOTRI", (ftnlen)32, (ftnlen)6);
00217         infoc_1.infot = 1;
00218         dpotri_("/", &c__0, a, &c__1, &info);
00219         chkxer_("DPOTRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00220                 infoc_1.ok);
00221         infoc_1.infot = 2;
00222         dpotri_("U", &c_n1, a, &c__1, &info);
00223         chkxer_("DPOTRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00224                 infoc_1.ok);
00225         infoc_1.infot = 4;
00226         dpotri_("U", &c__2, a, &c__1, &info);
00227         chkxer_("DPOTRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00228                 infoc_1.ok);
00229 
00230 /*        DPOTRS */
00231 
00232         s_copy(srnamc_1.srnamt, "DPOTRS", (ftnlen)32, (ftnlen)6);
00233         infoc_1.infot = 1;
00234         dpotrs_("/", &c__0, &c__0, a, &c__1, b, &c__1, &info);
00235         chkxer_("DPOTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00236                 infoc_1.ok);
00237         infoc_1.infot = 2;
00238         dpotrs_("U", &c_n1, &c__0, a, &c__1, b, &c__1, &info);
00239         chkxer_("DPOTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00240                 infoc_1.ok);
00241         infoc_1.infot = 3;
00242         dpotrs_("U", &c__0, &c_n1, a, &c__1, b, &c__1, &info);
00243         chkxer_("DPOTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00244                 infoc_1.ok);
00245         infoc_1.infot = 5;
00246         dpotrs_("U", &c__2, &c__1, a, &c__1, b, &c__2, &info);
00247         chkxer_("DPOTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00248                 infoc_1.ok);
00249         infoc_1.infot = 7;
00250         dpotrs_("U", &c__2, &c__1, a, &c__2, b, &c__1, &info);
00251         chkxer_("DPOTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00252                 infoc_1.ok);
00253 
00254 /*        DPORFS */
00255 
00256         s_copy(srnamc_1.srnamt, "DPORFS", (ftnlen)32, (ftnlen)6);
00257         infoc_1.infot = 1;
00258         dporfs_("/", &c__0, &c__0, a, &c__1, af, &c__1, b, &c__1, x, &c__1, 
00259                 r1, r2, w, iw, &info);
00260         chkxer_("DPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00261                 infoc_1.ok);
00262         infoc_1.infot = 2;
00263         dporfs_("U", &c_n1, &c__0, a, &c__1, af, &c__1, b, &c__1, x, &c__1, 
00264                 r1, r2, w, iw, &info);
00265         chkxer_("DPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00266                 infoc_1.ok);
00267         infoc_1.infot = 3;
00268         dporfs_("U", &c__0, &c_n1, a, &c__1, af, &c__1, b, &c__1, x, &c__1, 
00269                 r1, r2, w, iw, &info);
00270         chkxer_("DPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00271                 infoc_1.ok);
00272         infoc_1.infot = 5;
00273         dporfs_("U", &c__2, &c__1, a, &c__1, af, &c__2, b, &c__2, x, &c__2, 
00274                 r1, r2, w, iw, &info);
00275         chkxer_("DPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00276                 infoc_1.ok);
00277         infoc_1.infot = 7;
00278         dporfs_("U", &c__2, &c__1, a, &c__2, af, &c__1, b, &c__2, x, &c__2, 
00279                 r1, r2, w, iw, &info);
00280         chkxer_("DPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00281                 infoc_1.ok);
00282         infoc_1.infot = 9;
00283         dporfs_("U", &c__2, &c__1, a, &c__2, af, &c__2, b, &c__1, x, &c__2, 
00284                 r1, r2, w, iw, &info);
00285         chkxer_("DPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00286                 infoc_1.ok);
00287         infoc_1.infot = 11;
00288         dporfs_("U", &c__2, &c__1, a, &c__2, af, &c__2, b, &c__2, x, &c__1, 
00289                 r1, r2, w, iw, &info);
00290         chkxer_("DPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00291                 infoc_1.ok);
00292 
00293 /*        DPORFSX */
00294 
00295         n_err_bnds__ = 3;
00296         nparams = 0;
00297         s_copy(srnamc_1.srnamt, "DPORFSX", (ftnlen)32, (ftnlen)7);
00298         infoc_1.infot = 1;
00299         dporfsx_("/", eq, &c__0, &c__0, a, &c__1, af, &c__1, s, b, &c__1, x, &
00300                 c__1, &rcond, &berr, &n_err_bnds__, err_bnds_n__, 
00301                 err_bnds_c__, &nparams, &params, w, iw, &info);
00302         chkxer_("DPORFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00303                 infoc_1.ok);
00304         infoc_1.infot = 2;
00305         dporfsx_("U", eq, &c_n1, &c__0, a, &c__1, af, &c__1, s, b, &c__1, x, &
00306                 c__1, &rcond, &berr, &n_err_bnds__, err_bnds_n__, 
00307                 err_bnds_c__, &nparams, &params, w, iw, &info);
00308         chkxer_("DPORFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00309                 infoc_1.ok);
00310         *(unsigned char *)eq = 'N';
00311         infoc_1.infot = 3;
00312         dporfsx_("U", eq, &c_n1, &c__0, a, &c__1, af, &c__1, s, b, &c__1, x, &
00313                 c__1, &rcond, &berr, &n_err_bnds__, err_bnds_n__, 
00314                 err_bnds_c__, &nparams, &params, w, iw, &info);
00315         chkxer_("DPORFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00316                 infoc_1.ok);
00317         infoc_1.infot = 4;
00318         dporfsx_("U", eq, &c__0, &c_n1, a, &c__1, af, &c__1, s, b, &c__1, x, &
00319                 c__1, &rcond, &berr, &n_err_bnds__, err_bnds_n__, 
00320                 err_bnds_c__, &nparams, &params, w, iw, &info);
00321         chkxer_("DPORFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00322                 infoc_1.ok);
00323         infoc_1.infot = 6;
00324         dporfsx_("U", eq, &c__2, &c__1, a, &c__1, af, &c__2, s, b, &c__2, x, &
00325                 c__2, &rcond, &berr, &n_err_bnds__, err_bnds_n__, 
00326                 err_bnds_c__, &nparams, &params, w, iw, &info);
00327         chkxer_("DPORFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00328                 infoc_1.ok);
00329         infoc_1.infot = 8;
00330         dporfsx_("U", eq, &c__2, &c__1, a, &c__2, af, &c__1, s, b, &c__2, x, &
00331                 c__2, &rcond, &berr, &n_err_bnds__, err_bnds_n__, 
00332                 err_bnds_c__, &nparams, &params, w, iw, &info);
00333         chkxer_("DPORFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00334                 infoc_1.ok);
00335         infoc_1.infot = 11;
00336         dporfsx_("U", eq, &c__2, &c__1, a, &c__2, af, &c__2, s, b, &c__1, x, &
00337                 c__2, &rcond, &berr, &n_err_bnds__, err_bnds_n__, 
00338                 err_bnds_c__, &nparams, &params, w, iw, &info);
00339         chkxer_("DPORFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00340                 infoc_1.ok);
00341         infoc_1.infot = 13;
00342         dporfsx_("U", eq, &c__2, &c__1, a, &c__2, af, &c__2, s, b, &c__2, x, &
00343                 c__1, &rcond, &berr, &n_err_bnds__, err_bnds_n__, 
00344                 err_bnds_c__, &nparams, &params, w, iw, &info);
00345         chkxer_("DPORFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00346                 infoc_1.ok);
00347 
00348 /*        DPOCON */
00349 
00350         s_copy(srnamc_1.srnamt, "DPOCON", (ftnlen)32, (ftnlen)6);
00351         infoc_1.infot = 1;
00352         dpocon_("/", &c__0, a, &c__1, &anrm, &rcond, w, iw, &info);
00353         chkxer_("DPOCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00354                 infoc_1.ok);
00355         infoc_1.infot = 2;
00356         dpocon_("U", &c_n1, a, &c__1, &anrm, &rcond, w, iw, &info);
00357         chkxer_("DPOCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00358                 infoc_1.ok);
00359         infoc_1.infot = 4;
00360         dpocon_("U", &c__2, a, &c__1, &anrm, &rcond, w, iw, &info);
00361         chkxer_("DPOCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00362                 infoc_1.ok);
00363 
00364 /*        DPOEQU */
00365 
00366         s_copy(srnamc_1.srnamt, "DPOEQU", (ftnlen)32, (ftnlen)6);
00367         infoc_1.infot = 1;
00368         dpoequ_(&c_n1, a, &c__1, r1, &rcond, &anrm, &info);
00369         chkxer_("DPOEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00370                 infoc_1.ok);
00371         infoc_1.infot = 3;
00372         dpoequ_(&c__2, a, &c__1, r1, &rcond, &anrm, &info);
00373         chkxer_("DPOEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00374                 infoc_1.ok);
00375 
00376 /*        DPOEQUB */
00377 
00378         s_copy(srnamc_1.srnamt, "DPOEQUB", (ftnlen)32, (ftnlen)7);
00379         infoc_1.infot = 1;
00380         dpoequb_(&c_n1, a, &c__1, r1, &rcond, &anrm, &info);
00381         chkxer_("DPOEQUB", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00382                 infoc_1.ok);
00383         infoc_1.infot = 3;
00384         dpoequb_(&c__2, a, &c__1, r1, &rcond, &anrm, &info);
00385         chkxer_("DPOEQUB", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00386                 infoc_1.ok);
00387 
00388     } else if (lsamen_(&c__2, c2, "PP")) {
00389 
00390 /*        Test error exits of the routines that use the Cholesky */
00391 /*        decomposition of a symmetric positive definite packed matrix. */
00392 
00393 /*        DPPTRF */
00394 
00395         s_copy(srnamc_1.srnamt, "DPPTRF", (ftnlen)32, (ftnlen)6);
00396         infoc_1.infot = 1;
00397         dpptrf_("/", &c__0, a, &info);
00398         chkxer_("DPPTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00399                 infoc_1.ok);
00400         infoc_1.infot = 2;
00401         dpptrf_("U", &c_n1, a, &info);
00402         chkxer_("DPPTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00403                 infoc_1.ok);
00404 
00405 /*        DPPTRI */
00406 
00407         s_copy(srnamc_1.srnamt, "DPPTRI", (ftnlen)32, (ftnlen)6);
00408         infoc_1.infot = 1;
00409         dpptri_("/", &c__0, a, &info);
00410         chkxer_("DPPTRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00411                 infoc_1.ok);
00412         infoc_1.infot = 2;
00413         dpptri_("U", &c_n1, a, &info);
00414         chkxer_("DPPTRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00415                 infoc_1.ok);
00416 
00417 /*        DPPTRS */
00418 
00419         s_copy(srnamc_1.srnamt, "DPPTRS", (ftnlen)32, (ftnlen)6);
00420         infoc_1.infot = 1;
00421         dpptrs_("/", &c__0, &c__0, a, b, &c__1, &info);
00422         chkxer_("DPPTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00423                 infoc_1.ok);
00424         infoc_1.infot = 2;
00425         dpptrs_("U", &c_n1, &c__0, a, b, &c__1, &info);
00426         chkxer_("DPPTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00427                 infoc_1.ok);
00428         infoc_1.infot = 3;
00429         dpptrs_("U", &c__0, &c_n1, a, b, &c__1, &info);
00430         chkxer_("DPPTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00431                 infoc_1.ok);
00432         infoc_1.infot = 6;
00433         dpptrs_("U", &c__2, &c__1, a, b, &c__1, &info);
00434         chkxer_("DPPTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00435                 infoc_1.ok);
00436 
00437 /*        DPPRFS */
00438 
00439         s_copy(srnamc_1.srnamt, "DPPRFS", (ftnlen)32, (ftnlen)6);
00440         infoc_1.infot = 1;
00441         dpprfs_("/", &c__0, &c__0, a, af, b, &c__1, x, &c__1, r1, r2, w, iw, &
00442                 info);
00443         chkxer_("DPPRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00444                 infoc_1.ok);
00445         infoc_1.infot = 2;
00446         dpprfs_("U", &c_n1, &c__0, a, af, b, &c__1, x, &c__1, r1, r2, w, iw, &
00447                 info);
00448         chkxer_("DPPRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00449                 infoc_1.ok);
00450         infoc_1.infot = 3;
00451         dpprfs_("U", &c__0, &c_n1, a, af, b, &c__1, x, &c__1, r1, r2, w, iw, &
00452                 info);
00453         chkxer_("DPPRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00454                 infoc_1.ok);
00455         infoc_1.infot = 7;
00456         dpprfs_("U", &c__2, &c__1, a, af, b, &c__1, x, &c__2, r1, r2, w, iw, &
00457                 info);
00458         chkxer_("DPPRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00459                 infoc_1.ok);
00460         infoc_1.infot = 9;
00461         dpprfs_("U", &c__2, &c__1, a, af, b, &c__2, x, &c__1, r1, r2, w, iw, &
00462                 info);
00463         chkxer_("DPPRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00464                 infoc_1.ok);
00465 
00466 /*        DPPCON */
00467 
00468         s_copy(srnamc_1.srnamt, "DPPCON", (ftnlen)32, (ftnlen)6);
00469         infoc_1.infot = 1;
00470         dppcon_("/", &c__0, a, &anrm, &rcond, w, iw, &info);
00471         chkxer_("DPPCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00472                 infoc_1.ok);
00473         infoc_1.infot = 2;
00474         dppcon_("U", &c_n1, a, &anrm, &rcond, w, iw, &info);
00475         chkxer_("DPPCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00476                 infoc_1.ok);
00477 
00478 /*        DPPEQU */
00479 
00480         s_copy(srnamc_1.srnamt, "DPPEQU", (ftnlen)32, (ftnlen)6);
00481         infoc_1.infot = 1;
00482         dppequ_("/", &c__0, a, r1, &rcond, &anrm, &info);
00483         chkxer_("DPPEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00484                 infoc_1.ok);
00485         infoc_1.infot = 2;
00486         dppequ_("U", &c_n1, a, r1, &rcond, &anrm, &info);
00487         chkxer_("DPPEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00488                 infoc_1.ok);
00489 
00490     } else if (lsamen_(&c__2, c2, "PB")) {
00491 
00492 /*        Test error exits of the routines that use the Cholesky */
00493 /*        decomposition of a symmetric positive definite band matrix. */
00494 
00495 /*        DPBTRF */
00496 
00497         s_copy(srnamc_1.srnamt, "DPBTRF", (ftnlen)32, (ftnlen)6);
00498         infoc_1.infot = 1;
00499         dpbtrf_("/", &c__0, &c__0, a, &c__1, &info);
00500         chkxer_("DPBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00501                 infoc_1.ok);
00502         infoc_1.infot = 2;
00503         dpbtrf_("U", &c_n1, &c__0, a, &c__1, &info);
00504         chkxer_("DPBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00505                 infoc_1.ok);
00506         infoc_1.infot = 3;
00507         dpbtrf_("U", &c__1, &c_n1, a, &c__1, &info);
00508         chkxer_("DPBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00509                 infoc_1.ok);
00510         infoc_1.infot = 5;
00511         dpbtrf_("U", &c__2, &c__1, a, &c__1, &info);
00512         chkxer_("DPBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00513                 infoc_1.ok);
00514 
00515 /*        DPBTF2 */
00516 
00517         s_copy(srnamc_1.srnamt, "DPBTF2", (ftnlen)32, (ftnlen)6);
00518         infoc_1.infot = 1;
00519         dpbtf2_("/", &c__0, &c__0, a, &c__1, &info);
00520         chkxer_("DPBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00521                 infoc_1.ok);
00522         infoc_1.infot = 2;
00523         dpbtf2_("U", &c_n1, &c__0, a, &c__1, &info);
00524         chkxer_("DPBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00525                 infoc_1.ok);
00526         infoc_1.infot = 3;
00527         dpbtf2_("U", &c__1, &c_n1, a, &c__1, &info);
00528         chkxer_("DPBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00529                 infoc_1.ok);
00530         infoc_1.infot = 5;
00531         dpbtf2_("U", &c__2, &c__1, a, &c__1, &info);
00532         chkxer_("DPBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00533                 infoc_1.ok);
00534 
00535 /*        DPBTRS */
00536 
00537         s_copy(srnamc_1.srnamt, "DPBTRS", (ftnlen)32, (ftnlen)6);
00538         infoc_1.infot = 1;
00539         dpbtrs_("/", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, &info);
00540         chkxer_("DPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00541                 infoc_1.ok);
00542         infoc_1.infot = 2;
00543         dpbtrs_("U", &c_n1, &c__0, &c__0, a, &c__1, b, &c__1, &info);
00544         chkxer_("DPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00545                 infoc_1.ok);
00546         infoc_1.infot = 3;
00547         dpbtrs_("U", &c__1, &c_n1, &c__0, a, &c__1, b, &c__1, &info);
00548         chkxer_("DPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00549                 infoc_1.ok);
00550         infoc_1.infot = 4;
00551         dpbtrs_("U", &c__0, &c__0, &c_n1, a, &c__1, b, &c__1, &info);
00552         chkxer_("DPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00553                 infoc_1.ok);
00554         infoc_1.infot = 6;
00555         dpbtrs_("U", &c__2, &c__1, &c__1, a, &c__1, b, &c__1, &info);
00556         chkxer_("DPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00557                 infoc_1.ok);
00558         infoc_1.infot = 8;
00559         dpbtrs_("U", &c__2, &c__0, &c__1, a, &c__1, b, &c__1, &info);
00560         chkxer_("DPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00561                 infoc_1.ok);
00562 
00563 /*        DPBRFS */
00564 
00565         s_copy(srnamc_1.srnamt, "DPBRFS", (ftnlen)32, (ftnlen)6);
00566         infoc_1.infot = 1;
00567         dpbrfs_("/", &c__0, &c__0, &c__0, a, &c__1, af, &c__1, b, &c__1, x, &
00568                 c__1, r1, r2, w, iw, &info);
00569         chkxer_("DPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00570                 infoc_1.ok);
00571         infoc_1.infot = 2;
00572         dpbrfs_("U", &c_n1, &c__0, &c__0, a, &c__1, af, &c__1, b, &c__1, x, &
00573                 c__1, r1, r2, w, iw, &info);
00574         chkxer_("DPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00575                 infoc_1.ok);
00576         infoc_1.infot = 3;
00577         dpbrfs_("U", &c__1, &c_n1, &c__0, a, &c__1, af, &c__1, b, &c__1, x, &
00578                 c__1, r1, r2, w, iw, &info);
00579         chkxer_("DPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00580                 infoc_1.ok);
00581         infoc_1.infot = 4;
00582         dpbrfs_("U", &c__0, &c__0, &c_n1, a, &c__1, af, &c__1, b, &c__1, x, &
00583                 c__1, r1, r2, w, iw, &info);
00584         chkxer_("DPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00585                 infoc_1.ok);
00586         infoc_1.infot = 6;
00587         dpbrfs_("U", &c__2, &c__1, &c__1, a, &c__1, af, &c__2, b, &c__2, x, &
00588                 c__2, r1, r2, w, iw, &info);
00589         chkxer_("DPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00590                 infoc_1.ok);
00591         infoc_1.infot = 8;
00592         dpbrfs_("U", &c__2, &c__1, &c__1, a, &c__2, af, &c__1, b, &c__2, x, &
00593                 c__2, r1, r2, w, iw, &info);
00594         chkxer_("DPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00595                 infoc_1.ok);
00596         infoc_1.infot = 10;
00597         dpbrfs_("U", &c__2, &c__0, &c__1, a, &c__1, af, &c__1, b, &c__1, x, &
00598                 c__2, r1, r2, w, iw, &info);
00599         chkxer_("DPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00600                 infoc_1.ok);
00601         infoc_1.infot = 12;
00602         dpbrfs_("U", &c__2, &c__0, &c__1, a, &c__1, af, &c__1, b, &c__2, x, &
00603                 c__1, r1, r2, w, iw, &info);
00604         chkxer_("DPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00605                 infoc_1.ok);
00606 
00607 /*        DPBCON */
00608 
00609         s_copy(srnamc_1.srnamt, "DPBCON", (ftnlen)32, (ftnlen)6);
00610         infoc_1.infot = 1;
00611         dpbcon_("/", &c__0, &c__0, a, &c__1, &anrm, &rcond, w, iw, &info);
00612         chkxer_("DPBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00613                 infoc_1.ok);
00614         infoc_1.infot = 2;
00615         dpbcon_("U", &c_n1, &c__0, a, &c__1, &anrm, &rcond, w, iw, &info);
00616         chkxer_("DPBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00617                 infoc_1.ok);
00618         infoc_1.infot = 3;
00619         dpbcon_("U", &c__1, &c_n1, a, &c__1, &anrm, &rcond, w, iw, &info);
00620         chkxer_("DPBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00621                 infoc_1.ok);
00622         infoc_1.infot = 5;
00623         dpbcon_("U", &c__2, &c__1, a, &c__1, &anrm, &rcond, w, iw, &info);
00624         chkxer_("DPBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00625                 infoc_1.ok);
00626 
00627 /*        DPBEQU */
00628 
00629         s_copy(srnamc_1.srnamt, "DPBEQU", (ftnlen)32, (ftnlen)6);
00630         infoc_1.infot = 1;
00631         dpbequ_("/", &c__0, &c__0, a, &c__1, r1, &rcond, &anrm, &info);
00632         chkxer_("DPBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00633                 infoc_1.ok);
00634         infoc_1.infot = 2;
00635         dpbequ_("U", &c_n1, &c__0, a, &c__1, r1, &rcond, &anrm, &info);
00636         chkxer_("DPBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00637                 infoc_1.ok);
00638         infoc_1.infot = 3;
00639         dpbequ_("U", &c__1, &c_n1, a, &c__1, r1, &rcond, &anrm, &info);
00640         chkxer_("DPBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00641                 infoc_1.ok);
00642         infoc_1.infot = 5;
00643         dpbequ_("U", &c__2, &c__1, a, &c__1, r1, &rcond, &anrm, &info);
00644         chkxer_("DPBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00645                 infoc_1.ok);
00646     }
00647 
00648 /*     Print a summary line. */
00649 
00650     alaesm_(path, &infoc_1.ok, &infoc_1.nout);
00651 
00652     return 0;
00653 
00654 /*     End of DERRPO */
00655 
00656 } /* derrpo_ */