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 static integer c__3 = 3;
00038 static integer c__4 = 4;
00039
00040 int derrhs_(char *path, integer *nunit)
00041 {
00042
00043 static char fmt_9999[] = "(1x,a3,\002 routines passed the tests of the e"
00044 "rror exits\002,\002 (\002,i3,\002 tests done)\002)";
00045 static char fmt_9998[] = "(\002 *** \002,a3,\002 routines failed the tes"
00046 "ts of the error \002,\002exits ***\002)";
00047
00048
00049 integer s_wsle(cilist *), e_wsle(void);
00050 int s_copy(char *, char *, ftnlen, ftnlen);
00051 integer s_wsfe(cilist *), do_fio(integer *, char *, ftnlen), e_wsfe(void);
00052
00053
00054 doublereal a[9] , c__[9] ;
00055 integer i__, j, m;
00056 doublereal s[3], w[28];
00057 char c2[2];
00058 doublereal wi[3];
00059 integer nt;
00060 doublereal vl[9] , vr[9] , wr[3];
00061 integer ihi, ilo;
00062 logical sel[3];
00063 doublereal tau[3];
00064 integer info;
00065 extern int dgebak_(char *, char *, integer *, integer *,
00066 integer *, doublereal *, integer *, doublereal *, integer *,
00067 integer *), dgebal_(char *, integer *, doublereal
00068 *, integer *, integer *, integer *, doublereal *, integer *), dgehrd_(integer *, integer *, integer *, doublereal *,
00069 integer *, doublereal *, doublereal *, integer *, integer *);
00070 integer ifaill[3], ifailr[3];
00071 extern int dhsein_(char *, char *, char *, logical *,
00072 integer *, doublereal *, integer *, doublereal *, doublereal *,
00073 doublereal *, integer *, doublereal *, integer *, integer *,
00074 integer *, doublereal *, integer *, integer *, integer *);
00075 extern logical lsamen_(integer *, char *, char *);
00076 extern int chkxer_(char *, integer *, integer *, logical
00077 *, logical *), dorghr_(integer *, integer *, integer *,
00078 doublereal *, integer *, doublereal *, doublereal *, integer *,
00079 integer *), dhseqr_(char *, char *, integer *, integer *, integer
00080 *, doublereal *, integer *, doublereal *, doublereal *,
00081 doublereal *, integer *, doublereal *, integer *, integer *), dtrevc_(char *, char *, logical *, integer *,
00082 doublereal *, integer *, doublereal *, integer *, doublereal *,
00083 integer *, integer *, integer *, doublereal *, integer *), dormhr_(char *, char *, integer *, integer *, integer *,
00084 integer *, doublereal *, integer *, doublereal *, doublereal *,
00085 integer *, doublereal *, integer *, integer *);
00086
00087
00088 static cilist io___1 = { 0, 0, 0, 0, 0 };
00089 static cilist io___22 = { 0, 0, 0, fmt_9999, 0 };
00090 static cilist io___23 = { 0, 0, 0, fmt_9998, 0 };
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
00131
00132
00133
00134
00135
00136 infoc_1.nout = *nunit;
00137 io___1.ciunit = infoc_1.nout;
00138 s_wsle(&io___1);
00139 e_wsle();
00140 s_copy(c2, path + 1, (ftnlen)2, (ftnlen)2);
00141
00142
00143
00144 for (j = 1; j <= 3; ++j) {
00145 for (i__ = 1; i__ <= 3; ++i__) {
00146 a[i__ + j * 3 - 4] = 1. / (doublereal) (i__ + j);
00147
00148 }
00149 wi[j - 1] = (doublereal) j;
00150 sel[j - 1] = TRUE_;
00151
00152 }
00153 infoc_1.ok = TRUE_;
00154 nt = 0;
00155
00156
00157
00158 if (lsamen_(&c__2, c2, "HS")) {
00159
00160
00161
00162 s_copy(srnamc_1.srnamt, "DGEBAL", (ftnlen)32, (ftnlen)6);
00163 infoc_1.infot = 1;
00164 dgebal_("/", &c__0, a, &c__1, &ilo, &ihi, s, &info);
00165 chkxer_("DGEBAL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00166 infoc_1.ok);
00167 infoc_1.infot = 2;
00168 dgebal_("N", &c_n1, a, &c__1, &ilo, &ihi, s, &info);
00169 chkxer_("DGEBAL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00170 infoc_1.ok);
00171 infoc_1.infot = 4;
00172 dgebal_("N", &c__2, a, &c__1, &ilo, &ihi, s, &info);
00173 chkxer_("DGEBAL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00174 infoc_1.ok);
00175 nt += 3;
00176
00177
00178
00179 s_copy(srnamc_1.srnamt, "DGEBAK", (ftnlen)32, (ftnlen)6);
00180 infoc_1.infot = 1;
00181 dgebak_("/", "R", &c__0, &c__1, &c__0, s, &c__0, a, &c__1, &info);
00182 chkxer_("DGEBAK", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00183 infoc_1.ok);
00184 infoc_1.infot = 2;
00185 dgebak_("N", "/", &c__0, &c__1, &c__0, s, &c__0, a, &c__1, &info);
00186 chkxer_("DGEBAK", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00187 infoc_1.ok);
00188 infoc_1.infot = 3;
00189 dgebak_("N", "R", &c_n1, &c__1, &c__0, s, &c__0, a, &c__1, &info);
00190 chkxer_("DGEBAK", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00191 infoc_1.ok);
00192 infoc_1.infot = 4;
00193 dgebak_("N", "R", &c__0, &c__0, &c__0, s, &c__0, a, &c__1, &info);
00194 chkxer_("DGEBAK", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00195 infoc_1.ok);
00196 infoc_1.infot = 4;
00197 dgebak_("N", "R", &c__0, &c__2, &c__0, s, &c__0, a, &c__1, &info);
00198 chkxer_("DGEBAK", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00199 infoc_1.ok);
00200 infoc_1.infot = 5;
00201 dgebak_("N", "R", &c__2, &c__2, &c__1, s, &c__0, a, &c__2, &info);
00202 chkxer_("DGEBAK", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00203 infoc_1.ok);
00204 infoc_1.infot = 5;
00205 dgebak_("N", "R", &c__0, &c__1, &c__1, s, &c__0, a, &c__1, &info);
00206 chkxer_("DGEBAK", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00207 infoc_1.ok);
00208 infoc_1.infot = 7;
00209 dgebak_("N", "R", &c__0, &c__1, &c__0, s, &c_n1, a, &c__1, &info);
00210 chkxer_("DGEBAK", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00211 infoc_1.ok);
00212 infoc_1.infot = 9;
00213 dgebak_("N", "R", &c__2, &c__1, &c__2, s, &c__0, a, &c__1, &info);
00214 chkxer_("DGEBAK", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00215 infoc_1.ok);
00216 nt += 9;
00217
00218
00219
00220 s_copy(srnamc_1.srnamt, "DGEHRD", (ftnlen)32, (ftnlen)6);
00221 infoc_1.infot = 1;
00222 dgehrd_(&c_n1, &c__1, &c__1, a, &c__1, tau, w, &c__1, &info);
00223 chkxer_("DGEHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00224 infoc_1.ok);
00225 infoc_1.infot = 2;
00226 dgehrd_(&c__0, &c__0, &c__0, a, &c__1, tau, w, &c__1, &info);
00227 chkxer_("DGEHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00228 infoc_1.ok);
00229 infoc_1.infot = 2;
00230 dgehrd_(&c__0, &c__2, &c__0, a, &c__1, tau, w, &c__1, &info);
00231 chkxer_("DGEHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00232 infoc_1.ok);
00233 infoc_1.infot = 3;
00234 dgehrd_(&c__1, &c__1, &c__0, a, &c__1, tau, w, &c__1, &info);
00235 chkxer_("DGEHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00236 infoc_1.ok);
00237 infoc_1.infot = 3;
00238 dgehrd_(&c__0, &c__1, &c__1, a, &c__1, tau, w, &c__1, &info);
00239 chkxer_("DGEHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00240 infoc_1.ok);
00241 infoc_1.infot = 5;
00242 dgehrd_(&c__2, &c__1, &c__1, a, &c__1, tau, w, &c__2, &info);
00243 chkxer_("DGEHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00244 infoc_1.ok);
00245 infoc_1.infot = 8;
00246 dgehrd_(&c__2, &c__1, &c__2, a, &c__2, tau, w, &c__1, &info);
00247 chkxer_("DGEHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00248 infoc_1.ok);
00249 nt += 7;
00250
00251
00252
00253 s_copy(srnamc_1.srnamt, "DORGHR", (ftnlen)32, (ftnlen)6);
00254 infoc_1.infot = 1;
00255 dorghr_(&c_n1, &c__1, &c__1, a, &c__1, tau, w, &c__1, &info);
00256 chkxer_("DORGHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00257 infoc_1.ok);
00258 infoc_1.infot = 2;
00259 dorghr_(&c__0, &c__0, &c__0, a, &c__1, tau, w, &c__1, &info);
00260 chkxer_("DORGHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00261 infoc_1.ok);
00262 infoc_1.infot = 2;
00263 dorghr_(&c__0, &c__2, &c__0, a, &c__1, tau, w, &c__1, &info);
00264 chkxer_("DORGHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00265 infoc_1.ok);
00266 infoc_1.infot = 3;
00267 dorghr_(&c__1, &c__1, &c__0, a, &c__1, tau, w, &c__1, &info);
00268 chkxer_("DORGHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00269 infoc_1.ok);
00270 infoc_1.infot = 3;
00271 dorghr_(&c__0, &c__1, &c__1, a, &c__1, tau, w, &c__1, &info);
00272 chkxer_("DORGHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00273 infoc_1.ok);
00274 infoc_1.infot = 5;
00275 dorghr_(&c__2, &c__1, &c__1, a, &c__1, tau, w, &c__1, &info);
00276 chkxer_("DORGHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00277 infoc_1.ok);
00278 infoc_1.infot = 8;
00279 dorghr_(&c__3, &c__1, &c__3, a, &c__3, tau, w, &c__1, &info);
00280 chkxer_("DORGHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00281 infoc_1.ok);
00282 nt += 7;
00283
00284
00285
00286 s_copy(srnamc_1.srnamt, "DORMHR", (ftnlen)32, (ftnlen)6);
00287 infoc_1.infot = 1;
00288 dormhr_("/", "N", &c__0, &c__0, &c__1, &c__0, a, &c__1, tau, c__, &
00289 c__1, w, &c__1, &info);
00290 chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00291 infoc_1.ok);
00292 infoc_1.infot = 2;
00293 dormhr_("L", "/", &c__0, &c__0, &c__1, &c__0, a, &c__1, tau, c__, &
00294 c__1, w, &c__1, &info);
00295 chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00296 infoc_1.ok);
00297 infoc_1.infot = 3;
00298 dormhr_("L", "N", &c_n1, &c__0, &c__1, &c__0, a, &c__1, tau, c__, &
00299 c__1, w, &c__1, &info);
00300 chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00301 infoc_1.ok);
00302 infoc_1.infot = 4;
00303 dormhr_("L", "N", &c__0, &c_n1, &c__1, &c__0, a, &c__1, tau, c__, &
00304 c__1, w, &c__1, &info);
00305 chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00306 infoc_1.ok);
00307 infoc_1.infot = 5;
00308 dormhr_("L", "N", &c__0, &c__0, &c__0, &c__0, a, &c__1, tau, c__, &
00309 c__1, w, &c__1, &info);
00310 chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00311 infoc_1.ok);
00312 infoc_1.infot = 5;
00313 dormhr_("L", "N", &c__0, &c__0, &c__2, &c__0, a, &c__1, tau, c__, &
00314 c__1, w, &c__1, &info);
00315 chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00316 infoc_1.ok);
00317 infoc_1.infot = 5;
00318 dormhr_("L", "N", &c__1, &c__2, &c__2, &c__1, a, &c__1, tau, c__, &
00319 c__1, w, &c__2, &info);
00320 chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00321 infoc_1.ok);
00322 infoc_1.infot = 5;
00323 dormhr_("R", "N", &c__2, &c__1, &c__2, &c__1, a, &c__1, tau, c__, &
00324 c__2, w, &c__2, &info);
00325 chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00326 infoc_1.ok);
00327 infoc_1.infot = 6;
00328 dormhr_("L", "N", &c__1, &c__1, &c__1, &c__0, a, &c__1, tau, c__, &
00329 c__1, w, &c__1, &info);
00330 chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00331 infoc_1.ok);
00332 infoc_1.infot = 6;
00333 dormhr_("L", "N", &c__0, &c__1, &c__1, &c__1, a, &c__1, tau, c__, &
00334 c__1, w, &c__1, &info);
00335 chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00336 infoc_1.ok);
00337 infoc_1.infot = 6;
00338 dormhr_("R", "N", &c__1, &c__0, &c__1, &c__1, a, &c__1, tau, c__, &
00339 c__1, w, &c__1, &info);
00340 chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00341 infoc_1.ok);
00342 infoc_1.infot = 8;
00343 dormhr_("L", "N", &c__2, &c__1, &c__1, &c__1, a, &c__1, tau, c__, &
00344 c__2, w, &c__1, &info);
00345 chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00346 infoc_1.ok);
00347 infoc_1.infot = 8;
00348 dormhr_("R", "N", &c__1, &c__2, &c__1, &c__1, a, &c__1, tau, c__, &
00349 c__1, w, &c__1, &info);
00350 chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00351 infoc_1.ok);
00352 infoc_1.infot = 11;
00353 dormhr_("L", "N", &c__2, &c__1, &c__1, &c__1, a, &c__2, tau, c__, &
00354 c__1, w, &c__1, &info);
00355 chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00356 infoc_1.ok);
00357 infoc_1.infot = 13;
00358 dormhr_("L", "N", &c__1, &c__2, &c__1, &c__1, a, &c__1, tau, c__, &
00359 c__1, w, &c__1, &info);
00360 chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00361 infoc_1.ok);
00362 infoc_1.infot = 13;
00363 dormhr_("R", "N", &c__2, &c__1, &c__1, &c__1, a, &c__1, tau, c__, &
00364 c__2, w, &c__1, &info);
00365 chkxer_("DORMHR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00366 infoc_1.ok);
00367 nt += 16;
00368
00369
00370
00371 s_copy(srnamc_1.srnamt, "DHSEQR", (ftnlen)32, (ftnlen)6);
00372 infoc_1.infot = 1;
00373 dhseqr_("/", "N", &c__0, &c__1, &c__0, a, &c__1, wr, wi, c__, &c__1,
00374 w, &c__1, &info);
00375 chkxer_("DHSEQR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00376 infoc_1.ok);
00377 infoc_1.infot = 2;
00378 dhseqr_("E", "/", &c__0, &c__1, &c__0, a, &c__1, wr, wi, c__, &c__1,
00379 w, &c__1, &info);
00380 chkxer_("DHSEQR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00381 infoc_1.ok);
00382 infoc_1.infot = 3;
00383 dhseqr_("E", "N", &c_n1, &c__1, &c__0, a, &c__1, wr, wi, c__, &c__1,
00384 w, &c__1, &info);
00385 chkxer_("DHSEQR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00386 infoc_1.ok);
00387 infoc_1.infot = 4;
00388 dhseqr_("E", "N", &c__0, &c__0, &c__0, a, &c__1, wr, wi, c__, &c__1,
00389 w, &c__1, &info);
00390 chkxer_("DHSEQR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00391 infoc_1.ok);
00392 infoc_1.infot = 4;
00393 dhseqr_("E", "N", &c__0, &c__2, &c__0, a, &c__1, wr, wi, c__, &c__1,
00394 w, &c__1, &info);
00395 chkxer_("DHSEQR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00396 infoc_1.ok);
00397 infoc_1.infot = 5;
00398 dhseqr_("E", "N", &c__1, &c__1, &c__0, a, &c__1, wr, wi, c__, &c__1,
00399 w, &c__1, &info);
00400 chkxer_("DHSEQR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00401 infoc_1.ok);
00402 infoc_1.infot = 5;
00403 dhseqr_("E", "N", &c__1, &c__1, &c__2, a, &c__1, wr, wi, c__, &c__1,
00404 w, &c__1, &info);
00405 chkxer_("DHSEQR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00406 infoc_1.ok);
00407 infoc_1.infot = 7;
00408 dhseqr_("E", "N", &c__2, &c__1, &c__2, a, &c__1, wr, wi, c__, &c__2,
00409 w, &c__1, &info);
00410 chkxer_("DHSEQR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00411 infoc_1.ok);
00412 infoc_1.infot = 11;
00413 dhseqr_("E", "V", &c__2, &c__1, &c__2, a, &c__2, wr, wi, c__, &c__1,
00414 w, &c__1, &info);
00415 chkxer_("DHSEQR", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00416 infoc_1.ok);
00417 nt += 9;
00418
00419
00420
00421 s_copy(srnamc_1.srnamt, "DHSEIN", (ftnlen)32, (ftnlen)6);
00422 infoc_1.infot = 1;
00423 dhsein_("/", "N", "N", sel, &c__0, a, &c__1, wr, wi, vl, &c__1, vr, &
00424 c__1, &c__0, &m, w, ifaill, ifailr, &info);
00425 chkxer_("DHSEIN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00426 infoc_1.ok);
00427 infoc_1.infot = 2;
00428 dhsein_("R", "/", "N", sel, &c__0, a, &c__1, wr, wi, vl, &c__1, vr, &
00429 c__1, &c__0, &m, w, ifaill, ifailr, &info);
00430 chkxer_("DHSEIN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00431 infoc_1.ok);
00432 infoc_1.infot = 3;
00433 dhsein_("R", "N", "/", sel, &c__0, a, &c__1, wr, wi, vl, &c__1, vr, &
00434 c__1, &c__0, &m, w, ifaill, ifailr, &info);
00435 chkxer_("DHSEIN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00436 infoc_1.ok);
00437 infoc_1.infot = 5;
00438 dhsein_("R", "N", "N", sel, &c_n1, a, &c__1, wr, wi, vl, &c__1, vr, &
00439 c__1, &c__0, &m, w, ifaill, ifailr, &info);
00440 chkxer_("DHSEIN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00441 infoc_1.ok);
00442 infoc_1.infot = 7;
00443 dhsein_("R", "N", "N", sel, &c__2, a, &c__1, wr, wi, vl, &c__1, vr, &
00444 c__2, &c__4, &m, w, ifaill, ifailr, &info);
00445 chkxer_("DHSEIN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00446 infoc_1.ok);
00447 infoc_1.infot = 11;
00448 dhsein_("L", "N", "N", sel, &c__2, a, &c__2, wr, wi, vl, &c__1, vr, &
00449 c__1, &c__4, &m, w, ifaill, ifailr, &info);
00450 chkxer_("DHSEIN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00451 infoc_1.ok);
00452 infoc_1.infot = 13;
00453 dhsein_("R", "N", "N", sel, &c__2, a, &c__2, wr, wi, vl, &c__1, vr, &
00454 c__1, &c__4, &m, w, ifaill, ifailr, &info);
00455 chkxer_("DHSEIN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00456 infoc_1.ok);
00457 infoc_1.infot = 14;
00458 dhsein_("R", "N", "N", sel, &c__2, a, &c__2, wr, wi, vl, &c__1, vr, &
00459 c__2, &c__1, &m, w, ifaill, ifailr, &info);
00460 chkxer_("DHSEIN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00461 infoc_1.ok);
00462 nt += 8;
00463
00464
00465
00466 s_copy(srnamc_1.srnamt, "DTREVC", (ftnlen)32, (ftnlen)6);
00467 infoc_1.infot = 1;
00468 dtrevc_("/", "A", sel, &c__0, a, &c__1, vl, &c__1, vr, &c__1, &c__0, &
00469 m, w, &info);
00470 chkxer_("DTREVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00471 infoc_1.ok);
00472 infoc_1.infot = 2;
00473 dtrevc_("L", "/", sel, &c__0, a, &c__1, vl, &c__1, vr, &c__1, &c__0, &
00474 m, w, &info);
00475 chkxer_("DTREVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00476 infoc_1.ok);
00477 infoc_1.infot = 4;
00478 dtrevc_("L", "A", sel, &c_n1, a, &c__1, vl, &c__1, vr, &c__1, &c__0, &
00479 m, w, &info);
00480 chkxer_("DTREVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00481 infoc_1.ok);
00482 infoc_1.infot = 6;
00483 dtrevc_("L", "A", sel, &c__2, a, &c__1, vl, &c__2, vr, &c__1, &c__4, &
00484 m, w, &info);
00485 chkxer_("DTREVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00486 infoc_1.ok);
00487 infoc_1.infot = 8;
00488 dtrevc_("L", "A", sel, &c__2, a, &c__2, vl, &c__1, vr, &c__1, &c__4, &
00489 m, w, &info);
00490 chkxer_("DTREVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00491 infoc_1.ok);
00492 infoc_1.infot = 10;
00493 dtrevc_("R", "A", sel, &c__2, a, &c__2, vl, &c__1, vr, &c__1, &c__4, &
00494 m, w, &info);
00495 chkxer_("DTREVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00496 infoc_1.ok);
00497 infoc_1.infot = 11;
00498 dtrevc_("L", "A", sel, &c__2, a, &c__2, vl, &c__2, vr, &c__1, &c__1, &
00499 m, w, &info);
00500 chkxer_("DTREVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00501 infoc_1.ok);
00502 nt += 7;
00503 }
00504
00505
00506
00507 if (infoc_1.ok) {
00508 io___22.ciunit = infoc_1.nout;
00509 s_wsfe(&io___22);
00510 do_fio(&c__1, path, (ftnlen)3);
00511 do_fio(&c__1, (char *)&nt, (ftnlen)sizeof(integer));
00512 e_wsfe();
00513 } else {
00514 io___23.ciunit = infoc_1.nout;
00515 s_wsfe(&io___23);
00516 do_fio(&c__1, path, (ftnlen)3);
00517 e_wsfe();
00518 }
00519
00520
00521 return 0;
00522
00523
00524
00525 }