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