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