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
00040 int zerrge_(char *path, integer *nunit)
00041 {
00042
00043 integer i__1;
00044 doublereal d__1, d__2;
00045 doublecomplex z__1;
00046
00047
00048 integer s_wsle(cilist *), e_wsle(void);
00049 int s_copy(char *, char *, ftnlen, ftnlen);
00050
00051
00052 doublecomplex a[16] , b[4];
00053 integer i__, j;
00054 doublereal r__[4];
00055 doublecomplex w[8], x[4];
00056 char c2[2];
00057 doublereal r1[4], r2[4];
00058 doublecomplex af[16] ;
00059 integer ip[4], info;
00060 doublereal anrm, ccond, rcond;
00061 extern int zgbtf2_(integer *, integer *, integer *,
00062 integer *, doublecomplex *, integer *, integer *, integer *),
00063 zgetf2_(integer *, integer *, doublecomplex *, integer *, integer
00064 *, integer *), alaesm_(char *, logical *, integer *);
00065 extern logical lsamen_(integer *, char *, char *);
00066 extern int zgbcon_(char *, integer *, integer *, integer
00067 *, doublecomplex *, integer *, integer *, doublereal *,
00068 doublereal *, doublecomplex *, doublereal *, integer *),
00069 chkxer_(char *, integer *, integer *, logical *, logical *), zgecon_(char *, integer *, doublecomplex *, integer *,
00070 doublereal *, doublereal *, doublecomplex *, doublereal *,
00071 integer *), zgbequ_(integer *, integer *, integer *,
00072 integer *, doublecomplex *, integer *, doublereal *, doublereal *,
00073 doublereal *, doublereal *, doublereal *, integer *), zgbrfs_(
00074 char *, integer *, integer *, integer *, integer *, doublecomplex
00075 *, integer *, doublecomplex *, integer *, integer *,
00076 doublecomplex *, integer *, doublecomplex *, integer *,
00077 doublereal *, doublereal *, doublecomplex *, doublereal *,
00078 integer *), zgbtrf_(integer *, integer *, integer *,
00079 integer *, doublecomplex *, integer *, integer *, integer *),
00080 zgeequ_(integer *, integer *, doublecomplex *, integer *,
00081 doublereal *, doublereal *, doublereal *, doublereal *,
00082 doublereal *, integer *), zgerfs_(char *, integer *, integer *,
00083 doublecomplex *, integer *, doublecomplex *, integer *, integer *,
00084 doublecomplex *, integer *, doublecomplex *, integer *,
00085 doublereal *, doublereal *, doublecomplex *, doublereal *,
00086 integer *), zgetrf_(integer *, integer *, doublecomplex *,
00087 integer *, integer *, integer *), zgetri_(integer *,
00088 doublecomplex *, integer *, integer *, doublecomplex *, integer *,
00089 integer *), zgbtrs_(char *, integer *, integer *, integer *,
00090 integer *, doublecomplex *, integer *, integer *, doublecomplex *,
00091 integer *, integer *), zgetrs_(char *, integer *,
00092 integer *, doublecomplex *, integer *, integer *, doublecomplex *,
00093 integer *, integer *);
00094
00095
00096 static cilist io___1 = { 0, 0, 0, 0, 0 };
00097
00098
00099
00100
00101
00102
00103
00104
00105
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 infoc_1.nout = *nunit;
00143 io___1.ciunit = infoc_1.nout;
00144 s_wsle(&io___1);
00145 e_wsle();
00146 s_copy(c2, path + 1, (ftnlen)2, (ftnlen)2);
00147
00148
00149
00150 for (j = 1; j <= 4; ++j) {
00151 for (i__ = 1; i__ <= 4; ++i__) {
00152 i__1 = i__ + (j << 2) - 5;
00153 d__1 = 1. / (doublereal) (i__ + j);
00154 d__2 = -1. / (doublereal) (i__ + j);
00155 z__1.r = d__1, z__1.i = d__2;
00156 a[i__1].r = z__1.r, a[i__1].i = z__1.i;
00157 i__1 = i__ + (j << 2) - 5;
00158 d__1 = 1. / (doublereal) (i__ + j);
00159 d__2 = -1. / (doublereal) (i__ + j);
00160 z__1.r = d__1, z__1.i = d__2;
00161 af[i__1].r = z__1.r, af[i__1].i = z__1.i;
00162
00163 }
00164 i__1 = j - 1;
00165 b[i__1].r = 0., b[i__1].i = 0.;
00166 r1[j - 1] = 0.;
00167 r2[j - 1] = 0.;
00168 i__1 = j - 1;
00169 w[i__1].r = 0., w[i__1].i = 0.;
00170 i__1 = j - 1;
00171 x[i__1].r = 0., x[i__1].i = 0.;
00172 ip[j - 1] = j;
00173
00174 }
00175 infoc_1.ok = TRUE_;
00176
00177
00178
00179
00180 if (lsamen_(&c__2, c2, "GE")) {
00181
00182
00183
00184 s_copy(srnamc_1.srnamt, "ZGETRF", (ftnlen)32, (ftnlen)6);
00185 infoc_1.infot = 1;
00186 zgetrf_(&c_n1, &c__0, a, &c__1, ip, &info);
00187 chkxer_("ZGETRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00188 infoc_1.ok);
00189 infoc_1.infot = 2;
00190 zgetrf_(&c__0, &c_n1, a, &c__1, ip, &info);
00191 chkxer_("ZGETRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00192 infoc_1.ok);
00193 infoc_1.infot = 4;
00194 zgetrf_(&c__2, &c__1, a, &c__1, ip, &info);
00195 chkxer_("ZGETRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00196 infoc_1.ok);
00197
00198
00199
00200 s_copy(srnamc_1.srnamt, "ZGETF2", (ftnlen)32, (ftnlen)6);
00201 infoc_1.infot = 1;
00202 zgetf2_(&c_n1, &c__0, a, &c__1, ip, &info);
00203 chkxer_("ZGETF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00204 infoc_1.ok);
00205 infoc_1.infot = 2;
00206 zgetf2_(&c__0, &c_n1, a, &c__1, ip, &info);
00207 chkxer_("ZGETF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00208 infoc_1.ok);
00209 infoc_1.infot = 4;
00210 zgetf2_(&c__2, &c__1, a, &c__1, ip, &info);
00211 chkxer_("ZGETF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00212 infoc_1.ok);
00213
00214
00215
00216 s_copy(srnamc_1.srnamt, "ZGETRI", (ftnlen)32, (ftnlen)6);
00217 infoc_1.infot = 1;
00218 zgetri_(&c_n1, a, &c__1, ip, w, &c__1, &info);
00219 chkxer_("ZGETRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00220 infoc_1.ok);
00221 infoc_1.infot = 3;
00222 zgetri_(&c__2, a, &c__1, ip, w, &c__2, &info);
00223 chkxer_("ZGETRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00224 infoc_1.ok);
00225 infoc_1.infot = 6;
00226 zgetri_(&c__2, a, &c__2, ip, w, &c__1, &info);
00227 chkxer_("ZGETRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00228 infoc_1.ok);
00229
00230
00231
00232 s_copy(srnamc_1.srnamt, "ZGETRS", (ftnlen)32, (ftnlen)6);
00233 infoc_1.infot = 1;
00234 zgetrs_("/", &c__0, &c__0, a, &c__1, ip, b, &c__1, &info);
00235 chkxer_("ZGETRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00236 infoc_1.ok);
00237 infoc_1.infot = 2;
00238 zgetrs_("N", &c_n1, &c__0, a, &c__1, ip, b, &c__1, &info);
00239 chkxer_("ZGETRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00240 infoc_1.ok);
00241 infoc_1.infot = 3;
00242 zgetrs_("N", &c__0, &c_n1, a, &c__1, ip, b, &c__1, &info);
00243 chkxer_("ZGETRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00244 infoc_1.ok);
00245 infoc_1.infot = 5;
00246 zgetrs_("N", &c__2, &c__1, a, &c__1, ip, b, &c__2, &info);
00247 chkxer_("ZGETRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00248 infoc_1.ok);
00249 infoc_1.infot = 8;
00250 zgetrs_("N", &c__2, &c__1, a, &c__2, ip, b, &c__1, &info);
00251 chkxer_("ZGETRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00252 infoc_1.ok);
00253
00254
00255
00256 s_copy(srnamc_1.srnamt, "ZGERFS", (ftnlen)32, (ftnlen)6);
00257 infoc_1.infot = 1;
00258 zgerfs_("/", &c__0, &c__0, a, &c__1, af, &c__1, ip, b, &c__1, x, &
00259 c__1, r1, r2, w, r__, &info);
00260 chkxer_("ZGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00261 infoc_1.ok);
00262 infoc_1.infot = 2;
00263 zgerfs_("N", &c_n1, &c__0, a, &c__1, af, &c__1, ip, b, &c__1, x, &
00264 c__1, r1, r2, w, r__, &info);
00265 chkxer_("ZGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00266 infoc_1.ok);
00267 infoc_1.infot = 3;
00268 zgerfs_("N", &c__0, &c_n1, a, &c__1, af, &c__1, ip, b, &c__1, x, &
00269 c__1, r1, r2, w, r__, &info);
00270 chkxer_("ZGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00271 infoc_1.ok);
00272 infoc_1.infot = 5;
00273 zgerfs_("N", &c__2, &c__1, a, &c__1, af, &c__2, ip, b, &c__2, x, &
00274 c__2, r1, r2, w, r__, &info);
00275 chkxer_("ZGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00276 infoc_1.ok);
00277 infoc_1.infot = 7;
00278 zgerfs_("N", &c__2, &c__1, a, &c__2, af, &c__1, ip, b, &c__2, x, &
00279 c__2, r1, r2, w, r__, &info);
00280 chkxer_("ZGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00281 infoc_1.ok);
00282 infoc_1.infot = 10;
00283 zgerfs_("N", &c__2, &c__1, a, &c__2, af, &c__2, ip, b, &c__1, x, &
00284 c__2, r1, r2, w, r__, &info);
00285 chkxer_("ZGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00286 infoc_1.ok);
00287 infoc_1.infot = 12;
00288 zgerfs_("N", &c__2, &c__1, a, &c__2, af, &c__2, ip, b, &c__2, x, &
00289 c__1, r1, r2, w, r__, &info);
00290 chkxer_("ZGERFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00291 infoc_1.ok);
00292
00293
00294
00295 s_copy(srnamc_1.srnamt, "ZGECON", (ftnlen)32, (ftnlen)6);
00296 infoc_1.infot = 1;
00297 zgecon_("/", &c__0, a, &c__1, &anrm, &rcond, w, r__, &info)
00298 ;
00299 chkxer_("ZGECON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00300 infoc_1.ok);
00301 infoc_1.infot = 2;
00302 zgecon_("1", &c_n1, a, &c__1, &anrm, &rcond, w, r__, &info)
00303 ;
00304 chkxer_("ZGECON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00305 infoc_1.ok);
00306 infoc_1.infot = 4;
00307 zgecon_("1", &c__2, a, &c__1, &anrm, &rcond, w, r__, &info)
00308 ;
00309 chkxer_("ZGECON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00310 infoc_1.ok);
00311
00312
00313
00314 s_copy(srnamc_1.srnamt, "ZGEEQU", (ftnlen)32, (ftnlen)6);
00315 infoc_1.infot = 1;
00316 zgeequ_(&c_n1, &c__0, a, &c__1, r1, r2, &rcond, &ccond, &anrm, &info);
00317 chkxer_("ZGEEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00318 infoc_1.ok);
00319 infoc_1.infot = 2;
00320 zgeequ_(&c__0, &c_n1, a, &c__1, r1, r2, &rcond, &ccond, &anrm, &info);
00321 chkxer_("ZGEEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00322 infoc_1.ok);
00323 infoc_1.infot = 4;
00324 zgeequ_(&c__2, &c__2, a, &c__1, r1, r2, &rcond, &ccond, &anrm, &info);
00325 chkxer_("ZGEEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00326 infoc_1.ok);
00327
00328
00329
00330
00331 } else if (lsamen_(&c__2, c2, "GB")) {
00332
00333
00334
00335 s_copy(srnamc_1.srnamt, "ZGBTRF", (ftnlen)32, (ftnlen)6);
00336 infoc_1.infot = 1;
00337 zgbtrf_(&c_n1, &c__0, &c__0, &c__0, a, &c__1, ip, &info);
00338 chkxer_("ZGBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00339 infoc_1.ok);
00340 infoc_1.infot = 2;
00341 zgbtrf_(&c__0, &c_n1, &c__0, &c__0, a, &c__1, ip, &info);
00342 chkxer_("ZGBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00343 infoc_1.ok);
00344 infoc_1.infot = 3;
00345 zgbtrf_(&c__1, &c__1, &c_n1, &c__0, a, &c__1, ip, &info);
00346 chkxer_("ZGBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00347 infoc_1.ok);
00348 infoc_1.infot = 4;
00349 zgbtrf_(&c__1, &c__1, &c__0, &c_n1, a, &c__1, ip, &info);
00350 chkxer_("ZGBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00351 infoc_1.ok);
00352 infoc_1.infot = 6;
00353 zgbtrf_(&c__2, &c__2, &c__1, &c__1, a, &c__3, ip, &info);
00354 chkxer_("ZGBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00355 infoc_1.ok);
00356
00357
00358
00359 s_copy(srnamc_1.srnamt, "ZGBTF2", (ftnlen)32, (ftnlen)6);
00360 infoc_1.infot = 1;
00361 zgbtf2_(&c_n1, &c__0, &c__0, &c__0, a, &c__1, ip, &info);
00362 chkxer_("ZGBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00363 infoc_1.ok);
00364 infoc_1.infot = 2;
00365 zgbtf2_(&c__0, &c_n1, &c__0, &c__0, a, &c__1, ip, &info);
00366 chkxer_("ZGBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00367 infoc_1.ok);
00368 infoc_1.infot = 3;
00369 zgbtf2_(&c__1, &c__1, &c_n1, &c__0, a, &c__1, ip, &info);
00370 chkxer_("ZGBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00371 infoc_1.ok);
00372 infoc_1.infot = 4;
00373 zgbtf2_(&c__1, &c__1, &c__0, &c_n1, a, &c__1, ip, &info);
00374 chkxer_("ZGBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00375 infoc_1.ok);
00376 infoc_1.infot = 6;
00377 zgbtf2_(&c__2, &c__2, &c__1, &c__1, a, &c__3, ip, &info);
00378 chkxer_("ZGBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00379 infoc_1.ok);
00380
00381
00382
00383 s_copy(srnamc_1.srnamt, "ZGBTRS", (ftnlen)32, (ftnlen)6);
00384 infoc_1.infot = 1;
00385 zgbtrs_("/", &c__0, &c__0, &c__0, &c__1, a, &c__1, ip, b, &c__1, &
00386 info);
00387 chkxer_("ZGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00388 infoc_1.ok);
00389 infoc_1.infot = 2;
00390 zgbtrs_("N", &c_n1, &c__0, &c__0, &c__1, a, &c__1, ip, b, &c__1, &
00391 info);
00392 chkxer_("ZGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00393 infoc_1.ok);
00394 infoc_1.infot = 3;
00395 zgbtrs_("N", &c__1, &c_n1, &c__0, &c__1, a, &c__1, ip, b, &c__1, &
00396 info);
00397 chkxer_("ZGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00398 infoc_1.ok);
00399 infoc_1.infot = 4;
00400 zgbtrs_("N", &c__1, &c__0, &c_n1, &c__1, a, &c__1, ip, b, &c__1, &
00401 info);
00402 chkxer_("ZGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00403 infoc_1.ok);
00404 infoc_1.infot = 5;
00405 zgbtrs_("N", &c__1, &c__0, &c__0, &c_n1, a, &c__1, ip, b, &c__1, &
00406 info);
00407 chkxer_("ZGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00408 infoc_1.ok);
00409 infoc_1.infot = 7;
00410 zgbtrs_("N", &c__2, &c__1, &c__1, &c__1, a, &c__3, ip, b, &c__2, &
00411 info);
00412 chkxer_("ZGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00413 infoc_1.ok);
00414 infoc_1.infot = 10;
00415 zgbtrs_("N", &c__2, &c__0, &c__0, &c__1, a, &c__1, ip, b, &c__1, &
00416 info);
00417 chkxer_("ZGBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00418 infoc_1.ok);
00419
00420
00421
00422 s_copy(srnamc_1.srnamt, "ZGBRFS", (ftnlen)32, (ftnlen)6);
00423 infoc_1.infot = 1;
00424 zgbrfs_("/", &c__0, &c__0, &c__0, &c__0, a, &c__1, af, &c__1, ip, b, &
00425 c__1, x, &c__1, r1, r2, w, r__, &info);
00426 chkxer_("ZGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00427 infoc_1.ok);
00428 infoc_1.infot = 2;
00429 zgbrfs_("N", &c_n1, &c__0, &c__0, &c__0, a, &c__1, af, &c__1, ip, b, &
00430 c__1, x, &c__1, r1, r2, w, r__, &info);
00431 chkxer_("ZGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00432 infoc_1.ok);
00433 infoc_1.infot = 3;
00434 zgbrfs_("N", &c__1, &c_n1, &c__0, &c__0, a, &c__1, af, &c__1, ip, b, &
00435 c__1, x, &c__1, r1, r2, w, r__, &info);
00436 chkxer_("ZGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00437 infoc_1.ok);
00438 infoc_1.infot = 4;
00439 zgbrfs_("N", &c__1, &c__0, &c_n1, &c__0, a, &c__1, af, &c__1, ip, b, &
00440 c__1, x, &c__1, r1, r2, w, r__, &info);
00441 chkxer_("ZGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00442 infoc_1.ok);
00443 infoc_1.infot = 5;
00444 zgbrfs_("N", &c__1, &c__0, &c__0, &c_n1, a, &c__1, af, &c__1, ip, b, &
00445 c__1, x, &c__1, r1, r2, w, r__, &info);
00446 chkxer_("ZGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00447 infoc_1.ok);
00448 infoc_1.infot = 7;
00449 zgbrfs_("N", &c__2, &c__1, &c__1, &c__1, a, &c__2, af, &c__4, ip, b, &
00450 c__2, x, &c__2, r1, r2, w, r__, &info);
00451 chkxer_("ZGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00452 infoc_1.ok);
00453 infoc_1.infot = 9;
00454 zgbrfs_("N", &c__2, &c__1, &c__1, &c__1, a, &c__3, af, &c__3, ip, b, &
00455 c__2, x, &c__2, r1, r2, w, r__, &info);
00456 chkxer_("ZGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00457 infoc_1.ok);
00458 infoc_1.infot = 12;
00459 zgbrfs_("N", &c__2, &c__0, &c__0, &c__1, a, &c__1, af, &c__1, ip, b, &
00460 c__1, x, &c__2, r1, r2, w, r__, &info);
00461 chkxer_("ZGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00462 infoc_1.ok);
00463 infoc_1.infot = 14;
00464 zgbrfs_("N", &c__2, &c__0, &c__0, &c__1, a, &c__1, af, &c__1, ip, b, &
00465 c__2, x, &c__1, r1, r2, w, r__, &info);
00466 chkxer_("ZGBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00467 infoc_1.ok);
00468
00469
00470
00471 s_copy(srnamc_1.srnamt, "ZGBCON", (ftnlen)32, (ftnlen)6);
00472 infoc_1.infot = 1;
00473 zgbcon_("/", &c__0, &c__0, &c__0, a, &c__1, ip, &anrm, &rcond, w, r__,
00474 &info);
00475 chkxer_("ZGBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00476 infoc_1.ok);
00477 infoc_1.infot = 2;
00478 zgbcon_("1", &c_n1, &c__0, &c__0, a, &c__1, ip, &anrm, &rcond, w, r__,
00479 &info);
00480 chkxer_("ZGBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00481 infoc_1.ok);
00482 infoc_1.infot = 3;
00483 zgbcon_("1", &c__1, &c_n1, &c__0, a, &c__1, ip, &anrm, &rcond, w, r__,
00484 &info);
00485 chkxer_("ZGBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00486 infoc_1.ok);
00487 infoc_1.infot = 4;
00488 zgbcon_("1", &c__1, &c__0, &c_n1, a, &c__1, ip, &anrm, &rcond, w, r__,
00489 &info);
00490 chkxer_("ZGBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00491 infoc_1.ok);
00492 infoc_1.infot = 6;
00493 zgbcon_("1", &c__2, &c__1, &c__1, a, &c__3, ip, &anrm, &rcond, w, r__,
00494 &info);
00495 chkxer_("ZGBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00496 infoc_1.ok);
00497
00498
00499
00500 s_copy(srnamc_1.srnamt, "ZGBEQU", (ftnlen)32, (ftnlen)6);
00501 infoc_1.infot = 1;
00502 zgbequ_(&c_n1, &c__0, &c__0, &c__0, a, &c__1, r1, r2, &rcond, &ccond,
00503 &anrm, &info);
00504 chkxer_("ZGBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00505 infoc_1.ok);
00506 infoc_1.infot = 2;
00507 zgbequ_(&c__0, &c_n1, &c__0, &c__0, a, &c__1, r1, r2, &rcond, &ccond,
00508 &anrm, &info);
00509 chkxer_("ZGBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00510 infoc_1.ok);
00511 infoc_1.infot = 3;
00512 zgbequ_(&c__1, &c__1, &c_n1, &c__0, a, &c__1, r1, r2, &rcond, &ccond,
00513 &anrm, &info);
00514 chkxer_("ZGBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00515 infoc_1.ok);
00516 infoc_1.infot = 4;
00517 zgbequ_(&c__1, &c__1, &c__0, &c_n1, a, &c__1, r1, r2, &rcond, &ccond,
00518 &anrm, &info);
00519 chkxer_("ZGBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00520 infoc_1.ok);
00521 infoc_1.infot = 6;
00522 zgbequ_(&c__2, &c__2, &c__1, &c__1, a, &c__2, r1, r2, &rcond, &ccond,
00523 &anrm, &info);
00524 chkxer_("ZGBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00525 infoc_1.ok);
00526 }
00527
00528
00529
00530 alaesm_(path, &infoc_1.ok, &infoc_1.nout);
00531
00532 return 0;
00533
00534
00535
00536 }