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


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