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