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 nproc, nshift, maxb;
00020 } cenvir_;
00021
00022 #define cenvir_1 cenvir_
00023
00024 struct {
00025 integer infot, nunit;
00026 logical ok, lerr;
00027 } infoc_;
00028
00029 #define infoc_1 infoc_
00030
00031 struct {
00032 char srnamt[32];
00033 } srnamc_;
00034
00035 #define srnamc_1 srnamc_
00036
00037 struct {
00038 integer selopt, seldim;
00039 logical selval[20];
00040 doublereal selwr[20], selwi[20];
00041 } sslct_;
00042
00043 #define sslct_1 sslct_
00044
00045 struct {
00046 integer iparms[100];
00047 } zlaenv_;
00048
00049 #define zlaenv_1 zlaenv_
00050
00051
00052
00053 static integer c__1 = 1;
00054 static integer c__3 = 3;
00055 static integer c__5 = 5;
00056 static integer c__6 = 6;
00057 static integer c__20 = 20;
00058 static integer c__0 = 0;
00059 static integer c__132 = 132;
00060 static integer c__2 = 2;
00061 static integer c__12 = 12;
00062 static integer c__13 = 13;
00063 static integer c__14 = 14;
00064 static integer c__15 = 15;
00065 static integer c__16 = 16;
00066 static integer c__4 = 4;
00067 static integer c__8 = 8;
00068 static integer c__89760 = 89760;
00069 static integer c__9 = 9;
00070 static integer c__25 = 25;
00071 static integer c__20064 = 20064;
00072 static integer c__18 = 18;
00073 static integer c__400 = 400;
00074 static integer c__20062 = 20062;
00075 static integer c__264 = 264;
00076
00077 int MAIN__(void)
00078 {
00079
00080
00081 static char intstr[10] = "0123456789";
00082 static integer ioldsd[4] = { 0,0,0,1 };
00083
00084
00085 static char fmt_9987[] = "(\002 Tests of the Nonsymmetric Eigenvalue Pro"
00086 "blem routines\002)";
00087 static char fmt_9986[] = "(\002 Tests of the Hermitian Eigenvalue Proble"
00088 "m routines\002)";
00089 static char fmt_9985[] = "(\002 Tests of the Singular Value Decompositio"
00090 "n routines\002)";
00091 static char fmt_9979[] = "(/\002 Tests of the Nonsymmetric Eigenvalue Pr"
00092 "oblem Driver\002,/\002 ZGEEV (eigenvalues and eigevectors)"
00093 "\002)";
00094 static char fmt_9978[] = "(/\002 Tests of the Nonsymmetric Eigenvalue Pr"
00095 "oblem Driver\002,/\002 ZGEES (Schur form)\002)";
00096 static char fmt_9977[] = "(/\002 Tests of the Nonsymmetric Eigenvalue Pr"
00097 "oblem Expert\002,\002 Driver\002,/\002 ZGEEVX (eigenvalues, e"
00098 "igenvectors and\002,\002 condition numbers)\002)";
00099 static char fmt_9976[] = "(/\002 Tests of the Nonsymmetric Eigenvalue Pr"
00100 "oblem Expert\002,\002 Driver\002,/\002 ZGEESX (Schur form and"
00101 " condition\002,\002 numbers)\002)";
00102 static char fmt_9975[] = "(/\002 Tests of the Generalized Nonsymmetric E"
00103 "igenvalue \002,\002Problem routines\002)";
00104 static char fmt_9964[] = "(/\002 Tests of the Generalized Nonsymmetric E"
00105 "igenvalue \002,\002Problem Driver ZGGES\002)";
00106 static char fmt_9965[] = "(/\002 Tests of the Generalized Nonsymmetric E"
00107 "igenvalue \002,\002Problem Expert Driver ZGGESX\002)";
00108 static char fmt_9963[] = "(/\002 Tests of the Generalized Nonsymmetric E"
00109 "igenvalue \002,\002Problem Driver ZGGEV\002)";
00110 static char fmt_9962[] = "(/\002 Tests of the Generalized Nonsymmetric E"
00111 "igenvalue \002,\002Problem Expert Driver ZGGEVX\002)";
00112 static char fmt_9974[] = "(\002 Tests of ZHBTRD\002,/\002 (reduction of "
00113 "a Hermitian band \002,\002matrix to real tridiagonal form)\002)";
00114 static char fmt_9967[] = "(\002 Tests of ZGBBRD\002,/\002 (reduction of "
00115 "a general band \002,\002matrix to real bidiagonal form)\002)";
00116 static char fmt_9971[] = "(/\002 Tests of the Generalized Linear Regress"
00117 "ion Model \002,\002routines\002)";
00118 static char fmt_9970[] = "(/\002 Tests of the Generalized QR and RQ rout"
00119 "ines\002)";
00120 static char fmt_9969[] = "(/\002 Tests of the Generalized Singular Valu"
00121 "e\002,\002 Decomposition routines\002)";
00122 static char fmt_9968[] = "(/\002 Tests of the Linear Least Squares routi"
00123 "nes\002)";
00124 static char fmt_9992[] = "(1x,a3,\002: Unrecognized path name\002)";
00125 static char fmt_9972[] = "(/\002 LAPACK VERSION \002,i1,\002.\002,i1,"
00126 "\002.\002,i1)";
00127 static char fmt_9984[] = "(/\002 The following parameter values will be "
00128 "used:\002)";
00129 static char fmt_9989[] = "(\002 Invalid input value: \002,a,\002=\002,"
00130 "i6,\002; must be >=\002,i6)";
00131 static char fmt_9988[] = "(\002 Invalid input value: \002,a,\002=\002,"
00132 "i6,\002; must be <=\002,i6)";
00133 static char fmt_9983[] = "(4x,a,10i6,/10x,10i6)";
00134 static char fmt_9981[] = "(\002 Relative machine \002,a,\002 is taken to"
00135 " be\002,d16.6)";
00136 static char fmt_9982[] = "(/\002 Routines pass computational tests if te"
00137 "st ratio is \002,\002less than\002,f8.2,/)";
00138 static char fmt_9999[] = "(/\002 Execution not attempted due to input er"
00139 "rors\002)";
00140 static char fmt_9991[] = "(//\002 *** Invalid integer value in column"
00141 " \002,i2,\002 of input\002,\002 line:\002,/a79)";
00142 static char fmt_9990[] = "(//1x,a3,\002 routines were not tested\002)";
00143 static char fmt_9961[] = "(//1x,a3,\002: NB =\002,i4,\002, NBMIN =\002,"
00144 "i4,\002, NX =\002,i4,\002, INMIN=\002,i4,\002, INWIN =\002,i4"
00145 ",\002, INIBL =\002,i4,\002, ISHFTS =\002,i4,\002, IACC22 =\002,i"
00146 "4)";
00147 static char fmt_9980[] = "(\002 *** Error code from \002,a,\002 = \002,i"
00148 "4)";
00149 static char fmt_9997[] = "(//1x,a3,\002: NB =\002,i4,\002, NBMIN =\002,"
00150 "i4,\002, NX =\002,i4)";
00151 static char fmt_9995[] = "(//1x,a3,\002: NB =\002,i4,\002, NBMIN =\002,"
00152 "i4,\002, NX =\002,i4,\002, NRHS =\002,i4)";
00153 static char fmt_9973[] = "(/1x,71(\002-\002))";
00154 static char fmt_9996[] = "(//1x,a3,\002: NB =\002,i4,\002, NBMIN =\002,"
00155 "i4,\002, NS =\002,i4,\002, MAXB =\002,i4,\002, NBCOL =\002,i4)";
00156 static char fmt_9966[] = "(//1x,a3,\002: NRHS =\002,i4)";
00157 static char fmt_9994[] = "(//\002 End of tests\002)";
00158 static char fmt_9993[] = "(\002 Total time used = \002,f12.2,\002 seco"
00159 "nds\002,/)";
00160
00161
00162 integer i__1, i__2, i__3, i__4;
00163 doublereal d__1;
00164 cilist ci__1;
00165
00166
00167 integer s_rsfe(cilist *), do_fio(integer *, char *, ftnlen), e_rsfe(void);
00168 int s_copy(char *, char *, ftnlen, ftnlen);
00169 integer s_cmp(char *, char *, ftnlen, ftnlen), s_wsfe(cilist *), e_wsfe(
00170 void), s_rsle(cilist *), do_lio(integer *, integer *, char *,
00171 ftnlen), e_rsle(void), s_wsle(cilist *), e_wsle(void);
00172 int s_stop(char *, ftnlen);
00173 integer i_len(char *, ftnlen);
00174
00175
00176 doublecomplex a[243936] , b[87120]
00177 , c__[160000] ;
00178 integer i__, k;
00179 doublereal s[17424];
00180 doublecomplex x[660];
00181 char c1[1], c3[3];
00182 integer i1;
00183 doublereal s1, s2;
00184 doublecomplex dc[792] ;
00185 integer ic;
00186 doublereal dr[1584] ;
00187 integer nk, nn, vers_patch__, vers_major__, vers_minor__;
00188 logical zbb, glm, nep, lse, zbk, zbl, sep, gqr, zgg, zgk, zgl, svd, zhb,
00189 gsv;
00190 doublereal eps;
00191 logical zes, zgs, zev, zgv, zgx, zsx, zvx, zxv;
00192 doublereal beta[132];
00193 char line[80];
00194 doublecomplex taua[132];
00195 integer info;
00196 char path[3];
00197 integer kval[20], lenp, mval[20], nval[20];
00198 doublecomplex taub[132];
00199 integer pval[20], itmp, nrhs;
00200 doublecomplex work[89760];
00201 integer iacc22[20];
00202 doublereal alpha[132];
00203 logical fatal;
00204 integer iseed[4], nbcol[20], inibl[20], nbval[20], nbmin[20];
00205 char vname[32];
00206 integer inmin[20], newsd, nsval[20], inwin[20], nxval[20], iwork[20064];
00207 doublereal rwork[89760];
00208 extern doublereal dlamch_(char *), dsecnd_(void);
00209 extern int zchkbb_(integer *, integer *, integer *,
00210 integer *, integer *, integer *, logical *, integer *, integer *,
00211 doublereal *, integer *, doublecomplex *, integer *,
00212 doublecomplex *, integer *, doublereal *, doublereal *,
00213 doublecomplex *, integer *, doublecomplex *, integer *,
00214 doublecomplex *, integer *, doublecomplex *, doublecomplex *,
00215 integer *, doublereal *, doublereal *, integer *), alareq_(char *,
00216 integer *, logical *, integer *, integer *, integer *),
00217 zchkbd_(integer *, integer *, integer *, integer *, logical *,
00218 integer *, integer *, doublereal *, doublecomplex *, integer *,
00219 doublereal *, doublereal *, doublereal *, doublereal *,
00220 doublecomplex *, integer *, doublecomplex *, doublecomplex *,
00221 doublecomplex *, integer *, doublecomplex *, integer *,
00222 doublecomplex *, doublecomplex *, doublecomplex *, integer *,
00223 doublereal *, integer *, integer *), zchkec_(doublereal *,
00224 logical *, integer *, integer *), zchkhb_(integer *, integer *,
00225 integer *, integer *, integer *, logical *, integer *, doublereal
00226 *, integer *, doublecomplex *, integer *, doublereal *,
00227 doublereal *, doublecomplex *, integer *, doublecomplex *,
00228 integer *, doublereal *, doublereal *, integer *), zchkbk_(
00229 integer *, integer *), zchkbl_(integer *, integer *);
00230 extern logical lsamen_(integer *, char *, char *);
00231 extern int zchkgg_(integer *, integer *, integer *,
00232 logical *, integer *, doublereal *, logical *, doublereal *,
00233 integer *, doublecomplex *, integer *, doublecomplex *,
00234 doublecomplex *, doublecomplex *, doublecomplex *, doublecomplex *
00235 , doublecomplex *, doublecomplex *, doublecomplex *, integer *,
00236 doublecomplex *, doublecomplex *, doublecomplex *, doublecomplex *
00237 , doublecomplex *, doublecomplex *, doublecomplex *,
00238 doublecomplex *, doublecomplex *, doublecomplex *, integer *,
00239 doublereal *, logical *, doublereal *, integer *), zchkgk_(
00240 integer *, integer *), zchkgl_(integer *, integer *), ilaver_(
00241 integer *, integer *, integer *), zckglm_(integer *, integer *,
00242 integer *, integer *, integer *, integer *, doublereal *, integer
00243 *, doublecomplex *, doublecomplex *, doublecomplex *,
00244 doublecomplex *, doublecomplex *, doublecomplex *, doublereal *,
00245 integer *, integer *, integer *), zerrbd_(char *, integer *);
00246 integer mxbval[20];
00247 extern int zchkhs_(integer *, integer *, integer *,
00248 logical *, integer *, doublereal *, integer *, doublecomplex *,
00249 integer *, doublecomplex *, doublecomplex *, doublecomplex *,
00250 doublecomplex *, integer *, doublecomplex *, doublecomplex *,
00251 doublecomplex *, doublecomplex *, doublecomplex *, doublecomplex *
00252 , doublecomplex *, doublecomplex *, doublecomplex *,
00253 doublecomplex *, doublecomplex *, integer *, doublereal *,
00254 integer *, logical *, doublereal *, integer *), zcklse_(integer *,
00255 integer *, integer *, integer *, integer *, integer *,
00256 doublereal *, integer *, doublecomplex *, doublecomplex *,
00257 doublecomplex *, doublecomplex *, doublecomplex *, doublecomplex *
00258 , doublereal *, integer *, integer *, integer *), zdrvbd_(integer
00259 *, integer *, integer *, integer *, logical *, integer *,
00260 doublereal *, doublecomplex *, integer *, doublecomplex *,
00261 integer *, doublecomplex *, integer *, doublecomplex *,
00262 doublecomplex *, doublecomplex *, doublereal *, doublereal *,
00263 doublereal *, doublecomplex *, integer *, doublereal *, integer *,
00264 integer *, integer *);
00265 logical tstdif;
00266 doublereal thresh;
00267 extern int xlaenv_(integer *, integer *);
00268 logical tstchk;
00269 integer nparms, ishfts[20];
00270 extern int zckgqr_(integer *, integer *, integer *,
00271 integer *, integer *, integer *, integer *, integer *, doublereal
00272 *, integer *, doublecomplex *, doublecomplex *, doublecomplex *,
00273 doublecomplex *, doublecomplex *, doublecomplex *, doublecomplex *
00274 , doublecomplex *, doublecomplex *, doublecomplex *,
00275 doublecomplex *, doublecomplex *, doublereal *, integer *,
00276 integer *, integer *);
00277 logical dotype[30], logwrk[132];
00278 doublereal thrshn;
00279 extern int zchkst_(integer *, integer *, integer *,
00280 logical *, integer *, doublereal *, integer *, doublecomplex *,
00281 integer *, doublecomplex *, doublereal *, doublereal *,
00282 doublereal *, doublereal *, doublereal *, doublereal *,
00283 doublereal *, doublereal *, doublereal *, doublereal *,
00284 doublereal *, doublecomplex *, integer *, doublecomplex *,
00285 doublecomplex *, doublecomplex *, doublecomplex *, doublecomplex *
00286 , integer *, doublereal *, integer *, integer *, integer *,
00287 doublereal *, integer *), zckgsv_(integer *, integer *, integer *,
00288 integer *, integer *, integer *, doublereal *, integer *,
00289 doublecomplex *, doublecomplex *, doublecomplex *, doublecomplex *
00290 , doublecomplex *, doublecomplex *, doublecomplex *, doublereal *,
00291 doublereal *, doublecomplex *, integer *, doublecomplex *,
00292 doublereal *, integer *, integer *, integer *), zdrges_(integer *,
00293 integer *, integer *, logical *, integer *, doublereal *,
00294 integer *, doublecomplex *, integer *, doublecomplex *,
00295 doublecomplex *, doublecomplex *, doublecomplex *, integer *,
00296 doublecomplex *, doublecomplex *, doublecomplex *, doublecomplex *
00297 , integer *, doublereal *, doublereal *, logical *, integer *),
00298 zdrgev_(integer *, integer *, integer *, logical *, integer *,
00299 doublereal *, integer *, doublecomplex *, integer *,
00300 doublecomplex *, doublecomplex *, doublecomplex *, doublecomplex *
00301 , integer *, doublecomplex *, doublecomplex *, integer *,
00302 doublecomplex *, doublecomplex *, doublecomplex *, doublecomplex *
00303 , doublecomplex *, integer *, doublereal *, doublereal *, integer
00304 *), zdrvgg_(integer *, integer *, integer *, logical *, integer *,
00305 doublereal *, doublereal *, integer *, doublecomplex *, integer *
00306 , doublecomplex *, doublecomplex *, doublecomplex *,
00307 doublecomplex *, doublecomplex *, doublecomplex *, integer *,
00308 doublecomplex *, doublecomplex *, doublecomplex *, doublecomplex *
00309 , doublecomplex *, doublecomplex *, doublecomplex *,
00310 doublecomplex *, integer *, doublereal *, doublereal *, integer *)
00311 , zdrves_(integer *, integer *, integer *, logical *, integer *,
00312 doublereal *, integer *, doublecomplex *, integer *,
00313 doublecomplex *, doublecomplex *, doublecomplex *, doublecomplex *
00314 , doublecomplex *, integer *, doublereal *, doublecomplex *,
00315 integer *, doublereal *, integer *, logical *, integer *);
00316 doublereal result[500];
00317 extern int zdrvsg_(integer *, integer *, integer *,
00318 logical *, integer *, doublereal *, integer *, doublecomplex *,
00319 integer *, doublecomplex *, integer *, doublereal *,
00320 doublecomplex *, integer *, doublecomplex *, doublecomplex *,
00321 doublecomplex *, doublecomplex *, doublecomplex *, integer *,
00322 doublereal *, integer *, integer *, integer *, doublereal *,
00323 integer *), zdrvev_(integer *, integer *, integer *, logical *,
00324 integer *, doublereal *, integer *, doublecomplex *, integer *,
00325 doublecomplex *, doublecomplex *, doublecomplex *, doublecomplex *
00326 , integer *, doublecomplex *, integer *, doublecomplex *, integer
00327 *, doublereal *, doublecomplex *, integer *, doublereal *,
00328 integer *, integer *), zdrgsx_(integer *, integer *, doublereal *,
00329 integer *, integer *, doublecomplex *, integer *, doublecomplex *
00330 , doublecomplex *, doublecomplex *, doublecomplex *,
00331 doublecomplex *, doublecomplex *, doublecomplex *, doublecomplex *
00332 , integer *, doublereal *, doublecomplex *, integer *, doublereal
00333 *, integer *, integer *, logical *, integer *);
00334 integer maxtyp;
00335 logical tsterr;
00336 integer ntypes;
00337 logical tstdrv;
00338 extern int zdrgvx_(integer *, doublereal *, integer *,
00339 integer *, doublecomplex *, integer *, doublecomplex *,
00340 doublecomplex *, doublecomplex *, doublecomplex *, doublecomplex *
00341 , doublecomplex *, doublecomplex *, integer *, integer *,
00342 doublereal *, doublereal *, doublereal *, doublereal *,
00343 doublereal *, doublereal *, doublecomplex *, integer *,
00344 doublereal *, integer *, integer *, doublereal *, logical *,
00345 integer *), zerred_(char *, integer *), zerrgg_(char *,
00346 integer *), zerrhs_(char *, integer *), zerrst_(
00347 char *, integer *), zdrvst_(integer *, integer *, integer
00348 *, logical *, integer *, doublereal *, integer *, doublecomplex *,
00349 integer *, doublereal *, doublereal *, doublereal *, doublereal *
00350 , doublereal *, doublereal *, doublecomplex *, integer *,
00351 doublecomplex *, doublecomplex *, doublecomplex *, doublecomplex *
00352 , integer *, doublereal *, integer *, integer *, integer *,
00353 doublereal *, integer *), zdrvsx_(integer *, integer *, integer *,
00354 logical *, integer *, doublereal *, integer *, integer *,
00355 doublecomplex *, integer *, doublecomplex *, doublecomplex *,
00356 doublecomplex *, doublecomplex *, doublecomplex *, doublecomplex *
00357 , integer *, doublecomplex *, doublereal *, doublecomplex *,
00358 integer *, doublereal *, logical *, integer *), zdrvvx_(integer *,
00359 integer *, integer *, logical *, integer *, doublereal *,
00360 integer *, integer *, doublecomplex *, integer *, doublecomplex *,
00361 doublecomplex *, doublecomplex *, doublecomplex *, integer *,
00362 doublecomplex *, integer *, doublecomplex *, integer *,
00363 doublereal *, doublereal *, doublereal *, doublereal *,
00364 doublereal *, doublereal *, doublereal *, doublereal *,
00365 doublereal *, doublecomplex *, integer *, doublereal *, integer *)
00366 ;
00367
00368
00369 static cilist io___29 = { 0, 6, 0, fmt_9987, 0 };
00370 static cilist io___30 = { 0, 6, 0, fmt_9986, 0 };
00371 static cilist io___31 = { 0, 6, 0, fmt_9985, 0 };
00372 static cilist io___32 = { 0, 6, 0, fmt_9979, 0 };
00373 static cilist io___33 = { 0, 6, 0, fmt_9978, 0 };
00374 static cilist io___34 = { 0, 6, 0, fmt_9977, 0 };
00375 static cilist io___35 = { 0, 6, 0, fmt_9976, 0 };
00376 static cilist io___36 = { 0, 6, 0, fmt_9975, 0 };
00377 static cilist io___37 = { 0, 6, 0, fmt_9964, 0 };
00378 static cilist io___38 = { 0, 6, 0, fmt_9965, 0 };
00379 static cilist io___39 = { 0, 6, 0, fmt_9963, 0 };
00380 static cilist io___40 = { 0, 6, 0, fmt_9962, 0 };
00381 static cilist io___41 = { 0, 6, 0, fmt_9974, 0 };
00382 static cilist io___42 = { 0, 6, 0, fmt_9967, 0 };
00383 static cilist io___43 = { 0, 6, 0, fmt_9971, 0 };
00384 static cilist io___44 = { 0, 6, 0, fmt_9970, 0 };
00385 static cilist io___45 = { 0, 6, 0, fmt_9969, 0 };
00386 static cilist io___46 = { 0, 6, 0, fmt_9968, 0 };
00387 static cilist io___47 = { 0, 5, 0, 0, 0 };
00388 static cilist io___50 = { 0, 6, 0, fmt_9992, 0 };
00389 static cilist io___54 = { 0, 6, 0, fmt_9972, 0 };
00390 static cilist io___55 = { 0, 6, 0, fmt_9984, 0 };
00391 static cilist io___56 = { 0, 5, 0, 0, 0 };
00392 static cilist io___58 = { 0, 6, 0, fmt_9989, 0 };
00393 static cilist io___59 = { 0, 6, 0, fmt_9988, 0 };
00394 static cilist io___60 = { 0, 5, 0, 0, 0 };
00395 static cilist io___64 = { 0, 6, 0, fmt_9989, 0 };
00396 static cilist io___65 = { 0, 6, 0, fmt_9988, 0 };
00397 static cilist io___66 = { 0, 6, 0, fmt_9983, 0 };
00398 static cilist io___67 = { 0, 5, 0, 0, 0 };
00399 static cilist io___69 = { 0, 6, 0, fmt_9989, 0 };
00400 static cilist io___70 = { 0, 6, 0, fmt_9988, 0 };
00401 static cilist io___71 = { 0, 6, 0, fmt_9983, 0 };
00402 static cilist io___72 = { 0, 5, 0, 0, 0 };
00403 static cilist io___74 = { 0, 6, 0, fmt_9989, 0 };
00404 static cilist io___75 = { 0, 6, 0, fmt_9988, 0 };
00405 static cilist io___76 = { 0, 6, 0, fmt_9983, 0 };
00406 static cilist io___77 = { 0, 6, 0, fmt_9983, 0 };
00407 static cilist io___78 = { 0, 5, 0, 0, 0 };
00408 static cilist io___80 = { 0, 5, 0, 0, 0 };
00409 static cilist io___82 = { 0, 6, 0, fmt_9989, 0 };
00410 static cilist io___83 = { 0, 6, 0, fmt_9988, 0 };
00411 static cilist io___84 = { 0, 6, 0, fmt_9983, 0 };
00412 static cilist io___85 = { 0, 5, 0, 0, 0 };
00413 static cilist io___94 = { 0, 6, 0, fmt_9989, 0 };
00414 static cilist io___95 = { 0, 6, 0, fmt_9989, 0 };
00415 static cilist io___96 = { 0, 6, 0, fmt_9989, 0 };
00416 static cilist io___97 = { 0, 6, 0, fmt_9989, 0 };
00417 static cilist io___98 = { 0, 6, 0, fmt_9989, 0 };
00418 static cilist io___99 = { 0, 6, 0, fmt_9989, 0 };
00419 static cilist io___100 = { 0, 6, 0, fmt_9989, 0 };
00420 static cilist io___101 = { 0, 6, 0, fmt_9989, 0 };
00421 static cilist io___102 = { 0, 6, 0, fmt_9983, 0 };
00422 static cilist io___103 = { 0, 6, 0, fmt_9983, 0 };
00423 static cilist io___104 = { 0, 6, 0, fmt_9983, 0 };
00424 static cilist io___105 = { 0, 6, 0, fmt_9983, 0 };
00425 static cilist io___106 = { 0, 6, 0, fmt_9983, 0 };
00426 static cilist io___107 = { 0, 6, 0, fmt_9983, 0 };
00427 static cilist io___108 = { 0, 6, 0, fmt_9983, 0 };
00428 static cilist io___109 = { 0, 6, 0, fmt_9983, 0 };
00429 static cilist io___110 = { 0, 5, 0, 0, 0 };
00430 static cilist io___113 = { 0, 6, 0, fmt_9989, 0 };
00431 static cilist io___114 = { 0, 6, 0, fmt_9989, 0 };
00432 static cilist io___115 = { 0, 6, 0, fmt_9989, 0 };
00433 static cilist io___116 = { 0, 6, 0, fmt_9989, 0 };
00434 static cilist io___117 = { 0, 6, 0, fmt_9989, 0 };
00435 static cilist io___118 = { 0, 6, 0, fmt_9983, 0 };
00436 static cilist io___119 = { 0, 6, 0, fmt_9983, 0 };
00437 static cilist io___120 = { 0, 6, 0, fmt_9983, 0 };
00438 static cilist io___121 = { 0, 6, 0, fmt_9983, 0 };
00439 static cilist io___122 = { 0, 6, 0, fmt_9983, 0 };
00440 static cilist io___123 = { 0, 5, 0, 0, 0 };
00441 static cilist io___125 = { 0, 6, 0, fmt_9989, 0 };
00442 static cilist io___126 = { 0, 6, 0, fmt_9988, 0 };
00443 static cilist io___127 = { 0, 5, 0, 0, 0 };
00444 static cilist io___128 = { 0, 6, 0, fmt_9989, 0 };
00445 static cilist io___129 = { 0, 6, 0, fmt_9988, 0 };
00446 static cilist io___130 = { 0, 6, 0, fmt_9983, 0 };
00447 static cilist io___131 = { 0, 5, 0, 0, 0 };
00448 static cilist io___132 = { 0, 6, 0, fmt_9989, 0 };
00449 static cilist io___133 = { 0, 6, 0, fmt_9988, 0 };
00450 static cilist io___134 = { 0, 6, 0, fmt_9983, 0 };
00451 static cilist io___135 = { 0, 5, 0, 0, 0 };
00452 static cilist io___136 = { 0, 6, 0, fmt_9989, 0 };
00453 static cilist io___137 = { 0, 6, 0, fmt_9988, 0 };
00454 static cilist io___138 = { 0, 6, 0, fmt_9983, 0 };
00455 static cilist io___139 = { 0, 5, 0, 0, 0 };
00456 static cilist io___140 = { 0, 6, 0, fmt_9989, 0 };
00457 static cilist io___141 = { 0, 6, 0, fmt_9988, 0 };
00458 static cilist io___142 = { 0, 6, 0, fmt_9983, 0 };
00459 static cilist io___143 = { 0, 5, 0, 0, 0 };
00460 static cilist io___144 = { 0, 6, 0, fmt_9989, 0 };
00461 static cilist io___145 = { 0, 6, 0, fmt_9988, 0 };
00462 static cilist io___146 = { 0, 6, 0, fmt_9983, 0 };
00463 static cilist io___147 = { 0, 5, 0, 0, 0 };
00464 static cilist io___148 = { 0, 6, 0, fmt_9989, 0 };
00465 static cilist io___149 = { 0, 6, 0, fmt_9983, 0 };
00466 static cilist io___150 = { 0, 5, 0, 0, 0 };
00467 static cilist io___151 = { 0, 6, 0, fmt_9989, 0 };
00468 static cilist io___152 = { 0, 6, 0, fmt_9983, 0 };
00469 static cilist io___153 = { 0, 5, 0, 0, 0 };
00470 static cilist io___154 = { 0, 6, 0, fmt_9989, 0 };
00471 static cilist io___155 = { 0, 6, 0, fmt_9983, 0 };
00472 static cilist io___156 = { 0, 5, 0, 0, 0 };
00473 static cilist io___157 = { 0, 6, 0, fmt_9989, 0 };
00474 static cilist io___158 = { 0, 6, 0, fmt_9983, 0 };
00475 static cilist io___159 = { 0, 5, 0, 0, 0 };
00476 static cilist io___160 = { 0, 6, 0, fmt_9989, 0 };
00477 static cilist io___161 = { 0, 6, 0, fmt_9983, 0 };
00478 static cilist io___162 = { 0, 5, 0, 0, 0 };
00479 static cilist io___164 = { 0, 6, 0, fmt_9989, 0 };
00480 static cilist io___165 = { 0, 6, 0, fmt_9988, 0 };
00481 static cilist io___166 = { 0, 6, 0, fmt_9983, 0 };
00482 static cilist io___167 = { 0, 6, 0, 0, 0 };
00483 static cilist io___169 = { 0, 6, 0, fmt_9981, 0 };
00484 static cilist io___170 = { 0, 6, 0, fmt_9981, 0 };
00485 static cilist io___171 = { 0, 6, 0, fmt_9981, 0 };
00486 static cilist io___172 = { 0, 5, 0, 0, 0 };
00487 static cilist io___173 = { 0, 6, 0, fmt_9982, 0 };
00488 static cilist io___174 = { 0, 5, 0, 0, 0 };
00489 static cilist io___176 = { 0, 5, 0, 0, 0 };
00490 static cilist io___178 = { 0, 5, 0, 0, 0 };
00491 static cilist io___179 = { 0, 5, 0, 0, 0 };
00492 static cilist io___181 = { 0, 5, 0, 0, 0 };
00493 static cilist io___183 = { 0, 6, 0, fmt_9999, 0 };
00494 static cilist io___192 = { 0, 6, 0, fmt_9991, 0 };
00495 static cilist io___193 = { 0, 6, 0, fmt_9990, 0 };
00496 static cilist io___196 = { 0, 6, 0, fmt_9961, 0 };
00497 static cilist io___205 = { 0, 6, 0, fmt_9980, 0 };
00498 static cilist io___206 = { 0, 6, 0, fmt_9997, 0 };
00499 static cilist io___208 = { 0, 6, 0, fmt_9980, 0 };
00500 static cilist io___209 = { 0, 6, 0, fmt_9980, 0 };
00501 static cilist io___210 = { 0, 6, 0, fmt_9997, 0 };
00502 static cilist io___211 = { 0, 6, 0, fmt_9980, 0 };
00503 static cilist io___213 = { 0, 6, 0, fmt_9995, 0 };
00504 static cilist io___214 = { 0, 6, 0, fmt_9980, 0 };
00505 static cilist io___215 = { 0, 6, 0, fmt_9990, 0 };
00506 static cilist io___216 = { 0, 6, 0, fmt_9980, 0 };
00507 static cilist io___217 = { 0, 6, 0, fmt_9973, 0 };
00508 static cilist io___218 = { 0, 6, 0, fmt_9990, 0 };
00509 static cilist io___219 = { 0, 6, 0, fmt_9980, 0 };
00510 static cilist io___220 = { 0, 6, 0, fmt_9973, 0 };
00511 static cilist io___221 = { 0, 6, 0, fmt_9990, 0 };
00512 static cilist io___222 = { 0, 6, 0, fmt_9980, 0 };
00513 static cilist io___223 = { 0, 6, 0, fmt_9973, 0 };
00514 static cilist io___224 = { 0, 6, 0, fmt_9990, 0 };
00515 static cilist io___225 = { 0, 6, 0, fmt_9980, 0 };
00516 static cilist io___226 = { 0, 6, 0, fmt_9973, 0 };
00517 static cilist io___227 = { 0, 6, 0, fmt_9996, 0 };
00518 static cilist io___230 = { 0, 6, 0, fmt_9980, 0 };
00519 static cilist io___231 = { 0, 6, 0, fmt_9980, 0 };
00520 static cilist io___232 = { 0, 6, 0, fmt_9990, 0 };
00521 static cilist io___233 = { 0, 6, 0, fmt_9980, 0 };
00522 static cilist io___234 = { 0, 6, 0, fmt_9973, 0 };
00523 static cilist io___235 = { 0, 6, 0, fmt_9990, 0 };
00524 static cilist io___238 = { 0, 6, 0, fmt_9980, 0 };
00525 static cilist io___239 = { 0, 6, 0, fmt_9973, 0 };
00526 static cilist io___240 = { 0, 6, 0, fmt_9990, 0 };
00527 static cilist io___241 = { 0, 6, 0, fmt_9980, 0 };
00528 static cilist io___242 = { 0, 6, 0, fmt_9973, 0 };
00529 static cilist io___243 = { 0, 6, 0, fmt_9990, 0 };
00530 static cilist io___244 = { 0, 6, 0, fmt_9980, 0 };
00531 static cilist io___245 = { 0, 6, 0, fmt_9973, 0 };
00532 static cilist io___246 = { 0, 6, 0, fmt_9980, 0 };
00533 static cilist io___247 = { 0, 6, 0, fmt_9966, 0 };
00534 static cilist io___248 = { 0, 6, 0, fmt_9980, 0 };
00535 static cilist io___251 = { 0, 6, 0, fmt_9980, 0 };
00536 static cilist io___254 = { 0, 6, 0, fmt_9980, 0 };
00537 static cilist io___257 = { 0, 6, 0, fmt_9980, 0 };
00538 static cilist io___258 = { 0, 6, 0, fmt_9980, 0 };
00539 static cilist io___259 = { 0, 6, 0, 0, 0 };
00540 static cilist io___260 = { 0, 6, 0, 0, 0 };
00541 static cilist io___261 = { 0, 6, 0, fmt_9992, 0 };
00542 static cilist io___262 = { 0, 6, 0, fmt_9994, 0 };
00543 static cilist io___264 = { 0, 6, 0, fmt_9993, 0 };
00544
00545
00546
00547
00548
00549
00550
00551
00552
00553
00554
00555
00556
00557
00558
00559
00560
00561
00562
00563
00564
00565
00566
00567
00568
00569
00570
00571
00572
00573
00574
00575
00576
00577
00578
00579
00580
00581
00582
00583
00584
00585
00586
00587
00588
00589
00590
00591
00592
00593
00594
00595
00596
00597
00598
00599
00600
00601
00602
00603
00604
00605
00606
00607
00608
00609
00610
00611
00612
00613
00614
00615
00616
00617
00618
00619
00620
00621
00622
00623
00624
00625
00626
00627
00628
00629
00630
00631
00632
00633
00634
00635
00636
00637
00638
00639
00640
00641
00642
00643
00644
00645
00646
00647
00648
00649
00650
00651
00652
00653
00654
00655
00656
00657
00658
00659
00660
00661
00662
00663
00664
00665
00666
00667
00668
00669
00670
00671
00672
00673
00674
00675
00676
00677
00678
00679
00680
00681
00682
00683
00684
00685
00686
00687
00688
00689
00690
00691
00692
00693
00694
00695
00696
00697
00698
00699
00700
00701
00702
00703
00704
00705
00706
00707
00708
00709
00710
00711
00712
00713
00714
00715
00716
00717
00718
00719
00720
00721
00722
00723
00724
00725
00726
00727
00728
00729
00730
00731
00732
00733
00734
00735
00736
00737
00738
00739
00740
00741
00742
00743
00744
00745
00746
00747
00748
00749
00750
00751
00752
00753
00754
00755
00756
00757
00758
00759
00760
00761
00762
00763
00764
00765
00766
00767
00768
00769
00770
00771
00772
00773
00774
00775
00776
00777
00778
00779
00780
00781
00782
00783
00784
00785
00786
00787
00788
00789
00790
00791
00792
00793
00794
00795
00796
00797
00798
00799
00800
00801
00802
00803
00804
00805
00806
00807
00808
00809
00810
00811
00812
00813
00814
00815
00816
00817
00818
00819
00820
00821
00822
00823
00824
00825
00826
00827
00828
00829
00830
00831
00832
00833
00834
00835
00836
00837
00838
00839
00840
00841
00842
00843
00844
00845
00846
00847
00848
00849
00850
00851
00852
00853
00854
00855
00856
00857
00858
00859
00860
00861
00862
00863
00864
00865
00866
00867
00868
00869
00870
00871
00872
00873
00874
00875
00876
00877
00878
00879
00880
00881
00882
00883
00884
00885
00886
00887
00888
00889
00890
00891
00892
00893
00894
00895
00896
00897
00898
00899
00900
00901
00902
00903
00904
00905
00906
00907
00908
00909
00910
00911
00912
00913
00914
00915
00916
00917
00918
00919
00920
00921
00922
00923
00924
00925
00926
00927
00928
00929
00930
00931
00932
00933
00934
00935
00936
00937
00938
00939
00940
00941
00942
00943
00944
00945
00946
00947
00948
00949
00950
00951
00952
00953
00954
00955
00956
00957
00958
00959
00960
00961
00962
00963
00964
00965
00966
00967
00968
00969
00970
00971
00972
00973
00974
00975
00976
00977
00978
00979
00980
00981
00982
00983
00984
00985
00986
00987
00988
00989
00990
00991
00992
00993
00994
00995
00996
00997
00998
00999
01000
01001
01002
01003
01004
01005
01006
01007
01008
01009
01010
01011
01012
01013
01014
01015
01016
01017
01018
01019
01020
01021
01022
01023
01024
01025
01026
01027
01028
01029
01030
01031
01032
01033
01034
01035
01036
01037
01038
01039
01040
01041
01042
01043
01044
01045
01046
01047
01048
01049
01050
01051
01052
01053
01054
01055
01056
01057
01058
01059
01060
01061
01062
01063
01064
01065
01066
01067
01068
01069
01070
01071
01072
01073
01074
01075
01076
01077
01078
01079
01080
01081
01082
01083
01084
01085
01086
01087
01088
01089
01090
01091
01092
01093
01094
01095
01096
01097
01098
01099
01100
01101
01102
01103
01104
01105
01106
01107
01108
01109
01110
01111
01112
01113
01114
01115
01116
01117
01118
01119
01120
01121
01122
01123
01124
01125
01126
01127
01128
01129
01130
01131
01132
01133
01134
01135
01136
01137
01138
01139
01140
01141
01142
01143
01144
01145
01146
01147
01148
01149
01150
01151
01152
01153
01154
01155
01156
01157
01158
01159
01160
01161
01162
01163
01164
01165
01166
01167
01168
01169
01170
01171
01172
01173
01174
01175
01176
01177
01178
01179
01180
01181
01182
01183
01184
01185
01186
01187
01188
01189
01190
01191
01192
01193
01194
01195
01196
01197
01198
01199
01200
01201
01202
01203
01204
01205
01206
01207
01208
01209
01210
01211
01212
01213
01214
01215
01216
01217
01218
01219
01220
01221
01222
01223
01224
01225
01226
01227
01228
01229
01230
01231
01232
01233
01234
01235
01236
01237
01238
01239
01240
01241
01242
01243
01244
01245
01246
01247
01248
01249
01250
01251
01252
01253
01254
01255
01256
01257
01258
01259
01260
01261
01262
01263
01264
01265
01266
01267
01268
01269
01270
01271
01272
01273
01274
01275
01276
01277
01278
01279
01280
01281
01282
01283
01284
01285
01286
01287
01288
01289
01290
01291
01292
01293
01294
01295
01296
01297
01298
01299
01300
01301
01302
01303
01304
01305
01306
01307
01308
01309
01310
01311
01312
01313
01314
01315
01316
01317
01318
01319
01320
01321
01322
01323
01324
01325
01326
01327
01328
01329
01330
01331
01332
01333
01334
01335
01336
01337
01338
01339
01340
01341
01342
01343
01344
01345
01346
01347
01348
01349
01350
01351
01352
01353
01354
01355
01356
01357
01358
01359
01360
01361
01362
01363
01364
01365
01366
01367
01368
01369
01370
01371
01372
01373
01374
01375
01376
01377
01378
01379
01380
01381
01382
01383
01384
01385
01386
01387
01388
01389
01390
01391
01392
01393
01394
01395
01396
01397
01398
01399
01400
01401
01402
01403
01404
01405
01406
01407
01408
01409
01410
01411
01412
01413
01414
01415
01416
01417
01418
01419
01420
01421
01422
01423
01424
01425
01426
01427
01428
01429
01430
01431
01432
01433
01434
01435
01436
01437
01438
01439
01440
01441
01442
01443
01444
01445
01446
01447
01448
01449
01450
01451
01452
01453
01454
01455
01456
01457
01458
01459
01460
01461
01462
01463
01464
01465
01466
01467
01468
01469
01470
01471
01472
01473
01474
01475
01476
01477
01478
01479
01480
01481
01482
01483
01484
01485
01486
01487
01488
01489
01490
01491
01492
01493
01494
01495
01496
01497
01498
01499
01500
01501
01502
01503
01504
01505
01506
01507
01508
01509
01510
01511
01512
01513
01514
01515
01516
01517
01518
01519
01520
01521
01522
01523
01524
01525
01526
01527
01528
01529 s1 = dsecnd_();
01530 fatal = FALSE_;
01531 infoc_1.nunit = 6;
01532
01533
01534
01535 L10:
01536
01537
01538
01539 ci__1.cierr = 0;
01540 ci__1.ciend = 1;
01541 ci__1.ciunit = 5;
01542 ci__1.cifmt = "(A80)";
01543 i__1 = s_rsfe(&ci__1);
01544 if (i__1 != 0) {
01545 goto L380;
01546 }
01547 i__1 = do_fio(&c__1, line, (ftnlen)80);
01548 if (i__1 != 0) {
01549 goto L380;
01550 }
01551 i__1 = e_rsfe();
01552 if (i__1 != 0) {
01553 goto L380;
01554 }
01555 s_copy(path, line, (ftnlen)3, (ftnlen)3);
01556 nep = lsamen_(&c__3, path, "NEP") || lsamen_(&c__3,
01557 path, "ZHS");
01558 sep = lsamen_(&c__3, path, "SEP") || lsamen_(&c__3,
01559 path, "ZST") || lsamen_(&c__3, path, "ZSG");
01560 svd = lsamen_(&c__3, path, "SVD") || lsamen_(&c__3,
01561 path, "ZBD");
01562 zev = lsamen_(&c__3, path, "ZEV");
01563 zes = lsamen_(&c__3, path, "ZES");
01564 zvx = lsamen_(&c__3, path, "ZVX");
01565 zsx = lsamen_(&c__3, path, "ZSX");
01566 zgg = lsamen_(&c__3, path, "ZGG");
01567 zgs = lsamen_(&c__3, path, "ZGS");
01568 zgx = lsamen_(&c__3, path, "ZGX");
01569 zgv = lsamen_(&c__3, path, "ZGV");
01570 zxv = lsamen_(&c__3, path, "ZXV");
01571 zhb = lsamen_(&c__3, path, "ZHB");
01572 zbb = lsamen_(&c__3, path, "ZBB");
01573 glm = lsamen_(&c__3, path, "GLM");
01574 gqr = lsamen_(&c__3, path, "GQR") || lsamen_(&c__3,
01575 path, "GRQ");
01576 gsv = lsamen_(&c__3, path, "GSV");
01577 lse = lsamen_(&c__3, path, "LSE");
01578 zbl = lsamen_(&c__3, path, "ZBL");
01579 zbk = lsamen_(&c__3, path, "ZBK");
01580 zgl = lsamen_(&c__3, path, "ZGL");
01581 zgk = lsamen_(&c__3, path, "ZGK");
01582
01583
01584
01585 if (s_cmp(path, " ", (ftnlen)3, (ftnlen)3) == 0) {
01586 goto L10;
01587 } else if (nep) {
01588 s_wsfe(&io___29);
01589 e_wsfe();
01590 } else if (sep) {
01591 s_wsfe(&io___30);
01592 e_wsfe();
01593 } else if (svd) {
01594 s_wsfe(&io___31);
01595 e_wsfe();
01596 } else if (zev) {
01597 s_wsfe(&io___32);
01598 e_wsfe();
01599 } else if (zes) {
01600 s_wsfe(&io___33);
01601 e_wsfe();
01602 } else if (zvx) {
01603 s_wsfe(&io___34);
01604 e_wsfe();
01605 } else if (zsx) {
01606 s_wsfe(&io___35);
01607 e_wsfe();
01608 } else if (zgg) {
01609 s_wsfe(&io___36);
01610 e_wsfe();
01611 } else if (zgs) {
01612 s_wsfe(&io___37);
01613 e_wsfe();
01614 } else if (zgx) {
01615 s_wsfe(&io___38);
01616 e_wsfe();
01617 } else if (zgv) {
01618 s_wsfe(&io___39);
01619 e_wsfe();
01620 } else if (zxv) {
01621 s_wsfe(&io___40);
01622 e_wsfe();
01623 } else if (zhb) {
01624 s_wsfe(&io___41);
01625 e_wsfe();
01626 } else if (zbb) {
01627 s_wsfe(&io___42);
01628 e_wsfe();
01629 } else if (glm) {
01630 s_wsfe(&io___43);
01631 e_wsfe();
01632 } else if (gqr) {
01633 s_wsfe(&io___44);
01634 e_wsfe();
01635 } else if (gsv) {
01636 s_wsfe(&io___45);
01637 e_wsfe();
01638 } else if (lse) {
01639 s_wsfe(&io___46);
01640 e_wsfe();
01641 } else if (zbl) {
01642
01643
01644
01645 zchkbl_(&c__5, &c__6);
01646 goto L380;
01647 } else if (zbk) {
01648
01649
01650
01651 zchkbk_(&c__5, &c__6);
01652 goto L380;
01653 } else if (zgl) {
01654
01655
01656
01657 zchkgl_(&c__5, &c__6);
01658 goto L380;
01659 } else if (zgk) {
01660
01661
01662
01663 zchkgk_(&c__5, &c__6);
01664 goto L380;
01665 } else if (lsamen_(&c__3, path, "ZEC")) {
01666
01667
01668
01669 s_rsle(&io___47);
01670 do_lio(&c__5, &c__1, (char *)&thresh, (ftnlen)sizeof(doublereal));
01671 e_rsle();
01672 xlaenv_(&c__1, &c__1);
01673 tsterr = TRUE_;
01674 zchkec_(&thresh, &tsterr, &c__5, &c__6);
01675 goto L380;
01676 } else {
01677 s_wsfe(&io___50);
01678 do_fio(&c__1, path, (ftnlen)3);
01679 e_wsfe();
01680 goto L380;
01681 }
01682 ilaver_(&vers_major__, &vers_minor__, &vers_patch__);
01683 s_wsfe(&io___54);
01684 do_fio(&c__1, (char *)&vers_major__, (ftnlen)sizeof(integer));
01685 do_fio(&c__1, (char *)&vers_minor__, (ftnlen)sizeof(integer));
01686 do_fio(&c__1, (char *)&vers_patch__, (ftnlen)sizeof(integer));
01687 e_wsfe();
01688 s_wsfe(&io___55);
01689 e_wsfe();
01690
01691
01692
01693 s_rsle(&io___56);
01694 do_lio(&c__3, &c__1, (char *)&nn, (ftnlen)sizeof(integer));
01695 e_rsle();
01696 if (nn < 0) {
01697 s_wsfe(&io___58);
01698 do_fio(&c__1, " NN ", (ftnlen)6);
01699 do_fio(&c__1, (char *)&nn, (ftnlen)sizeof(integer));
01700 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
01701 e_wsfe();
01702 nn = 0;
01703 fatal = TRUE_;
01704 } else if (nn > 20) {
01705 s_wsfe(&io___59);
01706 do_fio(&c__1, " NN ", (ftnlen)6);
01707 do_fio(&c__1, (char *)&nn, (ftnlen)sizeof(integer));
01708 do_fio(&c__1, (char *)&c__20, (ftnlen)sizeof(integer));
01709 e_wsfe();
01710 nn = 0;
01711 fatal = TRUE_;
01712 }
01713
01714
01715
01716 if (! (zgx || zxv)) {
01717 s_rsle(&io___60);
01718 i__1 = nn;
01719 for (i__ = 1; i__ <= i__1; ++i__) {
01720 do_lio(&c__3, &c__1, (char *)&mval[i__ - 1], (ftnlen)sizeof(
01721 integer));
01722 }
01723 e_rsle();
01724 if (svd) {
01725 s_copy(vname, " M ", (ftnlen)32, (ftnlen)6);
01726 } else {
01727 s_copy(vname, " N ", (ftnlen)32, (ftnlen)6);
01728 }
01729 i__1 = nn;
01730 for (i__ = 1; i__ <= i__1; ++i__) {
01731 if (mval[i__ - 1] < 0) {
01732 s_wsfe(&io___64);
01733 do_fio(&c__1, vname, (ftnlen)32);
01734 do_fio(&c__1, (char *)&mval[i__ - 1], (ftnlen)sizeof(integer))
01735 ;
01736 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
01737 e_wsfe();
01738 fatal = TRUE_;
01739 } else if (mval[i__ - 1] > 132) {
01740 s_wsfe(&io___65);
01741 do_fio(&c__1, vname, (ftnlen)32);
01742 do_fio(&c__1, (char *)&mval[i__ - 1], (ftnlen)sizeof(integer))
01743 ;
01744 do_fio(&c__1, (char *)&c__132, (ftnlen)sizeof(integer));
01745 e_wsfe();
01746 fatal = TRUE_;
01747 }
01748
01749 }
01750 s_wsfe(&io___66);
01751 do_fio(&c__1, "M: ", (ftnlen)6);
01752 i__1 = nn;
01753 for (i__ = 1; i__ <= i__1; ++i__) {
01754 do_fio(&c__1, (char *)&mval[i__ - 1], (ftnlen)sizeof(integer));
01755 }
01756 e_wsfe();
01757 }
01758
01759
01760
01761 if (glm || gqr || gsv || lse) {
01762 s_rsle(&io___67);
01763 i__1 = nn;
01764 for (i__ = 1; i__ <= i__1; ++i__) {
01765 do_lio(&c__3, &c__1, (char *)&pval[i__ - 1], (ftnlen)sizeof(
01766 integer));
01767 }
01768 e_rsle();
01769 i__1 = nn;
01770 for (i__ = 1; i__ <= i__1; ++i__) {
01771 if (pval[i__ - 1] < 0) {
01772 s_wsfe(&io___69);
01773 do_fio(&c__1, " P ", (ftnlen)4);
01774 do_fio(&c__1, (char *)&pval[i__ - 1], (ftnlen)sizeof(integer))
01775 ;
01776 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
01777 e_wsfe();
01778 fatal = TRUE_;
01779 } else if (pval[i__ - 1] > 132) {
01780 s_wsfe(&io___70);
01781 do_fio(&c__1, " P ", (ftnlen)4);
01782 do_fio(&c__1, (char *)&pval[i__ - 1], (ftnlen)sizeof(integer))
01783 ;
01784 do_fio(&c__1, (char *)&c__132, (ftnlen)sizeof(integer));
01785 e_wsfe();
01786 fatal = TRUE_;
01787 }
01788
01789 }
01790 s_wsfe(&io___71);
01791 do_fio(&c__1, "P: ", (ftnlen)6);
01792 i__1 = nn;
01793 for (i__ = 1; i__ <= i__1; ++i__) {
01794 do_fio(&c__1, (char *)&pval[i__ - 1], (ftnlen)sizeof(integer));
01795 }
01796 e_wsfe();
01797 }
01798
01799
01800
01801 if (svd || zbb || glm || gqr || gsv || lse) {
01802 s_rsle(&io___72);
01803 i__1 = nn;
01804 for (i__ = 1; i__ <= i__1; ++i__) {
01805 do_lio(&c__3, &c__1, (char *)&nval[i__ - 1], (ftnlen)sizeof(
01806 integer));
01807 }
01808 e_rsle();
01809 i__1 = nn;
01810 for (i__ = 1; i__ <= i__1; ++i__) {
01811 if (nval[i__ - 1] < 0) {
01812 s_wsfe(&io___74);
01813 do_fio(&c__1, " N ", (ftnlen)6);
01814 do_fio(&c__1, (char *)&nval[i__ - 1], (ftnlen)sizeof(integer))
01815 ;
01816 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
01817 e_wsfe();
01818 fatal = TRUE_;
01819 } else if (nval[i__ - 1] > 132) {
01820 s_wsfe(&io___75);
01821 do_fio(&c__1, " N ", (ftnlen)6);
01822 do_fio(&c__1, (char *)&nval[i__ - 1], (ftnlen)sizeof(integer))
01823 ;
01824 do_fio(&c__1, (char *)&c__132, (ftnlen)sizeof(integer));
01825 e_wsfe();
01826 fatal = TRUE_;
01827 }
01828
01829 }
01830 } else {
01831 i__1 = nn;
01832 for (i__ = 1; i__ <= i__1; ++i__) {
01833 nval[i__ - 1] = mval[i__ - 1];
01834
01835 }
01836 }
01837 if (! (zgx || zxv)) {
01838 s_wsfe(&io___76);
01839 do_fio(&c__1, "N: ", (ftnlen)6);
01840 i__1 = nn;
01841 for (i__ = 1; i__ <= i__1; ++i__) {
01842 do_fio(&c__1, (char *)&nval[i__ - 1], (ftnlen)sizeof(integer));
01843 }
01844 e_wsfe();
01845 } else {
01846 s_wsfe(&io___77);
01847 do_fio(&c__1, "N: ", (ftnlen)6);
01848 do_fio(&c__1, (char *)&nn, (ftnlen)sizeof(integer));
01849 e_wsfe();
01850 }
01851
01852
01853
01854 if (zhb || zbb) {
01855 s_rsle(&io___78);
01856 do_lio(&c__3, &c__1, (char *)&nk, (ftnlen)sizeof(integer));
01857 e_rsle();
01858 s_rsle(&io___80);
01859 i__1 = nk;
01860 for (i__ = 1; i__ <= i__1; ++i__) {
01861 do_lio(&c__3, &c__1, (char *)&kval[i__ - 1], (ftnlen)sizeof(
01862 integer));
01863 }
01864 e_rsle();
01865 i__1 = nk;
01866 for (i__ = 1; i__ <= i__1; ++i__) {
01867 if (kval[i__ - 1] < 0) {
01868 s_wsfe(&io___82);
01869 do_fio(&c__1, " K ", (ftnlen)6);
01870 do_fio(&c__1, (char *)&kval[i__ - 1], (ftnlen)sizeof(integer))
01871 ;
01872 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
01873 e_wsfe();
01874 fatal = TRUE_;
01875 } else if (kval[i__ - 1] > 132) {
01876 s_wsfe(&io___83);
01877 do_fio(&c__1, " K ", (ftnlen)6);
01878 do_fio(&c__1, (char *)&kval[i__ - 1], (ftnlen)sizeof(integer))
01879 ;
01880 do_fio(&c__1, (char *)&c__132, (ftnlen)sizeof(integer));
01881 e_wsfe();
01882 fatal = TRUE_;
01883 }
01884
01885 }
01886 s_wsfe(&io___84);
01887 do_fio(&c__1, "K: ", (ftnlen)6);
01888 i__1 = nk;
01889 for (i__ = 1; i__ <= i__1; ++i__) {
01890 do_fio(&c__1, (char *)&kval[i__ - 1], (ftnlen)sizeof(integer));
01891 }
01892 e_wsfe();
01893 }
01894
01895 if (zev || zes || zvx || zsx) {
01896
01897
01898
01899
01900 s_rsle(&io___85);
01901 do_lio(&c__3, &c__1, (char *)&nbval[0], (ftnlen)sizeof(integer));
01902 do_lio(&c__3, &c__1, (char *)&nbmin[0], (ftnlen)sizeof(integer));
01903 do_lio(&c__3, &c__1, (char *)&nxval[0], (ftnlen)sizeof(integer));
01904 do_lio(&c__3, &c__1, (char *)&inmin[0], (ftnlen)sizeof(integer));
01905 do_lio(&c__3, &c__1, (char *)&inwin[0], (ftnlen)sizeof(integer));
01906 do_lio(&c__3, &c__1, (char *)&inibl[0], (ftnlen)sizeof(integer));
01907 do_lio(&c__3, &c__1, (char *)&ishfts[0], (ftnlen)sizeof(integer));
01908 do_lio(&c__3, &c__1, (char *)&iacc22[0], (ftnlen)sizeof(integer));
01909 e_rsle();
01910 if (nbval[0] < 1) {
01911 s_wsfe(&io___94);
01912 do_fio(&c__1, " NB ", (ftnlen)6);
01913 do_fio(&c__1, (char *)&nbval[0], (ftnlen)sizeof(integer));
01914 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
01915 e_wsfe();
01916 fatal = TRUE_;
01917 } else if (nbmin[0] < 1) {
01918 s_wsfe(&io___95);
01919 do_fio(&c__1, "NBMIN ", (ftnlen)6);
01920 do_fio(&c__1, (char *)&nbmin[0], (ftnlen)sizeof(integer));
01921 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
01922 e_wsfe();
01923 fatal = TRUE_;
01924 } else if (nxval[0] < 1) {
01925 s_wsfe(&io___96);
01926 do_fio(&c__1, " NX ", (ftnlen)6);
01927 do_fio(&c__1, (char *)&nxval[0], (ftnlen)sizeof(integer));
01928 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
01929 e_wsfe();
01930 fatal = TRUE_;
01931 } else if (inmin[0] < 1) {
01932 s_wsfe(&io___97);
01933 do_fio(&c__1, " INMIN ", (ftnlen)9);
01934 do_fio(&c__1, (char *)&inmin[0], (ftnlen)sizeof(integer));
01935 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
01936 e_wsfe();
01937 fatal = TRUE_;
01938 } else if (inwin[0] < 1) {
01939 s_wsfe(&io___98);
01940 do_fio(&c__1, " INWIN ", (ftnlen)9);
01941 do_fio(&c__1, (char *)&inwin[0], (ftnlen)sizeof(integer));
01942 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
01943 e_wsfe();
01944 fatal = TRUE_;
01945 } else if (inibl[0] < 1) {
01946 s_wsfe(&io___99);
01947 do_fio(&c__1, " INIBL ", (ftnlen)9);
01948 do_fio(&c__1, (char *)&inibl[0], (ftnlen)sizeof(integer));
01949 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
01950 e_wsfe();
01951 fatal = TRUE_;
01952 } else if (ishfts[0] < 1) {
01953 s_wsfe(&io___100);
01954 do_fio(&c__1, " ISHFTS ", (ftnlen)10);
01955 do_fio(&c__1, (char *)&ishfts[0], (ftnlen)sizeof(integer));
01956 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
01957 e_wsfe();
01958 fatal = TRUE_;
01959 } else if (iacc22[0] < 0) {
01960 s_wsfe(&io___101);
01961 do_fio(&c__1, " IACC22 ", (ftnlen)10);
01962 do_fio(&c__1, (char *)&iacc22[0], (ftnlen)sizeof(integer));
01963 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
01964 e_wsfe();
01965 fatal = TRUE_;
01966 }
01967 xlaenv_(&c__1, nbval);
01968 xlaenv_(&c__2, nbmin);
01969 xlaenv_(&c__3, nxval);
01970 i__1 = max(11,inmin[0]);
01971 xlaenv_(&c__12, &i__1);
01972 xlaenv_(&c__13, inwin);
01973 xlaenv_(&c__14, inibl);
01974 xlaenv_(&c__15, ishfts);
01975 xlaenv_(&c__16, iacc22);
01976 s_wsfe(&io___102);
01977 do_fio(&c__1, "NB: ", (ftnlen)6);
01978 do_fio(&c__1, (char *)&nbval[0], (ftnlen)sizeof(integer));
01979 e_wsfe();
01980 s_wsfe(&io___103);
01981 do_fio(&c__1, "NBMIN:", (ftnlen)6);
01982 do_fio(&c__1, (char *)&nbmin[0], (ftnlen)sizeof(integer));
01983 e_wsfe();
01984 s_wsfe(&io___104);
01985 do_fio(&c__1, "NX: ", (ftnlen)6);
01986 do_fio(&c__1, (char *)&nxval[0], (ftnlen)sizeof(integer));
01987 e_wsfe();
01988 s_wsfe(&io___105);
01989 do_fio(&c__1, "INMIN: ", (ftnlen)9);
01990 do_fio(&c__1, (char *)&inmin[0], (ftnlen)sizeof(integer));
01991 e_wsfe();
01992 s_wsfe(&io___106);
01993 do_fio(&c__1, "INWIN: ", (ftnlen)7);
01994 do_fio(&c__1, (char *)&inwin[0], (ftnlen)sizeof(integer));
01995 e_wsfe();
01996 s_wsfe(&io___107);
01997 do_fio(&c__1, "INIBL: ", (ftnlen)7);
01998 do_fio(&c__1, (char *)&inibl[0], (ftnlen)sizeof(integer));
01999 e_wsfe();
02000 s_wsfe(&io___108);
02001 do_fio(&c__1, "ISHFTS: ", (ftnlen)8);
02002 do_fio(&c__1, (char *)&ishfts[0], (ftnlen)sizeof(integer));
02003 e_wsfe();
02004 s_wsfe(&io___109);
02005 do_fio(&c__1, "IACC22: ", (ftnlen)8);
02006 do_fio(&c__1, (char *)&iacc22[0], (ftnlen)sizeof(integer));
02007 e_wsfe();
02008
02009 } else if (zgs || zgx || zgv || zxv) {
02010
02011
02012
02013
02014 s_rsle(&io___110);
02015 do_lio(&c__3, &c__1, (char *)&nbval[0], (ftnlen)sizeof(integer));
02016 do_lio(&c__3, &c__1, (char *)&nbmin[0], (ftnlen)sizeof(integer));
02017 do_lio(&c__3, &c__1, (char *)&nxval[0], (ftnlen)sizeof(integer));
02018 do_lio(&c__3, &c__1, (char *)&nsval[0], (ftnlen)sizeof(integer));
02019 do_lio(&c__3, &c__1, (char *)&mxbval[0], (ftnlen)sizeof(integer));
02020 e_rsle();
02021 if (nbval[0] < 1) {
02022 s_wsfe(&io___113);
02023 do_fio(&c__1, " NB ", (ftnlen)6);
02024 do_fio(&c__1, (char *)&nbval[0], (ftnlen)sizeof(integer));
02025 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
02026 e_wsfe();
02027 fatal = TRUE_;
02028 } else if (nbmin[0] < 1) {
02029 s_wsfe(&io___114);
02030 do_fio(&c__1, "NBMIN ", (ftnlen)6);
02031 do_fio(&c__1, (char *)&nbmin[0], (ftnlen)sizeof(integer));
02032 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
02033 e_wsfe();
02034 fatal = TRUE_;
02035 } else if (nxval[0] < 1) {
02036 s_wsfe(&io___115);
02037 do_fio(&c__1, " NX ", (ftnlen)6);
02038 do_fio(&c__1, (char *)&nxval[0], (ftnlen)sizeof(integer));
02039 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
02040 e_wsfe();
02041 fatal = TRUE_;
02042 } else if (nsval[0] < 2) {
02043 s_wsfe(&io___116);
02044 do_fio(&c__1, " NS ", (ftnlen)6);
02045 do_fio(&c__1, (char *)&nsval[0], (ftnlen)sizeof(integer));
02046 do_fio(&c__1, (char *)&c__2, (ftnlen)sizeof(integer));
02047 e_wsfe();
02048 fatal = TRUE_;
02049 } else if (mxbval[0] < 1) {
02050 s_wsfe(&io___117);
02051 do_fio(&c__1, " MAXB ", (ftnlen)6);
02052 do_fio(&c__1, (char *)&mxbval[0], (ftnlen)sizeof(integer));
02053 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
02054 e_wsfe();
02055 fatal = TRUE_;
02056 }
02057 xlaenv_(&c__1, nbval);
02058 xlaenv_(&c__2, nbmin);
02059 xlaenv_(&c__3, nxval);
02060 xlaenv_(&c__4, nsval);
02061 xlaenv_(&c__8, mxbval);
02062 s_wsfe(&io___118);
02063 do_fio(&c__1, "NB: ", (ftnlen)6);
02064 do_fio(&c__1, (char *)&nbval[0], (ftnlen)sizeof(integer));
02065 e_wsfe();
02066 s_wsfe(&io___119);
02067 do_fio(&c__1, "NBMIN:", (ftnlen)6);
02068 do_fio(&c__1, (char *)&nbmin[0], (ftnlen)sizeof(integer));
02069 e_wsfe();
02070 s_wsfe(&io___120);
02071 do_fio(&c__1, "NX: ", (ftnlen)6);
02072 do_fio(&c__1, (char *)&nxval[0], (ftnlen)sizeof(integer));
02073 e_wsfe();
02074 s_wsfe(&io___121);
02075 do_fio(&c__1, "NS: ", (ftnlen)6);
02076 do_fio(&c__1, (char *)&nsval[0], (ftnlen)sizeof(integer));
02077 e_wsfe();
02078 s_wsfe(&io___122);
02079 do_fio(&c__1, "MAXB: ", (ftnlen)6);
02080 do_fio(&c__1, (char *)&mxbval[0], (ftnlen)sizeof(integer));
02081 e_wsfe();
02082 } else if (! zhb && ! glm && ! gqr && ! gsv && ! lse) {
02083
02084
02085
02086
02087 s_rsle(&io___123);
02088 do_lio(&c__3, &c__1, (char *)&nparms, (ftnlen)sizeof(integer));
02089 e_rsle();
02090 if (nparms < 1) {
02091 s_wsfe(&io___125);
02092 do_fio(&c__1, "NPARMS", (ftnlen)6);
02093 do_fio(&c__1, (char *)&nparms, (ftnlen)sizeof(integer));
02094 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
02095 e_wsfe();
02096 nparms = 0;
02097 fatal = TRUE_;
02098 } else if (nparms > 20) {
02099 s_wsfe(&io___126);
02100 do_fio(&c__1, "NPARMS", (ftnlen)6);
02101 do_fio(&c__1, (char *)&nparms, (ftnlen)sizeof(integer));
02102 do_fio(&c__1, (char *)&c__20, (ftnlen)sizeof(integer));
02103 e_wsfe();
02104 nparms = 0;
02105 fatal = TRUE_;
02106 }
02107
02108
02109
02110 if (! zbb) {
02111 s_rsle(&io___127);
02112 i__1 = nparms;
02113 for (i__ = 1; i__ <= i__1; ++i__) {
02114 do_lio(&c__3, &c__1, (char *)&nbval[i__ - 1], (ftnlen)sizeof(
02115 integer));
02116 }
02117 e_rsle();
02118 i__1 = nparms;
02119 for (i__ = 1; i__ <= i__1; ++i__) {
02120 if (nbval[i__ - 1] < 0) {
02121 s_wsfe(&io___128);
02122 do_fio(&c__1, " NB ", (ftnlen)6);
02123 do_fio(&c__1, (char *)&nbval[i__ - 1], (ftnlen)sizeof(
02124 integer));
02125 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
02126 e_wsfe();
02127 fatal = TRUE_;
02128 } else if (nbval[i__ - 1] > 132) {
02129 s_wsfe(&io___129);
02130 do_fio(&c__1, " NB ", (ftnlen)6);
02131 do_fio(&c__1, (char *)&nbval[i__ - 1], (ftnlen)sizeof(
02132 integer));
02133 do_fio(&c__1, (char *)&c__132, (ftnlen)sizeof(integer));
02134 e_wsfe();
02135 fatal = TRUE_;
02136 }
02137
02138 }
02139 s_wsfe(&io___130);
02140 do_fio(&c__1, "NB: ", (ftnlen)6);
02141 i__1 = nparms;
02142 for (i__ = 1; i__ <= i__1; ++i__) {
02143 do_fio(&c__1, (char *)&nbval[i__ - 1], (ftnlen)sizeof(integer)
02144 );
02145 }
02146 e_wsfe();
02147 }
02148
02149
02150
02151 if (nep || sep || svd || zgg) {
02152 s_rsle(&io___131);
02153 i__1 = nparms;
02154 for (i__ = 1; i__ <= i__1; ++i__) {
02155 do_lio(&c__3, &c__1, (char *)&nbmin[i__ - 1], (ftnlen)sizeof(
02156 integer));
02157 }
02158 e_rsle();
02159 i__1 = nparms;
02160 for (i__ = 1; i__ <= i__1; ++i__) {
02161 if (nbmin[i__ - 1] < 0) {
02162 s_wsfe(&io___132);
02163 do_fio(&c__1, "NBMIN ", (ftnlen)6);
02164 do_fio(&c__1, (char *)&nbmin[i__ - 1], (ftnlen)sizeof(
02165 integer));
02166 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
02167 e_wsfe();
02168 fatal = TRUE_;
02169 } else if (nbmin[i__ - 1] > 132) {
02170 s_wsfe(&io___133);
02171 do_fio(&c__1, "NBMIN ", (ftnlen)6);
02172 do_fio(&c__1, (char *)&nbmin[i__ - 1], (ftnlen)sizeof(
02173 integer));
02174 do_fio(&c__1, (char *)&c__132, (ftnlen)sizeof(integer));
02175 e_wsfe();
02176 fatal = TRUE_;
02177 }
02178
02179 }
02180 s_wsfe(&io___134);
02181 do_fio(&c__1, "NBMIN:", (ftnlen)6);
02182 i__1 = nparms;
02183 for (i__ = 1; i__ <= i__1; ++i__) {
02184 do_fio(&c__1, (char *)&nbmin[i__ - 1], (ftnlen)sizeof(integer)
02185 );
02186 }
02187 e_wsfe();
02188 } else {
02189 i__1 = nparms;
02190 for (i__ = 1; i__ <= i__1; ++i__) {
02191 nbmin[i__ - 1] = 1;
02192
02193 }
02194 }
02195
02196
02197
02198 if (nep || sep || svd) {
02199 s_rsle(&io___135);
02200 i__1 = nparms;
02201 for (i__ = 1; i__ <= i__1; ++i__) {
02202 do_lio(&c__3, &c__1, (char *)&nxval[i__ - 1], (ftnlen)sizeof(
02203 integer));
02204 }
02205 e_rsle();
02206 i__1 = nparms;
02207 for (i__ = 1; i__ <= i__1; ++i__) {
02208 if (nxval[i__ - 1] < 0) {
02209 s_wsfe(&io___136);
02210 do_fio(&c__1, " NX ", (ftnlen)6);
02211 do_fio(&c__1, (char *)&nxval[i__ - 1], (ftnlen)sizeof(
02212 integer));
02213 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
02214 e_wsfe();
02215 fatal = TRUE_;
02216 } else if (nxval[i__ - 1] > 132) {
02217 s_wsfe(&io___137);
02218 do_fio(&c__1, " NX ", (ftnlen)6);
02219 do_fio(&c__1, (char *)&nxval[i__ - 1], (ftnlen)sizeof(
02220 integer));
02221 do_fio(&c__1, (char *)&c__132, (ftnlen)sizeof(integer));
02222 e_wsfe();
02223 fatal = TRUE_;
02224 }
02225
02226 }
02227 s_wsfe(&io___138);
02228 do_fio(&c__1, "NX: ", (ftnlen)6);
02229 i__1 = nparms;
02230 for (i__ = 1; i__ <= i__1; ++i__) {
02231 do_fio(&c__1, (char *)&nxval[i__ - 1], (ftnlen)sizeof(integer)
02232 );
02233 }
02234 e_wsfe();
02235 } else {
02236 i__1 = nparms;
02237 for (i__ = 1; i__ <= i__1; ++i__) {
02238 nxval[i__ - 1] = 1;
02239
02240 }
02241 }
02242
02243
02244
02245
02246 if (svd || zbb || zgg) {
02247 s_rsle(&io___139);
02248 i__1 = nparms;
02249 for (i__ = 1; i__ <= i__1; ++i__) {
02250 do_lio(&c__3, &c__1, (char *)&nsval[i__ - 1], (ftnlen)sizeof(
02251 integer));
02252 }
02253 e_rsle();
02254 i__1 = nparms;
02255 for (i__ = 1; i__ <= i__1; ++i__) {
02256 if (nsval[i__ - 1] < 0) {
02257 s_wsfe(&io___140);
02258 do_fio(&c__1, " NS ", (ftnlen)6);
02259 do_fio(&c__1, (char *)&nsval[i__ - 1], (ftnlen)sizeof(
02260 integer));
02261 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
02262 e_wsfe();
02263 fatal = TRUE_;
02264 } else if (nsval[i__ - 1] > 132) {
02265 s_wsfe(&io___141);
02266 do_fio(&c__1, " NS ", (ftnlen)6);
02267 do_fio(&c__1, (char *)&nsval[i__ - 1], (ftnlen)sizeof(
02268 integer));
02269 do_fio(&c__1, (char *)&c__132, (ftnlen)sizeof(integer));
02270 e_wsfe();
02271 fatal = TRUE_;
02272 }
02273
02274 }
02275 s_wsfe(&io___142);
02276 do_fio(&c__1, "NS: ", (ftnlen)6);
02277 i__1 = nparms;
02278 for (i__ = 1; i__ <= i__1; ++i__) {
02279 do_fio(&c__1, (char *)&nsval[i__ - 1], (ftnlen)sizeof(integer)
02280 );
02281 }
02282 e_wsfe();
02283 } else {
02284 i__1 = nparms;
02285 for (i__ = 1; i__ <= i__1; ++i__) {
02286 nsval[i__ - 1] = 1;
02287
02288 }
02289 }
02290
02291
02292
02293 if (zgg) {
02294 s_rsle(&io___143);
02295 i__1 = nparms;
02296 for (i__ = 1; i__ <= i__1; ++i__) {
02297 do_lio(&c__3, &c__1, (char *)&mxbval[i__ - 1], (ftnlen)sizeof(
02298 integer));
02299 }
02300 e_rsle();
02301 i__1 = nparms;
02302 for (i__ = 1; i__ <= i__1; ++i__) {
02303 if (mxbval[i__ - 1] < 0) {
02304 s_wsfe(&io___144);
02305 do_fio(&c__1, " MAXB ", (ftnlen)6);
02306 do_fio(&c__1, (char *)&mxbval[i__ - 1], (ftnlen)sizeof(
02307 integer));
02308 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
02309 e_wsfe();
02310 fatal = TRUE_;
02311 } else if (mxbval[i__ - 1] > 132) {
02312 s_wsfe(&io___145);
02313 do_fio(&c__1, " MAXB ", (ftnlen)6);
02314 do_fio(&c__1, (char *)&mxbval[i__ - 1], (ftnlen)sizeof(
02315 integer));
02316 do_fio(&c__1, (char *)&c__132, (ftnlen)sizeof(integer));
02317 e_wsfe();
02318 fatal = TRUE_;
02319 }
02320
02321 }
02322 s_wsfe(&io___146);
02323 do_fio(&c__1, "MAXB: ", (ftnlen)6);
02324 i__1 = nparms;
02325 for (i__ = 1; i__ <= i__1; ++i__) {
02326 do_fio(&c__1, (char *)&mxbval[i__ - 1], (ftnlen)sizeof(
02327 integer));
02328 }
02329 e_wsfe();
02330 } else {
02331 i__1 = nparms;
02332 for (i__ = 1; i__ <= i__1; ++i__) {
02333 mxbval[i__ - 1] = 1;
02334
02335 }
02336 }
02337
02338
02339
02340 if (nep) {
02341 s_rsle(&io___147);
02342 i__1 = nparms;
02343 for (i__ = 1; i__ <= i__1; ++i__) {
02344 do_lio(&c__3, &c__1, (char *)&inmin[i__ - 1], (ftnlen)sizeof(
02345 integer));
02346 }
02347 e_rsle();
02348 i__1 = nparms;
02349 for (i__ = 1; i__ <= i__1; ++i__) {
02350 if (inmin[i__ - 1] < 0) {
02351 s_wsfe(&io___148);
02352 do_fio(&c__1, " INMIN ", (ftnlen)7);
02353 do_fio(&c__1, (char *)&inmin[i__ - 1], (ftnlen)sizeof(
02354 integer));
02355 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
02356 e_wsfe();
02357 fatal = TRUE_;
02358 }
02359
02360 }
02361 s_wsfe(&io___149);
02362 do_fio(&c__1, "INMIN: ", (ftnlen)7);
02363 i__1 = nparms;
02364 for (i__ = 1; i__ <= i__1; ++i__) {
02365 do_fio(&c__1, (char *)&inmin[i__ - 1], (ftnlen)sizeof(integer)
02366 );
02367 }
02368 e_wsfe();
02369 } else {
02370 i__1 = nparms;
02371 for (i__ = 1; i__ <= i__1; ++i__) {
02372 inmin[i__ - 1] = 1;
02373
02374 }
02375 }
02376
02377
02378
02379 if (nep) {
02380 s_rsle(&io___150);
02381 i__1 = nparms;
02382 for (i__ = 1; i__ <= i__1; ++i__) {
02383 do_lio(&c__3, &c__1, (char *)&inwin[i__ - 1], (ftnlen)sizeof(
02384 integer));
02385 }
02386 e_rsle();
02387 i__1 = nparms;
02388 for (i__ = 1; i__ <= i__1; ++i__) {
02389 if (inwin[i__ - 1] < 0) {
02390 s_wsfe(&io___151);
02391 do_fio(&c__1, " INWIN ", (ftnlen)7);
02392 do_fio(&c__1, (char *)&inwin[i__ - 1], (ftnlen)sizeof(
02393 integer));
02394 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
02395 e_wsfe();
02396 fatal = TRUE_;
02397 }
02398
02399 }
02400 s_wsfe(&io___152);
02401 do_fio(&c__1, "INWIN: ", (ftnlen)7);
02402 i__1 = nparms;
02403 for (i__ = 1; i__ <= i__1; ++i__) {
02404 do_fio(&c__1, (char *)&inwin[i__ - 1], (ftnlen)sizeof(integer)
02405 );
02406 }
02407 e_wsfe();
02408 } else {
02409 i__1 = nparms;
02410 for (i__ = 1; i__ <= i__1; ++i__) {
02411 inwin[i__ - 1] = 1;
02412
02413 }
02414 }
02415
02416
02417
02418 if (nep) {
02419 s_rsle(&io___153);
02420 i__1 = nparms;
02421 for (i__ = 1; i__ <= i__1; ++i__) {
02422 do_lio(&c__3, &c__1, (char *)&inibl[i__ - 1], (ftnlen)sizeof(
02423 integer));
02424 }
02425 e_rsle();
02426 i__1 = nparms;
02427 for (i__ = 1; i__ <= i__1; ++i__) {
02428 if (inibl[i__ - 1] < 0) {
02429 s_wsfe(&io___154);
02430 do_fio(&c__1, " INIBL ", (ftnlen)7);
02431 do_fio(&c__1, (char *)&inibl[i__ - 1], (ftnlen)sizeof(
02432 integer));
02433 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
02434 e_wsfe();
02435 fatal = TRUE_;
02436 }
02437
02438 }
02439 s_wsfe(&io___155);
02440 do_fio(&c__1, "INIBL: ", (ftnlen)7);
02441 i__1 = nparms;
02442 for (i__ = 1; i__ <= i__1; ++i__) {
02443 do_fio(&c__1, (char *)&inibl[i__ - 1], (ftnlen)sizeof(integer)
02444 );
02445 }
02446 e_wsfe();
02447 } else {
02448 i__1 = nparms;
02449 for (i__ = 1; i__ <= i__1; ++i__) {
02450 inibl[i__ - 1] = 1;
02451
02452 }
02453 }
02454
02455
02456
02457 if (nep) {
02458 s_rsle(&io___156);
02459 i__1 = nparms;
02460 for (i__ = 1; i__ <= i__1; ++i__) {
02461 do_lio(&c__3, &c__1, (char *)&ishfts[i__ - 1], (ftnlen)sizeof(
02462 integer));
02463 }
02464 e_rsle();
02465 i__1 = nparms;
02466 for (i__ = 1; i__ <= i__1; ++i__) {
02467 if (ishfts[i__ - 1] < 0) {
02468 s_wsfe(&io___157);
02469 do_fio(&c__1, " ISHFTS ", (ftnlen)8);
02470 do_fio(&c__1, (char *)&ishfts[i__ - 1], (ftnlen)sizeof(
02471 integer));
02472 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
02473 e_wsfe();
02474 fatal = TRUE_;
02475 }
02476
02477 }
02478 s_wsfe(&io___158);
02479 do_fio(&c__1, "ISHFTS: ", (ftnlen)8);
02480 i__1 = nparms;
02481 for (i__ = 1; i__ <= i__1; ++i__) {
02482 do_fio(&c__1, (char *)&ishfts[i__ - 1], (ftnlen)sizeof(
02483 integer));
02484 }
02485 e_wsfe();
02486 } else {
02487 i__1 = nparms;
02488 for (i__ = 1; i__ <= i__1; ++i__) {
02489 ishfts[i__ - 1] = 1;
02490
02491 }
02492 }
02493
02494
02495
02496 if (nep) {
02497 s_rsle(&io___159);
02498 i__1 = nparms;
02499 for (i__ = 1; i__ <= i__1; ++i__) {
02500 do_lio(&c__3, &c__1, (char *)&iacc22[i__ - 1], (ftnlen)sizeof(
02501 integer));
02502 }
02503 e_rsle();
02504 i__1 = nparms;
02505 for (i__ = 1; i__ <= i__1; ++i__) {
02506 if (iacc22[i__ - 1] < 0) {
02507 s_wsfe(&io___160);
02508 do_fio(&c__1, " IACC22 ", (ftnlen)8);
02509 do_fio(&c__1, (char *)&iacc22[i__ - 1], (ftnlen)sizeof(
02510 integer));
02511 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
02512 e_wsfe();
02513 fatal = TRUE_;
02514 }
02515
02516 }
02517 s_wsfe(&io___161);
02518 do_fio(&c__1, "IACC22: ", (ftnlen)8);
02519 i__1 = nparms;
02520 for (i__ = 1; i__ <= i__1; ++i__) {
02521 do_fio(&c__1, (char *)&iacc22[i__ - 1], (ftnlen)sizeof(
02522 integer));
02523 }
02524 e_wsfe();
02525 } else {
02526 i__1 = nparms;
02527 for (i__ = 1; i__ <= i__1; ++i__) {
02528 iacc22[i__ - 1] = 1;
02529
02530 }
02531 }
02532
02533
02534
02535 if (zgg) {
02536 s_rsle(&io___162);
02537 i__1 = nparms;
02538 for (i__ = 1; i__ <= i__1; ++i__) {
02539 do_lio(&c__3, &c__1, (char *)&nbcol[i__ - 1], (ftnlen)sizeof(
02540 integer));
02541 }
02542 e_rsle();
02543 i__1 = nparms;
02544 for (i__ = 1; i__ <= i__1; ++i__) {
02545 if (nbcol[i__ - 1] < 0) {
02546 s_wsfe(&io___164);
02547 do_fio(&c__1, "NBCOL ", (ftnlen)6);
02548 do_fio(&c__1, (char *)&nbcol[i__ - 1], (ftnlen)sizeof(
02549 integer));
02550 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
02551 e_wsfe();
02552 fatal = TRUE_;
02553 } else if (nbcol[i__ - 1] > 132) {
02554 s_wsfe(&io___165);
02555 do_fio(&c__1, "NBCOL ", (ftnlen)6);
02556 do_fio(&c__1, (char *)&nbcol[i__ - 1], (ftnlen)sizeof(
02557 integer));
02558 do_fio(&c__1, (char *)&c__132, (ftnlen)sizeof(integer));
02559 e_wsfe();
02560 fatal = TRUE_;
02561 }
02562
02563 }
02564 s_wsfe(&io___166);
02565 do_fio(&c__1, "NBCOL:", (ftnlen)6);
02566 i__1 = nparms;
02567 for (i__ = 1; i__ <= i__1; ++i__) {
02568 do_fio(&c__1, (char *)&nbcol[i__ - 1], (ftnlen)sizeof(integer)
02569 );
02570 }
02571 e_wsfe();
02572 } else {
02573 i__1 = nparms;
02574 for (i__ = 1; i__ <= i__1; ++i__) {
02575 nbcol[i__ - 1] = 1;
02576
02577 }
02578 }
02579 }
02580
02581
02582
02583 s_wsle(&io___167);
02584 e_wsle();
02585 eps = dlamch_("Underflow threshold");
02586 s_wsfe(&io___169);
02587 do_fio(&c__1, "underflow", (ftnlen)9);
02588 do_fio(&c__1, (char *)&eps, (ftnlen)sizeof(doublereal));
02589 e_wsfe();
02590 eps = dlamch_("Overflow threshold");
02591 s_wsfe(&io___170);
02592 do_fio(&c__1, "overflow ", (ftnlen)9);
02593 do_fio(&c__1, (char *)&eps, (ftnlen)sizeof(doublereal));
02594 e_wsfe();
02595 eps = dlamch_("Epsilon");
02596 s_wsfe(&io___171);
02597 do_fio(&c__1, "precision", (ftnlen)9);
02598 do_fio(&c__1, (char *)&eps, (ftnlen)sizeof(doublereal));
02599 e_wsfe();
02600
02601
02602
02603 s_rsle(&io___172);
02604 do_lio(&c__5, &c__1, (char *)&thresh, (ftnlen)sizeof(doublereal));
02605 e_rsle();
02606 s_wsfe(&io___173);
02607 do_fio(&c__1, (char *)&thresh, (ftnlen)sizeof(doublereal));
02608 e_wsfe();
02609 if (sep || svd || zgg) {
02610
02611
02612
02613 s_rsle(&io___174);
02614 do_lio(&c__8, &c__1, (char *)&tstchk, (ftnlen)sizeof(logical));
02615 e_rsle();
02616
02617
02618
02619 s_rsle(&io___176);
02620 do_lio(&c__8, &c__1, (char *)&tstdrv, (ftnlen)sizeof(logical));
02621 e_rsle();
02622 }
02623
02624
02625
02626 s_rsle(&io___178);
02627 do_lio(&c__8, &c__1, (char *)&tsterr, (ftnlen)sizeof(logical));
02628 e_rsle();
02629
02630
02631
02632 s_rsle(&io___179);
02633 do_lio(&c__3, &c__1, (char *)&newsd, (ftnlen)sizeof(integer));
02634 e_rsle();
02635
02636
02637
02638 if (newsd == 2) {
02639 s_rsle(&io___181);
02640 for (i__ = 1; i__ <= 4; ++i__) {
02641 do_lio(&c__3, &c__1, (char *)&ioldsd[i__ - 1], (ftnlen)sizeof(
02642 integer));
02643 }
02644 e_rsle();
02645 }
02646
02647 for (i__ = 1; i__ <= 4; ++i__) {
02648 iseed[i__ - 1] = ioldsd[i__ - 1];
02649
02650 }
02651
02652 if (fatal) {
02653 s_wsfe(&io___183);
02654 e_wsfe();
02655 s_stop("", (ftnlen)0);
02656 }
02657
02658
02659
02660
02661
02662
02663 L190:
02664
02665 if (! (zgx || zxv)) {
02666
02667 L200:
02668 ci__1.cierr = 0;
02669 ci__1.ciend = 1;
02670 ci__1.ciunit = 5;
02671 ci__1.cifmt = "(A80)";
02672 i__1 = s_rsfe(&ci__1);
02673 if (i__1 != 0) {
02674 goto L380;
02675 }
02676 i__1 = do_fio(&c__1, line, (ftnlen)80);
02677 if (i__1 != 0) {
02678 goto L380;
02679 }
02680 i__1 = e_rsfe();
02681 if (i__1 != 0) {
02682 goto L380;
02683 }
02684 s_copy(c3, line, (ftnlen)3, (ftnlen)3);
02685 lenp = i_len(line, (ftnlen)80);
02686 i__ = 3;
02687 itmp = 0;
02688 i1 = 0;
02689 L210:
02690 ++i__;
02691 if (i__ > lenp) {
02692 if (i1 > 0) {
02693 goto L240;
02694 } else {
02695 ntypes = 30;
02696 goto L240;
02697 }
02698 }
02699 if (*(unsigned char *)&line[i__ - 1] != ' ' && *(unsigned char *)&
02700 line[i__ - 1] != ',') {
02701 i1 = i__;
02702 *(unsigned char *)c1 = *(unsigned char *)&line[i1 - 1];
02703
02704
02705
02706 for (k = 1; k <= 10; ++k) {
02707 if (*(unsigned char *)c1 == *(unsigned char *)&intstr[k - 1])
02708 {
02709 ic = k - 1;
02710 goto L230;
02711 }
02712
02713 }
02714 s_wsfe(&io___192);
02715 do_fio(&c__1, (char *)&i__, (ftnlen)sizeof(integer));
02716 do_fio(&c__1, line, (ftnlen)80);
02717 e_wsfe();
02718 goto L200;
02719 L230:
02720 itmp = itmp * 10 + ic;
02721 goto L210;
02722 } else if (i1 > 0) {
02723 goto L240;
02724 } else {
02725 goto L210;
02726 }
02727 L240:
02728 ntypes = itmp;
02729
02730
02731
02732 if (! (zev || zes || zvx || zsx || zgv || zgs) && ntypes <= 0) {
02733 s_wsfe(&io___193);
02734 do_fio(&c__1, c3, (ftnlen)3);
02735 e_wsfe();
02736 goto L200;
02737 }
02738
02739 } else {
02740 if (zgx) {
02741 s_copy(c3, "ZGX", (ftnlen)3, (ftnlen)3);
02742 }
02743 if (zxv) {
02744 s_copy(c3, "ZXV", (ftnlen)3, (ftnlen)3);
02745 }
02746 }
02747
02748
02749
02750 if (newsd == 0) {
02751 for (k = 1; k <= 4; ++k) {
02752 iseed[k - 1] = ioldsd[k - 1];
02753
02754 }
02755 }
02756
02757 if (lsamen_(&c__3, c3, "ZHS") || lsamen_(&c__3, c3,
02758 "NEP")) {
02759
02760
02761
02762
02763
02764
02765
02766
02767
02768
02769
02770 maxtyp = 21;
02771 ntypes = min(maxtyp,ntypes);
02772 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
02773 xlaenv_(&c__1, &c__1);
02774 if (tsterr) {
02775 zerrhs_("ZHSEQR", &c__6);
02776 }
02777 i__1 = nparms;
02778 for (i__ = 1; i__ <= i__1; ++i__) {
02779 xlaenv_(&c__1, &nbval[i__ - 1]);
02780 xlaenv_(&c__2, &nbmin[i__ - 1]);
02781 xlaenv_(&c__3, &nxval[i__ - 1]);
02782
02783 i__3 = 11, i__4 = inmin[i__ - 1];
02784 i__2 = max(i__3,i__4);
02785 xlaenv_(&c__12, &i__2);
02786 xlaenv_(&c__13, &inwin[i__ - 1]);
02787 xlaenv_(&c__14, &inibl[i__ - 1]);
02788 xlaenv_(&c__15, &ishfts[i__ - 1]);
02789 xlaenv_(&c__16, &iacc22[i__ - 1]);
02790
02791 if (newsd == 0) {
02792 for (k = 1; k <= 4; ++k) {
02793 iseed[k - 1] = ioldsd[k - 1];
02794
02795 }
02796 }
02797 s_wsfe(&io___196);
02798 do_fio(&c__1, c3, (ftnlen)3);
02799 do_fio(&c__1, (char *)&nbval[i__ - 1], (ftnlen)sizeof(integer));
02800 do_fio(&c__1, (char *)&nbmin[i__ - 1], (ftnlen)sizeof(integer));
02801 do_fio(&c__1, (char *)&nxval[i__ - 1], (ftnlen)sizeof(integer));
02802
02803 i__3 = 11, i__4 = inmin[i__ - 1];
02804 i__2 = max(i__3,i__4);
02805 do_fio(&c__1, (char *)&i__2, (ftnlen)sizeof(integer));
02806 do_fio(&c__1, (char *)&inwin[i__ - 1], (ftnlen)sizeof(integer));
02807 do_fio(&c__1, (char *)&inibl[i__ - 1], (ftnlen)sizeof(integer));
02808 do_fio(&c__1, (char *)&ishfts[i__ - 1], (ftnlen)sizeof(integer));
02809 do_fio(&c__1, (char *)&iacc22[i__ - 1], (ftnlen)sizeof(integer));
02810 e_wsfe();
02811 zchkhs_(&nn, nval, &maxtyp, dotype, iseed, &thresh, &c__6, a, &
02812 c__132, &a[17424], &a[34848], &a[52272], &a[69696], &
02813 c__132, &a[87120], &a[104544], dc, &dc[132], &a[121968], &
02814 a[139392], &a[156816], &a[174240], &a[191664], &dc[264],
02815 work, &c__89760, rwork, iwork, logwrk, result, &info);
02816 if (info != 0) {
02817 s_wsfe(&io___205);
02818 do_fio(&c__1, "ZCHKHS", (ftnlen)6);
02819 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
02820 e_wsfe();
02821 }
02822
02823 }
02824
02825 } else if (lsamen_(&c__3, c3, "ZST") || lsamen_(&
02826 c__3, c3, "SEP")) {
02827
02828
02829
02830
02831
02832
02833
02834
02835
02836 maxtyp = 21;
02837 ntypes = min(maxtyp,ntypes);
02838 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
02839 xlaenv_(&c__1, &c__1);
02840 xlaenv_(&c__9, &c__25);
02841 if (tsterr) {
02842 zerrst_("ZST", &c__6);
02843 }
02844 i__1 = nparms;
02845 for (i__ = 1; i__ <= i__1; ++i__) {
02846 xlaenv_(&c__1, &nbval[i__ - 1]);
02847 xlaenv_(&c__2, &nbmin[i__ - 1]);
02848 xlaenv_(&c__3, &nxval[i__ - 1]);
02849
02850 if (newsd == 0) {
02851 for (k = 1; k <= 4; ++k) {
02852 iseed[k - 1] = ioldsd[k - 1];
02853
02854 }
02855 }
02856 s_wsfe(&io___206);
02857 do_fio(&c__1, c3, (ftnlen)3);
02858 do_fio(&c__1, (char *)&nbval[i__ - 1], (ftnlen)sizeof(integer));
02859 do_fio(&c__1, (char *)&nbmin[i__ - 1], (ftnlen)sizeof(integer));
02860 do_fio(&c__1, (char *)&nxval[i__ - 1], (ftnlen)sizeof(integer));
02861 e_wsfe();
02862 if (tstchk) {
02863 zchkst_(&nn, nval, &maxtyp, dotype, iseed, &thresh, &c__6, a,
02864 &c__132, &a[17424], dr, &dr[132], &dr[264], &dr[396],
02865 &dr[528], &dr[660], &dr[792], &dr[924], &dr[1056], &
02866 dr[1188], &dr[1320], &a[34848], &c__132, &a[52272], &
02867 a[69696], dc, &a[87120], work, &c__89760, rwork, &
02868 c__89760, iwork, &c__20064, result, &info);
02869 if (info != 0) {
02870 s_wsfe(&io___208);
02871 do_fio(&c__1, "ZCHKST", (ftnlen)6);
02872 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
02873 e_wsfe();
02874 }
02875 }
02876 if (tstdrv) {
02877 zdrvst_(&nn, nval, &c__18, dotype, iseed, &thresh, &c__6, a, &
02878 c__132, &dr[264], &dr[396], &dr[528], &dr[924], &dr[
02879 1056], &dr[1188], &a[17424], &c__132, &a[34848], dc, &
02880 a[52272], work, &c__89760, rwork, &c__89760, iwork, &
02881 c__20064, result, &info);
02882 if (info != 0) {
02883 s_wsfe(&io___209);
02884 do_fio(&c__1, "ZDRVST", (ftnlen)6);
02885 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
02886 e_wsfe();
02887 }
02888 }
02889
02890 }
02891
02892 } else if (lsamen_(&c__3, c3, "ZSG")) {
02893
02894
02895
02896
02897
02898
02899
02900
02901
02902 maxtyp = 21;
02903 ntypes = min(maxtyp,ntypes);
02904 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
02905 xlaenv_(&c__9, &c__25);
02906 i__1 = nparms;
02907 for (i__ = 1; i__ <= i__1; ++i__) {
02908 xlaenv_(&c__1, &nbval[i__ - 1]);
02909 xlaenv_(&c__2, &nbmin[i__ - 1]);
02910 xlaenv_(&c__3, &nxval[i__ - 1]);
02911
02912 if (newsd == 0) {
02913 for (k = 1; k <= 4; ++k) {
02914 iseed[k - 1] = ioldsd[k - 1];
02915
02916 }
02917 }
02918 s_wsfe(&io___210);
02919 do_fio(&c__1, c3, (ftnlen)3);
02920 do_fio(&c__1, (char *)&nbval[i__ - 1], (ftnlen)sizeof(integer));
02921 do_fio(&c__1, (char *)&nbmin[i__ - 1], (ftnlen)sizeof(integer));
02922 do_fio(&c__1, (char *)&nxval[i__ - 1], (ftnlen)sizeof(integer));
02923 e_wsfe();
02924 if (tstchk) {
02925 zdrvsg_(&nn, nval, &maxtyp, dotype, iseed, &thresh, &c__6, a,
02926 &c__132, &a[17424], &c__132, &dr[264], &a[34848], &
02927 c__132, &a[52272], &a[69696], &a[87120], &a[104544],
02928 work, &c__89760, rwork, &c__89760, iwork, &c__20064,
02929 result, &info);
02930 if (info != 0) {
02931 s_wsfe(&io___211);
02932 do_fio(&c__1, "ZDRVSG", (ftnlen)6);
02933 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
02934 e_wsfe();
02935 }
02936 }
02937
02938 }
02939
02940 } else if (lsamen_(&c__3, c3, "ZBD") || lsamen_(&
02941 c__3, c3, "SVD")) {
02942
02943
02944
02945
02946
02947
02948
02949
02950
02951
02952 maxtyp = 16;
02953 ntypes = min(maxtyp,ntypes);
02954 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
02955 xlaenv_(&c__9, &c__25);
02956
02957
02958
02959 xlaenv_(&c__1, &c__1);
02960 if (tsterr && tstchk) {
02961 zerrbd_("ZBD", &c__6);
02962 }
02963 if (tsterr && tstdrv) {
02964 zerred_("ZBD", &c__6);
02965 }
02966
02967 i__1 = nparms;
02968 for (i__ = 1; i__ <= i__1; ++i__) {
02969 nrhs = nsval[i__ - 1];
02970 xlaenv_(&c__1, &nbval[i__ - 1]);
02971 xlaenv_(&c__2, &nbmin[i__ - 1]);
02972 xlaenv_(&c__3, &nxval[i__ - 1]);
02973 if (newsd == 0) {
02974 for (k = 1; k <= 4; ++k) {
02975 iseed[k - 1] = ioldsd[k - 1];
02976
02977 }
02978 }
02979 s_wsfe(&io___213);
02980 do_fio(&c__1, c3, (ftnlen)3);
02981 do_fio(&c__1, (char *)&nbval[i__ - 1], (ftnlen)sizeof(integer));
02982 do_fio(&c__1, (char *)&nbmin[i__ - 1], (ftnlen)sizeof(integer));
02983 do_fio(&c__1, (char *)&nxval[i__ - 1], (ftnlen)sizeof(integer));
02984 do_fio(&c__1, (char *)&nrhs, (ftnlen)sizeof(integer));
02985 e_wsfe();
02986 if (tstchk) {
02987 zchkbd_(&nn, mval, nval, &maxtyp, dotype, &nrhs, iseed, &
02988 thresh, a, &c__132, dr, &dr[132], &dr[264], &dr[396],
02989 &a[17424], &c__132, &a[34848], &a[52272], &a[69696], &
02990 c__132, &a[87120], &c__132, &a[104544], &a[121968],
02991 work, &c__89760, rwork, &c__6, &info);
02992 if (info != 0) {
02993 s_wsfe(&io___214);
02994 do_fio(&c__1, "ZCHKBD", (ftnlen)6);
02995 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
02996 e_wsfe();
02997 }
02998 }
02999 if (tstdrv) {
03000 zdrvbd_(&nn, mval, nval, &maxtyp, dotype, iseed, &thresh, a, &
03001 c__132, &a[17424], &c__132, &a[34848], &c__132, &a[
03002 52272], &a[69696], &a[87120], dr, &dr[132], &dr[264],
03003 work, &c__89760, rwork, iwork, &c__6, &info);
03004 }
03005
03006 }
03007
03008 } else if (lsamen_(&c__3, c3, "ZEV")) {
03009
03010
03011
03012
03013
03014
03015 maxtyp = 21;
03016 ntypes = min(maxtyp,ntypes);
03017 if (ntypes <= 0) {
03018 s_wsfe(&io___215);
03019 do_fio(&c__1, c3, (ftnlen)3);
03020 e_wsfe();
03021 } else {
03022 if (tsterr) {
03023 zerred_(c3, &c__6);
03024 }
03025 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
03026 zdrvev_(&nn, nval, &ntypes, dotype, iseed, &thresh, &c__6, a, &
03027 c__132, &a[17424], dc, &dc[132], &a[34848], &c__132, &a[
03028 52272], &c__132, &a[69696], &c__132, result, work, &
03029 c__89760, rwork, iwork, &info);
03030 if (info != 0) {
03031 s_wsfe(&io___216);
03032 do_fio(&c__1, "ZGEEV", (ftnlen)5);
03033 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03034 e_wsfe();
03035 }
03036 }
03037 s_wsfe(&io___217);
03038 e_wsfe();
03039 goto L10;
03040
03041 } else if (lsamen_(&c__3, c3, "ZES")) {
03042
03043
03044
03045
03046
03047
03048 maxtyp = 21;
03049 ntypes = min(maxtyp,ntypes);
03050 if (ntypes <= 0) {
03051 s_wsfe(&io___218);
03052 do_fio(&c__1, c3, (ftnlen)3);
03053 e_wsfe();
03054 } else {
03055 if (tsterr) {
03056 zerred_(c3, &c__6);
03057 }
03058 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
03059 zdrves_(&nn, nval, &ntypes, dotype, iseed, &thresh, &c__6, a, &
03060 c__132, &a[17424], &a[34848], dc, &dc[132], &a[52272], &
03061 c__132, result, work, &c__89760, rwork, iwork, logwrk, &
03062 info);
03063 if (info != 0) {
03064 s_wsfe(&io___219);
03065 do_fio(&c__1, "ZGEES", (ftnlen)5);
03066 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03067 e_wsfe();
03068 }
03069 }
03070 s_wsfe(&io___220);
03071 e_wsfe();
03072 goto L10;
03073
03074 } else if (lsamen_(&c__3, c3, "ZVX")) {
03075
03076
03077
03078
03079
03080
03081 maxtyp = 21;
03082 ntypes = min(maxtyp,ntypes);
03083 if (ntypes < 0) {
03084 s_wsfe(&io___221);
03085 do_fio(&c__1, c3, (ftnlen)3);
03086 e_wsfe();
03087 } else {
03088 if (tsterr) {
03089 zerred_(c3, &c__6);
03090 }
03091 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
03092 zdrvvx_(&nn, nval, &ntypes, dotype, iseed, &thresh, &c__5, &c__6,
03093 a, &c__132, &a[17424], dc, &dc[132], &a[34848], &c__132, &
03094 a[52272], &c__132, &a[69696], &c__132, dr, &dr[132], &dr[
03095 264], &dr[396], &dr[528], &dr[660], &dr[792], &dr[924],
03096 result, work, &c__89760, rwork, &info);
03097 if (info != 0) {
03098 s_wsfe(&io___222);
03099 do_fio(&c__1, "ZGEEVX", (ftnlen)6);
03100 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03101 e_wsfe();
03102 }
03103 }
03104 s_wsfe(&io___223);
03105 e_wsfe();
03106 goto L10;
03107
03108 } else if (lsamen_(&c__3, c3, "ZSX")) {
03109
03110
03111
03112
03113
03114
03115 maxtyp = 21;
03116 ntypes = min(maxtyp,ntypes);
03117 if (ntypes < 0) {
03118 s_wsfe(&io___224);
03119 do_fio(&c__1, c3, (ftnlen)3);
03120 e_wsfe();
03121 } else {
03122 if (tsterr) {
03123 zerred_(c3, &c__6);
03124 }
03125 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
03126 zdrvsx_(&nn, nval, &ntypes, dotype, iseed, &thresh, &c__5, &c__6,
03127 a, &c__132, &a[17424], &a[34848], dc, &dc[132], &dc[264],
03128 &a[52272], &c__132, &a[69696], result, work, &c__89760,
03129 rwork, logwrk, &info);
03130 if (info != 0) {
03131 s_wsfe(&io___225);
03132 do_fio(&c__1, "ZGEESX", (ftnlen)6);
03133 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03134 e_wsfe();
03135 }
03136 }
03137 s_wsfe(&io___226);
03138 e_wsfe();
03139 goto L10;
03140
03141 } else if (lsamen_(&c__3, c3, "ZGG")) {
03142
03143
03144
03145
03146
03147
03148
03149
03150
03151
03152
03153 maxtyp = 26;
03154 ntypes = min(maxtyp,ntypes);
03155 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
03156 if (tstchk && tsterr) {
03157 zerrgg_(c3, &c__6);
03158 }
03159 i__1 = nparms;
03160 for (i__ = 1; i__ <= i__1; ++i__) {
03161 xlaenv_(&c__1, &nbval[i__ - 1]);
03162 xlaenv_(&c__2, &nbmin[i__ - 1]);
03163 xlaenv_(&c__4, &nsval[i__ - 1]);
03164 xlaenv_(&c__8, &mxbval[i__ - 1]);
03165 xlaenv_(&c__5, &nbcol[i__ - 1]);
03166
03167 if (newsd == 0) {
03168 for (k = 1; k <= 4; ++k) {
03169 iseed[k - 1] = ioldsd[k - 1];
03170
03171 }
03172 }
03173 s_wsfe(&io___227);
03174 do_fio(&c__1, c3, (ftnlen)3);
03175 do_fio(&c__1, (char *)&nbval[i__ - 1], (ftnlen)sizeof(integer));
03176 do_fio(&c__1, (char *)&nbmin[i__ - 1], (ftnlen)sizeof(integer));
03177 do_fio(&c__1, (char *)&nsval[i__ - 1], (ftnlen)sizeof(integer));
03178 do_fio(&c__1, (char *)&mxbval[i__ - 1], (ftnlen)sizeof(integer));
03179 do_fio(&c__1, (char *)&nbcol[i__ - 1], (ftnlen)sizeof(integer));
03180 e_wsfe();
03181 tstdif = FALSE_;
03182 thrshn = 10.;
03183 if (tstchk) {
03184 zchkgg_(&nn, nval, &maxtyp, dotype, iseed, &thresh, &tstdif, &
03185 thrshn, &c__6, a, &c__132, &a[17424], &a[34848], &a[
03186 52272], &a[69696], &a[87120], &a[104544], &a[121968],
03187 &a[139392], &c__132, &a[156816], &a[174240], &a[
03188 191664], dc, &dc[132], &dc[264], &dc[396], &a[209088],
03189 &a[226512], work, &c__89760, rwork, logwrk, result, &
03190 info);
03191 if (info != 0) {
03192 s_wsfe(&io___230);
03193 do_fio(&c__1, "ZCHKGG", (ftnlen)6);
03194 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03195 e_wsfe();
03196 }
03197 }
03198 xlaenv_(&c__1, &c__1);
03199 if (tstdrv) {
03200 zdrvgg_(&nn, nval, &maxtyp, dotype, iseed, &thresh, &thrshn, &
03201 c__6, a, &c__132, &a[17424], &a[34848], &a[52272], &a[
03202 69696], &a[87120], &a[104544], &c__132, &a[121968],
03203 dc, &dc[132], &dc[264], &dc[396], &a[121968], &a[
03204 139392], work, &c__89760, rwork, result, &info);
03205 if (info != 0) {
03206 s_wsfe(&io___231);
03207 do_fio(&c__1, "ZDRVGG", (ftnlen)6);
03208 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03209 e_wsfe();
03210 }
03211 }
03212
03213 }
03214
03215 } else if (lsamen_(&c__3, c3, "ZGS")) {
03216
03217
03218
03219
03220
03221
03222 maxtyp = 26;
03223 ntypes = min(maxtyp,ntypes);
03224 if (ntypes <= 0) {
03225 s_wsfe(&io___232);
03226 do_fio(&c__1, c3, (ftnlen)3);
03227 e_wsfe();
03228 } else {
03229 if (tsterr) {
03230 zerrgg_(c3, &c__6);
03231 }
03232 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
03233 zdrges_(&nn, nval, &maxtyp, dotype, iseed, &thresh, &c__6, a, &
03234 c__132, &a[17424], &a[34848], &a[52272], &a[104544], &
03235 c__132, &a[121968], dc, &dc[132], work, &c__89760, rwork,
03236 result, logwrk, &info);
03237
03238 if (info != 0) {
03239 s_wsfe(&io___233);
03240 do_fio(&c__1, "ZDRGES", (ftnlen)6);
03241 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03242 e_wsfe();
03243 }
03244 }
03245 s_wsfe(&io___234);
03246 e_wsfe();
03247 goto L10;
03248
03249 } else if (zgx) {
03250
03251
03252
03253
03254
03255
03256 maxtyp = 5;
03257 ntypes = maxtyp;
03258 if (nn < 0) {
03259 s_wsfe(&io___235);
03260 do_fio(&c__1, c3, (ftnlen)3);
03261 e_wsfe();
03262 } else {
03263 if (tsterr) {
03264 zerrgg_(c3, &c__6);
03265 }
03266 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
03267 xlaenv_(&c__5, &c__2);
03268 zdrgsx_(&nn, &c__20, &thresh, &c__5, &c__6, a, &c__132, &a[17424],
03269 &a[34848], &a[52272], &a[69696], &a[87120], dc, &dc[132],
03270 c__, &c__400, s, work, &c__89760, rwork, iwork, &
03271 c__20064, logwrk, &info);
03272 if (info != 0) {
03273 s_wsfe(&io___238);
03274 do_fio(&c__1, "ZDRGSX", (ftnlen)6);
03275 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03276 e_wsfe();
03277 }
03278 }
03279 s_wsfe(&io___239);
03280 e_wsfe();
03281 goto L10;
03282
03283 } else if (lsamen_(&c__3, c3, "ZGV")) {
03284
03285
03286
03287
03288
03289
03290 maxtyp = 26;
03291 ntypes = min(maxtyp,ntypes);
03292 if (ntypes <= 0) {
03293 s_wsfe(&io___240);
03294 do_fio(&c__1, c3, (ftnlen)3);
03295 e_wsfe();
03296 } else {
03297 if (tsterr) {
03298 zerrgg_(c3, &c__6);
03299 }
03300 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
03301 zdrgev_(&nn, nval, &maxtyp, dotype, iseed, &thresh, &c__6, a, &
03302 c__132, &a[17424], &a[34848], &a[52272], &a[104544], &
03303 c__132, &a[121968], &a[139392], &c__132, dc, &dc[132], &
03304 dc[264], &dc[396], work, &c__89760, rwork, result, &info);
03305 if (info != 0) {
03306 s_wsfe(&io___241);
03307 do_fio(&c__1, "ZDRGEV", (ftnlen)6);
03308 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03309 e_wsfe();
03310 }
03311 }
03312 s_wsfe(&io___242);
03313 e_wsfe();
03314 goto L10;
03315
03316 } else if (zxv) {
03317
03318
03319
03320
03321
03322
03323 maxtyp = 2;
03324 ntypes = maxtyp;
03325 if (nn < 0) {
03326 s_wsfe(&io___243);
03327 do_fio(&c__1, c3, (ftnlen)3);
03328 e_wsfe();
03329 } else {
03330 if (tsterr) {
03331 zerrgg_(c3, &c__6);
03332 }
03333 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
03334 zdrgvx_(&nn, &thresh, &c__5, &c__6, a, &c__132, &a[17424], &a[
03335 34848], &a[52272], dc, &dc[132], &a[69696], &a[87120],
03336 iwork, &iwork[1], dr, &dr[132], &dr[264], &dr[396], &dr[
03337 528], &dr[660], work, &c__89760, rwork, &iwork[2], &
03338 c__20062, result, logwrk, &info);
03339
03340 if (info != 0) {
03341 s_wsfe(&io___244);
03342 do_fio(&c__1, "ZDRGVX", (ftnlen)6);
03343 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03344 e_wsfe();
03345 }
03346 }
03347 s_wsfe(&io___245);
03348 e_wsfe();
03349 goto L10;
03350
03351 } else if (lsamen_(&c__3, c3, "ZHB")) {
03352
03353
03354
03355
03356
03357 maxtyp = 15;
03358 ntypes = min(maxtyp,ntypes);
03359 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
03360 if (tsterr) {
03361 zerrst_("ZHB", &c__6);
03362 }
03363 zchkhb_(&nn, nval, &nk, kval, &maxtyp, dotype, iseed, &thresh, &c__6,
03364 a, &c__132, dr, &dr[132], &a[17424], &c__132, work, &c__89760,
03365 rwork, result, &info);
03366 if (info != 0) {
03367 s_wsfe(&io___246);
03368 do_fio(&c__1, "ZCHKHB", (ftnlen)6);
03369 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03370 e_wsfe();
03371 }
03372
03373 } else if (lsamen_(&c__3, c3, "ZBB")) {
03374
03375
03376
03377
03378
03379 maxtyp = 15;
03380 ntypes = min(maxtyp,ntypes);
03381 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
03382 i__1 = nparms;
03383 for (i__ = 1; i__ <= i__1; ++i__) {
03384 nrhs = nsval[i__ - 1];
03385
03386 if (newsd == 0) {
03387 for (k = 1; k <= 4; ++k) {
03388 iseed[k - 1] = ioldsd[k - 1];
03389
03390 }
03391 }
03392 s_wsfe(&io___247);
03393 do_fio(&c__1, c3, (ftnlen)3);
03394 do_fio(&c__1, (char *)&nrhs, (ftnlen)sizeof(integer));
03395 e_wsfe();
03396 zchkbb_(&nn, mval, nval, &nk, kval, &maxtyp, dotype, &nrhs, iseed,
03397 &thresh, &c__6, a, &c__132, &a[17424], &c__264, dr, &dr[
03398 132], &a[52272], &c__132, &a[69696], &c__132, &a[87120], &
03399 c__132, &a[104544], work, &c__89760, rwork, result, &info)
03400 ;
03401 if (info != 0) {
03402 s_wsfe(&io___248);
03403 do_fio(&c__1, "ZCHKBB", (ftnlen)6);
03404 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03405 e_wsfe();
03406 }
03407
03408 }
03409
03410 } else if (lsamen_(&c__3, c3, "GLM")) {
03411
03412
03413
03414
03415
03416 xlaenv_(&c__1, &c__1);
03417 if (tsterr) {
03418 zerrgg_("GLM", &c__6);
03419 }
03420 zckglm_(&nn, nval, mval, pval, &ntypes, iseed, &thresh, &c__132, a, &
03421 a[17424], b, &b[17424], x, work, dr, &c__5, &c__6, &info);
03422 if (info != 0) {
03423 s_wsfe(&io___251);
03424 do_fio(&c__1, "ZCKGLM", (ftnlen)6);
03425 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03426 e_wsfe();
03427 }
03428
03429 } else if (lsamen_(&c__3, c3, "GQR")) {
03430
03431
03432
03433
03434
03435 xlaenv_(&c__1, &c__1);
03436 if (tsterr) {
03437 zerrgg_("GQR", &c__6);
03438 }
03439 zckgqr_(&nn, mval, &nn, pval, &nn, nval, &ntypes, iseed, &thresh, &
03440 c__132, a, &a[17424], &a[34848], &a[52272], taua, b, &b[17424]
03441 , &b[34848], &b[52272], &b[69696], taub, work, dr, &c__5, &
03442 c__6, &info);
03443 if (info != 0) {
03444 s_wsfe(&io___254);
03445 do_fio(&c__1, "ZCKGQR", (ftnlen)6);
03446 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03447 e_wsfe();
03448 }
03449
03450 } else if (lsamen_(&c__3, c3, "GSV")) {
03451
03452
03453
03454
03455
03456 if (tsterr) {
03457 zerrgg_("GSV", &c__6);
03458 }
03459 zckgsv_(&nn, mval, pval, nval, &ntypes, iseed, &thresh, &c__132, a, &
03460 a[17424], b, &b[17424], &a[34848], &b[34848], &a[52272],
03461 alpha, beta, &b[52272], iwork, work, dr, &c__5, &c__6, &info);
03462 if (info != 0) {
03463 s_wsfe(&io___257);
03464 do_fio(&c__1, "ZCKGSV", (ftnlen)6);
03465 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03466 e_wsfe();
03467 }
03468
03469 } else if (lsamen_(&c__3, c3, "LSE")) {
03470
03471
03472
03473
03474
03475 xlaenv_(&c__1, &c__1);
03476 if (tsterr) {
03477 zerrgg_("LSE", &c__6);
03478 }
03479 zcklse_(&nn, mval, pval, nval, &ntypes, iseed, &thresh, &c__132, a, &
03480 a[17424], b, &b[17424], x, work, dr, &c__5, &c__6, &info);
03481 if (info != 0) {
03482 s_wsfe(&io___258);
03483 do_fio(&c__1, "ZCKLSE", (ftnlen)6);
03484 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03485 e_wsfe();
03486 }
03487 } else {
03488 s_wsle(&io___259);
03489 e_wsle();
03490 s_wsle(&io___260);
03491 e_wsle();
03492 s_wsfe(&io___261);
03493 do_fio(&c__1, c3, (ftnlen)3);
03494 e_wsfe();
03495 }
03496 if (! (zgx || zxv)) {
03497 goto L190;
03498 }
03499 L380:
03500 s_wsfe(&io___262);
03501 e_wsfe();
03502 s2 = dsecnd_();
03503 s_wsfe(&io___264);
03504 d__1 = s2 - s1;
03505 do_fio(&c__1, (char *)&d__1, (ftnlen)sizeof(doublereal));
03506 e_wsfe();
03507
03508
03509
03510
03511
03512 return 0;
03513 }
03514
03515 int zchkee_ () { MAIN__ (); return 0; }