swiftnav: serrpo.c Source File

Go to the documentation of this file.
00001 /* serrpo.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 serrpo_(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     real a[16]  /* was [4][4] */, b[4];
00046     integer i__, j;
00047     real w[12], x[4];
00048     char c2[2];
00049     real r1[4], r2[4], af[16]   /* was [4][4] */;
00050     integer iw[4], info;
00051     real anrm, rcond;
00052     extern /* Subroutine */ int spbtf2_(char *, integer *, integer *, real *, 
00053             integer *, integer *), spotf2_(char *, integer *, real *, 
00054             integer *, integer *), alaesm_(char *, logical *, integer 
00055             *);
00056     extern logical lsamen_(integer *, char *, char *);
00057     extern /* Subroutine */ int chkxer_(char *, integer *, integer *, logical 
00058             *, logical *), spbcon_(char *, integer *, integer *, real 
00059             *, integer *, real *, real *, real *, integer *, integer *), spbequ_(char *, integer *, integer *, real *, integer *, 
00060             real *, real *, real *, integer *), spbrfs_(char *, 
00061             integer *, integer *, integer *, real *, integer *, real *, 
00062             integer *, real *, integer *, real *, integer *, real *, real *, 
00063             real *, integer *, integer *), spbtrf_(char *, integer *, 
00064             integer *, real *, integer *, integer *), spocon_(char *, 
00065             integer *, real *, integer *, real *, real *, real *, integer *, 
00066             integer *), sppcon_(char *, integer *, real *, real *, 
00067             real *, real *, integer *, integer *), spoequ_(integer *, 
00068             real *, integer *, real *, real *, real *, integer *), spbtrs_(
00069             char *, integer *, integer *, integer *, real *, integer *, real *
00070 , integer *, integer *), sporfs_(char *, integer *, 
00071             integer *, real *, integer *, real *, integer *, real *, integer *
00072 , real *, integer *, real *, real *, real *, integer *, integer *), spotrf_(char *, integer *, real *, integer *, integer *), spotri_(char *, integer *, real *, integer *, integer *), sppequ_(char *, integer *, real *, real *, real *, real 
00073             *, integer *), spprfs_(char *, integer *, integer *, real 
00074             *, real *, real *, integer *, real *, integer *, real *, real *, 
00075             real *, integer *, integer *), spptrf_(char *, integer *, 
00076             real *, integer *), spptri_(char *, integer *, real *, 
00077             integer *), spotrs_(char *, integer *, integer *, real *, 
00078             integer *, real *, integer *, integer *), spptrs_(char *, 
00079             integer *, integer *, real *, real *, integer *, integer *);
00080 
00081     /* Fortran I/O blocks */
00082     static cilist io___1 = { 0, 0, 0, 0, 0 };
00083 
00084 
00085 
00086 /*  -- LAPACK test routine (version 3.1) -- */
00087 /*     Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. */
00088 /*     November 2006 */
00089 
00090 /*     .. Scalar Arguments .. */
00091 /*     .. */
00092 
00093 /*  Purpose */
00094 /*  ======= */
00095 
00096 /*  SERRPO tests the error exits for the REAL routines */
00097 /*  for symmetric positive definite matrices. */
00098 
00099 /*  Arguments */
00100 /*  ========= */
00101 
00102 /*  PATH    (input) CHARACTER*3 */
00103 /*          The LAPACK path name for the routines to be tested. */
00104 
00105 /*  NUNIT   (input) INTEGER */
00106 /*          The unit number for output. */
00107 
00108 /*  ===================================================================== */
00109 
00110 /*     .. Parameters .. */
00111 /*     .. */
00112 /*     .. Local Scalars .. */
00113 /*     .. */
00114 /*     .. Local Arrays .. */
00115 /*     .. */
00116 /*     .. External Functions .. */
00117 /*     .. */
00118 /*     .. External Subroutines .. */
00119 /*     .. */
00120 /*     .. Scalars in Common .. */
00121 /*     .. */
00122 /*     .. Common blocks .. */
00123 /*     .. */
00124 /*     .. Intrinsic Functions .. */
00125 /*     .. */
00126 /*     .. Executable Statements .. */
00127 
00128     infoc_1.nout = *nunit;
00129     io___1.ciunit = infoc_1.nout;
00130     s_wsle(&io___1);
00131     e_wsle();
00132     s_copy(c2, path + 1, (ftnlen)2, (ftnlen)2);
00133 
00134 /*     Set the variables to innocuous values. */
00135 
00136     for (j = 1; j <= 4; ++j) {
00137         for (i__ = 1; i__ <= 4; ++i__) {
00138             a[i__ + (j << 2) - 5] = 1.f / (real) (i__ + j);
00139             af[i__ + (j << 2) - 5] = 1.f / (real) (i__ + j);
00140 /* L10: */
00141         }
00142         b[j - 1] = 0.f;
00143         r1[j - 1] = 0.f;
00144         r2[j - 1] = 0.f;
00145         w[j - 1] = 0.f;
00146         x[j - 1] = 0.f;
00147         iw[j - 1] = j;
00148 /* L20: */
00149     }
00150     infoc_1.ok = TRUE_;
00151 
00152     if (lsamen_(&c__2, c2, "PO")) {
00153 
00154 /*        Test error exits of the routines that use the Cholesky */
00155 /*        decomposition of a symmetric positive definite matrix. */
00156 
00157 /*        SPOTRF */
00158 
00159         s_copy(srnamc_1.srnamt, "SPOTRF", (ftnlen)32, (ftnlen)6);
00160         infoc_1.infot = 1;
00161         spotrf_("/", &c__0, a, &c__1, &info);
00162         chkxer_("SPOTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00163                 infoc_1.ok);
00164         infoc_1.infot = 2;
00165         spotrf_("U", &c_n1, a, &c__1, &info);
00166         chkxer_("SPOTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00167                 infoc_1.ok);
00168         infoc_1.infot = 4;
00169         spotrf_("U", &c__2, a, &c__1, &info);
00170         chkxer_("SPOTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00171                 infoc_1.ok);
00172 
00173 /*        SPOTF2 */
00174 
00175         s_copy(srnamc_1.srnamt, "SPOTF2", (ftnlen)32, (ftnlen)6);
00176         infoc_1.infot = 1;
00177         spotf2_("/", &c__0, a, &c__1, &info);
00178         chkxer_("SPOTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00179                 infoc_1.ok);
00180         infoc_1.infot = 2;
00181         spotf2_("U", &c_n1, a, &c__1, &info);
00182         chkxer_("SPOTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00183                 infoc_1.ok);
00184         infoc_1.infot = 4;
00185         spotf2_("U", &c__2, a, &c__1, &info);
00186         chkxer_("SPOTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00187                 infoc_1.ok);
00188 
00189 /*        SPOTRI */
00190 
00191         s_copy(srnamc_1.srnamt, "SPOTRI", (ftnlen)32, (ftnlen)6);
00192         infoc_1.infot = 1;
00193         spotri_("/", &c__0, a, &c__1, &info);
00194         chkxer_("SPOTRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00195                 infoc_1.ok);
00196         infoc_1.infot = 2;
00197         spotri_("U", &c_n1, a, &c__1, &info);
00198         chkxer_("SPOTRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00199                 infoc_1.ok);
00200         infoc_1.infot = 4;
00201         spotri_("U", &c__2, a, &c__1, &info);
00202         chkxer_("SPOTRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00203                 infoc_1.ok);
00204 
00205 /*        SPOTRS */
00206 
00207         s_copy(srnamc_1.srnamt, "SPOTRS", (ftnlen)32, (ftnlen)6);
00208         infoc_1.infot = 1;
00209         spotrs_("/", &c__0, &c__0, a, &c__1, b, &c__1, &info);
00210         chkxer_("SPOTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00211                 infoc_1.ok);
00212         infoc_1.infot = 2;
00213         spotrs_("U", &c_n1, &c__0, a, &c__1, b, &c__1, &info);
00214         chkxer_("SPOTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00215                 infoc_1.ok);
00216         infoc_1.infot = 3;
00217         spotrs_("U", &c__0, &c_n1, a, &c__1, b, &c__1, &info);
00218         chkxer_("SPOTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00219                 infoc_1.ok);
00220         infoc_1.infot = 5;
00221         spotrs_("U", &c__2, &c__1, a, &c__1, b, &c__2, &info);
00222         chkxer_("SPOTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00223                 infoc_1.ok);
00224         infoc_1.infot = 7;
00225         spotrs_("U", &c__2, &c__1, a, &c__2, b, &c__1, &info);
00226         chkxer_("SPOTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00227                 infoc_1.ok);
00228 
00229 /*        SPORFS */
00230 
00231         s_copy(srnamc_1.srnamt, "SPORFS", (ftnlen)32, (ftnlen)6);
00232         infoc_1.infot = 1;
00233         sporfs_("/", &c__0, &c__0, a, &c__1, af, &c__1, b, &c__1, x, &c__1, 
00234                 r1, r2, w, iw, &info);
00235         chkxer_("SPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00236                 infoc_1.ok);
00237         infoc_1.infot = 2;
00238         sporfs_("U", &c_n1, &c__0, a, &c__1, af, &c__1, b, &c__1, x, &c__1, 
00239                 r1, r2, w, iw, &info);
00240         chkxer_("SPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00241                 infoc_1.ok);
00242         infoc_1.infot = 3;
00243         sporfs_("U", &c__0, &c_n1, a, &c__1, af, &c__1, b, &c__1, x, &c__1, 
00244                 r1, r2, w, iw, &info);
00245         chkxer_("SPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00246                 infoc_1.ok);
00247         infoc_1.infot = 5;
00248         sporfs_("U", &c__2, &c__1, a, &c__1, af, &c__2, b, &c__2, x, &c__2, 
00249                 r1, r2, w, iw, &info);
00250         chkxer_("SPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00251                 infoc_1.ok);
00252         infoc_1.infot = 7;
00253         sporfs_("U", &c__2, &c__1, a, &c__2, af, &c__1, b, &c__2, x, &c__2, 
00254                 r1, r2, w, iw, &info);
00255         chkxer_("SPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00256                 infoc_1.ok);
00257         infoc_1.infot = 9;
00258         sporfs_("U", &c__2, &c__1, a, &c__2, af, &c__2, b, &c__1, x, &c__2, 
00259                 r1, r2, w, iw, &info);
00260         chkxer_("SPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00261                 infoc_1.ok);
00262         infoc_1.infot = 11;
00263         sporfs_("U", &c__2, &c__1, a, &c__2, af, &c__2, b, &c__2, x, &c__1, 
00264                 r1, r2, w, iw, &info);
00265         chkxer_("SPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00266                 infoc_1.ok);
00267 
00268 /*        SPOCON */
00269 
00270         s_copy(srnamc_1.srnamt, "SPOCON", (ftnlen)32, (ftnlen)6);
00271         infoc_1.infot = 1;
00272         spocon_("/", &c__0, a, &c__1, &anrm, &rcond, w, iw, &info);
00273         chkxer_("SPOCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00274                 infoc_1.ok);
00275         infoc_1.infot = 2;
00276         spocon_("U", &c_n1, a, &c__1, &anrm, &rcond, w, iw, &info);
00277         chkxer_("SPOCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00278                 infoc_1.ok);
00279         infoc_1.infot = 4;
00280         spocon_("U", &c__2, a, &c__1, &anrm, &rcond, w, iw, &info);
00281         chkxer_("SPOCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00282                 infoc_1.ok);
00283 
00284 /*        SPOEQU */
00285 
00286         s_copy(srnamc_1.srnamt, "SPOEQU", (ftnlen)32, (ftnlen)6);
00287         infoc_1.infot = 1;
00288         spoequ_(&c_n1, a, &c__1, r1, &rcond, &anrm, &info);
00289         chkxer_("SPOEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00290                 infoc_1.ok);
00291         infoc_1.infot = 3;
00292         spoequ_(&c__2, a, &c__1, r1, &rcond, &anrm, &info);
00293         chkxer_("SPOEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00294                 infoc_1.ok);
00295 
00296     } else if (lsamen_(&c__2, c2, "PP")) {
00297 
00298 /*        Test error exits of the routines that use the Cholesky */
00299 /*        decomposition of a symmetric positive definite packed matrix. */
00300 
00301 /*        SPPTRF */
00302 
00303         s_copy(srnamc_1.srnamt, "SPPTRF", (ftnlen)32, (ftnlen)6);
00304         infoc_1.infot = 1;
00305         spptrf_("/", &c__0, a, &info);
00306         chkxer_("SPPTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00307                 infoc_1.ok);
00308         infoc_1.infot = 2;
00309         spptrf_("U", &c_n1, a, &info);
00310         chkxer_("SPPTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00311                 infoc_1.ok);
00312 
00313 /*        SPPTRI */
00314 
00315         s_copy(srnamc_1.srnamt, "SPPTRI", (ftnlen)32, (ftnlen)6);
00316         infoc_1.infot = 1;
00317         spptri_("/", &c__0, a, &info);
00318         chkxer_("SPPTRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00319                 infoc_1.ok);
00320         infoc_1.infot = 2;
00321         spptri_("U", &c_n1, a, &info);
00322         chkxer_("SPPTRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00323                 infoc_1.ok);
00324 
00325 /*        SPPTRS */
00326 
00327         s_copy(srnamc_1.srnamt, "SPPTRS", (ftnlen)32, (ftnlen)6);
00328         infoc_1.infot = 1;
00329         spptrs_("/", &c__0, &c__0, a, b, &c__1, &info);
00330         chkxer_("SPPTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00331                 infoc_1.ok);
00332         infoc_1.infot = 2;
00333         spptrs_("U", &c_n1, &c__0, a, b, &c__1, &info);
00334         chkxer_("SPPTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00335                 infoc_1.ok);
00336         infoc_1.infot = 3;
00337         spptrs_("U", &c__0, &c_n1, a, b, &c__1, &info);
00338         chkxer_("SPPTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00339                 infoc_1.ok);
00340         infoc_1.infot = 6;
00341         spptrs_("U", &c__2, &c__1, a, b, &c__1, &info);
00342         chkxer_("SPPTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00343                 infoc_1.ok);
00344 
00345 /*        SPPRFS */
00346 
00347         s_copy(srnamc_1.srnamt, "SPPRFS", (ftnlen)32, (ftnlen)6);
00348         infoc_1.infot = 1;
00349         spprfs_("/", &c__0, &c__0, a, af, b, &c__1, x, &c__1, r1, r2, w, iw, &
00350                 info);
00351         chkxer_("SPPRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00352                 infoc_1.ok);
00353         infoc_1.infot = 2;
00354         spprfs_("U", &c_n1, &c__0, a, af, b, &c__1, x, &c__1, r1, r2, w, iw, &
00355                 info);
00356         chkxer_("SPPRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00357                 infoc_1.ok);
00358         infoc_1.infot = 3;
00359         spprfs_("U", &c__0, &c_n1, a, af, b, &c__1, x, &c__1, r1, r2, w, iw, &
00360                 info);
00361         chkxer_("SPPRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00362                 infoc_1.ok);
00363         infoc_1.infot = 7;
00364         spprfs_("U", &c__2, &c__1, a, af, b, &c__1, x, &c__2, r1, r2, w, iw, &
00365                 info);
00366         chkxer_("SPPRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00367                 infoc_1.ok);
00368         infoc_1.infot = 9;
00369         spprfs_("U", &c__2, &c__1, a, af, b, &c__2, x, &c__1, r1, r2, w, iw, &
00370                 info);
00371         chkxer_("SPPRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00372                 infoc_1.ok);
00373 
00374 /*        SPPCON */
00375 
00376         s_copy(srnamc_1.srnamt, "SPPCON", (ftnlen)32, (ftnlen)6);
00377         infoc_1.infot = 1;
00378         sppcon_("/", &c__0, a, &anrm, &rcond, w, iw, &info);
00379         chkxer_("SPPCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00380                 infoc_1.ok);
00381         infoc_1.infot = 2;
00382         sppcon_("U", &c_n1, a, &anrm, &rcond, w, iw, &info);
00383         chkxer_("SPPCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00384                 infoc_1.ok);
00385 
00386 /*        SPPEQU */
00387 
00388         s_copy(srnamc_1.srnamt, "SPPEQU", (ftnlen)32, (ftnlen)6);
00389         infoc_1.infot = 1;
00390         sppequ_("/", &c__0, a, r1, &rcond, &anrm, &info);
00391         chkxer_("SPPEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00392                 infoc_1.ok);
00393         infoc_1.infot = 2;
00394         sppequ_("U", &c_n1, a, r1, &rcond, &anrm, &info);
00395         chkxer_("SPPEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00396                 infoc_1.ok);
00397 
00398     } else if (lsamen_(&c__2, c2, "PB")) {
00399 
00400 /*        Test error exits of the routines that use the Cholesky */
00401 /*        decomposition of a symmetric positive definite band matrix. */
00402 
00403 /*        SPBTRF */
00404 
00405         s_copy(srnamc_1.srnamt, "SPBTRF", (ftnlen)32, (ftnlen)6);
00406         infoc_1.infot = 1;
00407         spbtrf_("/", &c__0, &c__0, a, &c__1, &info);
00408         chkxer_("SPBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00409                 infoc_1.ok);
00410         infoc_1.infot = 2;
00411         spbtrf_("U", &c_n1, &c__0, a, &c__1, &info);
00412         chkxer_("SPBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00413                 infoc_1.ok);
00414         infoc_1.infot = 3;
00415         spbtrf_("U", &c__1, &c_n1, a, &c__1, &info);
00416         chkxer_("SPBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00417                 infoc_1.ok);
00418         infoc_1.infot = 5;
00419         spbtrf_("U", &c__2, &c__1, a, &c__1, &info);
00420         chkxer_("SPBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00421                 infoc_1.ok);
00422 
00423 /*        SPBTF2 */
00424 
00425         s_copy(srnamc_1.srnamt, "SPBTF2", (ftnlen)32, (ftnlen)6);
00426         infoc_1.infot = 1;
00427         spbtf2_("/", &c__0, &c__0, a, &c__1, &info);
00428         chkxer_("SPBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00429                 infoc_1.ok);
00430         infoc_1.infot = 2;
00431         spbtf2_("U", &c_n1, &c__0, a, &c__1, &info);
00432         chkxer_("SPBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00433                 infoc_1.ok);
00434         infoc_1.infot = 3;
00435         spbtf2_("U", &c__1, &c_n1, a, &c__1, &info);
00436         chkxer_("SPBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00437                 infoc_1.ok);
00438         infoc_1.infot = 5;
00439         spbtf2_("U", &c__2, &c__1, a, &c__1, &info);
00440         chkxer_("SPBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00441                 infoc_1.ok);
00442 
00443 /*        SPBTRS */
00444 
00445         s_copy(srnamc_1.srnamt, "SPBTRS", (ftnlen)32, (ftnlen)6);
00446         infoc_1.infot = 1;
00447         spbtrs_("/", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, &info);
00448         chkxer_("SPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00449                 infoc_1.ok);
00450         infoc_1.infot = 2;
00451         spbtrs_("U", &c_n1, &c__0, &c__0, a, &c__1, b, &c__1, &info);
00452         chkxer_("SPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00453                 infoc_1.ok);
00454         infoc_1.infot = 3;
00455         spbtrs_("U", &c__1, &c_n1, &c__0, a, &c__1, b, &c__1, &info);
00456         chkxer_("SPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00457                 infoc_1.ok);
00458         infoc_1.infot = 4;
00459         spbtrs_("U", &c__0, &c__0, &c_n1, a, &c__1, b, &c__1, &info);
00460         chkxer_("SPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00461                 infoc_1.ok);
00462         infoc_1.infot = 6;
00463         spbtrs_("U", &c__2, &c__1, &c__1, a, &c__1, b, &c__1, &info);
00464         chkxer_("SPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00465                 infoc_1.ok);
00466         infoc_1.infot = 8;
00467         spbtrs_("U", &c__2, &c__0, &c__1, a, &c__1, b, &c__1, &info);
00468         chkxer_("SPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00469                 infoc_1.ok);
00470 
00471 /*        SPBRFS */
00472 
00473         s_copy(srnamc_1.srnamt, "SPBRFS", (ftnlen)32, (ftnlen)6);
00474         infoc_1.infot = 1;
00475         spbrfs_("/", &c__0, &c__0, &c__0, a, &c__1, af, &c__1, b, &c__1, x, &
00476                 c__1, r1, r2, w, iw, &info);
00477         chkxer_("SPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00478                 infoc_1.ok);
00479         infoc_1.infot = 2;
00480         spbrfs_("U", &c_n1, &c__0, &c__0, a, &c__1, af, &c__1, b, &c__1, x, &
00481                 c__1, r1, r2, w, iw, &info);
00482         chkxer_("SPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00483                 infoc_1.ok);
00484         infoc_1.infot = 3;
00485         spbrfs_("U", &c__1, &c_n1, &c__0, a, &c__1, af, &c__1, b, &c__1, x, &
00486                 c__1, r1, r2, w, iw, &info);
00487         chkxer_("SPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00488                 infoc_1.ok);
00489         infoc_1.infot = 4;
00490         spbrfs_("U", &c__0, &c__0, &c_n1, a, &c__1, af, &c__1, b, &c__1, x, &
00491                 c__1, r1, r2, w, iw, &info);
00492         chkxer_("SPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00493                 infoc_1.ok);
00494         infoc_1.infot = 6;
00495         spbrfs_("U", &c__2, &c__1, &c__1, a, &c__1, af, &c__2, b, &c__2, x, &
00496                 c__2, r1, r2, w, iw, &info);
00497         chkxer_("SPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00498                 infoc_1.ok);
00499         infoc_1.infot = 8;
00500         spbrfs_("U", &c__2, &c__1, &c__1, a, &c__2, af, &c__1, b, &c__2, x, &
00501                 c__2, r1, r2, w, iw, &info);
00502         chkxer_("SPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00503                 infoc_1.ok);
00504         infoc_1.infot = 10;
00505         spbrfs_("U", &c__2, &c__0, &c__1, a, &c__1, af, &c__1, b, &c__1, x, &
00506                 c__2, r1, r2, w, iw, &info);
00507         chkxer_("SPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00508                 infoc_1.ok);
00509         infoc_1.infot = 12;
00510         spbrfs_("U", &c__2, &c__0, &c__1, a, &c__1, af, &c__1, b, &c__2, x, &
00511                 c__1, r1, r2, w, iw, &info);
00512         chkxer_("SPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00513                 infoc_1.ok);
00514 
00515 /*        SPBCON */
00516 
00517         s_copy(srnamc_1.srnamt, "SPBCON", (ftnlen)32, (ftnlen)6);
00518         infoc_1.infot = 1;
00519         spbcon_("/", &c__0, &c__0, a, &c__1, &anrm, &rcond, w, iw, &info);
00520         chkxer_("SPBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00521                 infoc_1.ok);
00522         infoc_1.infot = 2;
00523         spbcon_("U", &c_n1, &c__0, a, &c__1, &anrm, &rcond, w, iw, &info);
00524         chkxer_("SPBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00525                 infoc_1.ok);
00526         infoc_1.infot = 3;
00527         spbcon_("U", &c__1, &c_n1, a, &c__1, &anrm, &rcond, w, iw, &info);
00528         chkxer_("SPBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00529                 infoc_1.ok);
00530         infoc_1.infot = 5;
00531         spbcon_("U", &c__2, &c__1, a, &c__1, &anrm, &rcond, w, iw, &info);
00532         chkxer_("SPBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00533                 infoc_1.ok);
00534 
00535 /*        SPBEQU */
00536 
00537         s_copy(srnamc_1.srnamt, "SPBEQU", (ftnlen)32, (ftnlen)6);
00538         infoc_1.infot = 1;
00539         spbequ_("/", &c__0, &c__0, a, &c__1, r1, &rcond, &anrm, &info);
00540         chkxer_("SPBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00541                 infoc_1.ok);
00542         infoc_1.infot = 2;
00543         spbequ_("U", &c_n1, &c__0, a, &c__1, r1, &rcond, &anrm, &info);
00544         chkxer_("SPBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00545                 infoc_1.ok);
00546         infoc_1.infot = 3;
00547         spbequ_("U", &c__1, &c_n1, a, &c__1, r1, &rcond, &anrm, &info);
00548         chkxer_("SPBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00549                 infoc_1.ok);
00550         infoc_1.infot = 5;
00551         spbequ_("U", &c__2, &c__1, a, &c__1, r1, &rcond, &anrm, &info);
00552         chkxer_("SPBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00553                 infoc_1.ok);
00554     }
00555 
00556 /*     Print a summary line. */
00557 
00558     alaesm_(path, &infoc_1.ok, &infoc_1.nout);
00559 
00560     return 0;
00561 
00562 /*     End of SERRPO */
00563 
00564 } /* serrpo_ */