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