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