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