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


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