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__0 = 0;
00035 static integer c__1 = 1;
00036 static integer c_n1 = -1;
00037
00038 int cerrpo_(char *path, integer *nunit)
00039 {
00040
00041 integer i__1;
00042 real r__1, r__2;
00043 complex q__1;
00044
00045
00046 integer s_wsle(cilist *), e_wsle(void);
00047 int s_copy(char *, char *, ftnlen, ftnlen);
00048
00049
00050 complex a[16] , b[4];
00051 integer i__, j;
00052 real r__[4];
00053 complex w[8], x[4];
00054 char c2[2];
00055 real r1[4], r2[4];
00056 complex af[16] ;
00057 integer info;
00058 real anrm, rcond;
00059 extern int cpbtf2_(char *, integer *, integer *, complex
00060 *, integer *, integer *), cpotf2_(char *, integer *,
00061 complex *, integer *, integer *), alaesm_(char *, logical
00062 *, integer *), cpbcon_(char *, integer *, integer *,
00063 complex *, integer *, real *, real *, complex *, real *, integer *
00064 );
00065 extern logical lsamen_(integer *, char *, char *);
00066 extern int cpbequ_(char *, integer *, integer *, complex
00067 *, integer *, real *, real *, real *, integer *), cpbrfs_(
00068 char *, integer *, integer *, integer *, complex *, integer *,
00069 complex *, integer *, complex *, integer *, complex *, integer *,
00070 real *, real *, complex *, real *, integer *), cpbtrf_(
00071 char *, integer *, integer *, complex *, integer *, integer *), cpocon_(char *, integer *, complex *, integer *, real *,
00072 real *, complex *, real *, integer *), chkxer_(char *,
00073 integer *, integer *, logical *, logical *), cppcon_(char
00074 *, integer *, complex *, real *, real *, complex *, real *,
00075 integer *), cpoequ_(integer *, complex *, integer *, real
00076 *, real *, real *, integer *), cpbtrs_(char *, integer *, integer
00077 *, integer *, complex *, integer *, complex *, integer *, integer
00078 *), cporfs_(char *, integer *, integer *, complex *,
00079 integer *, complex *, integer *, complex *, integer *, complex *,
00080 integer *, real *, real *, complex *, real *, integer *),
00081 cpotrf_(char *, integer *, complex *, integer *, integer *), cpotri_(char *, integer *, complex *, integer *, integer
00082 *), cppequ_(char *, integer *, complex *, real *, real *,
00083 real *, integer *), cpprfs_(char *, integer *, integer *,
00084 complex *, complex *, complex *, integer *, complex *, integer *,
00085 real *, real *, complex *, real *, integer *), cpptrf_(
00086 char *, integer *, complex *, integer *), cpptri_(char *,
00087 integer *, complex *, integer *), cpotrs_(char *, integer
00088 *, integer *, complex *, integer *, complex *, integer *, integer
00089 *), cpptrs_(char *, integer *, integer *, complex *,
00090 complex *, integer *, integer *);
00091
00092
00093 static cilist io___1 = { 0, 0, 0, 0, 0 };
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
00139 infoc_1.nout = *nunit;
00140 io___1.ciunit = infoc_1.nout;
00141 s_wsle(&io___1);
00142 e_wsle();
00143 s_copy(c2, path + 1, (ftnlen)2, (ftnlen)2);
00144
00145
00146
00147 for (j = 1; j <= 4; ++j) {
00148 for (i__ = 1; i__ <= 4; ++i__) {
00149 i__1 = i__ + (j << 2) - 5;
00150 r__1 = 1.f / (real) (i__ + j);
00151 r__2 = -1.f / (real) (i__ + j);
00152 q__1.r = r__1, q__1.i = r__2;
00153 a[i__1].r = q__1.r, a[i__1].i = q__1.i;
00154 i__1 = i__ + (j << 2) - 5;
00155 r__1 = 1.f / (real) (i__ + j);
00156 r__2 = -1.f / (real) (i__ + j);
00157 q__1.r = r__1, q__1.i = r__2;
00158 af[i__1].r = q__1.r, af[i__1].i = q__1.i;
00159
00160 }
00161 i__1 = j - 1;
00162 b[i__1].r = 0.f, b[i__1].i = 0.f;
00163 r1[j - 1] = 0.f;
00164 r2[j - 1] = 0.f;
00165 i__1 = j - 1;
00166 w[i__1].r = 0.f, w[i__1].i = 0.f;
00167 i__1 = j - 1;
00168 x[i__1].r = 0.f, x[i__1].i = 0.f;
00169
00170 }
00171 anrm = 1.f;
00172 infoc_1.ok = TRUE_;
00173
00174
00175
00176
00177 if (lsamen_(&c__2, c2, "PO")) {
00178
00179
00180
00181 s_copy(srnamc_1.srnamt, "CPOTRF", (ftnlen)32, (ftnlen)6);
00182 infoc_1.infot = 1;
00183 cpotrf_("/", &c__0, a, &c__1, &info);
00184 chkxer_("CPOTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00185 infoc_1.ok);
00186 infoc_1.infot = 2;
00187 cpotrf_("U", &c_n1, a, &c__1, &info);
00188 chkxer_("CPOTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00189 infoc_1.ok);
00190 infoc_1.infot = 4;
00191 cpotrf_("U", &c__2, a, &c__1, &info);
00192 chkxer_("CPOTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00193 infoc_1.ok);
00194
00195
00196
00197 s_copy(srnamc_1.srnamt, "CPOTF2", (ftnlen)32, (ftnlen)6);
00198 infoc_1.infot = 1;
00199 cpotf2_("/", &c__0, a, &c__1, &info);
00200 chkxer_("CPOTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00201 infoc_1.ok);
00202 infoc_1.infot = 2;
00203 cpotf2_("U", &c_n1, a, &c__1, &info);
00204 chkxer_("CPOTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00205 infoc_1.ok);
00206 infoc_1.infot = 4;
00207 cpotf2_("U", &c__2, a, &c__1, &info);
00208 chkxer_("CPOTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00209 infoc_1.ok);
00210
00211
00212
00213 s_copy(srnamc_1.srnamt, "CPOTRI", (ftnlen)32, (ftnlen)6);
00214 infoc_1.infot = 1;
00215 cpotri_("/", &c__0, a, &c__1, &info);
00216 chkxer_("CPOTRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00217 infoc_1.ok);
00218 infoc_1.infot = 2;
00219 cpotri_("U", &c_n1, a, &c__1, &info);
00220 chkxer_("CPOTRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00221 infoc_1.ok);
00222 infoc_1.infot = 4;
00223 cpotri_("U", &c__2, a, &c__1, &info);
00224 chkxer_("CPOTRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00225 infoc_1.ok);
00226
00227
00228
00229 s_copy(srnamc_1.srnamt, "CPOTRS", (ftnlen)32, (ftnlen)6);
00230 infoc_1.infot = 1;
00231 cpotrs_("/", &c__0, &c__0, a, &c__1, b, &c__1, &info);
00232 chkxer_("CPOTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00233 infoc_1.ok);
00234 infoc_1.infot = 2;
00235 cpotrs_("U", &c_n1, &c__0, a, &c__1, b, &c__1, &info);
00236 chkxer_("CPOTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00237 infoc_1.ok);
00238 infoc_1.infot = 3;
00239 cpotrs_("U", &c__0, &c_n1, a, &c__1, b, &c__1, &info);
00240 chkxer_("CPOTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00241 infoc_1.ok);
00242 infoc_1.infot = 5;
00243 cpotrs_("U", &c__2, &c__1, a, &c__1, b, &c__2, &info);
00244 chkxer_("CPOTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00245 infoc_1.ok);
00246 infoc_1.infot = 7;
00247 cpotrs_("U", &c__2, &c__1, a, &c__2, b, &c__1, &info);
00248 chkxer_("CPOTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00249 infoc_1.ok);
00250
00251
00252
00253 s_copy(srnamc_1.srnamt, "CPORFS", (ftnlen)32, (ftnlen)6);
00254 infoc_1.infot = 1;
00255 cporfs_("/", &c__0, &c__0, a, &c__1, af, &c__1, b, &c__1, x, &c__1,
00256 r1, r2, w, r__, &info);
00257 chkxer_("CPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00258 infoc_1.ok);
00259 infoc_1.infot = 2;
00260 cporfs_("U", &c_n1, &c__0, a, &c__1, af, &c__1, b, &c__1, x, &c__1,
00261 r1, r2, w, r__, &info);
00262 chkxer_("CPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00263 infoc_1.ok);
00264 infoc_1.infot = 3;
00265 cporfs_("U", &c__0, &c_n1, a, &c__1, af, &c__1, b, &c__1, x, &c__1,
00266 r1, r2, w, r__, &info);
00267 chkxer_("CPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00268 infoc_1.ok);
00269 infoc_1.infot = 5;
00270 cporfs_("U", &c__2, &c__1, a, &c__1, af, &c__2, b, &c__2, x, &c__2,
00271 r1, r2, w, r__, &info);
00272 chkxer_("CPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00273 infoc_1.ok);
00274 infoc_1.infot = 7;
00275 cporfs_("U", &c__2, &c__1, a, &c__2, af, &c__1, b, &c__2, x, &c__2,
00276 r1, r2, w, r__, &info);
00277 chkxer_("CPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00278 infoc_1.ok);
00279 infoc_1.infot = 9;
00280 cporfs_("U", &c__2, &c__1, a, &c__2, af, &c__2, b, &c__1, x, &c__2,
00281 r1, r2, w, r__, &info);
00282 chkxer_("CPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00283 infoc_1.ok);
00284 infoc_1.infot = 11;
00285 cporfs_("U", &c__2, &c__1, a, &c__2, af, &c__2, b, &c__2, x, &c__1,
00286 r1, r2, w, r__, &info);
00287 chkxer_("CPORFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00288 infoc_1.ok);
00289
00290
00291
00292 s_copy(srnamc_1.srnamt, "CPOCON", (ftnlen)32, (ftnlen)6);
00293 infoc_1.infot = 1;
00294 cpocon_("/", &c__0, a, &c__1, &anrm, &rcond, w, r__, &info)
00295 ;
00296 chkxer_("CPOCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00297 infoc_1.ok);
00298 infoc_1.infot = 2;
00299 cpocon_("U", &c_n1, a, &c__1, &anrm, &rcond, w, r__, &info)
00300 ;
00301 chkxer_("CPOCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00302 infoc_1.ok);
00303 infoc_1.infot = 4;
00304 cpocon_("U", &c__2, a, &c__1, &anrm, &rcond, w, r__, &info)
00305 ;
00306 chkxer_("CPOCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00307 infoc_1.ok);
00308 infoc_1.infot = 5;
00309 r__1 = -anrm;
00310 cpocon_("U", &c__1, a, &c__1, &r__1, &rcond, w, r__, &info)
00311 ;
00312 chkxer_("CPOCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00313 infoc_1.ok);
00314
00315
00316
00317 s_copy(srnamc_1.srnamt, "CPOEQU", (ftnlen)32, (ftnlen)6);
00318 infoc_1.infot = 1;
00319 cpoequ_(&c_n1, a, &c__1, r1, &rcond, &anrm, &info);
00320 chkxer_("CPOEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00321 infoc_1.ok);
00322 infoc_1.infot = 3;
00323 cpoequ_(&c__2, a, &c__1, r1, &rcond, &anrm, &info);
00324 chkxer_("CPOEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00325 infoc_1.ok);
00326
00327
00328
00329
00330 } else if (lsamen_(&c__2, c2, "PP")) {
00331
00332
00333
00334 s_copy(srnamc_1.srnamt, "CPPTRF", (ftnlen)32, (ftnlen)6);
00335 infoc_1.infot = 1;
00336 cpptrf_("/", &c__0, a, &info);
00337 chkxer_("CPPTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00338 infoc_1.ok);
00339 infoc_1.infot = 2;
00340 cpptrf_("U", &c_n1, a, &info);
00341 chkxer_("CPPTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00342 infoc_1.ok);
00343
00344
00345
00346 s_copy(srnamc_1.srnamt, "CPPTRI", (ftnlen)32, (ftnlen)6);
00347 infoc_1.infot = 1;
00348 cpptri_("/", &c__0, a, &info);
00349 chkxer_("CPPTRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00350 infoc_1.ok);
00351 infoc_1.infot = 2;
00352 cpptri_("U", &c_n1, a, &info);
00353 chkxer_("CPPTRI", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00354 infoc_1.ok);
00355
00356
00357
00358 s_copy(srnamc_1.srnamt, "CPPTRS", (ftnlen)32, (ftnlen)6);
00359 infoc_1.infot = 1;
00360 cpptrs_("/", &c__0, &c__0, a, b, &c__1, &info);
00361 chkxer_("CPPTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00362 infoc_1.ok);
00363 infoc_1.infot = 2;
00364 cpptrs_("U", &c_n1, &c__0, a, b, &c__1, &info);
00365 chkxer_("CPPTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00366 infoc_1.ok);
00367 infoc_1.infot = 3;
00368 cpptrs_("U", &c__0, &c_n1, a, b, &c__1, &info);
00369 chkxer_("CPPTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00370 infoc_1.ok);
00371 infoc_1.infot = 6;
00372 cpptrs_("U", &c__2, &c__1, a, b, &c__1, &info);
00373 chkxer_("CPPTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00374 infoc_1.ok);
00375
00376
00377
00378 s_copy(srnamc_1.srnamt, "CPPRFS", (ftnlen)32, (ftnlen)6);
00379 infoc_1.infot = 1;
00380 cpprfs_("/", &c__0, &c__0, a, af, b, &c__1, x, &c__1, r1, r2, w, r__,
00381 &info);
00382 chkxer_("CPPRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00383 infoc_1.ok);
00384 infoc_1.infot = 2;
00385 cpprfs_("U", &c_n1, &c__0, a, af, b, &c__1, x, &c__1, r1, r2, w, r__,
00386 &info);
00387 chkxer_("CPPRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00388 infoc_1.ok);
00389 infoc_1.infot = 3;
00390 cpprfs_("U", &c__0, &c_n1, a, af, b, &c__1, x, &c__1, r1, r2, w, r__,
00391 &info);
00392 chkxer_("CPPRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00393 infoc_1.ok);
00394 infoc_1.infot = 7;
00395 cpprfs_("U", &c__2, &c__1, a, af, b, &c__1, x, &c__2, r1, r2, w, r__,
00396 &info);
00397 chkxer_("CPPRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00398 infoc_1.ok);
00399 infoc_1.infot = 9;
00400 cpprfs_("U", &c__2, &c__1, a, af, b, &c__2, x, &c__1, r1, r2, w, r__,
00401 &info);
00402 chkxer_("CPPRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00403 infoc_1.ok);
00404
00405
00406
00407 s_copy(srnamc_1.srnamt, "CPPCON", (ftnlen)32, (ftnlen)6);
00408 infoc_1.infot = 1;
00409 cppcon_("/", &c__0, a, &anrm, &rcond, w, r__, &info);
00410 chkxer_("CPPCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00411 infoc_1.ok);
00412 infoc_1.infot = 2;
00413 cppcon_("U", &c_n1, a, &anrm, &rcond, w, r__, &info);
00414 chkxer_("CPPCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00415 infoc_1.ok);
00416 infoc_1.infot = 4;
00417 r__1 = -anrm;
00418 cppcon_("U", &c__1, a, &r__1, &rcond, w, r__, &info);
00419 chkxer_("CPPCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00420 infoc_1.ok);
00421
00422
00423
00424 s_copy(srnamc_1.srnamt, "CPPEQU", (ftnlen)32, (ftnlen)6);
00425 infoc_1.infot = 1;
00426 cppequ_("/", &c__0, a, r1, &rcond, &anrm, &info);
00427 chkxer_("CPPEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00428 infoc_1.ok);
00429 infoc_1.infot = 2;
00430 cppequ_("U", &c_n1, a, r1, &rcond, &anrm, &info);
00431 chkxer_("CPPEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00432 infoc_1.ok);
00433
00434
00435
00436
00437 } else if (lsamen_(&c__2, c2, "PB")) {
00438
00439
00440
00441 s_copy(srnamc_1.srnamt, "CPBTRF", (ftnlen)32, (ftnlen)6);
00442 infoc_1.infot = 1;
00443 cpbtrf_("/", &c__0, &c__0, a, &c__1, &info);
00444 chkxer_("CPBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00445 infoc_1.ok);
00446 infoc_1.infot = 2;
00447 cpbtrf_("U", &c_n1, &c__0, a, &c__1, &info);
00448 chkxer_("CPBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00449 infoc_1.ok);
00450 infoc_1.infot = 3;
00451 cpbtrf_("U", &c__1, &c_n1, a, &c__1, &info);
00452 chkxer_("CPBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00453 infoc_1.ok);
00454 infoc_1.infot = 5;
00455 cpbtrf_("U", &c__2, &c__1, a, &c__1, &info);
00456 chkxer_("CPBTRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00457 infoc_1.ok);
00458
00459
00460
00461 s_copy(srnamc_1.srnamt, "CPBTF2", (ftnlen)32, (ftnlen)6);
00462 infoc_1.infot = 1;
00463 cpbtf2_("/", &c__0, &c__0, a, &c__1, &info);
00464 chkxer_("CPBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00465 infoc_1.ok);
00466 infoc_1.infot = 2;
00467 cpbtf2_("U", &c_n1, &c__0, a, &c__1, &info);
00468 chkxer_("CPBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00469 infoc_1.ok);
00470 infoc_1.infot = 3;
00471 cpbtf2_("U", &c__1, &c_n1, a, &c__1, &info);
00472 chkxer_("CPBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00473 infoc_1.ok);
00474 infoc_1.infot = 5;
00475 cpbtf2_("U", &c__2, &c__1, a, &c__1, &info);
00476 chkxer_("CPBTF2", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00477 infoc_1.ok);
00478
00479
00480
00481 s_copy(srnamc_1.srnamt, "CPBTRS", (ftnlen)32, (ftnlen)6);
00482 infoc_1.infot = 1;
00483 cpbtrs_("/", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, &info);
00484 chkxer_("CPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00485 infoc_1.ok);
00486 infoc_1.infot = 2;
00487 cpbtrs_("U", &c_n1, &c__0, &c__0, a, &c__1, b, &c__1, &info);
00488 chkxer_("CPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00489 infoc_1.ok);
00490 infoc_1.infot = 3;
00491 cpbtrs_("U", &c__1, &c_n1, &c__0, a, &c__1, b, &c__1, &info);
00492 chkxer_("CPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00493 infoc_1.ok);
00494 infoc_1.infot = 4;
00495 cpbtrs_("U", &c__0, &c__0, &c_n1, a, &c__1, b, &c__1, &info);
00496 chkxer_("CPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00497 infoc_1.ok);
00498 infoc_1.infot = 6;
00499 cpbtrs_("U", &c__2, &c__1, &c__1, a, &c__1, b, &c__1, &info);
00500 chkxer_("CPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00501 infoc_1.ok);
00502 infoc_1.infot = 8;
00503 cpbtrs_("U", &c__2, &c__0, &c__1, a, &c__1, b, &c__1, &info);
00504 chkxer_("CPBTRS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00505 infoc_1.ok);
00506
00507
00508
00509 s_copy(srnamc_1.srnamt, "CPBRFS", (ftnlen)32, (ftnlen)6);
00510 infoc_1.infot = 1;
00511 cpbrfs_("/", &c__0, &c__0, &c__0, a, &c__1, af, &c__1, b, &c__1, x, &
00512 c__1, r1, r2, w, r__, &info);
00513 chkxer_("CPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00514 infoc_1.ok);
00515 infoc_1.infot = 2;
00516 cpbrfs_("U", &c_n1, &c__0, &c__0, a, &c__1, af, &c__1, b, &c__1, x, &
00517 c__1, r1, r2, w, r__, &info);
00518 chkxer_("CPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00519 infoc_1.ok);
00520 infoc_1.infot = 3;
00521 cpbrfs_("U", &c__1, &c_n1, &c__0, a, &c__1, af, &c__1, b, &c__1, x, &
00522 c__1, r1, r2, w, r__, &info);
00523 chkxer_("CPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00524 infoc_1.ok);
00525 infoc_1.infot = 4;
00526 cpbrfs_("U", &c__0, &c__0, &c_n1, a, &c__1, af, &c__1, b, &c__1, x, &
00527 c__1, r1, r2, w, r__, &info);
00528 chkxer_("CPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00529 infoc_1.ok);
00530 infoc_1.infot = 6;
00531 cpbrfs_("U", &c__2, &c__1, &c__1, a, &c__1, af, &c__2, b, &c__2, x, &
00532 c__2, r1, r2, w, r__, &info);
00533 chkxer_("CPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00534 infoc_1.ok);
00535 infoc_1.infot = 8;
00536 cpbrfs_("U", &c__2, &c__1, &c__1, a, &c__2, af, &c__1, b, &c__2, x, &
00537 c__2, r1, r2, w, r__, &info);
00538 chkxer_("CPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00539 infoc_1.ok);
00540 infoc_1.infot = 10;
00541 cpbrfs_("U", &c__2, &c__0, &c__1, a, &c__1, af, &c__1, b, &c__1, x, &
00542 c__2, r1, r2, w, r__, &info);
00543 chkxer_("CPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00544 infoc_1.ok);
00545 infoc_1.infot = 12;
00546 cpbrfs_("U", &c__2, &c__0, &c__1, a, &c__1, af, &c__1, b, &c__2, x, &
00547 c__1, r1, r2, w, r__, &info);
00548 chkxer_("CPBRFS", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00549 infoc_1.ok);
00550
00551
00552
00553 s_copy(srnamc_1.srnamt, "CPBCON", (ftnlen)32, (ftnlen)6);
00554 infoc_1.infot = 1;
00555 cpbcon_("/", &c__0, &c__0, a, &c__1, &anrm, &rcond, w, r__, &info);
00556 chkxer_("CPBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00557 infoc_1.ok);
00558 infoc_1.infot = 2;
00559 cpbcon_("U", &c_n1, &c__0, a, &c__1, &anrm, &rcond, w, r__, &info);
00560 chkxer_("CPBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00561 infoc_1.ok);
00562 infoc_1.infot = 3;
00563 cpbcon_("U", &c__1, &c_n1, a, &c__1, &anrm, &rcond, w, r__, &info);
00564 chkxer_("CPBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00565 infoc_1.ok);
00566 infoc_1.infot = 5;
00567 cpbcon_("U", &c__2, &c__1, a, &c__1, &anrm, &rcond, w, r__, &info);
00568 chkxer_("CPBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00569 infoc_1.ok);
00570 infoc_1.infot = 6;
00571 r__1 = -anrm;
00572 cpbcon_("U", &c__1, &c__0, a, &c__1, &r__1, &rcond, w, r__, &info);
00573 chkxer_("CPBCON", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00574 infoc_1.ok);
00575
00576
00577
00578 s_copy(srnamc_1.srnamt, "CPBEQU", (ftnlen)32, (ftnlen)6);
00579 infoc_1.infot = 1;
00580 cpbequ_("/", &c__0, &c__0, a, &c__1, r1, &rcond, &anrm, &info);
00581 chkxer_("CPBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00582 infoc_1.ok);
00583 infoc_1.infot = 2;
00584 cpbequ_("U", &c_n1, &c__0, a, &c__1, r1, &rcond, &anrm, &info);
00585 chkxer_("CPBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00586 infoc_1.ok);
00587 infoc_1.infot = 3;
00588 cpbequ_("U", &c__1, &c_n1, a, &c__1, r1, &rcond, &anrm, &info);
00589 chkxer_("CPBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00590 infoc_1.ok);
00591 infoc_1.infot = 5;
00592 cpbequ_("U", &c__2, &c__1, a, &c__1, r1, &rcond, &anrm, &info);
00593 chkxer_("CPBEQU", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00594 infoc_1.ok);
00595 }
00596
00597
00598
00599 alaesm_(path, &infoc_1.ok, &infoc_1.nout);
00600
00601 return 0;
00602
00603
00604
00605 }