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__3 = 3;
00038 static integer c__4 = 4;
00039 static integer c__5 = 5;
00040
00041 int cerrge_(char *path, integer *nunit)
00042 {
00043
00044 integer i__1;
00045 real r__1, r__2;
00046 complex q__1;
00047
00048
00049 integer s_wsle(cilist *), e_wsle(void);
00050 int s_copy(char *, char *, ftnlen, ftnlen);
00051
00052
00053 complex a[16] , b[4];
00054 integer i__, j;
00055 real r__[4];
00056 complex w[8], x[4];
00057 char c2[2];
00058 real r1[4], r2[4];
00059 complex af[16] ;
00060 char eq[1];
00061 real cs[4];
00062 integer ip[4];
00063 real rs[4];
00064 complex err_bnds_c__[12] ;
00065 integer n_err_bnds__;
00066 complex err_bnds_n__[12] ;
00067 real berr;
00068 integer info;
00069 real anrm, ccond, rcond;
00070 extern int cgbtf2_(integer *, integer *, integer *,
00071 integer *, complex *, integer *, integer *, integer *), cgetf2_(
00072 integer *, integer *, complex *, integer *, integer *, integer *),
00073 cgbcon_(char *, integer *, integer *, integer *, complex *,
00074 integer *, integer *, real *, real *, complex *, real *, integer *
00075 ), cgecon_(char *, integer *, complex *, integer *, real *
00076 , real *, complex *, real *, integer *), alaesm_(char *,
00077 logical *, integer *), cgbequ_(integer *, integer *,
00078 integer *, integer *, complex *, integer *, real *, real *, real *
00079 , real *, real *, integer *), cgbrfs_(char *, integer *, integer *
00080 , integer *, integer *, complex *, integer *, complex *, integer *
00081 , integer *, complex *, integer *, complex *, integer *, real *,
00082 real *, complex *, real *, integer *), cgbtrf_(integer *,
00083 integer *, integer *, integer *, complex *, integer *, integer *,
00084 integer *), cgeequ_(integer *, integer *, complex *, integer *,
00085 real *, real *, real *, real *, real *, integer *), cgerfs_(char *
00086 , integer *, integer *, complex *, integer *, complex *, integer *
00087 , integer *, complex *, integer *, complex *, integer *, real *,
00088 real *, complex *, real *, integer *), cgetrf_(integer *,
00089 integer *, complex *, integer *, integer *, integer *), cgetri_(
00090 integer *, complex *, integer *, integer *, complex *, integer *,
00091 integer *);
00092 extern logical lsamen_(integer *, char *, char *);
00093 complex params[1];
00094 extern int cgbtrs_(char *, integer *, integer *, integer
00095 *, integer *, complex *, integer *, integer *, complex *, integer
00096 *, integer *), chkxer_(char *, integer *, integer *,
00097 logical *, logical *), cgetrs_(char *, integer *, integer
00098 *, complex *, integer *, integer *, complex *, integer *, integer
00099 *), cgbequb_(integer *, integer *, integer *, integer *,
00100 complex *, integer *, real *, real *, real *, real *, real *,
00101 integer *), cgeequb_(integer *, integer *, complex *, integer *,
00102 real *, real *, real *, real *, real *, integer *), cgbrfsx_(char
00103 *, char *, integer *, integer *, integer *, integer *, complex *,
00104 integer *, complex *, integer *, integer *, real *, real *,
00105 complex *, integer *, complex *, integer *, real *, real *,
00106 integer *, complex *, complex *, integer *, complex *, complex *,
00107 real *, integer *);
00108 integer nparams;
00109 extern int cgerfsx_(char *, char *, integer *, integer *,
00110 complex *, integer *, complex *, integer *, integer *, real *,
00111 real *, complex *, integer *, complex *, integer *, real *, real *
00112 , integer *, complex *, complex *, integer *, complex *, complex *
00113 , real *, integer *);
00114
00115
00116 static cilist io___1 = { 0, 0, 0, 0, 0 };
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
00155
00156
00157
00158
00159
00160
00161
00162
00163
00164
00165 infoc_1.nout = *nunit;
00166 io___1.ciunit = infoc_1.nout;
00167 s_wsle(&io___1);
00168 e_wsle();
00169 s_copy(c2, path + 1, (ftnlen)2, (ftnlen)2);
00170
00171
00172
00173 for (j = 1; j <= 4; ++j) {
00174 for (i__ = 1; i__ <= 4; ++i__) {
00175 i__1 = i__ + (j << 2) - 5;
00176 r__1 = 1.f / (real) (i__ + j);
00177 r__2 = -1.f / (real) (i__ + j);
00178 q__1.r = r__1, q__1.i = r__2;
00179 a[i__1].r = q__1.r, a[i__1].i = q__1.i;
00180 i__1 = i__ + (j << 2) - 5;
00181 r__1 = 1.f / (real) (i__ + j);
00182 r__2 = -1.f / (real) (i__ + j);
00183 q__1.r = r__1, q__1.i = r__2;
00184 af[i__1].r = q__1.r, af[i__1].i = q__1.i;
00185
00186 }
00187 i__1 = j - 1;
00188 b[i__1].r = 0.f, b[i__1].i = 0.f;
00189 r1[j - 1] = 0.f;
00190 r2[j - 1] = 0.f;
00191 i__1 = j - 1;
00192 w[i__1].r = 0.f, w[i__1].i = 0.f;
00193 i__1 = j - 1;
00194 x[i__1].r = 0.f, x[i__1].i = 0.f;
00195 cs[j - 1] = 0.f;
00196 rs[j - 1] = 0.f;
00197 ip[j - 1] = j;
00198
00199 }
00200 infoc_1.ok = TRUE_;
00201
00202
00203
00204
00205 if (lsamen_(&c__2, c2, "GE")) {
00206
00207
00208
00209 s_copy(srnamc_1.srnamt, "CGETRF", (ftnlen)32, (ftnlen)6);
00210 infoc_1.infot = 1;
00211 cgetrf_(&c_n1, &c__0, a, &c__1, ip, &info);
00212 chkxer_("CGETRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00213 infoc_1.ok);
00214 infoc_1.infot = 2;
00215 cgetrf_(&c__0, &c_n1, a, &c__1, ip, &info);
00216 chkxer_("CGETRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00217 infoc_1.ok);
00218 infoc_1.infot = 4;
00219 cgetrf_(&c__2, &c__1, a, &c__1, ip, &info);
00220 chkxer_("CGETRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00221 infoc_1.ok);
00222
00223
00224
00225 s_copy(srnamc_1.srnamt, "CGETF2", (ftnlen)32, (ftnlen)6);
00226 infoc_1.infot = 1;
00227 cgetf2_(&c_n1, &c__0, a, &c__1, ip, &info);
00228 chkxer_("CGETF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00229 infoc_1.ok);
00230 infoc_1.infot = 2;
00231 cgetf2_(&c__0, &c_n1, a, &c__1, ip, &info);
00232 chkxer_("CGETF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00233 infoc_1.ok);
00234 infoc_1.infot = 4;
00235 cgetf2_(&c__2, &c__1, a, &c__1, ip, &info);
00236 chkxer_("CGETF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00237 infoc_1.ok);
00238
00239
00240
00241 s_copy(srnamc_1.srnamt, "CGETRI", (ftnlen)32, (ftnlen)6);
00242 infoc_1.infot = 1;
00243 cgetri_(&c_n1, a, &c__1, ip, w, &c__1, &info);
00244 chkxer_("CGETRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00245 infoc_1.ok);
00246 infoc_1.infot = 3;
00247 cgetri_(&c__2, a, &c__1, ip, w, &c__2, &info);
00248 chkxer_("CGETRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00249 infoc_1.ok);
00250 infoc_1.infot = 6;
00251 cgetri_(&c__2, a, &c__2, ip, w, &c__1, &info);
00252 chkxer_("CGETRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00253 infoc_1.ok);
00254
00255
00256
00257 s_copy(srnamc_1.srnamt, "CGETRS", (ftnlen)32, (ftnlen)6);
00258 infoc_1.infot = 1;
00259 cgetrs_("/", &c__0, &c__0, a, &c__1, ip, b, &c__1, &info);
00260 chkxer_("CGETRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00261 infoc_1.ok);
00262 infoc_1.infot = 2;
00263 cgetrs_("N", &c_n1, &c__0, a, &c__1, ip, b, &c__1, &info);
00264 chkxer_("CGETRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00265 infoc_1.ok);
00266 infoc_1.infot = 3;
00267 cgetrs_("N", &c__0, &c_n1, a, &c__1, ip, b, &c__1, &info);
00268 chkxer_("CGETRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00269 infoc_1.ok);
00270 infoc_1.infot = 5;
00271 cgetrs_("N", &c__2, &c__1, a, &c__1, ip, b, &c__2, &info);
00272 chkxer_("CGETRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00273 infoc_1.ok);
00274 infoc_1.infot = 8;
00275 cgetrs_("N", &c__2, &c__1, a, &c__2, ip, b, &c__1, &info);
00276 chkxer_("CGETRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00277 infoc_1.ok);
00278
00279
00280
00281 s_copy(srnamc_1.srnamt, "CGERFS", (ftnlen)32, (ftnlen)6);
00282 infoc_1.infot = 1;
00283 cgerfs_("/", &c__0, &c__0, a, &c__1, af, &c__1, ip, b, &c__1, x, &
00284 c__1, r1, r2, w, r__, &info);
00285 chkxer_("CGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00286 infoc_1.ok);
00287 infoc_1.infot = 2;
00288 cgerfs_("N", &c_n1, &c__0, a, &c__1, af, &c__1, ip, b, &c__1, x, &
00289 c__1, r1, r2, w, r__, &info);
00290 chkxer_("CGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00291 infoc_1.ok);
00292 infoc_1.infot = 3;
00293 cgerfs_("N", &c__0, &c_n1, a, &c__1, af, &c__1, ip, b, &c__1, x, &
00294 c__1, r1, r2, w, r__, &info);
00295 chkxer_("CGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00296 infoc_1.ok);
00297 infoc_1.infot = 5;
00298 cgerfs_("N", &c__2, &c__1, a, &c__1, af, &c__2, ip, b, &c__2, x, &
00299 c__2, r1, r2, w, r__, &info);
00300 chkxer_("CGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00301 infoc_1.ok);
00302 infoc_1.infot = 7;
00303 cgerfs_("N", &c__2, &c__1, a, &c__2, af, &c__1, ip, b, &c__2, x, &
00304 c__2, r1, r2, w, r__, &info);
00305 chkxer_("CGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00306 infoc_1.ok);
00307 infoc_1.infot = 10;
00308 cgerfs_("N", &c__2, &c__1, a, &c__2, af, &c__2, ip, b, &c__1, x, &
00309 c__2, r1, r2, w, r__, &info);
00310 chkxer_("CGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00311 infoc_1.ok);
00312 infoc_1.infot = 12;
00313 cgerfs_("N", &c__2, &c__1, a, &c__2, af, &c__2, ip, b, &c__2, x, &
00314 c__1, r1, r2, w, r__, &info);
00315 chkxer_("CGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00316 infoc_1.ok);
00317
00318
00319
00320 n_err_bnds__ = 3;
00321 nparams = 0;
00322 s_copy(srnamc_1.srnamt, "CGERFSX", (ftnlen)32, (ftnlen)7);
00323 infoc_1.infot = 1;
00324 cgerfsx_("/", eq, &c__0, &c__0, a, &c__1, af, &c__1, ip, rs, cs, b, &
00325 c__1, x, &c__1, &rcond, &berr, &n_err_bnds__, err_bnds_n__,
00326 err_bnds_c__, &nparams, params, w, r__, &info);
00327 chkxer_("CGERFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00328 infoc_1.ok);
00329 infoc_1.infot = 2;
00330 *(unsigned char *)eq = '/';
00331 cgerfsx_("N", eq, &c__2, &c__1, a, &c__1, af, &c__2, ip, rs, cs, b, &
00332 c__2, x, &c__2, &rcond, &berr, &n_err_bnds__, err_bnds_n__,
00333 err_bnds_c__, &nparams, params, w, r__, &info);
00334 chkxer_("CGERFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00335 infoc_1.ok);
00336 infoc_1.infot = 3;
00337 *(unsigned char *)eq = 'R';
00338 cgerfsx_("N", eq, &c_n1, &c__0, a, &c__1, af, &c__1, ip, rs, cs, b, &
00339 c__1, x, &c__1, &rcond, &berr, &n_err_bnds__, err_bnds_n__,
00340 err_bnds_c__, &nparams, params, w, r__, &info);
00341 chkxer_("CGERFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00342 infoc_1.ok);
00343 infoc_1.infot = 4;
00344 cgerfsx_("N", eq, &c__0, &c_n1, a, &c__1, af, &c__1, ip, rs, cs, b, &
00345 c__1, x, &c__1, &rcond, &berr, &n_err_bnds__, err_bnds_n__,
00346 err_bnds_c__, &nparams, params, w, r__, &info);
00347 chkxer_("CGERFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00348 infoc_1.ok);
00349 infoc_1.infot = 6;
00350 cgerfsx_("N", eq, &c__2, &c__1, a, &c__1, af, &c__2, ip, rs, cs, b, &
00351 c__2, x, &c__2, &rcond, &berr, &n_err_bnds__, err_bnds_n__,
00352 err_bnds_c__, &nparams, params, w, r__, &info);
00353 chkxer_("CGERFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00354 infoc_1.ok);
00355 infoc_1.infot = 8;
00356 cgerfsx_("N", eq, &c__2, &c__1, a, &c__2, af, &c__1, ip, rs, cs, b, &
00357 c__2, x, &c__2, &rcond, &berr, &n_err_bnds__, err_bnds_n__,
00358 err_bnds_c__, &nparams, params, w, r__, &info);
00359 chkxer_("CGERFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00360 infoc_1.ok);
00361 infoc_1.infot = 13;
00362 *(unsigned char *)eq = 'C';
00363 cgerfsx_("N", eq, &c__2, &c__1, a, &c__2, af, &c__2, ip, rs, cs, b, &
00364 c__1, x, &c__2, &rcond, &berr, &n_err_bnds__, err_bnds_n__,
00365 err_bnds_c__, &nparams, params, w, r__, &info);
00366 chkxer_("CGERFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00367 infoc_1.ok);
00368 infoc_1.infot = 15;
00369 cgerfsx_("N", eq, &c__2, &c__1, a, &c__2, af, &c__2, ip, rs, cs, b, &
00370 c__2, x, &c__1, &rcond, &berr, &n_err_bnds__, err_bnds_n__,
00371 err_bnds_c__, &nparams, params, w, r__, &info);
00372 chkxer_("CGERFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00373 infoc_1.ok);
00374
00375
00376
00377 s_copy(srnamc_1.srnamt, "CGECON", (ftnlen)32, (ftnlen)6);
00378 infoc_1.infot = 1;
00379 cgecon_("/", &c__0, a, &c__1, &anrm, &rcond, w, r__, &info)
00380 ;
00381 chkxer_("CGECON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00382 infoc_1.ok);
00383 infoc_1.infot = 2;
00384 cgecon_("1", &c_n1, a, &c__1, &anrm, &rcond, w, r__, &info)
00385 ;
00386 chkxer_("CGECON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00387 infoc_1.ok);
00388 infoc_1.infot = 4;
00389 cgecon_("1", &c__2, a, &c__1, &anrm, &rcond, w, r__, &info)
00390 ;
00391 chkxer_("CGECON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00392 infoc_1.ok);
00393
00394
00395
00396 s_copy(srnamc_1.srnamt, "CGEEQU", (ftnlen)32, (ftnlen)6);
00397 infoc_1.infot = 1;
00398 cgeequ_(&c_n1, &c__0, a, &c__1, r1, r2, &rcond, &ccond, &anrm, &info);
00399 chkxer_("CGEEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00400 infoc_1.ok);
00401 infoc_1.infot = 2;
00402 cgeequ_(&c__0, &c_n1, a, &c__1, r1, r2, &rcond, &ccond, &anrm, &info);
00403 chkxer_("CGEEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00404 infoc_1.ok);
00405 infoc_1.infot = 4;
00406 cgeequ_(&c__2, &c__2, a, &c__1, r1, r2, &rcond, &ccond, &anrm, &info);
00407 chkxer_("CGEEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00408 infoc_1.ok);
00409
00410
00411
00412 s_copy(srnamc_1.srnamt, "CGEEQUB", (ftnlen)32, (ftnlen)7);
00413 infoc_1.infot = 1;
00414 cgeequb_(&c_n1, &c__0, a, &c__1, r1, r2, &rcond, &ccond, &anrm, &info)
00415 ;
00416 chkxer_("CGEEQUB", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00417 infoc_1.ok);
00418 infoc_1.infot = 2;
00419 cgeequb_(&c__0, &c_n1, a, &c__1, r1, r2, &rcond, &ccond, &anrm, &info)
00420 ;
00421 chkxer_("CGEEQUB", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00422 infoc_1.ok);
00423 infoc_1.infot = 4;
00424 cgeequb_(&c__2, &c__2, a, &c__1, r1, r2, &rcond, &ccond, &anrm, &info)
00425 ;
00426 chkxer_("CGEEQUB", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00427 infoc_1.ok);
00428
00429
00430
00431
00432 } else if (lsamen_(&c__2, c2, "GB")) {
00433
00434
00435
00436 s_copy(srnamc_1.srnamt, "CGBTRF", (ftnlen)32, (ftnlen)6);
00437 infoc_1.infot = 1;
00438 cgbtrf_(&c_n1, &c__0, &c__0, &c__0, a, &c__1, ip, &info);
00439 chkxer_("CGBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00440 infoc_1.ok);
00441 infoc_1.infot = 2;
00442 cgbtrf_(&c__0, &c_n1, &c__0, &c__0, a, &c__1, ip, &info);
00443 chkxer_("CGBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00444 infoc_1.ok);
00445 infoc_1.infot = 3;
00446 cgbtrf_(&c__1, &c__1, &c_n1, &c__0, a, &c__1, ip, &info);
00447 chkxer_("CGBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00448 infoc_1.ok);
00449 infoc_1.infot = 4;
00450 cgbtrf_(&c__1, &c__1, &c__0, &c_n1, a, &c__1, ip, &info);
00451 chkxer_("CGBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00452 infoc_1.ok);
00453 infoc_1.infot = 6;
00454 cgbtrf_(&c__2, &c__2, &c__1, &c__1, a, &c__3, ip, &info);
00455 chkxer_("CGBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00456 infoc_1.ok);
00457
00458
00459
00460 s_copy(srnamc_1.srnamt, "CGBTF2", (ftnlen)32, (ftnlen)6);
00461 infoc_1.infot = 1;
00462 cgbtf2_(&c_n1, &c__0, &c__0, &c__0, a, &c__1, ip, &info);
00463 chkxer_("CGBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00464 infoc_1.ok);
00465 infoc_1.infot = 2;
00466 cgbtf2_(&c__0, &c_n1, &c__0, &c__0, a, &c__1, ip, &info);
00467 chkxer_("CGBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00468 infoc_1.ok);
00469 infoc_1.infot = 3;
00470 cgbtf2_(&c__1, &c__1, &c_n1, &c__0, a, &c__1, ip, &info);
00471 chkxer_("CGBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00472 infoc_1.ok);
00473 infoc_1.infot = 4;
00474 cgbtf2_(&c__1, &c__1, &c__0, &c_n1, a, &c__1, ip, &info);
00475 chkxer_("CGBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00476 infoc_1.ok);
00477 infoc_1.infot = 6;
00478 cgbtf2_(&c__2, &c__2, &c__1, &c__1, a, &c__3, ip, &info);
00479 chkxer_("CGBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00480 infoc_1.ok);
00481
00482
00483
00484 s_copy(srnamc_1.srnamt, "CGBTRS", (ftnlen)32, (ftnlen)6);
00485 infoc_1.infot = 1;
00486 cgbtrs_("/", &c__0, &c__0, &c__0, &c__1, a, &c__1, ip, b, &c__1, &
00487 info);
00488 chkxer_("CGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00489 infoc_1.ok);
00490 infoc_1.infot = 2;
00491 cgbtrs_("N", &c_n1, &c__0, &c__0, &c__1, a, &c__1, ip, b, &c__1, &
00492 info);
00493 chkxer_("CGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00494 infoc_1.ok);
00495 infoc_1.infot = 3;
00496 cgbtrs_("N", &c__1, &c_n1, &c__0, &c__1, a, &c__1, ip, b, &c__1, &
00497 info);
00498 chkxer_("CGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00499 infoc_1.ok);
00500 infoc_1.infot = 4;
00501 cgbtrs_("N", &c__1, &c__0, &c_n1, &c__1, a, &c__1, ip, b, &c__1, &
00502 info);
00503 chkxer_("CGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00504 infoc_1.ok);
00505 infoc_1.infot = 5;
00506 cgbtrs_("N", &c__1, &c__0, &c__0, &c_n1, a, &c__1, ip, b, &c__1, &
00507 info);
00508 chkxer_("CGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00509 infoc_1.ok);
00510 infoc_1.infot = 7;
00511 cgbtrs_("N", &c__2, &c__1, &c__1, &c__1, a, &c__3, ip, b, &c__2, &
00512 info);
00513 chkxer_("CGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00514 infoc_1.ok);
00515 infoc_1.infot = 10;
00516 cgbtrs_("N", &c__2, &c__0, &c__0, &c__1, a, &c__1, ip, b, &c__1, &
00517 info);
00518 chkxer_("CGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00519 infoc_1.ok);
00520
00521
00522
00523 s_copy(srnamc_1.srnamt, "CGBRFS", (ftnlen)32, (ftnlen)6);
00524 infoc_1.infot = 1;
00525 cgbrfs_("/", &c__0, &c__0, &c__0, &c__0, a, &c__1, af, &c__1, ip, b, &
00526 c__1, x, &c__1, r1, r2, w, r__, &info);
00527 chkxer_("CGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00528 infoc_1.ok);
00529 infoc_1.infot = 2;
00530 cgbrfs_("N", &c_n1, &c__0, &c__0, &c__0, a, &c__1, af, &c__1, ip, b, &
00531 c__1, x, &c__1, r1, r2, w, r__, &info);
00532 chkxer_("CGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00533 infoc_1.ok);
00534 infoc_1.infot = 3;
00535 cgbrfs_("N", &c__1, &c_n1, &c__0, &c__0, a, &c__1, af, &c__1, ip, b, &
00536 c__1, x, &c__1, r1, r2, w, r__, &info);
00537 chkxer_("CGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00538 infoc_1.ok);
00539 infoc_1.infot = 4;
00540 cgbrfs_("N", &c__1, &c__0, &c_n1, &c__0, a, &c__1, af, &c__1, ip, b, &
00541 c__1, x, &c__1, r1, r2, w, r__, &info);
00542 chkxer_("CGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00543 infoc_1.ok);
00544 infoc_1.infot = 5;
00545 cgbrfs_("N", &c__1, &c__0, &c__0, &c_n1, a, &c__1, af, &c__1, ip, b, &
00546 c__1, x, &c__1, r1, r2, w, r__, &info);
00547 chkxer_("CGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00548 infoc_1.ok);
00549 infoc_1.infot = 7;
00550 cgbrfs_("N", &c__2, &c__1, &c__1, &c__1, a, &c__2, af, &c__4, ip, b, &
00551 c__2, x, &c__2, r1, r2, w, r__, &info);
00552 chkxer_("CGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00553 infoc_1.ok);
00554 infoc_1.infot = 9;
00555 cgbrfs_("N", &c__2, &c__1, &c__1, &c__1, a, &c__3, af, &c__3, ip, b, &
00556 c__2, x, &c__2, r1, r2, w, r__, &info);
00557 chkxer_("CGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00558 infoc_1.ok);
00559 infoc_1.infot = 12;
00560 cgbrfs_("N", &c__2, &c__0, &c__0, &c__1, a, &c__1, af, &c__1, ip, b, &
00561 c__1, x, &c__2, r1, r2, w, r__, &info);
00562 chkxer_("CGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00563 infoc_1.ok);
00564 infoc_1.infot = 14;
00565 cgbrfs_("N", &c__2, &c__0, &c__0, &c__1, a, &c__1, af, &c__1, ip, b, &
00566 c__2, x, &c__1, r1, r2, w, r__, &info);
00567 chkxer_("CGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00568 infoc_1.ok);
00569
00570
00571
00572 n_err_bnds__ = 3;
00573 nparams = 0;
00574 s_copy(srnamc_1.srnamt, "CGBRFSX", (ftnlen)32, (ftnlen)7);
00575 infoc_1.infot = 1;
00576 cgbrfsx_("/", eq, &c__0, &c__0, &c__0, &c__0, a, &c__1, af, &c__1, ip,
00577 rs, cs, b, &c__1, x, &c__1, &rcond, &berr, &n_err_bnds__,
00578 err_bnds_n__, err_bnds_c__, &nparams, params, w, r__, &info);
00579 chkxer_("CGBRFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00580 infoc_1.ok);
00581 infoc_1.infot = 2;
00582 *(unsigned char *)eq = '/';
00583 cgbrfsx_("N", eq, &c__2, &c__1, &c__1, &c__1, a, &c__1, af, &c__2, ip,
00584 rs, cs, b, &c__2, x, &c__2, &rcond, &berr, &n_err_bnds__,
00585 err_bnds_n__, err_bnds_c__, &nparams, params, w, r__, &info);
00586 chkxer_("CGBRFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00587 infoc_1.ok);
00588 infoc_1.infot = 3;
00589 *(unsigned char *)eq = 'R';
00590 cgbrfsx_("N", eq, &c_n1, &c__1, &c__1, &c__0, a, &c__1, af, &c__1, ip,
00591 rs, cs, b, &c__1, x, &c__1, &rcond, &berr, &n_err_bnds__,
00592 err_bnds_n__, err_bnds_c__, &nparams, params, w, r__, &info);
00593 chkxer_("CGBRFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00594 infoc_1.ok);
00595 infoc_1.infot = 4;
00596 *(unsigned char *)eq = 'R';
00597 cgbrfsx_("N", eq, &c__2, &c_n1, &c__1, &c__1, a, &c__3, af, &c__4, ip,
00598 rs, cs, b, &c__1, x, &c__1, &rcond, &berr, &n_err_bnds__,
00599 err_bnds_n__, err_bnds_c__, &nparams, params, w, r__, &info);
00600 chkxer_("CGBRFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00601 infoc_1.ok);
00602 infoc_1.infot = 5;
00603 *(unsigned char *)eq = 'R';
00604 cgbrfsx_("N", eq, &c__2, &c__1, &c_n1, &c__1, a, &c__3, af, &c__4, ip,
00605 rs, cs, b, &c__1, x, &c__1, &rcond, &berr, &n_err_bnds__,
00606 err_bnds_n__, err_bnds_c__, &nparams, params, w, r__, &info);
00607 chkxer_("CGBRFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00608 infoc_1.ok);
00609 infoc_1.infot = 6;
00610 cgbrfsx_("N", eq, &c__0, &c__0, &c__0, &c_n1, a, &c__1, af, &c__1, ip,
00611 rs, cs, b, &c__1, x, &c__1, &rcond, &berr, &n_err_bnds__,
00612 err_bnds_n__, err_bnds_c__, &nparams, params, w, r__, &info);
00613 chkxer_("CGBRFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00614 infoc_1.ok);
00615 infoc_1.infot = 8;
00616 cgbrfsx_("N", eq, &c__2, &c__1, &c__1, &c__1, a, &c__1, af, &c__2, ip,
00617 rs, cs, b, &c__2, x, &c__2, &rcond, &berr, &n_err_bnds__,
00618 err_bnds_n__, err_bnds_c__, &nparams, params, w, r__, &info);
00619 chkxer_("CGBRFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00620 infoc_1.ok);
00621 infoc_1.infot = 10;
00622 cgbrfsx_("N", eq, &c__2, &c__1, &c__1, &c__1, a, &c__3, af, &c__3, ip,
00623 rs, cs, b, &c__2, x, &c__2, &rcond, &berr, &n_err_bnds__,
00624 err_bnds_n__, err_bnds_c__, &nparams, params, w, r__, &info);
00625 chkxer_("CGBRFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00626 infoc_1.ok);
00627 infoc_1.infot = 13;
00628 *(unsigned char *)eq = 'C';
00629 cgbrfsx_("N", eq, &c__2, &c__1, &c__1, &c__1, a, &c__3, af, &c__5, ip,
00630 rs, cs, b, &c__1, x, &c__2, &rcond, &berr, &n_err_bnds__,
00631 err_bnds_n__, err_bnds_c__, &nparams, params, w, r__, &info);
00632 chkxer_("CGBRFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00633 infoc_1.ok);
00634 infoc_1.infot = 15;
00635 cgbrfsx_("N", eq, &c__2, &c__1, &c__1, &c__1, a, &c__3, af, &c__5, ip,
00636 rs, cs, b, &c__2, x, &c__1, &rcond, &berr, &n_err_bnds__,
00637 err_bnds_n__, err_bnds_c__, &nparams, params, w, r__, &info);
00638 chkxer_("CGBRFSX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00639 infoc_1.ok);
00640
00641
00642
00643 s_copy(srnamc_1.srnamt, "CGBCON", (ftnlen)32, (ftnlen)6);
00644 infoc_1.infot = 1;
00645 cgbcon_("/", &c__0, &c__0, &c__0, a, &c__1, ip, &anrm, &rcond, w, r__,
00646 &info);
00647 chkxer_("CGBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00648 infoc_1.ok);
00649 infoc_1.infot = 2;
00650 cgbcon_("1", &c_n1, &c__0, &c__0, a, &c__1, ip, &anrm, &rcond, w, r__,
00651 &info);
00652 chkxer_("CGBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00653 infoc_1.ok);
00654 infoc_1.infot = 3;
00655 cgbcon_("1", &c__1, &c_n1, &c__0, a, &c__1, ip, &anrm, &rcond, w, r__,
00656 &info);
00657 chkxer_("CGBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00658 infoc_1.ok);
00659 infoc_1.infot = 4;
00660 cgbcon_("1", &c__1, &c__0, &c_n1, a, &c__1, ip, &anrm, &rcond, w, r__,
00661 &info);
00662 chkxer_("CGBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00663 infoc_1.ok);
00664 infoc_1.infot = 6;
00665 cgbcon_("1", &c__2, &c__1, &c__1, a, &c__3, ip, &anrm, &rcond, w, r__,
00666 &info);
00667 chkxer_("CGBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00668 infoc_1.ok);
00669
00670
00671
00672 s_copy(srnamc_1.srnamt, "CGBEQU", (ftnlen)32, (ftnlen)6);
00673 infoc_1.infot = 1;
00674 cgbequ_(&c_n1, &c__0, &c__0, &c__0, a, &c__1, r1, r2, &rcond, &ccond,
00675 &anrm, &info);
00676 chkxer_("CGBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00677 infoc_1.ok);
00678 infoc_1.infot = 2;
00679 cgbequ_(&c__0, &c_n1, &c__0, &c__0, a, &c__1, r1, r2, &rcond, &ccond,
00680 &anrm, &info);
00681 chkxer_("CGBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00682 infoc_1.ok);
00683 infoc_1.infot = 3;
00684 cgbequ_(&c__1, &c__1, &c_n1, &c__0, a, &c__1, r1, r2, &rcond, &ccond,
00685 &anrm, &info);
00686 chkxer_("CGBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00687 infoc_1.ok);
00688 infoc_1.infot = 4;
00689 cgbequ_(&c__1, &c__1, &c__0, &c_n1, a, &c__1, r1, r2, &rcond, &ccond,
00690 &anrm, &info);
00691 chkxer_("CGBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00692 infoc_1.ok);
00693 infoc_1.infot = 6;
00694 cgbequ_(&c__2, &c__2, &c__1, &c__1, a, &c__2, r1, r2, &rcond, &ccond,
00695 &anrm, &info);
00696 chkxer_("CGBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00697 infoc_1.ok);
00698
00699
00700
00701 s_copy(srnamc_1.srnamt, "CGBEQUB", (ftnlen)32, (ftnlen)7);
00702 infoc_1.infot = 1;
00703 cgbequb_(&c_n1, &c__0, &c__0, &c__0, a, &c__1, r1, r2, &rcond, &ccond,
00704 &anrm, &info);
00705 chkxer_("CGBEQUB", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00706 infoc_1.ok);
00707 infoc_1.infot = 2;
00708 cgbequb_(&c__0, &c_n1, &c__0, &c__0, a, &c__1, r1, r2, &rcond, &ccond,
00709 &anrm, &info);
00710 chkxer_("CGBEQUB", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00711 infoc_1.ok);
00712 infoc_1.infot = 3;
00713 cgbequb_(&c__1, &c__1, &c_n1, &c__0, a, &c__1, r1, r2, &rcond, &ccond,
00714 &anrm, &info);
00715 chkxer_("CGBEQUB", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00716 infoc_1.ok);
00717 infoc_1.infot = 4;
00718 cgbequb_(&c__1, &c__1, &c__0, &c_n1, a, &c__1, r1, r2, &rcond, &ccond,
00719 &anrm, &info);
00720 chkxer_("CGBEQUB", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00721 infoc_1.ok);
00722 infoc_1.infot = 6;
00723 cgbequb_(&c__2, &c__2, &c__1, &c__1, a, &c__2, r1, r2, &rcond, &ccond,
00724 &anrm, &info);
00725 chkxer_("CGBEQUB", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00726 infoc_1.ok);
00727 }
00728
00729
00730
00731 alaesm_(path, &infoc_1.ok, &infoc_1.nout);
00732
00733 return 0;
00734
00735
00736
00737 }