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


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