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