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