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


swiftnav
Author(s):
autogenerated on Sat Jun 8 2019 18:55:41