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


swiftnav
Author(s):
autogenerated on Sat Jun 8 2019 18:56:04