00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013 #include "f2c.h"
00014 #include "blaswrap.h"
00015
00016
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
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 int serrge_(char *path, integer *nunit)
00043 {
00044
00045 integer s_wsle(cilist *), e_wsle(void);
00046 int s_copy(char *, char *, ftnlen, ftnlen);
00047
00048
00049 real a[16] , 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] ;
00054 char eq[1];
00055 integer ip[4], iw[4];
00056 real err_bnds_c__[12] ;
00057 integer n_err_bnds__;
00058 real err_bnds_n__[12] , berr;
00059 integer info;
00060 real anrm, ccond, rcond;
00061 extern 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 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 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
00105 static cilist io___1 = { 0, 0, 0, 0, 0 };
00106
00107
00108
00109
00110
00111
00112
00113
00114
00115
00116
00117
00118
00119
00120
00121
00122
00123
00124
00125
00126
00127
00128
00129
00130
00131
00132
00133
00134
00135
00136
00137
00138
00139
00140
00141
00142
00143
00144
00145
00146
00147
00148
00149
00150
00151
00152
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
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
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
00178 }
00179 infoc_1.ok = TRUE_;
00180
00181 if (lsamen_(&c__2, c2, "GE")) {
00182
00183
00184
00185
00186
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
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
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
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
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
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
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
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
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
00404
00405
00406
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
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
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
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
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
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
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
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
00702
00703 alaesm_(path, &infoc_1.ok, &infoc_1.nout);
00704
00705 return 0;
00706
00707
00708
00709 }