zerrgex.c
Go to the documentation of this file.
00001 /* zerrgex.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_n1 = -1;
00035 static integer c__0 = 0;
00036 static integer c__1 = 1;
00037 static integer c__3 = 3;
00038 static integer c__4 = 4;
00039 static integer c__5 = 5;
00040 
00041 /* Subroutine */ int zerrge_(char *path, integer *nunit)
00042 {
00043     /* System generated locals */
00044     integer i__1;
00045     doublereal d__1, d__2;
00046     doublecomplex z__1;
00047 
00048     /* Builtin functions */
00049     integer s_wsle(cilist *), e_wsle(void);
00050     /* Subroutine */ int s_copy(char *, char *, ftnlen, ftnlen);
00051 
00052     /* Local variables */
00053     doublecomplex a[16] /* was [4][4] */, b[4];
00054     integer i__, j;
00055     doublereal r__[4];
00056     doublecomplex w[8], x[4];
00057     char c2[2];
00058     doublereal r1[4], r2[4];
00059     doublecomplex af[16]        /* was [4][4] */;
00060     char eq[1];
00061     doublereal cs[4];
00062     integer ip[4];
00063     doublereal rs[4];
00064     doublecomplex err_bnds_c__[12]      /* was [4][3] */;
00065     integer n_err_bnds__;
00066     doublecomplex err_bnds_n__[12]      /* was [4][3] */;
00067     doublereal berr;
00068     integer info;
00069     doublereal anrm, ccond, rcond;
00070     extern /* Subroutine */ int zgbtf2_(integer *, integer *, integer *, 
00071             integer *, doublecomplex *, integer *, integer *, integer *), 
00072             zgetf2_(integer *, integer *, doublecomplex *, integer *, integer 
00073             *, integer *), alaesm_(char *, logical *, integer *);
00074     extern logical lsamen_(integer *, char *, char *);
00075     extern /* Subroutine */ int zgbcon_(char *, integer *, integer *, integer 
00076             *, doublecomplex *, integer *, integer *, doublereal *, 
00077             doublereal *, doublecomplex *, doublereal *, integer *);
00078     doublecomplex params;
00079     extern /* Subroutine */ int chkxer_(char *, integer *, integer *, logical 
00080             *, logical *), zgecon_(char *, integer *, doublecomplex *, 
00081              integer *, doublereal *, doublereal *, doublecomplex *, 
00082             doublereal *, integer *), zgbequ_(integer *, integer *, 
00083             integer *, integer *, doublecomplex *, integer *, doublereal *, 
00084             doublereal *, doublereal *, doublereal *, doublereal *, integer *)
00085             , zgbrfs_(char *, integer *, integer *, integer *, integer *, 
00086             doublecomplex *, integer *, doublecomplex *, integer *, integer *, 
00087              doublecomplex *, integer *, doublecomplex *, integer *, 
00088             doublereal *, doublereal *, doublecomplex *, doublereal *, 
00089             integer *), zgbtrf_(integer *, integer *, integer *, 
00090             integer *, doublecomplex *, integer *, integer *, integer *), 
00091             zgeequ_(integer *, integer *, doublecomplex *, integer *, 
00092             doublereal *, doublereal *, doublereal *, doublereal *, 
00093             doublereal *, integer *), zgerfs_(char *, integer *, integer *, 
00094             doublecomplex *, integer *, doublecomplex *, integer *, integer *, 
00095              doublecomplex *, integer *, doublecomplex *, integer *, 
00096             doublereal *, doublereal *, doublecomplex *, doublereal *, 
00097             integer *), zgetrf_(integer *, integer *, doublecomplex *, 
00098              integer *, integer *, integer *), zgetri_(integer *, 
00099             doublecomplex *, integer *, integer *, doublecomplex *, integer *, 
00100              integer *), zgbtrs_(char *, integer *, integer *, integer *, 
00101             integer *, doublecomplex *, integer *, integer *, doublecomplex *, 
00102              integer *, integer *), zgetrs_(char *, integer *, 
00103             integer *, doublecomplex *, integer *, integer *, doublecomplex *, 
00104              integer *, integer *), zgbequb_(integer *, integer *, 
00105             integer *, integer *, doublecomplex *, integer *, doublereal *, 
00106             doublereal *, doublereal *, doublereal *, doublereal *, integer *)
00107             ;
00108     integer nparams;
00109     extern /* Subroutine */ int zgeequb_(integer *, integer *, doublecomplex *
00110 , integer *, doublereal *, doublereal *, doublereal *, doublereal 
00111             *, doublereal *, integer *), zgbrfsx_(char *, char *, integer *, 
00112             integer *, integer *, integer *, doublecomplex *, integer *, 
00113             doublecomplex *, integer *, integer *, doublereal *, doublereal *, 
00114              doublecomplex *, integer *, doublecomplex *, integer *, 
00115             doublereal *, doublereal *, integer *, doublecomplex *, 
00116             doublecomplex *, integer *, doublecomplex *, doublecomplex *, 
00117             doublereal *, integer *), zgerfsx_(char *, char *, 
00118              integer *, integer *, doublecomplex *, integer *, doublecomplex *
00119 , integer *, integer *, doublereal *, doublereal *, doublecomplex 
00120             *, integer *, doublecomplex *, integer *, doublereal *, 
00121             doublereal *, integer *, doublecomplex *, doublecomplex *, 
00122             integer *, doublecomplex *, doublecomplex *, doublereal *, 
00123             integer *);
00124 
00125     /* Fortran I/O blocks */
00126     static cilist io___1 = { 0, 0, 0, 0, 0 };
00127 
00128 
00129 
00130 /*  -- LAPACK test routine (version 3.1) -- */
00131 /*     Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. */
00132 /*     November 2006 */
00133 
00134 /*     .. Scalar Arguments .. */
00135 /*     .. */
00136 
00137 /*  Purpose */
00138 /*  ======= */
00139 
00140 /*  ZERRGE tests the error exits for the COMPLEX*16 routines */
00141 /*  for general matrices. */
00142 
00143 /*  Note that this file is used only when the XBLAS are available, */
00144 /*  otherwise zerrge.f defines this subroutine. */
00145 
00146 /*  Arguments */
00147 /*  ========= */
00148 
00149 /*  PATH    (input) CHARACTER*3 */
00150 /*          The LAPACK path name for the routines to be tested. */
00151 
00152 /*  NUNIT   (input) INTEGER */
00153 /*          The unit number for output. */
00154 
00155 /*  ===================================================================== */
00156 
00157 /*     .. Parameters .. */
00158 /*     .. */
00159 /*     .. Local Scalars .. */
00160 /*     .. */
00161 /*     .. Local Arrays .. */
00162 /*     .. */
00163 /*     .. External Functions .. */
00164 /*     .. */
00165 /*     .. External Subroutines .. */
00166 /*     .. */
00167 /*     .. Scalars in Common .. */
00168 /*     .. */
00169 /*     .. Common blocks .. */
00170 /*     .. */
00171 /*     .. Intrinsic Functions .. */
00172 /*     .. */
00173 /*     .. Executable Statements .. */
00174 
00175     infoc_1.nout = *nunit;
00176     io___1.ciunit = infoc_1.nout;
00177     s_wsle(&io___1);
00178     e_wsle();
00179     s_copy(c2, path + 1, (ftnlen)2, (ftnlen)2);
00180 
00181 /*     Set the variables to innocuous values. */
00182 
00183     for (j = 1; j <= 4; ++j) {
00184         for (i__ = 1; i__ <= 4; ++i__) {
00185             i__1 = i__ + (j << 2) - 5;
00186             d__1 = 1. / (doublereal) (i__ + j);
00187             d__2 = -1. / (doublereal) (i__ + j);
00188             z__1.r = d__1, z__1.i = d__2;
00189             a[i__1].r = z__1.r, a[i__1].i = z__1.i;
00190             i__1 = i__ + (j << 2) - 5;
00191             d__1 = 1. / (doublereal) (i__ + j);
00192             d__2 = -1. / (doublereal) (i__ + j);
00193             z__1.r = d__1, z__1.i = d__2;
00194             af[i__1].r = z__1.r, af[i__1].i = z__1.i;
00195 /* L10: */
00196         }
00197         i__1 = j - 1;
00198         b[i__1].r = 0., b[i__1].i = 0.;
00199         r1[j - 1] = 0.;
00200         r2[j - 1] = 0.;
00201         i__1 = j - 1;
00202         w[i__1].r = 0., w[i__1].i = 0.;
00203         i__1 = j - 1;
00204         x[i__1].r = 0., x[i__1].i = 0.;
00205         cs[j - 1] = 0.;
00206         rs[j - 1] = 0.;
00207         ip[j - 1] = j;
00208 /* L20: */
00209     }
00210     infoc_1.ok = TRUE_;
00211 
00212 /*     Test error exits of the routines that use the LU decomposition */
00213 /*     of a general matrix. */
00214 
00215     if (lsamen_(&c__2, c2, "GE")) {
00216 
00217 /*        ZGETRF */
00218 
00219         s_copy(srnamc_1.srnamt, "ZGETRF", (ftnlen)32, (ftnlen)6);
00220         infoc_1.infot = 1;
00221         zgetrf_(&c_n1, &c__0, a, &c__1, ip, &info);
00222         chkxer_("ZGETRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00223                 infoc_1.ok);
00224         infoc_1.infot = 2;
00225         zgetrf_(&c__0, &c_n1, a, &c__1, ip, &info);
00226         chkxer_("ZGETRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00227                 infoc_1.ok);
00228         infoc_1.infot = 4;
00229         zgetrf_(&c__2, &c__1, a, &c__1, ip, &info);
00230         chkxer_("ZGETRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00231                 infoc_1.ok);
00232 
00233 /*        ZGETF2 */
00234 
00235         s_copy(srnamc_1.srnamt, "ZGETF2", (ftnlen)32, (ftnlen)6);
00236         infoc_1.infot = 1;
00237         zgetf2_(&c_n1, &c__0, a, &c__1, ip, &info);
00238         chkxer_("ZGETF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00239                 infoc_1.ok);
00240         infoc_1.infot = 2;
00241         zgetf2_(&c__0, &c_n1, a, &c__1, ip, &info);
00242         chkxer_("ZGETF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00243                 infoc_1.ok);
00244         infoc_1.infot = 4;
00245         zgetf2_(&c__2, &c__1, a, &c__1, ip, &info);
00246         chkxer_("ZGETF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00247                 infoc_1.ok);
00248 
00249 /*        ZGETRI */
00250 
00251         s_copy(srnamc_1.srnamt, "ZGETRI", (ftnlen)32, (ftnlen)6);
00252         infoc_1.infot = 1;
00253         zgetri_(&c_n1, a, &c__1, ip, w, &c__1, &info);
00254         chkxer_("ZGETRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00255                 infoc_1.ok);
00256         infoc_1.infot = 3;
00257         zgetri_(&c__2, a, &c__1, ip, w, &c__2, &info);
00258         chkxer_("ZGETRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00259                 infoc_1.ok);
00260         infoc_1.infot = 6;
00261         zgetri_(&c__2, a, &c__2, ip, w, &c__1, &info);
00262         chkxer_("ZGETRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00263                 infoc_1.ok);
00264 
00265 /*        ZGETRS */
00266 
00267         s_copy(srnamc_1.srnamt, "ZGETRS", (ftnlen)32, (ftnlen)6);
00268         infoc_1.infot = 1;
00269         zgetrs_("/", &c__0, &c__0, a, &c__1, ip, b, &c__1, &info);
00270         chkxer_("ZGETRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00271                 infoc_1.ok);
00272         infoc_1.infot = 2;
00273         zgetrs_("N", &c_n1, &c__0, a, &c__1, ip, b, &c__1, &info);
00274         chkxer_("ZGETRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00275                 infoc_1.ok);
00276         infoc_1.infot = 3;
00277         zgetrs_("N", &c__0, &c_n1, a, &c__1, ip, b, &c__1, &info);
00278         chkxer_("ZGETRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00279                 infoc_1.ok);
00280         infoc_1.infot = 5;
00281         zgetrs_("N", &c__2, &c__1, a, &c__1, ip, b, &c__2, &info);
00282         chkxer_("ZGETRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00283                 infoc_1.ok);
00284         infoc_1.infot = 8;
00285         zgetrs_("N", &c__2, &c__1, a, &c__2, ip, b, &c__1, &info);
00286         chkxer_("ZGETRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00287                 infoc_1.ok);
00288 
00289 /*        ZGERFS */
00290 
00291         s_copy(srnamc_1.srnamt, "ZGERFS", (ftnlen)32, (ftnlen)6);
00292         infoc_1.infot = 1;
00293         zgerfs_("/", &c__0, &c__0, a, &c__1, af, &c__1, ip, b, &c__1, x, &
00294                 c__1, r1, r2, w, r__, &info);
00295         chkxer_("ZGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00296                 infoc_1.ok);
00297         infoc_1.infot = 2;
00298         zgerfs_("N", &c_n1, &c__0, a, &c__1, af, &c__1, ip, b, &c__1, x, &
00299                 c__1, r1, r2, w, r__, &info);
00300         chkxer_("ZGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00301                 infoc_1.ok);
00302         infoc_1.infot = 3;
00303         zgerfs_("N", &c__0, &c_n1, a, &c__1, af, &c__1, ip, b, &c__1, x, &
00304                 c__1, r1, r2, w, r__, &info);
00305         chkxer_("ZGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00306                 infoc_1.ok);
00307         infoc_1.infot = 5;
00308         zgerfs_("N", &c__2, &c__1, a, &c__1, af, &c__2, ip, b, &c__2, x, &
00309                 c__2, r1, r2, w, r__, &info);
00310         chkxer_("ZGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00311                 infoc_1.ok);
00312         infoc_1.infot = 7;
00313         zgerfs_("N", &c__2, &c__1, a, &c__2, af, &c__1, ip, b, &c__2, x, &
00314                 c__2, r1, r2, w, r__, &info);
00315         chkxer_("ZGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00316                 infoc_1.ok);
00317         infoc_1.infot = 10;
00318         zgerfs_("N", &c__2, &c__1, a, &c__2, af, &c__2, ip, b, &c__1, x, &
00319                 c__2, r1, r2, w, r__, &info);
00320         chkxer_("ZGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00321                 infoc_1.ok);
00322         infoc_1.infot = 12;
00323         zgerfs_("N", &c__2, &c__1, a, &c__2, af, &c__2, ip, b, &c__2, x, &
00324                 c__1, r1, r2, w, r__, &info);
00325         chkxer_("ZGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00326                 infoc_1.ok);
00327 
00328 /*        ZGERFSX */
00329 
00330         n_err_bnds__ = 3;
00331         nparams = 0;
00332         s_copy(srnamc_1.srnamt, "ZGERFSX", (ftnlen)32, (ftnlen)7);
00333         infoc_1.infot = 1;
00334         zgerfsx_("/", eq, &c__0, &c__0, a, &c__1, af, &c__1, ip, rs, cs, b, &
00335                 c__1, x, &c__1, &rcond, &berr, &n_err_bnds__, err_bnds_n__, 
00336                 err_bnds_c__, &nparams, &params, w, r__, &info);
00337         chkxer_("ZGERFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00338                 infoc_1.ok);
00339         infoc_1.infot = 2;
00340         *(unsigned char *)eq = '/';
00341         zgerfsx_("N", eq, &c__2, &c__1, a, &c__1, af, &c__2, ip, rs, cs, b, &
00342                 c__2, x, &c__2, &rcond, &berr, &n_err_bnds__, err_bnds_n__, 
00343                 err_bnds_c__, &nparams, &params, w, r__, &info);
00344         chkxer_("ZGERFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00345                 infoc_1.ok);
00346         infoc_1.infot = 3;
00347         *(unsigned char *)eq = 'R';
00348         zgerfsx_("N", eq, &c_n1, &c__0, a, &c__1, af, &c__1, ip, rs, cs, b, &
00349                 c__1, x, &c__1, &rcond, &berr, &n_err_bnds__, err_bnds_n__, 
00350                 err_bnds_c__, &nparams, &params, w, r__, &info);
00351         chkxer_("ZGERFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00352                 infoc_1.ok);
00353         infoc_1.infot = 4;
00354         zgerfsx_("N", eq, &c__0, &c_n1, a, &c__1, af, &c__1, ip, rs, cs, b, &
00355                 c__1, x, &c__1, &rcond, &berr, &n_err_bnds__, err_bnds_n__, 
00356                 err_bnds_c__, &nparams, &params, w, r__, &info);
00357         chkxer_("ZGERFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00358                 infoc_1.ok);
00359         infoc_1.infot = 6;
00360         zgerfsx_("N", eq, &c__2, &c__1, a, &c__1, af, &c__2, ip, rs, cs, b, &
00361                 c__2, x, &c__2, &rcond, &berr, &n_err_bnds__, err_bnds_n__, 
00362                 err_bnds_c__, &nparams, &params, w, r__, &info);
00363         chkxer_("ZGERFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00364                 infoc_1.ok);
00365         infoc_1.infot = 8;
00366         zgerfsx_("N", eq, &c__2, &c__1, a, &c__2, af, &c__1, ip, rs, cs, b, &
00367                 c__2, x, &c__2, &rcond, &berr, &n_err_bnds__, err_bnds_n__, 
00368                 err_bnds_c__, &nparams, &params, w, r__, &info);
00369         chkxer_("ZGERFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00370                 infoc_1.ok);
00371         infoc_1.infot = 13;
00372         *(unsigned char *)eq = 'C';
00373         zgerfsx_("N", eq, &c__2, &c__1, a, &c__2, af, &c__2, ip, rs, cs, b, &
00374                 c__1, x, &c__2, &rcond, &berr, &n_err_bnds__, err_bnds_n__, 
00375                 err_bnds_c__, &nparams, &params, w, r__, &info);
00376         chkxer_("ZGERFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00377                 infoc_1.ok);
00378         infoc_1.infot = 15;
00379         zgerfsx_("N", eq, &c__2, &c__1, a, &c__2, af, &c__2, ip, rs, cs, b, &
00380                 c__2, x, &c__1, &rcond, &berr, &n_err_bnds__, err_bnds_n__, 
00381                 err_bnds_c__, &nparams, &params, w, r__, &info);
00382         chkxer_("ZGERFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00383                 infoc_1.ok);
00384 
00385 /*        ZGECON */
00386 
00387         s_copy(srnamc_1.srnamt, "ZGECON", (ftnlen)32, (ftnlen)6);
00388         infoc_1.infot = 1;
00389         zgecon_("/", &c__0, a, &c__1, &anrm, &rcond, w, r__, &info)
00390                 ;
00391         chkxer_("ZGECON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00392                 infoc_1.ok);
00393         infoc_1.infot = 2;
00394         zgecon_("1", &c_n1, a, &c__1, &anrm, &rcond, w, r__, &info)
00395                 ;
00396         chkxer_("ZGECON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00397                 infoc_1.ok);
00398         infoc_1.infot = 4;
00399         zgecon_("1", &c__2, a, &c__1, &anrm, &rcond, w, r__, &info)
00400                 ;
00401         chkxer_("ZGECON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00402                 infoc_1.ok);
00403 
00404 /*        ZGEEQU */
00405 
00406         s_copy(srnamc_1.srnamt, "ZGEEQU", (ftnlen)32, (ftnlen)6);
00407         infoc_1.infot = 1;
00408         zgeequ_(&c_n1, &c__0, a, &c__1, r1, r2, &rcond, &ccond, &anrm, &info);
00409         chkxer_("ZGEEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00410                 infoc_1.ok);
00411         infoc_1.infot = 2;
00412         zgeequ_(&c__0, &c_n1, a, &c__1, r1, r2, &rcond, &ccond, &anrm, &info);
00413         chkxer_("ZGEEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00414                 infoc_1.ok);
00415         infoc_1.infot = 4;
00416         zgeequ_(&c__2, &c__2, a, &c__1, r1, r2, &rcond, &ccond, &anrm, &info);
00417         chkxer_("ZGEEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00418                 infoc_1.ok);
00419 
00420 /*        ZGEEQUB */
00421 
00422         s_copy(srnamc_1.srnamt, "ZGEEQUB", (ftnlen)32, (ftnlen)7);
00423         infoc_1.infot = 1;
00424         zgeequb_(&c_n1, &c__0, a, &c__1, r1, r2, &rcond, &ccond, &anrm, &info)
00425                 ;
00426         chkxer_("ZGEEQUB", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00427                 infoc_1.ok);
00428         infoc_1.infot = 2;
00429         zgeequb_(&c__0, &c_n1, a, &c__1, r1, r2, &rcond, &ccond, &anrm, &info)
00430                 ;
00431         chkxer_("ZGEEQUB", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00432                 infoc_1.ok);
00433         infoc_1.infot = 4;
00434         zgeequb_(&c__2, &c__2, a, &c__1, r1, r2, &rcond, &ccond, &anrm, &info)
00435                 ;
00436         chkxer_("ZGEEQUB", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00437                 infoc_1.ok);
00438 
00439 /*     Test error exits of the routines that use the LU decomposition */
00440 /*     of a general band matrix. */
00441 
00442     } else if (lsamen_(&c__2, c2, "GB")) {
00443 
00444 /*        ZGBTRF */
00445 
00446         s_copy(srnamc_1.srnamt, "ZGBTRF", (ftnlen)32, (ftnlen)6);
00447         infoc_1.infot = 1;
00448         zgbtrf_(&c_n1, &c__0, &c__0, &c__0, a, &c__1, ip, &info);
00449         chkxer_("ZGBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00450                 infoc_1.ok);
00451         infoc_1.infot = 2;
00452         zgbtrf_(&c__0, &c_n1, &c__0, &c__0, a, &c__1, ip, &info);
00453         chkxer_("ZGBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00454                 infoc_1.ok);
00455         infoc_1.infot = 3;
00456         zgbtrf_(&c__1, &c__1, &c_n1, &c__0, a, &c__1, ip, &info);
00457         chkxer_("ZGBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00458                 infoc_1.ok);
00459         infoc_1.infot = 4;
00460         zgbtrf_(&c__1, &c__1, &c__0, &c_n1, a, &c__1, ip, &info);
00461         chkxer_("ZGBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00462                 infoc_1.ok);
00463         infoc_1.infot = 6;
00464         zgbtrf_(&c__2, &c__2, &c__1, &c__1, a, &c__3, ip, &info);
00465         chkxer_("ZGBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00466                 infoc_1.ok);
00467 
00468 /*        ZGBTF2 */
00469 
00470         s_copy(srnamc_1.srnamt, "ZGBTF2", (ftnlen)32, (ftnlen)6);
00471         infoc_1.infot = 1;
00472         zgbtf2_(&c_n1, &c__0, &c__0, &c__0, a, &c__1, ip, &info);
00473         chkxer_("ZGBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00474                 infoc_1.ok);
00475         infoc_1.infot = 2;
00476         zgbtf2_(&c__0, &c_n1, &c__0, &c__0, a, &c__1, ip, &info);
00477         chkxer_("ZGBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00478                 infoc_1.ok);
00479         infoc_1.infot = 3;
00480         zgbtf2_(&c__1, &c__1, &c_n1, &c__0, a, &c__1, ip, &info);
00481         chkxer_("ZGBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00482                 infoc_1.ok);
00483         infoc_1.infot = 4;
00484         zgbtf2_(&c__1, &c__1, &c__0, &c_n1, a, &c__1, ip, &info);
00485         chkxer_("ZGBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00486                 infoc_1.ok);
00487         infoc_1.infot = 6;
00488         zgbtf2_(&c__2, &c__2, &c__1, &c__1, a, &c__3, ip, &info);
00489         chkxer_("ZGBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00490                 infoc_1.ok);
00491 
00492 /*        ZGBTRS */
00493 
00494         s_copy(srnamc_1.srnamt, "ZGBTRS", (ftnlen)32, (ftnlen)6);
00495         infoc_1.infot = 1;
00496         zgbtrs_("/", &c__0, &c__0, &c__0, &c__1, a, &c__1, ip, b, &c__1, &
00497                 info);
00498         chkxer_("ZGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00499                 infoc_1.ok);
00500         infoc_1.infot = 2;
00501         zgbtrs_("N", &c_n1, &c__0, &c__0, &c__1, a, &c__1, ip, b, &c__1, &
00502                 info);
00503         chkxer_("ZGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00504                 infoc_1.ok);
00505         infoc_1.infot = 3;
00506         zgbtrs_("N", &c__1, &c_n1, &c__0, &c__1, a, &c__1, ip, b, &c__1, &
00507                 info);
00508         chkxer_("ZGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00509                 infoc_1.ok);
00510         infoc_1.infot = 4;
00511         zgbtrs_("N", &c__1, &c__0, &c_n1, &c__1, a, &c__1, ip, b, &c__1, &
00512                 info);
00513         chkxer_("ZGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00514                 infoc_1.ok);
00515         infoc_1.infot = 5;
00516         zgbtrs_("N", &c__1, &c__0, &c__0, &c_n1, a, &c__1, ip, b, &c__1, &
00517                 info);
00518         chkxer_("ZGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00519                 infoc_1.ok);
00520         infoc_1.infot = 7;
00521         zgbtrs_("N", &c__2, &c__1, &c__1, &c__1, a, &c__3, ip, b, &c__2, &
00522                 info);
00523         chkxer_("ZGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00524                 infoc_1.ok);
00525         infoc_1.infot = 10;
00526         zgbtrs_("N", &c__2, &c__0, &c__0, &c__1, a, &c__1, ip, b, &c__1, &
00527                 info);
00528         chkxer_("ZGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00529                 infoc_1.ok);
00530 
00531 /*        ZGBRFS */
00532 
00533         s_copy(srnamc_1.srnamt, "ZGBRFS", (ftnlen)32, (ftnlen)6);
00534         infoc_1.infot = 1;
00535         zgbrfs_("/", &c__0, &c__0, &c__0, &c__0, a, &c__1, af, &c__1, ip, b, &
00536                 c__1, x, &c__1, r1, r2, w, r__, &info);
00537         chkxer_("ZGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00538                 infoc_1.ok);
00539         infoc_1.infot = 2;
00540         zgbrfs_("N", &c_n1, &c__0, &c__0, &c__0, a, &c__1, af, &c__1, ip, b, &
00541                 c__1, x, &c__1, r1, r2, w, r__, &info);
00542         chkxer_("ZGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00543                 infoc_1.ok);
00544         infoc_1.infot = 3;
00545         zgbrfs_("N", &c__1, &c_n1, &c__0, &c__0, a, &c__1, af, &c__1, ip, b, &
00546                 c__1, x, &c__1, r1, r2, w, r__, &info);
00547         chkxer_("ZGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00548                 infoc_1.ok);
00549         infoc_1.infot = 4;
00550         zgbrfs_("N", &c__1, &c__0, &c_n1, &c__0, a, &c__1, af, &c__1, ip, b, &
00551                 c__1, x, &c__1, r1, r2, w, r__, &info);
00552         chkxer_("ZGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00553                 infoc_1.ok);
00554         infoc_1.infot = 5;
00555         zgbrfs_("N", &c__1, &c__0, &c__0, &c_n1, a, &c__1, af, &c__1, ip, b, &
00556                 c__1, x, &c__1, r1, r2, w, r__, &info);
00557         chkxer_("ZGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00558                 infoc_1.ok);
00559         infoc_1.infot = 7;
00560         zgbrfs_("N", &c__2, &c__1, &c__1, &c__1, a, &c__2, af, &c__4, ip, b, &
00561                 c__2, x, &c__2, r1, r2, w, r__, &info);
00562         chkxer_("ZGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00563                 infoc_1.ok);
00564         infoc_1.infot = 9;
00565         zgbrfs_("N", &c__2, &c__1, &c__1, &c__1, a, &c__3, af, &c__3, ip, b, &
00566                 c__2, x, &c__2, r1, r2, w, r__, &info);
00567         chkxer_("ZGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00568                 infoc_1.ok);
00569         infoc_1.infot = 12;
00570         zgbrfs_("N", &c__2, &c__0, &c__0, &c__1, a, &c__1, af, &c__1, ip, b, &
00571                 c__1, x, &c__2, r1, r2, w, r__, &info);
00572         chkxer_("ZGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00573                 infoc_1.ok);
00574         infoc_1.infot = 14;
00575         zgbrfs_("N", &c__2, &c__0, &c__0, &c__1, a, &c__1, af, &c__1, ip, b, &
00576                 c__2, x, &c__1, r1, r2, w, r__, &info);
00577         chkxer_("ZGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00578                 infoc_1.ok);
00579 
00580 /*        ZGBRFSX */
00581 
00582         n_err_bnds__ = 3;
00583         nparams = 0;
00584         s_copy(srnamc_1.srnamt, "ZGBRFSX", (ftnlen)32, (ftnlen)7);
00585         infoc_1.infot = 1;
00586         zgbrfsx_("/", eq, &c__0, &c__0, &c__0, &c__0, a, &c__1, af, &c__1, ip, 
00587                  rs, cs, b, &c__1, x, &c__1, &rcond, &berr, &n_err_bnds__, 
00588                 err_bnds_n__, err_bnds_c__, &nparams, &params, w, r__, &info);
00589         chkxer_("ZGBRFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00590                 infoc_1.ok);
00591         infoc_1.infot = 2;
00592         *(unsigned char *)eq = '/';
00593         zgbrfsx_("N", eq, &c__2, &c__1, &c__1, &c__1, a, &c__1, af, &c__2, ip, 
00594                  rs, cs, b, &c__2, x, &c__2, &rcond, &berr, &n_err_bnds__, 
00595                 err_bnds_n__, err_bnds_c__, &nparams, &params, w, r__, &info);
00596         chkxer_("ZGBRFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00597                 infoc_1.ok);
00598         infoc_1.infot = 3;
00599         *(unsigned char *)eq = 'R';
00600         zgbrfsx_("N", eq, &c_n1, &c__1, &c__1, &c__0, a, &c__1, af, &c__1, ip, 
00601                  rs, cs, b, &c__1, x, &c__1, &rcond, &berr, &n_err_bnds__, 
00602                 err_bnds_n__, err_bnds_c__, &nparams, &params, w, r__, &info);
00603         chkxer_("ZGBRFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00604                 infoc_1.ok);
00605         infoc_1.infot = 4;
00606         *(unsigned char *)eq = 'R';
00607         zgbrfsx_("N", eq, &c__2, &c_n1, &c__1, &c__1, a, &c__3, af, &c__4, ip, 
00608                  rs, cs, b, &c__1, x, &c__1, &rcond, &berr, &n_err_bnds__, 
00609                 err_bnds_n__, err_bnds_c__, &nparams, &params, w, r__, &info);
00610         chkxer_("ZGBRFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00611                 infoc_1.ok);
00612         infoc_1.infot = 5;
00613         *(unsigned char *)eq = 'R';
00614         zgbrfsx_("N", eq, &c__2, &c__1, &c_n1, &c__1, a, &c__3, af, &c__4, ip, 
00615                  rs, cs, b, &c__1, x, &c__1, &rcond, &berr, &n_err_bnds__, 
00616                 err_bnds_n__, err_bnds_c__, &nparams, &params, w, r__, &info);
00617         chkxer_("ZGBRFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00618                 infoc_1.ok);
00619         infoc_1.infot = 6;
00620         zgbrfsx_("N", eq, &c__0, &c__0, &c__0, &c_n1, a, &c__1, af, &c__1, ip, 
00621                  rs, cs, b, &c__1, x, &c__1, &rcond, &berr, &n_err_bnds__, 
00622                 err_bnds_n__, err_bnds_c__, &nparams, &params, w, r__, &info);
00623         chkxer_("ZGBRFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00624                 infoc_1.ok);
00625         infoc_1.infot = 8;
00626         zgbrfsx_("N", eq, &c__2, &c__1, &c__1, &c__1, a, &c__1, af, &c__2, ip, 
00627                  rs, cs, b, &c__2, x, &c__2, &rcond, &berr, &n_err_bnds__, 
00628                 err_bnds_n__, err_bnds_c__, &nparams, &params, w, r__, &info);
00629         chkxer_("ZGBRFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00630                 infoc_1.ok);
00631         infoc_1.infot = 10;
00632         zgbrfsx_("N", eq, &c__2, &c__1, &c__1, &c__1, a, &c__3, af, &c__3, ip, 
00633                  rs, cs, b, &c__2, x, &c__2, &rcond, &berr, &n_err_bnds__, 
00634                 err_bnds_n__, err_bnds_c__, &nparams, &params, w, r__, &info);
00635         chkxer_("ZGBRFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00636                 infoc_1.ok);
00637         infoc_1.infot = 13;
00638         *(unsigned char *)eq = 'C';
00639         zgbrfsx_("N", eq, &c__2, &c__1, &c__1, &c__1, a, &c__3, af, &c__5, ip, 
00640                  rs, cs, b, &c__1, x, &c__2, &rcond, &berr, &n_err_bnds__, 
00641                 err_bnds_n__, err_bnds_c__, &nparams, &params, w, r__, &info);
00642         chkxer_("ZGBRFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00643                 infoc_1.ok);
00644         infoc_1.infot = 15;
00645         zgbrfsx_("N", eq, &c__2, &c__1, &c__1, &c__1, a, &c__3, af, &c__5, ip, 
00646                  rs, cs, b, &c__2, x, &c__1, &rcond, &berr, &n_err_bnds__, 
00647                 err_bnds_n__, err_bnds_c__, &nparams, &params, w, r__, &info);
00648         chkxer_("ZGBRFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00649                 infoc_1.ok);
00650 
00651 /*        ZGBCON */
00652 
00653         s_copy(srnamc_1.srnamt, "ZGBCON", (ftnlen)32, (ftnlen)6);
00654         infoc_1.infot = 1;
00655         zgbcon_("/", &c__0, &c__0, &c__0, a, &c__1, ip, &anrm, &rcond, w, r__, 
00656                  &info);
00657         chkxer_("ZGBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00658                 infoc_1.ok);
00659         infoc_1.infot = 2;
00660         zgbcon_("1", &c_n1, &c__0, &c__0, a, &c__1, ip, &anrm, &rcond, w, r__, 
00661                  &info);
00662         chkxer_("ZGBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00663                 infoc_1.ok);
00664         infoc_1.infot = 3;
00665         zgbcon_("1", &c__1, &c_n1, &c__0, a, &c__1, ip, &anrm, &rcond, w, r__, 
00666                  &info);
00667         chkxer_("ZGBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00668                 infoc_1.ok);
00669         infoc_1.infot = 4;
00670         zgbcon_("1", &c__1, &c__0, &c_n1, a, &c__1, ip, &anrm, &rcond, w, r__, 
00671                  &info);
00672         chkxer_("ZGBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00673                 infoc_1.ok);
00674         infoc_1.infot = 6;
00675         zgbcon_("1", &c__2, &c__1, &c__1, a, &c__3, ip, &anrm, &rcond, w, r__, 
00676                  &info);
00677         chkxer_("ZGBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00678                 infoc_1.ok);
00679 
00680 /*        ZGBEQU */
00681 
00682         s_copy(srnamc_1.srnamt, "ZGBEQU", (ftnlen)32, (ftnlen)6);
00683         infoc_1.infot = 1;
00684         zgbequ_(&c_n1, &c__0, &c__0, &c__0, a, &c__1, r1, r2, &rcond, &ccond, 
00685                 &anrm, &info);
00686         chkxer_("ZGBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00687                 infoc_1.ok);
00688         infoc_1.infot = 2;
00689         zgbequ_(&c__0, &c_n1, &c__0, &c__0, a, &c__1, r1, r2, &rcond, &ccond, 
00690                 &anrm, &info);
00691         chkxer_("ZGBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00692                 infoc_1.ok);
00693         infoc_1.infot = 3;
00694         zgbequ_(&c__1, &c__1, &c_n1, &c__0, a, &c__1, r1, r2, &rcond, &ccond, 
00695                 &anrm, &info);
00696         chkxer_("ZGBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00697                 infoc_1.ok);
00698         infoc_1.infot = 4;
00699         zgbequ_(&c__1, &c__1, &c__0, &c_n1, a, &c__1, r1, r2, &rcond, &ccond, 
00700                 &anrm, &info);
00701         chkxer_("ZGBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00702                 infoc_1.ok);
00703         infoc_1.infot = 6;
00704         zgbequ_(&c__2, &c__2, &c__1, &c__1, a, &c__2, r1, r2, &rcond, &ccond, 
00705                 &anrm, &info);
00706         chkxer_("ZGBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00707                 infoc_1.ok);
00708 
00709 /*        ZGBEQUB */
00710 
00711         s_copy(srnamc_1.srnamt, "ZGBEQUB", (ftnlen)32, (ftnlen)7);
00712         infoc_1.infot = 1;
00713         zgbequb_(&c_n1, &c__0, &c__0, &c__0, a, &c__1, r1, r2, &rcond, &ccond, 
00714                  &anrm, &info);
00715         chkxer_("ZGBEQUB", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00716                 infoc_1.ok);
00717         infoc_1.infot = 2;
00718         zgbequb_(&c__0, &c_n1, &c__0, &c__0, a, &c__1, r1, r2, &rcond, &ccond, 
00719                  &anrm, &info);
00720         chkxer_("ZGBEQUB", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00721                 infoc_1.ok);
00722         infoc_1.infot = 3;
00723         zgbequb_(&c__1, &c__1, &c_n1, &c__0, a, &c__1, r1, r2, &rcond, &ccond, 
00724                  &anrm, &info);
00725         chkxer_("ZGBEQUB", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00726                 infoc_1.ok);
00727         infoc_1.infot = 4;
00728         zgbequb_(&c__1, &c__1, &c__0, &c_n1, a, &c__1, r1, r2, &rcond, &ccond, 
00729                  &anrm, &info);
00730         chkxer_("ZGBEQUB", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00731                 infoc_1.ok);
00732         infoc_1.infot = 6;
00733         zgbequb_(&c__2, &c__2, &c__1, &c__1, a, &c__2, r1, r2, &rcond, &ccond, 
00734                  &anrm, &info);
00735         chkxer_("ZGBEQUB", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00736                 infoc_1.ok);
00737     }
00738 
00739 /*     Print a summary line. */
00740 
00741     alaesm_(path, &infoc_1.ok, &infoc_1.nout);
00742 
00743     return 0;
00744 
00745 /*     End of ZERRGE */
00746 
00747 } /* zerrge_ */


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