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