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