derrge.c
Go to the documentation of this file.
00001 /* derrge.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 
00041 /* Subroutine */ int derrge_(char *path, integer *nunit)
00042 {
00043     /* Builtin functions */
00044     integer s_wsle(cilist *), e_wsle(void);
00045     /* Subroutine */ int s_copy(char *, char *, ftnlen, ftnlen);
00046 
00047     /* Local variables */
00048     doublereal a[16]    /* was [4][4] */, b[4];
00049     integer i__, j;
00050     doublereal w[12], x[4];
00051     char c2[2];
00052     doublereal r1[4], r2[4], af[16]     /* was [4][4] */;
00053     integer ip[4], iw[4], info;
00054     doublereal anrm, ccond, rcond;
00055     extern /* Subroutine */ int dgbtf2_(integer *, integer *, integer *, 
00056             integer *, doublereal *, integer *, integer *, integer *), 
00057             dgetf2_(integer *, integer *, doublereal *, integer *, integer *, 
00058             integer *), dgbcon_(char *, integer *, integer *, integer *, 
00059             doublereal *, integer *, integer *, doublereal *, doublereal *, 
00060             doublereal *, integer *, integer *), dgecon_(char *, 
00061             integer *, doublereal *, integer *, doublereal *, doublereal *, 
00062             doublereal *, integer *, integer *), alaesm_(char *, 
00063             logical *, integer *), dgbequ_(integer *, integer *, 
00064             integer *, integer *, doublereal *, integer *, doublereal *, 
00065             doublereal *, doublereal *, doublereal *, doublereal *, integer *)
00066             , dgbrfs_(char *, integer *, integer *, integer *, integer *, 
00067             doublereal *, integer *, doublereal *, integer *, integer *, 
00068             doublereal *, integer *, doublereal *, integer *, doublereal *, 
00069             doublereal *, doublereal *, integer *, integer *), 
00070             dgbtrf_(integer *, integer *, integer *, integer *, doublereal *, 
00071             integer *, integer *, integer *), dgeequ_(integer *, integer *, 
00072             doublereal *, integer *, doublereal *, doublereal *, doublereal *, 
00073              doublereal *, doublereal *, integer *), dgerfs_(char *, integer *
00074 , integer *, doublereal *, integer *, doublereal *, integer *, 
00075             integer *, doublereal *, integer *, doublereal *, integer *, 
00076             doublereal *, doublereal *, doublereal *, integer *, integer *), dgetrf_(integer *, integer *, doublereal *, integer *, 
00077             integer *, integer *), dgetri_(integer *, doublereal *, integer *, 
00078              integer *, doublereal *, integer *, integer *);
00079     extern logical lsamen_(integer *, char *, char *);
00080     extern /* Subroutine */ int chkxer_(char *, integer *, integer *, logical 
00081             *, logical *), dgbtrs_(char *, integer *, integer *, 
00082             integer *, integer *, doublereal *, integer *, integer *, 
00083             doublereal *, integer *, integer *), dgetrs_(char *, 
00084             integer *, integer *, doublereal *, integer *, integer *, 
00085             doublereal *, integer *, integer *);
00086 
00087     /* Fortran I/O blocks */
00088     static cilist io___1 = { 0, 0, 0, 0, 0 };
00089 
00090 
00091 
00092 /*  -- LAPACK test routine (version 3.1) -- */
00093 /*     Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. */
00094 /*     November 2006 */
00095 
00096 /*     .. Scalar Arguments .. */
00097 /*     .. */
00098 
00099 /*  Purpose */
00100 /*  ======= */
00101 
00102 /*  DERRGE tests the error exits for the DOUBLE PRECISION routines */
00103 /*  for general matrices. */
00104 
00105 /*  Arguments */
00106 /*  ========= */
00107 
00108 /*  PATH    (input) CHARACTER*3 */
00109 /*          The LAPACK path name for the routines to be tested. */
00110 
00111 /*  NUNIT   (input) INTEGER */
00112 /*          The unit number for output. */
00113 
00114 /*  ===================================================================== */
00115 
00116 /*     .. Parameters .. */
00117 /*     .. */
00118 /*     .. Local Scalars .. */
00119 /*     .. */
00120 /*     .. Local Arrays .. */
00121 /*     .. */
00122 /*     .. External Functions .. */
00123 /*     .. */
00124 /*     .. External Subroutines .. */
00125 /*     .. */
00126 /*     .. Scalars in Common .. */
00127 /*     .. */
00128 /*     .. Common blocks .. */
00129 /*     .. */
00130 /*     .. Intrinsic Functions .. */
00131 /*     .. */
00132 /*     .. Executable Statements .. */
00133 
00134     infoc_1.nout = *nunit;
00135     io___1.ciunit = infoc_1.nout;
00136     s_wsle(&io___1);
00137     e_wsle();
00138     s_copy(c2, path + 1, (ftnlen)2, (ftnlen)2);
00139 
00140 /*     Set the variables to innocuous values. */
00141 
00142     for (j = 1; j <= 4; ++j) {
00143         for (i__ = 1; i__ <= 4; ++i__) {
00144             a[i__ + (j << 2) - 5] = 1. / (doublereal) (i__ + j);
00145             af[i__ + (j << 2) - 5] = 1. / (doublereal) (i__ + j);
00146 /* L10: */
00147         }
00148         b[j - 1] = 0.;
00149         r1[j - 1] = 0.;
00150         r2[j - 1] = 0.;
00151         w[j - 1] = 0.;
00152         x[j - 1] = 0.;
00153         ip[j - 1] = j;
00154         iw[j - 1] = j;
00155 /* L20: */
00156     }
00157     infoc_1.ok = TRUE_;
00158 
00159     if (lsamen_(&c__2, c2, "GE")) {
00160 
00161 /*        Test error exits of the routines that use the LU decomposition */
00162 /*        of a general matrix. */
00163 
00164 /*        DGETRF */
00165 
00166         s_copy(srnamc_1.srnamt, "DGETRF", (ftnlen)32, (ftnlen)6);
00167         infoc_1.infot = 1;
00168         dgetrf_(&c_n1, &c__0, a, &c__1, ip, &info);
00169         chkxer_("DGETRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00170                 infoc_1.ok);
00171         infoc_1.infot = 2;
00172         dgetrf_(&c__0, &c_n1, a, &c__1, ip, &info);
00173         chkxer_("DGETRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00174                 infoc_1.ok);
00175         infoc_1.infot = 4;
00176         dgetrf_(&c__2, &c__1, a, &c__1, ip, &info);
00177         chkxer_("DGETRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00178                 infoc_1.ok);
00179 
00180 /*        DGETF2 */
00181 
00182         s_copy(srnamc_1.srnamt, "DGETF2", (ftnlen)32, (ftnlen)6);
00183         infoc_1.infot = 1;
00184         dgetf2_(&c_n1, &c__0, a, &c__1, ip, &info);
00185         chkxer_("DGETF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00186                 infoc_1.ok);
00187         infoc_1.infot = 2;
00188         dgetf2_(&c__0, &c_n1, a, &c__1, ip, &info);
00189         chkxer_("DGETF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00190                 infoc_1.ok);
00191         infoc_1.infot = 4;
00192         dgetf2_(&c__2, &c__1, a, &c__1, ip, &info);
00193         chkxer_("DGETF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00194                 infoc_1.ok);
00195 
00196 /*        DGETRI */
00197 
00198         s_copy(srnamc_1.srnamt, "DGETRI", (ftnlen)32, (ftnlen)6);
00199         infoc_1.infot = 1;
00200         dgetri_(&c_n1, a, &c__1, ip, w, &c__12, &info);
00201         chkxer_("DGETRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00202                 infoc_1.ok);
00203         infoc_1.infot = 3;
00204         dgetri_(&c__2, a, &c__1, ip, w, &c__12, &info);
00205         chkxer_("DGETRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00206                 infoc_1.ok);
00207 
00208 /*        DGETRS */
00209 
00210         s_copy(srnamc_1.srnamt, "DGETRS", (ftnlen)32, (ftnlen)6);
00211         infoc_1.infot = 1;
00212         dgetrs_("/", &c__0, &c__0, a, &c__1, ip, b, &c__1, &info);
00213         chkxer_("DGETRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00214                 infoc_1.ok);
00215         infoc_1.infot = 2;
00216         dgetrs_("N", &c_n1, &c__0, a, &c__1, ip, b, &c__1, &info);
00217         chkxer_("DGETRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00218                 infoc_1.ok);
00219         infoc_1.infot = 3;
00220         dgetrs_("N", &c__0, &c_n1, a, &c__1, ip, b, &c__1, &info);
00221         chkxer_("DGETRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00222                 infoc_1.ok);
00223         infoc_1.infot = 5;
00224         dgetrs_("N", &c__2, &c__1, a, &c__1, ip, b, &c__2, &info);
00225         chkxer_("DGETRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00226                 infoc_1.ok);
00227         infoc_1.infot = 8;
00228         dgetrs_("N", &c__2, &c__1, a, &c__2, ip, b, &c__1, &info);
00229         chkxer_("DGETRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00230                 infoc_1.ok);
00231 
00232 /*        DGERFS */
00233 
00234         s_copy(srnamc_1.srnamt, "DGERFS", (ftnlen)32, (ftnlen)6);
00235         infoc_1.infot = 1;
00236         dgerfs_("/", &c__0, &c__0, a, &c__1, af, &c__1, ip, b, &c__1, x, &
00237                 c__1, r1, r2, w, iw, &info);
00238         chkxer_("DGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00239                 infoc_1.ok);
00240         infoc_1.infot = 2;
00241         dgerfs_("N", &c_n1, &c__0, a, &c__1, af, &c__1, ip, b, &c__1, x, &
00242                 c__1, r1, r2, w, iw, &info);
00243         chkxer_("DGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00244                 infoc_1.ok);
00245         infoc_1.infot = 3;
00246         dgerfs_("N", &c__0, &c_n1, a, &c__1, af, &c__1, ip, b, &c__1, x, &
00247                 c__1, r1, r2, w, iw, &info);
00248         chkxer_("DGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00249                 infoc_1.ok);
00250         infoc_1.infot = 5;
00251         dgerfs_("N", &c__2, &c__1, a, &c__1, af, &c__2, ip, b, &c__2, x, &
00252                 c__2, r1, r2, w, iw, &info);
00253         chkxer_("DGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00254                 infoc_1.ok);
00255         infoc_1.infot = 7;
00256         dgerfs_("N", &c__2, &c__1, a, &c__2, af, &c__1, ip, b, &c__2, x, &
00257                 c__2, r1, r2, w, iw, &info);
00258         chkxer_("DGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00259                 infoc_1.ok);
00260         infoc_1.infot = 10;
00261         dgerfs_("N", &c__2, &c__1, a, &c__2, af, &c__2, ip, b, &c__1, x, &
00262                 c__2, r1, r2, w, iw, &info);
00263         chkxer_("DGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00264                 infoc_1.ok);
00265         infoc_1.infot = 12;
00266         dgerfs_("N", &c__2, &c__1, a, &c__2, af, &c__2, ip, b, &c__2, x, &
00267                 c__1, r1, r2, w, iw, &info);
00268         chkxer_("DGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00269                 infoc_1.ok);
00270 
00271 /*        DGECON */
00272 
00273         s_copy(srnamc_1.srnamt, "DGECON", (ftnlen)32, (ftnlen)6);
00274         infoc_1.infot = 1;
00275         dgecon_("/", &c__0, a, &c__1, &anrm, &rcond, w, iw, &info);
00276         chkxer_("DGECON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00277                 infoc_1.ok);
00278         infoc_1.infot = 2;
00279         dgecon_("1", &c_n1, a, &c__1, &anrm, &rcond, w, iw, &info);
00280         chkxer_("DGECON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00281                 infoc_1.ok);
00282         infoc_1.infot = 4;
00283         dgecon_("1", &c__2, a, &c__1, &anrm, &rcond, w, iw, &info);
00284         chkxer_("DGECON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00285                 infoc_1.ok);
00286 
00287 /*        DGEEQU */
00288 
00289         s_copy(srnamc_1.srnamt, "DGEEQU", (ftnlen)32, (ftnlen)6);
00290         infoc_1.infot = 1;
00291         dgeequ_(&c_n1, &c__0, a, &c__1, r1, r2, &rcond, &ccond, &anrm, &info);
00292         chkxer_("DGEEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00293                 infoc_1.ok);
00294         infoc_1.infot = 2;
00295         dgeequ_(&c__0, &c_n1, a, &c__1, r1, r2, &rcond, &ccond, &anrm, &info);
00296         chkxer_("DGEEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00297                 infoc_1.ok);
00298         infoc_1.infot = 4;
00299         dgeequ_(&c__2, &c__2, a, &c__1, r1, r2, &rcond, &ccond, &anrm, &info);
00300         chkxer_("DGEEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00301                 infoc_1.ok);
00302 
00303     } else if (lsamen_(&c__2, c2, "GB")) {
00304 
00305 /*        Test error exits of the routines that use the LU decomposition */
00306 /*        of a general band matrix. */
00307 
00308 /*        DGBTRF */
00309 
00310         s_copy(srnamc_1.srnamt, "DGBTRF", (ftnlen)32, (ftnlen)6);
00311         infoc_1.infot = 1;
00312         dgbtrf_(&c_n1, &c__0, &c__0, &c__0, a, &c__1, ip, &info);
00313         chkxer_("DGBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00314                 infoc_1.ok);
00315         infoc_1.infot = 2;
00316         dgbtrf_(&c__0, &c_n1, &c__0, &c__0, a, &c__1, ip, &info);
00317         chkxer_("DGBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00318                 infoc_1.ok);
00319         infoc_1.infot = 3;
00320         dgbtrf_(&c__1, &c__1, &c_n1, &c__0, a, &c__1, ip, &info);
00321         chkxer_("DGBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00322                 infoc_1.ok);
00323         infoc_1.infot = 4;
00324         dgbtrf_(&c__1, &c__1, &c__0, &c_n1, a, &c__1, ip, &info);
00325         chkxer_("DGBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00326                 infoc_1.ok);
00327         infoc_1.infot = 6;
00328         dgbtrf_(&c__2, &c__2, &c__1, &c__1, a, &c__3, ip, &info);
00329         chkxer_("DGBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00330                 infoc_1.ok);
00331 
00332 /*        DGBTF2 */
00333 
00334         s_copy(srnamc_1.srnamt, "DGBTF2", (ftnlen)32, (ftnlen)6);
00335         infoc_1.infot = 1;
00336         dgbtf2_(&c_n1, &c__0, &c__0, &c__0, a, &c__1, ip, &info);
00337         chkxer_("DGBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00338                 infoc_1.ok);
00339         infoc_1.infot = 2;
00340         dgbtf2_(&c__0, &c_n1, &c__0, &c__0, a, &c__1, ip, &info);
00341         chkxer_("DGBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00342                 infoc_1.ok);
00343         infoc_1.infot = 3;
00344         dgbtf2_(&c__1, &c__1, &c_n1, &c__0, a, &c__1, ip, &info);
00345         chkxer_("DGBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00346                 infoc_1.ok);
00347         infoc_1.infot = 4;
00348         dgbtf2_(&c__1, &c__1, &c__0, &c_n1, a, &c__1, ip, &info);
00349         chkxer_("DGBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00350                 infoc_1.ok);
00351         infoc_1.infot = 6;
00352         dgbtf2_(&c__2, &c__2, &c__1, &c__1, a, &c__3, ip, &info);
00353         chkxer_("DGBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00354                 infoc_1.ok);
00355 
00356 /*        DGBTRS */
00357 
00358         s_copy(srnamc_1.srnamt, "DGBTRS", (ftnlen)32, (ftnlen)6);
00359         infoc_1.infot = 1;
00360         dgbtrs_("/", &c__0, &c__0, &c__0, &c__1, a, &c__1, ip, b, &c__1, &
00361                 info);
00362         chkxer_("DGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00363                 infoc_1.ok);
00364         infoc_1.infot = 2;
00365         dgbtrs_("N", &c_n1, &c__0, &c__0, &c__1, a, &c__1, ip, b, &c__1, &
00366                 info);
00367         chkxer_("DGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00368                 infoc_1.ok);
00369         infoc_1.infot = 3;
00370         dgbtrs_("N", &c__1, &c_n1, &c__0, &c__1, a, &c__1, ip, b, &c__1, &
00371                 info);
00372         chkxer_("DGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00373                 infoc_1.ok);
00374         infoc_1.infot = 4;
00375         dgbtrs_("N", &c__1, &c__0, &c_n1, &c__1, a, &c__1, ip, b, &c__1, &
00376                 info);
00377         chkxer_("DGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00378                 infoc_1.ok);
00379         infoc_1.infot = 5;
00380         dgbtrs_("N", &c__1, &c__0, &c__0, &c_n1, a, &c__1, ip, b, &c__1, &
00381                 info);
00382         chkxer_("DGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00383                 infoc_1.ok);
00384         infoc_1.infot = 7;
00385         dgbtrs_("N", &c__2, &c__1, &c__1, &c__1, a, &c__3, ip, b, &c__2, &
00386                 info);
00387         chkxer_("DGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00388                 infoc_1.ok);
00389         infoc_1.infot = 10;
00390         dgbtrs_("N", &c__2, &c__0, &c__0, &c__1, a, &c__1, ip, b, &c__1, &
00391                 info);
00392         chkxer_("DGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00393                 infoc_1.ok);
00394 
00395 /*        DGBRFS */
00396 
00397         s_copy(srnamc_1.srnamt, "DGBRFS", (ftnlen)32, (ftnlen)6);
00398         infoc_1.infot = 1;
00399         dgbrfs_("/", &c__0, &c__0, &c__0, &c__0, a, &c__1, af, &c__1, ip, b, &
00400                 c__1, x, &c__1, r1, r2, w, iw, &info);
00401         chkxer_("DGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00402                 infoc_1.ok);
00403         infoc_1.infot = 2;
00404         dgbrfs_("N", &c_n1, &c__0, &c__0, &c__0, a, &c__1, af, &c__1, ip, b, &
00405                 c__1, x, &c__1, r1, r2, w, iw, &info);
00406         chkxer_("DGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00407                 infoc_1.ok);
00408         infoc_1.infot = 3;
00409         dgbrfs_("N", &c__1, &c_n1, &c__0, &c__0, a, &c__1, af, &c__1, ip, b, &
00410                 c__1, x, &c__1, r1, r2, w, iw, &info);
00411         chkxer_("DGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00412                 infoc_1.ok);
00413         infoc_1.infot = 4;
00414         dgbrfs_("N", &c__1, &c__0, &c_n1, &c__0, a, &c__1, af, &c__1, ip, b, &
00415                 c__1, x, &c__1, r1, r2, w, iw, &info);
00416         chkxer_("DGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00417                 infoc_1.ok);
00418         infoc_1.infot = 5;
00419         dgbrfs_("N", &c__1, &c__0, &c__0, &c_n1, a, &c__1, af, &c__1, ip, b, &
00420                 c__1, x, &c__1, r1, r2, w, iw, &info);
00421         chkxer_("DGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00422                 infoc_1.ok);
00423         infoc_1.infot = 7;
00424         dgbrfs_("N", &c__2, &c__1, &c__1, &c__1, a, &c__2, af, &c__4, ip, b, &
00425                 c__2, x, &c__2, r1, r2, w, iw, &info);
00426         chkxer_("DGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00427                 infoc_1.ok);
00428         infoc_1.infot = 9;
00429         dgbrfs_("N", &c__2, &c__1, &c__1, &c__1, a, &c__3, af, &c__3, ip, b, &
00430                 c__2, x, &c__2, r1, r2, w, iw, &info);
00431         chkxer_("DGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00432                 infoc_1.ok);
00433         infoc_1.infot = 12;
00434         dgbrfs_("N", &c__2, &c__0, &c__0, &c__1, a, &c__1, af, &c__1, ip, b, &
00435                 c__1, x, &c__2, r1, r2, w, iw, &info);
00436         chkxer_("DGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00437                 infoc_1.ok);
00438         infoc_1.infot = 14;
00439         dgbrfs_("N", &c__2, &c__0, &c__0, &c__1, a, &c__1, af, &c__1, ip, b, &
00440                 c__2, x, &c__1, r1, r2, w, iw, &info);
00441         chkxer_("DGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00442                 infoc_1.ok);
00443 
00444 /*        DGBCON */
00445 
00446         s_copy(srnamc_1.srnamt, "DGBCON", (ftnlen)32, (ftnlen)6);
00447         infoc_1.infot = 1;
00448         dgbcon_("/", &c__0, &c__0, &c__0, a, &c__1, ip, &anrm, &rcond, w, iw, 
00449                 &info);
00450         chkxer_("DGBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00451                 infoc_1.ok);
00452         infoc_1.infot = 2;
00453         dgbcon_("1", &c_n1, &c__0, &c__0, a, &c__1, ip, &anrm, &rcond, w, iw, 
00454                 &info);
00455         chkxer_("DGBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00456                 infoc_1.ok);
00457         infoc_1.infot = 3;
00458         dgbcon_("1", &c__1, &c_n1, &c__0, a, &c__1, ip, &anrm, &rcond, w, iw, 
00459                 &info);
00460         chkxer_("DGBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00461                 infoc_1.ok);
00462         infoc_1.infot = 4;
00463         dgbcon_("1", &c__1, &c__0, &c_n1, a, &c__1, ip, &anrm, &rcond, w, iw, 
00464                 &info);
00465         chkxer_("DGBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00466                 infoc_1.ok);
00467         infoc_1.infot = 6;
00468         dgbcon_("1", &c__2, &c__1, &c__1, a, &c__3, ip, &anrm, &rcond, w, iw, 
00469                 &info);
00470         chkxer_("DGBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00471                 infoc_1.ok);
00472 
00473 /*        DGBEQU */
00474 
00475         s_copy(srnamc_1.srnamt, "DGBEQU", (ftnlen)32, (ftnlen)6);
00476         infoc_1.infot = 1;
00477         dgbequ_(&c_n1, &c__0, &c__0, &c__0, a, &c__1, r1, r2, &rcond, &ccond, 
00478                 &anrm, &info);
00479         chkxer_("DGBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00480                 infoc_1.ok);
00481         infoc_1.infot = 2;
00482         dgbequ_(&c__0, &c_n1, &c__0, &c__0, a, &c__1, r1, r2, &rcond, &ccond, 
00483                 &anrm, &info);
00484         chkxer_("DGBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00485                 infoc_1.ok);
00486         infoc_1.infot = 3;
00487         dgbequ_(&c__1, &c__1, &c_n1, &c__0, a, &c__1, r1, r2, &rcond, &ccond, 
00488                 &anrm, &info);
00489         chkxer_("DGBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00490                 infoc_1.ok);
00491         infoc_1.infot = 4;
00492         dgbequ_(&c__1, &c__1, &c__0, &c_n1, a, &c__1, r1, r2, &rcond, &ccond, 
00493                 &anrm, &info);
00494         chkxer_("DGBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00495                 infoc_1.ok);
00496         infoc_1.infot = 6;
00497         dgbequ_(&c__2, &c__2, &c__1, &c__1, a, &c__2, r1, r2, &rcond, &ccond, 
00498                 &anrm, &info);
00499         chkxer_("DGBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00500                 infoc_1.ok);
00501     }
00502 
00503 /*     Print a summary line. */
00504 
00505     alaesm_(path, &infoc_1.ok, &infoc_1.nout);
00506 
00507     return 0;
00508 
00509 /*     End of DERRGE */
00510 
00511 } /* derrge_ */


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