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


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