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 iparms[100];
00026 } claenv_;
00027
00028 #define claenv_1 claenv_
00029
00030 struct {
00031 integer infot, nunit;
00032 logical ok, lerr;
00033 } infoc_;
00034
00035 #define infoc_1 infoc_
00036
00037 struct {
00038 char srnamt[32];
00039 } srnamc_;
00040
00041 #define srnamc_1 srnamc_
00042
00043 struct {
00044 integer selopt, seldim;
00045 logical selval[20];
00046 real selwr[20], selwi[20];
00047 } sslct_;
00048
00049 #define sslct_1 sslct_
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__4 = 4;
00058 static integer c__12 = 12;
00059 static integer c__11 = 11;
00060 static integer c__13 = 13;
00061 static integer c__2 = 2;
00062 static integer c__14 = 14;
00063 static integer c__0 = 0;
00064 static integer c__15 = 15;
00065 static integer c__16 = 16;
00066 static integer c__20 = 20;
00067 static integer c__132 = 132;
00068 static integer c__8 = 8;
00069 static integer c__87781 = 87781;
00070 static integer c__9 = 9;
00071 static integer c__25 = 25;
00072 static integer c__89760 = 89760;
00073 static integer c__18 = 18;
00074 static integer c__400 = 400;
00075 static integer c__89758 = 89758;
00076 static integer c__264 = 264;
00077
00078 int MAIN__(void)
00079 {
00080
00081
00082 static char intstr[10] = "0123456789";
00083 static integer ioldsd[4] = { 0,0,0,1 };
00084
00085
00086 static char fmt_9987[] = "(\002 Tests of the Nonsymmetric Eigenvalue Pro"
00087 "blem routines\002)";
00088 static char fmt_9986[] = "(\002 Tests of the Symmetric Eigenvalue Proble"
00089 "m routines\002)";
00090 static char fmt_9985[] = "(\002 Tests of the Singular Value Decompositio"
00091 "n routines\002)";
00092 static char fmt_9979[] = "(/\002 Tests of the Nonsymmetric Eigenvalue Pr"
00093 "oblem Driver\002,/\002 SGEEV (eigenvalues and eigevectors)"
00094 "\002)";
00095 static char fmt_9978[] = "(/\002 Tests of the Nonsymmetric Eigenvalue Pr"
00096 "oblem Driver\002,/\002 SGEES (Schur form)\002)";
00097 static char fmt_9977[] = "(/\002 Tests of the Nonsymmetric Eigenvalue Pr"
00098 "oblem Expert\002,\002 Driver\002,/\002 SGEEVX (eigenvalues, e"
00099 "igenvectors and\002,\002 condition numbers)\002)";
00100 static char fmt_9976[] = "(/\002 Tests of the Nonsymmetric Eigenvalue Pr"
00101 "oblem Expert\002,\002 Driver\002,/\002 SGEESX (Schur form and"
00102 " condition\002,\002 numbers)\002)";
00103 static char fmt_9975[] = "(/\002 Tests of the Generalized Nonsymmetric E"
00104 "igenvalue \002,\002Problem routines\002)";
00105 static char fmt_9964[] = "(/\002 Tests of the Generalized Nonsymmetric E"
00106 "igenvalue \002,\002Problem Driver SGGES\002)";
00107 static char fmt_9965[] = "(/\002 Tests of the Generalized Nonsymmetric E"
00108 "igenvalue \002,\002Problem Expert Driver SGGESX\002)";
00109 static char fmt_9963[] = "(/\002 Tests of the Generalized Nonsymmetric E"
00110 "igenvalue \002,\002Problem Driver SGGEV\002)";
00111 static char fmt_9962[] = "(/\002 Tests of the Generalized Nonsymmetric E"
00112 "igenvalue \002,\002Problem Expert Driver SGGEVX\002)";
00113 static char fmt_9974[] = "(\002 Tests of SSBTRD\002,/\002 (reduction of "
00114 "a symmetric band \002,\002matrix to tridiagonal form)\002)";
00115 static char fmt_9967[] = "(\002 Tests of SGBBRD\002,/\002 (reduction of "
00116 "a general band \002,\002matrix to real bidiagonal form)\002)";
00117 static char fmt_9971[] = "(/\002 Tests of the Generalized Linear Regress"
00118 "ion Model \002,\002routines\002)";
00119 static char fmt_9970[] = "(/\002 Tests of the Generalized QR and RQ rout"
00120 "ines\002)";
00121 static char fmt_9969[] = "(/\002 Tests of the Generalized Singular Valu"
00122 "e\002,\002 Decomposition routines\002)";
00123 static char fmt_9968[] = "(/\002 Tests of the Linear Least Squares routi"
00124 "nes\002)";
00125 static char fmt_9992[] = "(1x,a3,\002: Unrecognized path name\002)";
00126 static char fmt_9972[] = "(/\002 LAPACK VERSION \002,i1,\002.\002,i1,"
00127 "\002.\002,i1)";
00128 static char fmt_9984[] = "(/\002 The following parameter values will be "
00129 "used:\002)";
00130 static char fmt_9989[] = "(\002 Invalid input value: \002,a,\002=\002,"
00131 "i6,\002; must be >=\002,i6)";
00132 static char fmt_9988[] = "(\002 Invalid input value: \002,a,\002=\002,"
00133 "i6,\002; must be <=\002,i6)";
00134 static char fmt_9983[] = "(4x,a,10i6,/10x,10i6)";
00135 static char fmt_9981[] = "(\002 Relative machine \002,a,\002 is taken to"
00136 " be\002,e16.6)";
00137 static char fmt_9982[] = "(/\002 Routines pass computational tests if te"
00138 "st ratio is \002,\002less than\002,f8.2,/)";
00139 static char fmt_9999[] = "(/\002 Execution not attempted due to input er"
00140 "rors\002)";
00141 static char fmt_9991[] = "(//\002 *** Invalid integer value in column"
00142 " \002,i2,\002 of input\002,\002 line:\002,/a79)";
00143 static char fmt_9990[] = "(//1x,a3,\002 routines were not tested\002)";
00144 static char fmt_9961[] = "(//1x,a3,\002: NB =\002,i4,\002, NBMIN =\002,"
00145 "i4,\002, NX =\002,i4,\002, INMIN=\002,i4,\002, INWIN =\002,i4"
00146 ",\002, INIBL =\002,i4,\002, ISHFTS =\002,i4,\002, IACC22 =\002,i"
00147 "4)";
00148 static char fmt_9980[] = "(\002 *** Error code from \002,a,\002 = \002,i"
00149 "4)";
00150 static char fmt_9997[] = "(//1x,a3,\002: NB =\002,i4,\002, NBMIN =\002,"
00151 "i4,\002, NX =\002,i4)";
00152 static char fmt_9995[] = "(//1x,a3,\002: NB =\002,i4,\002, NBMIN =\002,"
00153 "i4,\002, NX =\002,i4,\002, NRHS =\002,i4)";
00154 static char fmt_9973[] = "(/1x,71(\002-\002))";
00155 static char fmt_9996[] = "(//1x,a3,\002: NB =\002,i4,\002, NBMIN =\002,"
00156 "i4,\002, NS =\002,i4,\002, MAXB =\002,i4,\002, NBCOL =\002,i4)";
00157 static char fmt_9966[] = "(//1x,a3,\002: NRHS =\002,i4)";
00158 static char fmt_9994[] = "(//\002 End of tests\002)";
00159 static char fmt_9993[] = "(\002 Total time used = \002,f12.2,\002 seco"
00160 "nds\002,/)";
00161
00162
00163 integer i__1, i__2, i__3, i__4;
00164 real r__1;
00165 cilist ci__1;
00166
00167
00168 integer s_rsfe(cilist *), do_fio(integer *, char *, ftnlen), e_rsfe(void);
00169 int s_copy(char *, char *, ftnlen, ftnlen);
00170 integer s_cmp(char *, char *, ftnlen, ftnlen), s_wsfe(cilist *), e_wsfe(
00171 void), s_rsle(cilist *), do_lio(integer *, integer *, char *,
00172 ftnlen), e_rsle(void), s_wsle(cilist *), e_wsle(void);
00173 int s_stop(char *, ftnlen);
00174 integer i_len(char *, ftnlen);
00175
00176
00177 real a[243936] , b[87120] ,
00178 c__[160000] , d__[1584] ;
00179 integer i__, k;
00180 real x[660];
00181 char c1[1], c3[3];
00182 integer i1;
00183 real s1, s2;
00184 integer ic, nk, nn, vers_patch__, vers_major__, vers_minor__;
00185 logical sbb, glm, sbk, sbl, nep, lse, sgg, sep, sgk, gqr, ses, sgl, sgs,
00186 sev, ssb, gsv, sgv, sgx, svd;
00187 real eps;
00188 logical ssx, svx, sxv;
00189 char line[80];
00190 real taua[132];
00191 integer info;
00192 char path[3];
00193 integer kval[20], lenp, mval[20], nval[20];
00194 real taub[132];
00195 integer pval[20], itmp, nrhs;
00196 real work[87781];
00197 integer iacc22[20];
00198 logical fatal;
00199 integer iseed[4], nbcol[20], inibl[20], nbval[20], nbmin[20];
00200 char vname[32];
00201 integer inmin[20], newsd, nsval[20], inwin[20], nxval[20], iwork[89760];
00202 extern int schkbb_(integer *, integer *, integer *,
00203 integer *, integer *, integer *, logical *, integer *, integer *,
00204 real *, integer *, real *, integer *, real *, integer *, real *,
00205 real *, real *, integer *, real *, integer *, real *, integer *,
00206 real *, real *, integer *, real *, integer *), schkbd_(integer *,
00207 integer *, integer *, integer *, logical *, integer *, integer *,
00208 real *, real *, integer *, real *, real *, real *, real *, real *,
00209 integer *, real *, real *, real *, integer *, real *, integer *,
00210 real *, real *, real *, integer *, integer *, integer *, integer *
00211 ), schkec_(real *, logical *, integer *, integer *), alareq_(char
00212 *, integer *, logical *, integer *, integer *, integer *),
00213 schkbk_(integer *, integer *);
00214 extern doublereal slamch_(char *);
00215 extern int schkbl_(integer *, integer *), schkgg_(
00216 integer *, integer *, integer *, logical *, integer *, real *,
00217 logical *, real *, integer *, real *, integer *, real *, real *,
00218 real *, real *, real *, real *, real *, real *, integer *, real *,
00219 real *, real *, real *, real *, real *, real *, real *, real *,
00220 real *, real *, real *, integer *, logical *, real *, integer *),
00221 schkgk_(integer *, integer *);
00222 extern doublereal second_(void);
00223 extern int schkgl_(integer *, integer *), schksb_(
00224 integer *, integer *, integer *, integer *, integer *, logical *,
00225 integer *, real *, integer *, real *, integer *, real *, real *,
00226 real *, integer *, real *, integer *, real *, integer *);
00227 extern logical lsamen_(integer *, char *, char *);
00228 extern int sckglm_(integer *, integer *, integer *,
00229 integer *, integer *, integer *, real *, integer *, real *, real *
00230 , real *, real *, real *, real *, real *, integer *, integer *,
00231 integer *), serrbd_(char *, integer *), ilaver_(integer *,
00232 integer *, integer *), schkhs_(integer *, integer *, integer *,
00233 logical *, integer *, real *, integer *, real *, integer *, real *
00234 , real *, real *, real *, integer *, real *, real *, real *, real
00235 *, real *, real *, real *, real *, real *, real *, real *, real *,
00236 real *, integer *, integer *, logical *, real *, integer *),
00237 scklse_(integer *, integer *, integer *, integer *, integer *,
00238 integer *, real *, integer *, real *, real *, real *, real *,
00239 real *, real *, real *, integer *, integer *, integer *), sdrvbd_(
00240 integer *, integer *, integer *, integer *, logical *, integer *,
00241 real *, real *, integer *, real *, integer *, real *, integer *,
00242 real *, real *, real *, real *, real *, real *, real *, integer *,
00243 integer *, integer *, integer *), serred_(char *, integer *), sdrges_(integer *, integer *, integer *, logical *,
00244 integer *, real *, integer *, real *, integer *, real *, real *,
00245 real *, real *, integer *, real *, real *, real *, real *, real *,
00246 integer *, real *, logical *, integer *);
00247 integer mxbval[20];
00248 extern int sckgqr_(integer *, integer *, integer *,
00249 integer *, integer *, integer *, integer *, integer *, real *,
00250 integer *, real *, real *, real *, real *, real *, real *, real *,
00251 real *, real *, real *, real *, real *, real *, integer *,
00252 integer *, integer *), sdrgev_(integer *, integer *, integer *,
00253 logical *, integer *, real *, integer *, real *, integer *, real *
00254 , real *, real *, real *, integer *, real *, real *, integer *,
00255 real *, real *, real *, real *, real *, real *, real *, integer *,
00256 real *, integer *), sdrvgg_(integer *, integer *, integer *,
00257 logical *, integer *, real *, real *, integer *, real *, integer *
00258 , real *, real *, real *, real *, real *, real *, integer *, real
00259 *, real *, real *, real *, real *, real *, real *, real *, real *,
00260 real *, integer *, real *, integer *);
00261 logical tstdif;
00262 real thresh;
00263 extern int schkst_(integer *, integer *, integer *,
00264 logical *, integer *, real *, integer *, real *, integer *, real *
00265 , real *, real *, real *, real *, real *, real *, real *, real *,
00266 real *, real *, real *, real *, integer *, real *, real *, real *,
00267 real *, real *, integer *, integer *, integer *, real *, integer
00268 *);
00269 logical tstchk;
00270 integer nparms, ishfts[20];
00271 extern int sckgsv_(integer *, integer *, integer *,
00272 integer *, integer *, integer *, real *, integer *, real *, real *
00273 , real *, real *, real *, real *, real *, real *, real *, real *,
00274 integer *, real *, real *, integer *, integer *, integer *);
00275 logical dotype[30], logwrk[132];
00276 real thrshn;
00277 extern int sdrves_(integer *, integer *, integer *,
00278 logical *, integer *, real *, integer *, real *, integer *, real *
00279 , real *, real *, real *, real *, real *, real *, integer *, real
00280 *, real *, integer *, integer *, logical *, integer *), sdrvev_(
00281 integer *, integer *, integer *, logical *, integer *, real *,
00282 integer *, real *, integer *, real *, real *, real *, real *,
00283 real *, real *, integer *, real *, integer *, real *, integer *,
00284 real *, real *, integer *, integer *, integer *), sdrgsx_(integer
00285 *, integer *, real *, integer *, integer *, real *, integer *,
00286 real *, real *, real *, real *, real *, real *, real *, real *,
00287 real *, integer *, real *, real *, integer *, integer *, integer *
00288 , logical *, integer *), sdrvsg_(integer *, integer *, integer *,
00289 logical *, integer *, real *, integer *, real *, integer *, real *
00290 , integer *, real *, real *, integer *, real *, real *, real *,
00291 real *, real *, integer *, integer *, integer *, real *, integer *
00292 ), serrgg_(char *, integer *), sdrgvx_(integer *, real *,
00293 integer *, integer *, real *, integer *, real *, real *, real *,
00294 real *, real *, real *, real *, real *, integer *, integer *,
00295 real *, real *, real *, real *, real *, real *, real *, integer *,
00296 integer *, integer *, real *, logical *, integer *);
00297 real result[500];
00298 extern int serrhs_(char *, integer *), xlaenv_(
00299 integer *, integer *);
00300 integer maxtyp;
00301 logical tsterr;
00302 integer ntypes;
00303 logical tstdrv;
00304 extern int sdrvst_(integer *, integer *, integer *,
00305 logical *, integer *, real *, integer *, real *, integer *, real *
00306 , real *, real *, real *, real *, real *, real *, real *, real *,
00307 integer *, real *, real *, real *, real *, integer *, integer *,
00308 integer *, real *, integer *), serrst_(char *, integer *),
00309 sdrvsx_(integer *, integer *, integer *, logical *, integer *,
00310 real *, integer *, integer *, real *, integer *, real *, real *,
00311 real *, real *, real *, real *, real *, real *, real *, integer *,
00312 real *, real *, real *, integer *, integer *, logical *, integer
00313 *), sdrvvx_(integer *, integer *, integer *, logical *, integer *,
00314 real *, integer *, integer *, real *, integer *, real *, real *,
00315 real *, real *, real *, real *, integer *, real *, integer *,
00316 real *, integer *, real *, real *, real *, real *, real *, real *,
00317 real *, real *, real *, real *, integer *, integer *, integer *);
00318
00319
00320 static cilist io___29 = { 0, 6, 0, fmt_9987, 0 };
00321 static cilist io___30 = { 0, 6, 0, fmt_9986, 0 };
00322 static cilist io___31 = { 0, 6, 0, fmt_9985, 0 };
00323 static cilist io___32 = { 0, 6, 0, fmt_9979, 0 };
00324 static cilist io___33 = { 0, 6, 0, fmt_9978, 0 };
00325 static cilist io___34 = { 0, 6, 0, fmt_9977, 0 };
00326 static cilist io___35 = { 0, 6, 0, fmt_9976, 0 };
00327 static cilist io___36 = { 0, 6, 0, fmt_9975, 0 };
00328 static cilist io___37 = { 0, 6, 0, fmt_9964, 0 };
00329 static cilist io___38 = { 0, 6, 0, fmt_9965, 0 };
00330 static cilist io___39 = { 0, 6, 0, fmt_9963, 0 };
00331 static cilist io___40 = { 0, 6, 0, fmt_9962, 0 };
00332 static cilist io___41 = { 0, 6, 0, fmt_9974, 0 };
00333 static cilist io___42 = { 0, 6, 0, fmt_9967, 0 };
00334 static cilist io___43 = { 0, 6, 0, fmt_9971, 0 };
00335 static cilist io___44 = { 0, 6, 0, fmt_9970, 0 };
00336 static cilist io___45 = { 0, 6, 0, fmt_9969, 0 };
00337 static cilist io___46 = { 0, 6, 0, fmt_9968, 0 };
00338 static cilist io___47 = { 0, 5, 0, 0, 0 };
00339 static cilist io___50 = { 0, 6, 0, fmt_9992, 0 };
00340 static cilist io___54 = { 0, 6, 0, fmt_9972, 0 };
00341 static cilist io___55 = { 0, 6, 0, fmt_9984, 0 };
00342 static cilist io___56 = { 0, 5, 0, 0, 0 };
00343 static cilist io___58 = { 0, 6, 0, fmt_9989, 0 };
00344 static cilist io___59 = { 0, 6, 0, fmt_9988, 0 };
00345 static cilist io___60 = { 0, 5, 0, 0, 0 };
00346 static cilist io___64 = { 0, 6, 0, fmt_9989, 0 };
00347 static cilist io___65 = { 0, 6, 0, fmt_9988, 0 };
00348 static cilist io___66 = { 0, 6, 0, fmt_9983, 0 };
00349 static cilist io___67 = { 0, 5, 0, 0, 0 };
00350 static cilist io___69 = { 0, 6, 0, fmt_9989, 0 };
00351 static cilist io___70 = { 0, 6, 0, fmt_9988, 0 };
00352 static cilist io___71 = { 0, 6, 0, fmt_9983, 0 };
00353 static cilist io___72 = { 0, 5, 0, 0, 0 };
00354 static cilist io___74 = { 0, 6, 0, fmt_9989, 0 };
00355 static cilist io___75 = { 0, 6, 0, fmt_9988, 0 };
00356 static cilist io___76 = { 0, 6, 0, fmt_9983, 0 };
00357 static cilist io___77 = { 0, 6, 0, fmt_9983, 0 };
00358 static cilist io___78 = { 0, 5, 0, 0, 0 };
00359 static cilist io___80 = { 0, 5, 0, 0, 0 };
00360 static cilist io___82 = { 0, 6, 0, fmt_9989, 0 };
00361 static cilist io___83 = { 0, 6, 0, fmt_9988, 0 };
00362 static cilist io___84 = { 0, 6, 0, fmt_9983, 0 };
00363 static cilist io___85 = { 0, 5, 0, 0, 0 };
00364 static cilist io___94 = { 0, 6, 0, fmt_9989, 0 };
00365 static cilist io___95 = { 0, 6, 0, fmt_9989, 0 };
00366 static cilist io___96 = { 0, 6, 0, fmt_9989, 0 };
00367 static cilist io___97 = { 0, 6, 0, fmt_9989, 0 };
00368 static cilist io___98 = { 0, 6, 0, fmt_9989, 0 };
00369 static cilist io___99 = { 0, 6, 0, fmt_9989, 0 };
00370 static cilist io___100 = { 0, 6, 0, fmt_9989, 0 };
00371 static cilist io___101 = { 0, 6, 0, fmt_9989, 0 };
00372 static cilist io___102 = { 0, 6, 0, fmt_9983, 0 };
00373 static cilist io___103 = { 0, 6, 0, fmt_9983, 0 };
00374 static cilist io___104 = { 0, 6, 0, fmt_9983, 0 };
00375 static cilist io___105 = { 0, 6, 0, fmt_9983, 0 };
00376 static cilist io___106 = { 0, 6, 0, fmt_9983, 0 };
00377 static cilist io___107 = { 0, 6, 0, fmt_9983, 0 };
00378 static cilist io___108 = { 0, 6, 0, fmt_9983, 0 };
00379 static cilist io___109 = { 0, 6, 0, fmt_9983, 0 };
00380 static cilist io___110 = { 0, 5, 0, 0, 0 };
00381 static cilist io___113 = { 0, 6, 0, fmt_9989, 0 };
00382 static cilist io___114 = { 0, 6, 0, fmt_9989, 0 };
00383 static cilist io___115 = { 0, 6, 0, fmt_9989, 0 };
00384 static cilist io___116 = { 0, 6, 0, fmt_9989, 0 };
00385 static cilist io___117 = { 0, 6, 0, fmt_9989, 0 };
00386 static cilist io___118 = { 0, 6, 0, fmt_9983, 0 };
00387 static cilist io___119 = { 0, 6, 0, fmt_9983, 0 };
00388 static cilist io___120 = { 0, 6, 0, fmt_9983, 0 };
00389 static cilist io___121 = { 0, 6, 0, fmt_9983, 0 };
00390 static cilist io___122 = { 0, 6, 0, fmt_9983, 0 };
00391 static cilist io___123 = { 0, 5, 0, 0, 0 };
00392 static cilist io___125 = { 0, 6, 0, fmt_9989, 0 };
00393 static cilist io___126 = { 0, 6, 0, fmt_9988, 0 };
00394 static cilist io___127 = { 0, 5, 0, 0, 0 };
00395 static cilist io___128 = { 0, 6, 0, fmt_9989, 0 };
00396 static cilist io___129 = { 0, 6, 0, fmt_9988, 0 };
00397 static cilist io___130 = { 0, 6, 0, fmt_9983, 0 };
00398 static cilist io___131 = { 0, 5, 0, 0, 0 };
00399 static cilist io___132 = { 0, 6, 0, fmt_9989, 0 };
00400 static cilist io___133 = { 0, 6, 0, fmt_9988, 0 };
00401 static cilist io___134 = { 0, 6, 0, fmt_9983, 0 };
00402 static cilist io___135 = { 0, 5, 0, 0, 0 };
00403 static cilist io___136 = { 0, 6, 0, fmt_9989, 0 };
00404 static cilist io___137 = { 0, 6, 0, fmt_9988, 0 };
00405 static cilist io___138 = { 0, 6, 0, fmt_9983, 0 };
00406 static cilist io___139 = { 0, 5, 0, 0, 0 };
00407 static cilist io___140 = { 0, 6, 0, fmt_9989, 0 };
00408 static cilist io___141 = { 0, 6, 0, fmt_9988, 0 };
00409 static cilist io___142 = { 0, 6, 0, fmt_9983, 0 };
00410 static cilist io___143 = { 0, 5, 0, 0, 0 };
00411 static cilist io___144 = { 0, 6, 0, fmt_9989, 0 };
00412 static cilist io___145 = { 0, 6, 0, fmt_9988, 0 };
00413 static cilist io___146 = { 0, 6, 0, fmt_9983, 0 };
00414 static cilist io___147 = { 0, 5, 0, 0, 0 };
00415 static cilist io___148 = { 0, 6, 0, fmt_9989, 0 };
00416 static cilist io___149 = { 0, 6, 0, fmt_9983, 0 };
00417 static cilist io___150 = { 0, 5, 0, 0, 0 };
00418 static cilist io___151 = { 0, 6, 0, fmt_9989, 0 };
00419 static cilist io___152 = { 0, 6, 0, fmt_9983, 0 };
00420 static cilist io___153 = { 0, 5, 0, 0, 0 };
00421 static cilist io___154 = { 0, 6, 0, fmt_9989, 0 };
00422 static cilist io___155 = { 0, 6, 0, fmt_9983, 0 };
00423 static cilist io___156 = { 0, 5, 0, 0, 0 };
00424 static cilist io___157 = { 0, 6, 0, fmt_9989, 0 };
00425 static cilist io___158 = { 0, 6, 0, fmt_9983, 0 };
00426 static cilist io___159 = { 0, 5, 0, 0, 0 };
00427 static cilist io___160 = { 0, 6, 0, fmt_9989, 0 };
00428 static cilist io___161 = { 0, 6, 0, fmt_9983, 0 };
00429 static cilist io___162 = { 0, 5, 0, 0, 0 };
00430 static cilist io___164 = { 0, 6, 0, fmt_9989, 0 };
00431 static cilist io___165 = { 0, 6, 0, fmt_9988, 0 };
00432 static cilist io___166 = { 0, 6, 0, fmt_9983, 0 };
00433 static cilist io___167 = { 0, 6, 0, 0, 0 };
00434 static cilist io___169 = { 0, 6, 0, fmt_9981, 0 };
00435 static cilist io___170 = { 0, 6, 0, fmt_9981, 0 };
00436 static cilist io___171 = { 0, 6, 0, fmt_9981, 0 };
00437 static cilist io___172 = { 0, 5, 0, 0, 0 };
00438 static cilist io___173 = { 0, 6, 0, fmt_9982, 0 };
00439 static cilist io___174 = { 0, 5, 0, 0, 0 };
00440 static cilist io___176 = { 0, 5, 0, 0, 0 };
00441 static cilist io___178 = { 0, 5, 0, 0, 0 };
00442 static cilist io___179 = { 0, 5, 0, 0, 0 };
00443 static cilist io___181 = { 0, 5, 0, 0, 0 };
00444 static cilist io___183 = { 0, 6, 0, fmt_9999, 0 };
00445 static cilist io___192 = { 0, 6, 0, fmt_9991, 0 };
00446 static cilist io___193 = { 0, 6, 0, fmt_9990, 0 };
00447 static cilist io___196 = { 0, 6, 0, fmt_9961, 0 };
00448 static cilist io___204 = { 0, 6, 0, fmt_9980, 0 };
00449 static cilist io___205 = { 0, 6, 0, fmt_9997, 0 };
00450 static cilist io___206 = { 0, 6, 0, fmt_9980, 0 };
00451 static cilist io___207 = { 0, 6, 0, fmt_9980, 0 };
00452 static cilist io___208 = { 0, 6, 0, fmt_9997, 0 };
00453 static cilist io___209 = { 0, 6, 0, fmt_9980, 0 };
00454 static cilist io___211 = { 0, 6, 0, fmt_9995, 0 };
00455 static cilist io___212 = { 0, 6, 0, fmt_9980, 0 };
00456 static cilist io___213 = { 0, 6, 0, fmt_9990, 0 };
00457 static cilist io___214 = { 0, 6, 0, fmt_9980, 0 };
00458 static cilist io___215 = { 0, 6, 0, fmt_9973, 0 };
00459 static cilist io___216 = { 0, 6, 0, fmt_9990, 0 };
00460 static cilist io___217 = { 0, 6, 0, fmt_9980, 0 };
00461 static cilist io___218 = { 0, 6, 0, fmt_9973, 0 };
00462 static cilist io___219 = { 0, 6, 0, fmt_9990, 0 };
00463 static cilist io___220 = { 0, 6, 0, fmt_9980, 0 };
00464 static cilist io___221 = { 0, 6, 0, fmt_9973, 0 };
00465 static cilist io___222 = { 0, 6, 0, fmt_9990, 0 };
00466 static cilist io___223 = { 0, 6, 0, fmt_9980, 0 };
00467 static cilist io___224 = { 0, 6, 0, fmt_9973, 0 };
00468 static cilist io___225 = { 0, 6, 0, fmt_9996, 0 };
00469 static cilist io___228 = { 0, 6, 0, fmt_9980, 0 };
00470 static cilist io___229 = { 0, 6, 0, fmt_9980, 0 };
00471 static cilist io___230 = { 0, 6, 0, fmt_9990, 0 };
00472 static cilist io___231 = { 0, 6, 0, fmt_9980, 0 };
00473 static cilist io___232 = { 0, 6, 0, fmt_9973, 0 };
00474 static cilist io___233 = { 0, 6, 0, fmt_9990, 0 };
00475 static cilist io___235 = { 0, 6, 0, fmt_9980, 0 };
00476 static cilist io___236 = { 0, 6, 0, fmt_9973, 0 };
00477 static cilist io___237 = { 0, 6, 0, fmt_9990, 0 };
00478 static cilist io___238 = { 0, 6, 0, fmt_9980, 0 };
00479 static cilist io___239 = { 0, 6, 0, fmt_9973, 0 };
00480 static cilist io___240 = { 0, 6, 0, fmt_9990, 0 };
00481 static cilist io___241 = { 0, 6, 0, fmt_9980, 0 };
00482 static cilist io___242 = { 0, 6, 0, fmt_9973, 0 };
00483 static cilist io___243 = { 0, 6, 0, fmt_9980, 0 };
00484 static cilist io___244 = { 0, 6, 0, fmt_9966, 0 };
00485 static cilist io___245 = { 0, 6, 0, fmt_9980, 0 };
00486 static cilist io___248 = { 0, 6, 0, fmt_9980, 0 };
00487 static cilist io___251 = { 0, 6, 0, fmt_9980, 0 };
00488 static cilist io___252 = { 0, 6, 0, fmt_9980, 0 };
00489 static cilist io___253 = { 0, 6, 0, fmt_9980, 0 };
00490 static cilist io___254 = { 0, 6, 0, 0, 0 };
00491 static cilist io___255 = { 0, 6, 0, 0, 0 };
00492 static cilist io___256 = { 0, 6, 0, fmt_9992, 0 };
00493 static cilist io___257 = { 0, 6, 0, fmt_9994, 0 };
00494 static cilist io___259 = { 0, 6, 0, fmt_9993, 0 };
00495
00496
00497
00498
00499
00500
00501
00502
00503
00504
00505
00506
00507
00508
00509
00510
00511
00512
00513
00514
00515
00516
00517
00518
00519
00520
00521
00522
00523
00524
00525
00526
00527
00528
00529
00530
00531
00532
00533
00534
00535
00536
00537
00538
00539
00540
00541
00542
00543
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 s1 = second_();
01487 fatal = FALSE_;
01488 infoc_1.nunit = 6;
01489
01490
01491
01492 L10:
01493
01494
01495
01496 ci__1.cierr = 0;
01497 ci__1.ciend = 1;
01498 ci__1.ciunit = 5;
01499 ci__1.cifmt = "(A80)";
01500 i__1 = s_rsfe(&ci__1);
01501 if (i__1 != 0) {
01502 goto L380;
01503 }
01504 i__1 = do_fio(&c__1, line, (ftnlen)80);
01505 if (i__1 != 0) {
01506 goto L380;
01507 }
01508 i__1 = e_rsfe();
01509 if (i__1 != 0) {
01510 goto L380;
01511 }
01512 s_copy(path, line, (ftnlen)3, (ftnlen)3);
01513 nep = lsamen_(&c__3, path, "NEP") || lsamen_(&c__3,
01514 path, "SHS");
01515 sep = lsamen_(&c__3, path, "SEP") || lsamen_(&c__3,
01516 path, "SST") || lsamen_(&c__3, path, "SSG");
01517 svd = lsamen_(&c__3, path, "SVD") || lsamen_(&c__3,
01518 path, "SBD");
01519 sev = lsamen_(&c__3, path, "SEV");
01520 ses = lsamen_(&c__3, path, "SES");
01521 svx = lsamen_(&c__3, path, "SVX");
01522 ssx = lsamen_(&c__3, path, "SSX");
01523 sgg = lsamen_(&c__3, path, "SGG");
01524 sgs = lsamen_(&c__3, path, "SGS");
01525 sgx = lsamen_(&c__3, path, "SGX");
01526 sgv = lsamen_(&c__3, path, "SGV");
01527 sxv = lsamen_(&c__3, path, "SXV");
01528 ssb = lsamen_(&c__3, path, "SSB");
01529 sbb = lsamen_(&c__3, path, "SBB");
01530 glm = lsamen_(&c__3, path, "GLM");
01531 gqr = lsamen_(&c__3, path, "GQR") || lsamen_(&c__3,
01532 path, "GRQ");
01533 gsv = lsamen_(&c__3, path, "GSV");
01534 lse = lsamen_(&c__3, path, "LSE");
01535 sbl = lsamen_(&c__3, path, "SBL");
01536 sbk = lsamen_(&c__3, path, "SBK");
01537 sgl = lsamen_(&c__3, path, "SGL");
01538 sgk = lsamen_(&c__3, path, "SGK");
01539
01540
01541
01542 if (s_cmp(path, " ", (ftnlen)3, (ftnlen)3) == 0) {
01543 goto L10;
01544 } else if (nep) {
01545 s_wsfe(&io___29);
01546 e_wsfe();
01547 } else if (sep) {
01548 s_wsfe(&io___30);
01549 e_wsfe();
01550 } else if (svd) {
01551 s_wsfe(&io___31);
01552 e_wsfe();
01553 } else if (sev) {
01554 s_wsfe(&io___32);
01555 e_wsfe();
01556 } else if (ses) {
01557 s_wsfe(&io___33);
01558 e_wsfe();
01559 } else if (svx) {
01560 s_wsfe(&io___34);
01561 e_wsfe();
01562 } else if (ssx) {
01563 s_wsfe(&io___35);
01564 e_wsfe();
01565 } else if (sgg) {
01566 s_wsfe(&io___36);
01567 e_wsfe();
01568 } else if (sgs) {
01569 s_wsfe(&io___37);
01570 e_wsfe();
01571 } else if (sgx) {
01572 s_wsfe(&io___38);
01573 e_wsfe();
01574 } else if (sgv) {
01575 s_wsfe(&io___39);
01576 e_wsfe();
01577 } else if (sxv) {
01578 s_wsfe(&io___40);
01579 e_wsfe();
01580 } else if (ssb) {
01581 s_wsfe(&io___41);
01582 e_wsfe();
01583 } else if (sbb) {
01584 s_wsfe(&io___42);
01585 e_wsfe();
01586 } else if (glm) {
01587 s_wsfe(&io___43);
01588 e_wsfe();
01589 } else if (gqr) {
01590 s_wsfe(&io___44);
01591 e_wsfe();
01592 } else if (gsv) {
01593 s_wsfe(&io___45);
01594 e_wsfe();
01595 } else if (lse) {
01596 s_wsfe(&io___46);
01597 e_wsfe();
01598 } else if (sbl) {
01599
01600
01601
01602 schkbl_(&c__5, &c__6);
01603 goto L10;
01604 } else if (sbk) {
01605
01606
01607
01608 schkbk_(&c__5, &c__6);
01609 goto L10;
01610 } else if (sgl) {
01611
01612
01613
01614 schkgl_(&c__5, &c__6);
01615 goto L10;
01616 } else if (sgk) {
01617
01618
01619
01620 schkgk_(&c__5, &c__6);
01621 goto L10;
01622 } else if (lsamen_(&c__3, path, "SEC")) {
01623
01624
01625
01626 s_rsle(&io___47);
01627 do_lio(&c__4, &c__1, (char *)&thresh, (ftnlen)sizeof(real));
01628 e_rsle();
01629 xlaenv_(&c__1, &c__1);
01630 xlaenv_(&c__12, &c__11);
01631 xlaenv_(&c__13, &c__2);
01632 xlaenv_(&c__14, &c__0);
01633 xlaenv_(&c__15, &c__2);
01634 xlaenv_(&c__16, &c__2);
01635 tsterr = TRUE_;
01636 schkec_(&thresh, &tsterr, &c__5, &c__6);
01637 goto L10;
01638 } else {
01639 s_wsfe(&io___50);
01640 do_fio(&c__1, path, (ftnlen)3);
01641 e_wsfe();
01642 goto L10;
01643 }
01644 ilaver_(&vers_major__, &vers_minor__, &vers_patch__);
01645 s_wsfe(&io___54);
01646 do_fio(&c__1, (char *)&vers_major__, (ftnlen)sizeof(integer));
01647 do_fio(&c__1, (char *)&vers_minor__, (ftnlen)sizeof(integer));
01648 do_fio(&c__1, (char *)&vers_patch__, (ftnlen)sizeof(integer));
01649 e_wsfe();
01650 s_wsfe(&io___55);
01651 e_wsfe();
01652
01653
01654
01655 s_rsle(&io___56);
01656 do_lio(&c__3, &c__1, (char *)&nn, (ftnlen)sizeof(integer));
01657 e_rsle();
01658 if (nn < 0) {
01659 s_wsfe(&io___58);
01660 do_fio(&c__1, " NN ", (ftnlen)6);
01661 do_fio(&c__1, (char *)&nn, (ftnlen)sizeof(integer));
01662 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
01663 e_wsfe();
01664 nn = 0;
01665 fatal = TRUE_;
01666 } else if (nn > 20) {
01667 s_wsfe(&io___59);
01668 do_fio(&c__1, " NN ", (ftnlen)6);
01669 do_fio(&c__1, (char *)&nn, (ftnlen)sizeof(integer));
01670 do_fio(&c__1, (char *)&c__20, (ftnlen)sizeof(integer));
01671 e_wsfe();
01672 nn = 0;
01673 fatal = TRUE_;
01674 }
01675
01676
01677
01678 if (! (sgx || sxv)) {
01679 s_rsle(&io___60);
01680 i__1 = nn;
01681 for (i__ = 1; i__ <= i__1; ++i__) {
01682 do_lio(&c__3, &c__1, (char *)&mval[i__ - 1], (ftnlen)sizeof(
01683 integer));
01684 }
01685 e_rsle();
01686 if (svd) {
01687 s_copy(vname, " M ", (ftnlen)32, (ftnlen)6);
01688 } else {
01689 s_copy(vname, " N ", (ftnlen)32, (ftnlen)6);
01690 }
01691 i__1 = nn;
01692 for (i__ = 1; i__ <= i__1; ++i__) {
01693 if (mval[i__ - 1] < 0) {
01694 s_wsfe(&io___64);
01695 do_fio(&c__1, vname, (ftnlen)32);
01696 do_fio(&c__1, (char *)&mval[i__ - 1], (ftnlen)sizeof(integer))
01697 ;
01698 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
01699 e_wsfe();
01700 fatal = TRUE_;
01701 } else if (mval[i__ - 1] > 132) {
01702 s_wsfe(&io___65);
01703 do_fio(&c__1, vname, (ftnlen)32);
01704 do_fio(&c__1, (char *)&mval[i__ - 1], (ftnlen)sizeof(integer))
01705 ;
01706 do_fio(&c__1, (char *)&c__132, (ftnlen)sizeof(integer));
01707 e_wsfe();
01708 fatal = TRUE_;
01709 }
01710
01711 }
01712 s_wsfe(&io___66);
01713 do_fio(&c__1, "M: ", (ftnlen)6);
01714 i__1 = nn;
01715 for (i__ = 1; i__ <= i__1; ++i__) {
01716 do_fio(&c__1, (char *)&mval[i__ - 1], (ftnlen)sizeof(integer));
01717 }
01718 e_wsfe();
01719 }
01720
01721
01722
01723 if (glm || gqr || gsv || lse) {
01724 s_rsle(&io___67);
01725 i__1 = nn;
01726 for (i__ = 1; i__ <= i__1; ++i__) {
01727 do_lio(&c__3, &c__1, (char *)&pval[i__ - 1], (ftnlen)sizeof(
01728 integer));
01729 }
01730 e_rsle();
01731 i__1 = nn;
01732 for (i__ = 1; i__ <= i__1; ++i__) {
01733 if (pval[i__ - 1] < 0) {
01734 s_wsfe(&io___69);
01735 do_fio(&c__1, " P ", (ftnlen)4);
01736 do_fio(&c__1, (char *)&pval[i__ - 1], (ftnlen)sizeof(integer))
01737 ;
01738 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
01739 e_wsfe();
01740 fatal = TRUE_;
01741 } else if (pval[i__ - 1] > 132) {
01742 s_wsfe(&io___70);
01743 do_fio(&c__1, " P ", (ftnlen)4);
01744 do_fio(&c__1, (char *)&pval[i__ - 1], (ftnlen)sizeof(integer))
01745 ;
01746 do_fio(&c__1, (char *)&c__132, (ftnlen)sizeof(integer));
01747 e_wsfe();
01748 fatal = TRUE_;
01749 }
01750
01751 }
01752 s_wsfe(&io___71);
01753 do_fio(&c__1, "P: ", (ftnlen)6);
01754 i__1 = nn;
01755 for (i__ = 1; i__ <= i__1; ++i__) {
01756 do_fio(&c__1, (char *)&pval[i__ - 1], (ftnlen)sizeof(integer));
01757 }
01758 e_wsfe();
01759 }
01760
01761
01762
01763 if (svd || sbb || glm || gqr || gsv || lse) {
01764 s_rsle(&io___72);
01765 i__1 = nn;
01766 for (i__ = 1; i__ <= i__1; ++i__) {
01767 do_lio(&c__3, &c__1, (char *)&nval[i__ - 1], (ftnlen)sizeof(
01768 integer));
01769 }
01770 e_rsle();
01771 i__1 = nn;
01772 for (i__ = 1; i__ <= i__1; ++i__) {
01773 if (nval[i__ - 1] < 0) {
01774 s_wsfe(&io___74);
01775 do_fio(&c__1, " N ", (ftnlen)6);
01776 do_fio(&c__1, (char *)&nval[i__ - 1], (ftnlen)sizeof(integer))
01777 ;
01778 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
01779 e_wsfe();
01780 fatal = TRUE_;
01781 } else if (nval[i__ - 1] > 132) {
01782 s_wsfe(&io___75);
01783 do_fio(&c__1, " N ", (ftnlen)6);
01784 do_fio(&c__1, (char *)&nval[i__ - 1], (ftnlen)sizeof(integer))
01785 ;
01786 do_fio(&c__1, (char *)&c__132, (ftnlen)sizeof(integer));
01787 e_wsfe();
01788 fatal = TRUE_;
01789 }
01790
01791 }
01792 } else {
01793 i__1 = nn;
01794 for (i__ = 1; i__ <= i__1; ++i__) {
01795 nval[i__ - 1] = mval[i__ - 1];
01796
01797 }
01798 }
01799 if (! (sgx || sxv)) {
01800 s_wsfe(&io___76);
01801 do_fio(&c__1, "N: ", (ftnlen)6);
01802 i__1 = nn;
01803 for (i__ = 1; i__ <= i__1; ++i__) {
01804 do_fio(&c__1, (char *)&nval[i__ - 1], (ftnlen)sizeof(integer));
01805 }
01806 e_wsfe();
01807 } else {
01808 s_wsfe(&io___77);
01809 do_fio(&c__1, "N: ", (ftnlen)6);
01810 do_fio(&c__1, (char *)&nn, (ftnlen)sizeof(integer));
01811 e_wsfe();
01812 }
01813
01814
01815
01816 if (ssb || sbb) {
01817 s_rsle(&io___78);
01818 do_lio(&c__3, &c__1, (char *)&nk, (ftnlen)sizeof(integer));
01819 e_rsle();
01820 s_rsle(&io___80);
01821 i__1 = nk;
01822 for (i__ = 1; i__ <= i__1; ++i__) {
01823 do_lio(&c__3, &c__1, (char *)&kval[i__ - 1], (ftnlen)sizeof(
01824 integer));
01825 }
01826 e_rsle();
01827 i__1 = nk;
01828 for (i__ = 1; i__ <= i__1; ++i__) {
01829 if (kval[i__ - 1] < 0) {
01830 s_wsfe(&io___82);
01831 do_fio(&c__1, " K ", (ftnlen)6);
01832 do_fio(&c__1, (char *)&kval[i__ - 1], (ftnlen)sizeof(integer))
01833 ;
01834 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
01835 e_wsfe();
01836 fatal = TRUE_;
01837 } else if (kval[i__ - 1] > 132) {
01838 s_wsfe(&io___83);
01839 do_fio(&c__1, " K ", (ftnlen)6);
01840 do_fio(&c__1, (char *)&kval[i__ - 1], (ftnlen)sizeof(integer))
01841 ;
01842 do_fio(&c__1, (char *)&c__132, (ftnlen)sizeof(integer));
01843 e_wsfe();
01844 fatal = TRUE_;
01845 }
01846
01847 }
01848 s_wsfe(&io___84);
01849 do_fio(&c__1, "K: ", (ftnlen)6);
01850 i__1 = nk;
01851 for (i__ = 1; i__ <= i__1; ++i__) {
01852 do_fio(&c__1, (char *)&kval[i__ - 1], (ftnlen)sizeof(integer));
01853 }
01854 e_wsfe();
01855 }
01856
01857 if (sev || ses || svx || ssx) {
01858
01859
01860
01861
01862 s_rsle(&io___85);
01863 do_lio(&c__3, &c__1, (char *)&nbval[0], (ftnlen)sizeof(integer));
01864 do_lio(&c__3, &c__1, (char *)&nbmin[0], (ftnlen)sizeof(integer));
01865 do_lio(&c__3, &c__1, (char *)&nxval[0], (ftnlen)sizeof(integer));
01866 do_lio(&c__3, &c__1, (char *)&inmin[0], (ftnlen)sizeof(integer));
01867 do_lio(&c__3, &c__1, (char *)&inwin[0], (ftnlen)sizeof(integer));
01868 do_lio(&c__3, &c__1, (char *)&inibl[0], (ftnlen)sizeof(integer));
01869 do_lio(&c__3, &c__1, (char *)&ishfts[0], (ftnlen)sizeof(integer));
01870 do_lio(&c__3, &c__1, (char *)&iacc22[0], (ftnlen)sizeof(integer));
01871 e_rsle();
01872 if (nbval[0] < 1) {
01873 s_wsfe(&io___94);
01874 do_fio(&c__1, " NB ", (ftnlen)6);
01875 do_fio(&c__1, (char *)&nbval[0], (ftnlen)sizeof(integer));
01876 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
01877 e_wsfe();
01878 fatal = TRUE_;
01879 } else if (nbmin[0] < 1) {
01880 s_wsfe(&io___95);
01881 do_fio(&c__1, "NBMIN ", (ftnlen)6);
01882 do_fio(&c__1, (char *)&nbmin[0], (ftnlen)sizeof(integer));
01883 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
01884 e_wsfe();
01885 fatal = TRUE_;
01886 } else if (nxval[0] < 1) {
01887 s_wsfe(&io___96);
01888 do_fio(&c__1, " NX ", (ftnlen)6);
01889 do_fio(&c__1, (char *)&nxval[0], (ftnlen)sizeof(integer));
01890 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
01891 e_wsfe();
01892 fatal = TRUE_;
01893 } else if (inmin[0] < 1) {
01894 s_wsfe(&io___97);
01895 do_fio(&c__1, " INMIN ", (ftnlen)9);
01896 do_fio(&c__1, (char *)&inmin[0], (ftnlen)sizeof(integer));
01897 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
01898 e_wsfe();
01899 fatal = TRUE_;
01900 } else if (inwin[0] < 1) {
01901 s_wsfe(&io___98);
01902 do_fio(&c__1, " INWIN ", (ftnlen)9);
01903 do_fio(&c__1, (char *)&inwin[0], (ftnlen)sizeof(integer));
01904 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
01905 e_wsfe();
01906 fatal = TRUE_;
01907 } else if (inibl[0] < 1) {
01908 s_wsfe(&io___99);
01909 do_fio(&c__1, " INIBL ", (ftnlen)9);
01910 do_fio(&c__1, (char *)&inibl[0], (ftnlen)sizeof(integer));
01911 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
01912 e_wsfe();
01913 fatal = TRUE_;
01914 } else if (ishfts[0] < 1) {
01915 s_wsfe(&io___100);
01916 do_fio(&c__1, " ISHFTS ", (ftnlen)10);
01917 do_fio(&c__1, (char *)&ishfts[0], (ftnlen)sizeof(integer));
01918 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
01919 e_wsfe();
01920 fatal = TRUE_;
01921 } else if (iacc22[0] < 0) {
01922 s_wsfe(&io___101);
01923 do_fio(&c__1, " IACC22 ", (ftnlen)10);
01924 do_fio(&c__1, (char *)&iacc22[0], (ftnlen)sizeof(integer));
01925 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
01926 e_wsfe();
01927 fatal = TRUE_;
01928 }
01929 xlaenv_(&c__1, nbval);
01930 xlaenv_(&c__2, nbmin);
01931 xlaenv_(&c__3, nxval);
01932 i__1 = max(11,inmin[0]);
01933 xlaenv_(&c__12, &i__1);
01934 xlaenv_(&c__13, inwin);
01935 xlaenv_(&c__14, inibl);
01936 xlaenv_(&c__15, ishfts);
01937 xlaenv_(&c__16, iacc22);
01938 s_wsfe(&io___102);
01939 do_fio(&c__1, "NB: ", (ftnlen)6);
01940 do_fio(&c__1, (char *)&nbval[0], (ftnlen)sizeof(integer));
01941 e_wsfe();
01942 s_wsfe(&io___103);
01943 do_fio(&c__1, "NBMIN:", (ftnlen)6);
01944 do_fio(&c__1, (char *)&nbmin[0], (ftnlen)sizeof(integer));
01945 e_wsfe();
01946 s_wsfe(&io___104);
01947 do_fio(&c__1, "NX: ", (ftnlen)6);
01948 do_fio(&c__1, (char *)&nxval[0], (ftnlen)sizeof(integer));
01949 e_wsfe();
01950 s_wsfe(&io___105);
01951 do_fio(&c__1, "INMIN: ", (ftnlen)9);
01952 do_fio(&c__1, (char *)&inmin[0], (ftnlen)sizeof(integer));
01953 e_wsfe();
01954 s_wsfe(&io___106);
01955 do_fio(&c__1, "INWIN: ", (ftnlen)7);
01956 do_fio(&c__1, (char *)&inwin[0], (ftnlen)sizeof(integer));
01957 e_wsfe();
01958 s_wsfe(&io___107);
01959 do_fio(&c__1, "INIBL: ", (ftnlen)7);
01960 do_fio(&c__1, (char *)&inibl[0], (ftnlen)sizeof(integer));
01961 e_wsfe();
01962 s_wsfe(&io___108);
01963 do_fio(&c__1, "ISHFTS: ", (ftnlen)8);
01964 do_fio(&c__1, (char *)&ishfts[0], (ftnlen)sizeof(integer));
01965 e_wsfe();
01966 s_wsfe(&io___109);
01967 do_fio(&c__1, "IACC22: ", (ftnlen)8);
01968 do_fio(&c__1, (char *)&iacc22[0], (ftnlen)sizeof(integer));
01969 e_wsfe();
01970
01971 } else if (sgs || sgx || sgv || sxv) {
01972
01973
01974
01975
01976 s_rsle(&io___110);
01977 do_lio(&c__3, &c__1, (char *)&nbval[0], (ftnlen)sizeof(integer));
01978 do_lio(&c__3, &c__1, (char *)&nbmin[0], (ftnlen)sizeof(integer));
01979 do_lio(&c__3, &c__1, (char *)&nxval[0], (ftnlen)sizeof(integer));
01980 do_lio(&c__3, &c__1, (char *)&nsval[0], (ftnlen)sizeof(integer));
01981 do_lio(&c__3, &c__1, (char *)&mxbval[0], (ftnlen)sizeof(integer));
01982 e_rsle();
01983 if (nbval[0] < 1) {
01984 s_wsfe(&io___113);
01985 do_fio(&c__1, " NB ", (ftnlen)6);
01986 do_fio(&c__1, (char *)&nbval[0], (ftnlen)sizeof(integer));
01987 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
01988 e_wsfe();
01989 fatal = TRUE_;
01990 } else if (nbmin[0] < 1) {
01991 s_wsfe(&io___114);
01992 do_fio(&c__1, "NBMIN ", (ftnlen)6);
01993 do_fio(&c__1, (char *)&nbmin[0], (ftnlen)sizeof(integer));
01994 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
01995 e_wsfe();
01996 fatal = TRUE_;
01997 } else if (nxval[0] < 1) {
01998 s_wsfe(&io___115);
01999 do_fio(&c__1, " NX ", (ftnlen)6);
02000 do_fio(&c__1, (char *)&nxval[0], (ftnlen)sizeof(integer));
02001 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
02002 e_wsfe();
02003 fatal = TRUE_;
02004 } else if (nsval[0] < 2) {
02005 s_wsfe(&io___116);
02006 do_fio(&c__1, " NS ", (ftnlen)6);
02007 do_fio(&c__1, (char *)&nsval[0], (ftnlen)sizeof(integer));
02008 do_fio(&c__1, (char *)&c__2, (ftnlen)sizeof(integer));
02009 e_wsfe();
02010 fatal = TRUE_;
02011 } else if (mxbval[0] < 1) {
02012 s_wsfe(&io___117);
02013 do_fio(&c__1, " MAXB ", (ftnlen)6);
02014 do_fio(&c__1, (char *)&mxbval[0], (ftnlen)sizeof(integer));
02015 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
02016 e_wsfe();
02017 fatal = TRUE_;
02018 }
02019 xlaenv_(&c__1, nbval);
02020 xlaenv_(&c__2, nbmin);
02021 xlaenv_(&c__3, nxval);
02022 xlaenv_(&c__4, nsval);
02023 xlaenv_(&c__8, mxbval);
02024 s_wsfe(&io___118);
02025 do_fio(&c__1, "NB: ", (ftnlen)6);
02026 do_fio(&c__1, (char *)&nbval[0], (ftnlen)sizeof(integer));
02027 e_wsfe();
02028 s_wsfe(&io___119);
02029 do_fio(&c__1, "NBMIN:", (ftnlen)6);
02030 do_fio(&c__1, (char *)&nbmin[0], (ftnlen)sizeof(integer));
02031 e_wsfe();
02032 s_wsfe(&io___120);
02033 do_fio(&c__1, "NX: ", (ftnlen)6);
02034 do_fio(&c__1, (char *)&nxval[0], (ftnlen)sizeof(integer));
02035 e_wsfe();
02036 s_wsfe(&io___121);
02037 do_fio(&c__1, "NS: ", (ftnlen)6);
02038 do_fio(&c__1, (char *)&nsval[0], (ftnlen)sizeof(integer));
02039 e_wsfe();
02040 s_wsfe(&io___122);
02041 do_fio(&c__1, "MAXB: ", (ftnlen)6);
02042 do_fio(&c__1, (char *)&mxbval[0], (ftnlen)sizeof(integer));
02043 e_wsfe();
02044
02045 } else if (! ssb && ! glm && ! gqr && ! gsv && ! lse) {
02046
02047
02048
02049
02050 s_rsle(&io___123);
02051 do_lio(&c__3, &c__1, (char *)&nparms, (ftnlen)sizeof(integer));
02052 e_rsle();
02053 if (nparms < 1) {
02054 s_wsfe(&io___125);
02055 do_fio(&c__1, "NPARMS", (ftnlen)6);
02056 do_fio(&c__1, (char *)&nparms, (ftnlen)sizeof(integer));
02057 do_fio(&c__1, (char *)&c__1, (ftnlen)sizeof(integer));
02058 e_wsfe();
02059 nparms = 0;
02060 fatal = TRUE_;
02061 } else if (nparms > 20) {
02062 s_wsfe(&io___126);
02063 do_fio(&c__1, "NPARMS", (ftnlen)6);
02064 do_fio(&c__1, (char *)&nparms, (ftnlen)sizeof(integer));
02065 do_fio(&c__1, (char *)&c__20, (ftnlen)sizeof(integer));
02066 e_wsfe();
02067 nparms = 0;
02068 fatal = TRUE_;
02069 }
02070
02071
02072
02073 if (! sbb) {
02074 s_rsle(&io___127);
02075 i__1 = nparms;
02076 for (i__ = 1; i__ <= i__1; ++i__) {
02077 do_lio(&c__3, &c__1, (char *)&nbval[i__ - 1], (ftnlen)sizeof(
02078 integer));
02079 }
02080 e_rsle();
02081 i__1 = nparms;
02082 for (i__ = 1; i__ <= i__1; ++i__) {
02083 if (nbval[i__ - 1] < 0) {
02084 s_wsfe(&io___128);
02085 do_fio(&c__1, " NB ", (ftnlen)6);
02086 do_fio(&c__1, (char *)&nbval[i__ - 1], (ftnlen)sizeof(
02087 integer));
02088 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
02089 e_wsfe();
02090 fatal = TRUE_;
02091 } else if (nbval[i__ - 1] > 132) {
02092 s_wsfe(&io___129);
02093 do_fio(&c__1, " NB ", (ftnlen)6);
02094 do_fio(&c__1, (char *)&nbval[i__ - 1], (ftnlen)sizeof(
02095 integer));
02096 do_fio(&c__1, (char *)&c__132, (ftnlen)sizeof(integer));
02097 e_wsfe();
02098 fatal = TRUE_;
02099 }
02100
02101 }
02102 s_wsfe(&io___130);
02103 do_fio(&c__1, "NB: ", (ftnlen)6);
02104 i__1 = nparms;
02105 for (i__ = 1; i__ <= i__1; ++i__) {
02106 do_fio(&c__1, (char *)&nbval[i__ - 1], (ftnlen)sizeof(integer)
02107 );
02108 }
02109 e_wsfe();
02110 }
02111
02112
02113
02114 if (nep || sep || svd || sgg) {
02115 s_rsle(&io___131);
02116 i__1 = nparms;
02117 for (i__ = 1; i__ <= i__1; ++i__) {
02118 do_lio(&c__3, &c__1, (char *)&nbmin[i__ - 1], (ftnlen)sizeof(
02119 integer));
02120 }
02121 e_rsle();
02122 i__1 = nparms;
02123 for (i__ = 1; i__ <= i__1; ++i__) {
02124 if (nbmin[i__ - 1] < 0) {
02125 s_wsfe(&io___132);
02126 do_fio(&c__1, "NBMIN ", (ftnlen)6);
02127 do_fio(&c__1, (char *)&nbmin[i__ - 1], (ftnlen)sizeof(
02128 integer));
02129 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
02130 e_wsfe();
02131 fatal = TRUE_;
02132 } else if (nbmin[i__ - 1] > 132) {
02133 s_wsfe(&io___133);
02134 do_fio(&c__1, "NBMIN ", (ftnlen)6);
02135 do_fio(&c__1, (char *)&nbmin[i__ - 1], (ftnlen)sizeof(
02136 integer));
02137 do_fio(&c__1, (char *)&c__132, (ftnlen)sizeof(integer));
02138 e_wsfe();
02139 fatal = TRUE_;
02140 }
02141
02142 }
02143 s_wsfe(&io___134);
02144 do_fio(&c__1, "NBMIN:", (ftnlen)6);
02145 i__1 = nparms;
02146 for (i__ = 1; i__ <= i__1; ++i__) {
02147 do_fio(&c__1, (char *)&nbmin[i__ - 1], (ftnlen)sizeof(integer)
02148 );
02149 }
02150 e_wsfe();
02151 } else {
02152 i__1 = nparms;
02153 for (i__ = 1; i__ <= i__1; ++i__) {
02154 nbmin[i__ - 1] = 1;
02155
02156 }
02157 }
02158
02159
02160
02161 if (nep || sep || svd) {
02162 s_rsle(&io___135);
02163 i__1 = nparms;
02164 for (i__ = 1; i__ <= i__1; ++i__) {
02165 do_lio(&c__3, &c__1, (char *)&nxval[i__ - 1], (ftnlen)sizeof(
02166 integer));
02167 }
02168 e_rsle();
02169 i__1 = nparms;
02170 for (i__ = 1; i__ <= i__1; ++i__) {
02171 if (nxval[i__ - 1] < 0) {
02172 s_wsfe(&io___136);
02173 do_fio(&c__1, " NX ", (ftnlen)6);
02174 do_fio(&c__1, (char *)&nxval[i__ - 1], (ftnlen)sizeof(
02175 integer));
02176 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
02177 e_wsfe();
02178 fatal = TRUE_;
02179 } else if (nxval[i__ - 1] > 132) {
02180 s_wsfe(&io___137);
02181 do_fio(&c__1, " NX ", (ftnlen)6);
02182 do_fio(&c__1, (char *)&nxval[i__ - 1], (ftnlen)sizeof(
02183 integer));
02184 do_fio(&c__1, (char *)&c__132, (ftnlen)sizeof(integer));
02185 e_wsfe();
02186 fatal = TRUE_;
02187 }
02188
02189 }
02190 s_wsfe(&io___138);
02191 do_fio(&c__1, "NX: ", (ftnlen)6);
02192 i__1 = nparms;
02193 for (i__ = 1; i__ <= i__1; ++i__) {
02194 do_fio(&c__1, (char *)&nxval[i__ - 1], (ftnlen)sizeof(integer)
02195 );
02196 }
02197 e_wsfe();
02198 } else {
02199 i__1 = nparms;
02200 for (i__ = 1; i__ <= i__1; ++i__) {
02201 nxval[i__ - 1] = 1;
02202
02203 }
02204 }
02205
02206
02207
02208
02209 if (svd || sbb || sgg) {
02210 s_rsle(&io___139);
02211 i__1 = nparms;
02212 for (i__ = 1; i__ <= i__1; ++i__) {
02213 do_lio(&c__3, &c__1, (char *)&nsval[i__ - 1], (ftnlen)sizeof(
02214 integer));
02215 }
02216 e_rsle();
02217 i__1 = nparms;
02218 for (i__ = 1; i__ <= i__1; ++i__) {
02219 if (nsval[i__ - 1] < 0) {
02220 s_wsfe(&io___140);
02221 do_fio(&c__1, " NS ", (ftnlen)6);
02222 do_fio(&c__1, (char *)&nsval[i__ - 1], (ftnlen)sizeof(
02223 integer));
02224 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
02225 e_wsfe();
02226 fatal = TRUE_;
02227 } else if (nsval[i__ - 1] > 132) {
02228 s_wsfe(&io___141);
02229 do_fio(&c__1, " NS ", (ftnlen)6);
02230 do_fio(&c__1, (char *)&nsval[i__ - 1], (ftnlen)sizeof(
02231 integer));
02232 do_fio(&c__1, (char *)&c__132, (ftnlen)sizeof(integer));
02233 e_wsfe();
02234 fatal = TRUE_;
02235 }
02236
02237 }
02238 s_wsfe(&io___142);
02239 do_fio(&c__1, "NS: ", (ftnlen)6);
02240 i__1 = nparms;
02241 for (i__ = 1; i__ <= i__1; ++i__) {
02242 do_fio(&c__1, (char *)&nsval[i__ - 1], (ftnlen)sizeof(integer)
02243 );
02244 }
02245 e_wsfe();
02246 } else {
02247 i__1 = nparms;
02248 for (i__ = 1; i__ <= i__1; ++i__) {
02249 nsval[i__ - 1] = 1;
02250
02251 }
02252 }
02253
02254
02255
02256 if (sgg) {
02257 s_rsle(&io___143);
02258 i__1 = nparms;
02259 for (i__ = 1; i__ <= i__1; ++i__) {
02260 do_lio(&c__3, &c__1, (char *)&mxbval[i__ - 1], (ftnlen)sizeof(
02261 integer));
02262 }
02263 e_rsle();
02264 i__1 = nparms;
02265 for (i__ = 1; i__ <= i__1; ++i__) {
02266 if (mxbval[i__ - 1] < 0) {
02267 s_wsfe(&io___144);
02268 do_fio(&c__1, " MAXB ", (ftnlen)6);
02269 do_fio(&c__1, (char *)&mxbval[i__ - 1], (ftnlen)sizeof(
02270 integer));
02271 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
02272 e_wsfe();
02273 fatal = TRUE_;
02274 } else if (mxbval[i__ - 1] > 132) {
02275 s_wsfe(&io___145);
02276 do_fio(&c__1, " MAXB ", (ftnlen)6);
02277 do_fio(&c__1, (char *)&mxbval[i__ - 1], (ftnlen)sizeof(
02278 integer));
02279 do_fio(&c__1, (char *)&c__132, (ftnlen)sizeof(integer));
02280 e_wsfe();
02281 fatal = TRUE_;
02282 }
02283
02284 }
02285 s_wsfe(&io___146);
02286 do_fio(&c__1, "MAXB: ", (ftnlen)6);
02287 i__1 = nparms;
02288 for (i__ = 1; i__ <= i__1; ++i__) {
02289 do_fio(&c__1, (char *)&mxbval[i__ - 1], (ftnlen)sizeof(
02290 integer));
02291 }
02292 e_wsfe();
02293 } else {
02294 i__1 = nparms;
02295 for (i__ = 1; i__ <= i__1; ++i__) {
02296 mxbval[i__ - 1] = 1;
02297
02298 }
02299 }
02300
02301
02302
02303 if (nep) {
02304 s_rsle(&io___147);
02305 i__1 = nparms;
02306 for (i__ = 1; i__ <= i__1; ++i__) {
02307 do_lio(&c__3, &c__1, (char *)&inmin[i__ - 1], (ftnlen)sizeof(
02308 integer));
02309 }
02310 e_rsle();
02311 i__1 = nparms;
02312 for (i__ = 1; i__ <= i__1; ++i__) {
02313 if (inmin[i__ - 1] < 0) {
02314 s_wsfe(&io___148);
02315 do_fio(&c__1, " INMIN ", (ftnlen)7);
02316 do_fio(&c__1, (char *)&inmin[i__ - 1], (ftnlen)sizeof(
02317 integer));
02318 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
02319 e_wsfe();
02320 fatal = TRUE_;
02321 }
02322
02323 }
02324 s_wsfe(&io___149);
02325 do_fio(&c__1, "INMIN: ", (ftnlen)7);
02326 i__1 = nparms;
02327 for (i__ = 1; i__ <= i__1; ++i__) {
02328 do_fio(&c__1, (char *)&inmin[i__ - 1], (ftnlen)sizeof(integer)
02329 );
02330 }
02331 e_wsfe();
02332 } else {
02333 i__1 = nparms;
02334 for (i__ = 1; i__ <= i__1; ++i__) {
02335 inmin[i__ - 1] = 1;
02336
02337 }
02338 }
02339
02340
02341
02342 if (nep) {
02343 s_rsle(&io___150);
02344 i__1 = nparms;
02345 for (i__ = 1; i__ <= i__1; ++i__) {
02346 do_lio(&c__3, &c__1, (char *)&inwin[i__ - 1], (ftnlen)sizeof(
02347 integer));
02348 }
02349 e_rsle();
02350 i__1 = nparms;
02351 for (i__ = 1; i__ <= i__1; ++i__) {
02352 if (inwin[i__ - 1] < 0) {
02353 s_wsfe(&io___151);
02354 do_fio(&c__1, " INWIN ", (ftnlen)7);
02355 do_fio(&c__1, (char *)&inwin[i__ - 1], (ftnlen)sizeof(
02356 integer));
02357 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
02358 e_wsfe();
02359 fatal = TRUE_;
02360 }
02361
02362 }
02363 s_wsfe(&io___152);
02364 do_fio(&c__1, "INWIN: ", (ftnlen)7);
02365 i__1 = nparms;
02366 for (i__ = 1; i__ <= i__1; ++i__) {
02367 do_fio(&c__1, (char *)&inwin[i__ - 1], (ftnlen)sizeof(integer)
02368 );
02369 }
02370 e_wsfe();
02371 } else {
02372 i__1 = nparms;
02373 for (i__ = 1; i__ <= i__1; ++i__) {
02374 inwin[i__ - 1] = 1;
02375
02376 }
02377 }
02378
02379
02380
02381 if (nep) {
02382 s_rsle(&io___153);
02383 i__1 = nparms;
02384 for (i__ = 1; i__ <= i__1; ++i__) {
02385 do_lio(&c__3, &c__1, (char *)&inibl[i__ - 1], (ftnlen)sizeof(
02386 integer));
02387 }
02388 e_rsle();
02389 i__1 = nparms;
02390 for (i__ = 1; i__ <= i__1; ++i__) {
02391 if (inibl[i__ - 1] < 0) {
02392 s_wsfe(&io___154);
02393 do_fio(&c__1, " INIBL ", (ftnlen)7);
02394 do_fio(&c__1, (char *)&inibl[i__ - 1], (ftnlen)sizeof(
02395 integer));
02396 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
02397 e_wsfe();
02398 fatal = TRUE_;
02399 }
02400
02401 }
02402 s_wsfe(&io___155);
02403 do_fio(&c__1, "INIBL: ", (ftnlen)7);
02404 i__1 = nparms;
02405 for (i__ = 1; i__ <= i__1; ++i__) {
02406 do_fio(&c__1, (char *)&inibl[i__ - 1], (ftnlen)sizeof(integer)
02407 );
02408 }
02409 e_wsfe();
02410 } else {
02411 i__1 = nparms;
02412 for (i__ = 1; i__ <= i__1; ++i__) {
02413 inibl[i__ - 1] = 1;
02414
02415 }
02416 }
02417
02418
02419
02420 if (nep) {
02421 s_rsle(&io___156);
02422 i__1 = nparms;
02423 for (i__ = 1; i__ <= i__1; ++i__) {
02424 do_lio(&c__3, &c__1, (char *)&ishfts[i__ - 1], (ftnlen)sizeof(
02425 integer));
02426 }
02427 e_rsle();
02428 i__1 = nparms;
02429 for (i__ = 1; i__ <= i__1; ++i__) {
02430 if (ishfts[i__ - 1] < 0) {
02431 s_wsfe(&io___157);
02432 do_fio(&c__1, " ISHFTS ", (ftnlen)8);
02433 do_fio(&c__1, (char *)&ishfts[i__ - 1], (ftnlen)sizeof(
02434 integer));
02435 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
02436 e_wsfe();
02437 fatal = TRUE_;
02438 }
02439
02440 }
02441 s_wsfe(&io___158);
02442 do_fio(&c__1, "ISHFTS: ", (ftnlen)8);
02443 i__1 = nparms;
02444 for (i__ = 1; i__ <= i__1; ++i__) {
02445 do_fio(&c__1, (char *)&ishfts[i__ - 1], (ftnlen)sizeof(
02446 integer));
02447 }
02448 e_wsfe();
02449 } else {
02450 i__1 = nparms;
02451 for (i__ = 1; i__ <= i__1; ++i__) {
02452 ishfts[i__ - 1] = 1;
02453
02454 }
02455 }
02456
02457
02458
02459 if (nep) {
02460 s_rsle(&io___159);
02461 i__1 = nparms;
02462 for (i__ = 1; i__ <= i__1; ++i__) {
02463 do_lio(&c__3, &c__1, (char *)&iacc22[i__ - 1], (ftnlen)sizeof(
02464 integer));
02465 }
02466 e_rsle();
02467 i__1 = nparms;
02468 for (i__ = 1; i__ <= i__1; ++i__) {
02469 if (iacc22[i__ - 1] < 0) {
02470 s_wsfe(&io___160);
02471 do_fio(&c__1, " IACC22 ", (ftnlen)8);
02472 do_fio(&c__1, (char *)&iacc22[i__ - 1], (ftnlen)sizeof(
02473 integer));
02474 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
02475 e_wsfe();
02476 fatal = TRUE_;
02477 }
02478
02479 }
02480 s_wsfe(&io___161);
02481 do_fio(&c__1, "IACC22: ", (ftnlen)8);
02482 i__1 = nparms;
02483 for (i__ = 1; i__ <= i__1; ++i__) {
02484 do_fio(&c__1, (char *)&iacc22[i__ - 1], (ftnlen)sizeof(
02485 integer));
02486 }
02487 e_wsfe();
02488 } else {
02489 i__1 = nparms;
02490 for (i__ = 1; i__ <= i__1; ++i__) {
02491 iacc22[i__ - 1] = 1;
02492
02493 }
02494 }
02495
02496
02497
02498 if (sgg) {
02499 s_rsle(&io___162);
02500 i__1 = nparms;
02501 for (i__ = 1; i__ <= i__1; ++i__) {
02502 do_lio(&c__3, &c__1, (char *)&nbcol[i__ - 1], (ftnlen)sizeof(
02503 integer));
02504 }
02505 e_rsle();
02506 i__1 = nparms;
02507 for (i__ = 1; i__ <= i__1; ++i__) {
02508 if (nbcol[i__ - 1] < 0) {
02509 s_wsfe(&io___164);
02510 do_fio(&c__1, "NBCOL ", (ftnlen)6);
02511 do_fio(&c__1, (char *)&nbcol[i__ - 1], (ftnlen)sizeof(
02512 integer));
02513 do_fio(&c__1, (char *)&c__0, (ftnlen)sizeof(integer));
02514 e_wsfe();
02515 fatal = TRUE_;
02516 } else if (nbcol[i__ - 1] > 132) {
02517 s_wsfe(&io___165);
02518 do_fio(&c__1, "NBCOL ", (ftnlen)6);
02519 do_fio(&c__1, (char *)&nbcol[i__ - 1], (ftnlen)sizeof(
02520 integer));
02521 do_fio(&c__1, (char *)&c__132, (ftnlen)sizeof(integer));
02522 e_wsfe();
02523 fatal = TRUE_;
02524 }
02525
02526 }
02527 s_wsfe(&io___166);
02528 do_fio(&c__1, "NBCOL:", (ftnlen)6);
02529 i__1 = nparms;
02530 for (i__ = 1; i__ <= i__1; ++i__) {
02531 do_fio(&c__1, (char *)&nbcol[i__ - 1], (ftnlen)sizeof(integer)
02532 );
02533 }
02534 e_wsfe();
02535 } else {
02536 i__1 = nparms;
02537 for (i__ = 1; i__ <= i__1; ++i__) {
02538 nbcol[i__ - 1] = 1;
02539
02540 }
02541 }
02542 }
02543
02544
02545
02546 s_wsle(&io___167);
02547 e_wsle();
02548 eps = slamch_("Underflow threshold");
02549 s_wsfe(&io___169);
02550 do_fio(&c__1, "underflow", (ftnlen)9);
02551 do_fio(&c__1, (char *)&eps, (ftnlen)sizeof(real));
02552 e_wsfe();
02553 eps = slamch_("Overflow threshold");
02554 s_wsfe(&io___170);
02555 do_fio(&c__1, "overflow ", (ftnlen)9);
02556 do_fio(&c__1, (char *)&eps, (ftnlen)sizeof(real));
02557 e_wsfe();
02558 eps = slamch_("Epsilon");
02559 s_wsfe(&io___171);
02560 do_fio(&c__1, "precision", (ftnlen)9);
02561 do_fio(&c__1, (char *)&eps, (ftnlen)sizeof(real));
02562 e_wsfe();
02563
02564
02565
02566 s_rsle(&io___172);
02567 do_lio(&c__4, &c__1, (char *)&thresh, (ftnlen)sizeof(real));
02568 e_rsle();
02569 s_wsfe(&io___173);
02570 do_fio(&c__1, (char *)&thresh, (ftnlen)sizeof(real));
02571 e_wsfe();
02572 if (sep || svd || sgg) {
02573
02574
02575
02576 s_rsle(&io___174);
02577 do_lio(&c__8, &c__1, (char *)&tstchk, (ftnlen)sizeof(logical));
02578 e_rsle();
02579
02580
02581
02582 s_rsle(&io___176);
02583 do_lio(&c__8, &c__1, (char *)&tstdrv, (ftnlen)sizeof(logical));
02584 e_rsle();
02585 }
02586
02587
02588
02589 s_rsle(&io___178);
02590 do_lio(&c__8, &c__1, (char *)&tsterr, (ftnlen)sizeof(logical));
02591 e_rsle();
02592
02593
02594
02595 s_rsle(&io___179);
02596 do_lio(&c__3, &c__1, (char *)&newsd, (ftnlen)sizeof(integer));
02597 e_rsle();
02598
02599
02600
02601 if (newsd == 2) {
02602 s_rsle(&io___181);
02603 for (i__ = 1; i__ <= 4; ++i__) {
02604 do_lio(&c__3, &c__1, (char *)&ioldsd[i__ - 1], (ftnlen)sizeof(
02605 integer));
02606 }
02607 e_rsle();
02608 }
02609
02610 for (i__ = 1; i__ <= 4; ++i__) {
02611 iseed[i__ - 1] = ioldsd[i__ - 1];
02612
02613 }
02614
02615 if (fatal) {
02616 s_wsfe(&io___183);
02617 e_wsfe();
02618 s_stop("", (ftnlen)0);
02619 }
02620
02621
02622
02623
02624
02625
02626 L190:
02627
02628 if (! (sgx || sxv)) {
02629
02630 L200:
02631 ci__1.cierr = 0;
02632 ci__1.ciend = 1;
02633 ci__1.ciunit = 5;
02634 ci__1.cifmt = "(A80)";
02635 i__1 = s_rsfe(&ci__1);
02636 if (i__1 != 0) {
02637 goto L380;
02638 }
02639 i__1 = do_fio(&c__1, line, (ftnlen)80);
02640 if (i__1 != 0) {
02641 goto L380;
02642 }
02643 i__1 = e_rsfe();
02644 if (i__1 != 0) {
02645 goto L380;
02646 }
02647 s_copy(c3, line, (ftnlen)3, (ftnlen)3);
02648 lenp = i_len(line, (ftnlen)80);
02649 i__ = 3;
02650 itmp = 0;
02651 i1 = 0;
02652 L210:
02653 ++i__;
02654 if (i__ > lenp) {
02655 if (i1 > 0) {
02656 goto L240;
02657 } else {
02658 ntypes = 30;
02659 goto L240;
02660 }
02661 }
02662 if (*(unsigned char *)&line[i__ - 1] != ' ' && *(unsigned char *)&
02663 line[i__ - 1] != ',') {
02664 i1 = i__;
02665 *(unsigned char *)c1 = *(unsigned char *)&line[i1 - 1];
02666
02667
02668
02669 for (k = 1; k <= 10; ++k) {
02670 if (*(unsigned char *)c1 == *(unsigned char *)&intstr[k - 1])
02671 {
02672 ic = k - 1;
02673 goto L230;
02674 }
02675
02676 }
02677 s_wsfe(&io___192);
02678 do_fio(&c__1, (char *)&i__, (ftnlen)sizeof(integer));
02679 do_fio(&c__1, line, (ftnlen)80);
02680 e_wsfe();
02681 goto L200;
02682 L230:
02683 itmp = itmp * 10 + ic;
02684 goto L210;
02685 } else if (i1 > 0) {
02686 goto L240;
02687 } else {
02688 goto L210;
02689 }
02690 L240:
02691 ntypes = itmp;
02692
02693
02694
02695 if (! (sev || ses || svx || ssx || sgv || sgs) && ntypes <= 0) {
02696 s_wsfe(&io___193);
02697 do_fio(&c__1, c3, (ftnlen)3);
02698 e_wsfe();
02699 goto L200;
02700 }
02701
02702 } else {
02703 if (sxv) {
02704 s_copy(c3, "SXV", (ftnlen)3, (ftnlen)3);
02705 }
02706 if (sgx) {
02707 s_copy(c3, "SGX", (ftnlen)3, (ftnlen)3);
02708 }
02709 }
02710
02711
02712
02713 if (newsd == 0) {
02714 for (k = 1; k <= 4; ++k) {
02715 iseed[k - 1] = ioldsd[k - 1];
02716
02717 }
02718 }
02719
02720 if (lsamen_(&c__3, c3, "SHS") || lsamen_(&c__3, c3,
02721 "NEP")) {
02722
02723
02724
02725
02726
02727
02728
02729
02730
02731
02732
02733 maxtyp = 21;
02734 ntypes = min(maxtyp,ntypes);
02735 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
02736 xlaenv_(&c__1, &c__1);
02737 if (tsterr) {
02738 serrhs_("SHSEQR", &c__6);
02739 }
02740 i__1 = nparms;
02741 for (i__ = 1; i__ <= i__1; ++i__) {
02742 xlaenv_(&c__1, &nbval[i__ - 1]);
02743 xlaenv_(&c__2, &nbmin[i__ - 1]);
02744 xlaenv_(&c__3, &nxval[i__ - 1]);
02745
02746 i__3 = 11, i__4 = inmin[i__ - 1];
02747 i__2 = max(i__3,i__4);
02748 xlaenv_(&c__12, &i__2);
02749 xlaenv_(&c__13, &inwin[i__ - 1]);
02750 xlaenv_(&c__14, &inibl[i__ - 1]);
02751 xlaenv_(&c__15, &ishfts[i__ - 1]);
02752 xlaenv_(&c__16, &iacc22[i__ - 1]);
02753
02754 if (newsd == 0) {
02755 for (k = 1; k <= 4; ++k) {
02756 iseed[k - 1] = ioldsd[k - 1];
02757
02758 }
02759 }
02760 s_wsfe(&io___196);
02761 do_fio(&c__1, c3, (ftnlen)3);
02762 do_fio(&c__1, (char *)&nbval[i__ - 1], (ftnlen)sizeof(integer));
02763 do_fio(&c__1, (char *)&nbmin[i__ - 1], (ftnlen)sizeof(integer));
02764 do_fio(&c__1, (char *)&nxval[i__ - 1], (ftnlen)sizeof(integer));
02765
02766 i__3 = 11, i__4 = inmin[i__ - 1];
02767 i__2 = max(i__3,i__4);
02768 do_fio(&c__1, (char *)&i__2, (ftnlen)sizeof(integer));
02769 do_fio(&c__1, (char *)&inwin[i__ - 1], (ftnlen)sizeof(integer));
02770 do_fio(&c__1, (char *)&inibl[i__ - 1], (ftnlen)sizeof(integer));
02771 do_fio(&c__1, (char *)&ishfts[i__ - 1], (ftnlen)sizeof(integer));
02772 do_fio(&c__1, (char *)&iacc22[i__ - 1], (ftnlen)sizeof(integer));
02773 e_wsfe();
02774 schkhs_(&nn, nval, &maxtyp, dotype, iseed, &thresh, &c__6, a, &
02775 c__132, &a[17424], &a[34848], &a[52272], &a[69696], &
02776 c__132, &a[87120], &a[104544], d__, &d__[132], &d__[264],
02777 &d__[396], &a[121968], &a[139392], &a[156816], &a[174240],
02778 &a[191664], &d__[528], work, &c__87781, iwork, logwrk,
02779 result, &info);
02780 if (info != 0) {
02781 s_wsfe(&io___204);
02782 do_fio(&c__1, "SCHKHS", (ftnlen)6);
02783 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
02784 e_wsfe();
02785 }
02786
02787 }
02788
02789 } else if (lsamen_(&c__3, c3, "SST") || lsamen_(&
02790 c__3, c3, "SEP")) {
02791
02792
02793
02794
02795
02796
02797
02798
02799
02800 maxtyp = 21;
02801 ntypes = min(maxtyp,ntypes);
02802 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
02803 xlaenv_(&c__1, &c__1);
02804 xlaenv_(&c__9, &c__25);
02805 if (tsterr) {
02806 serrst_("SST", &c__6);
02807 }
02808 i__1 = nparms;
02809 for (i__ = 1; i__ <= i__1; ++i__) {
02810 xlaenv_(&c__1, &nbval[i__ - 1]);
02811 xlaenv_(&c__2, &nbmin[i__ - 1]);
02812 xlaenv_(&c__3, &nxval[i__ - 1]);
02813
02814 if (newsd == 0) {
02815 for (k = 1; k <= 4; ++k) {
02816 iseed[k - 1] = ioldsd[k - 1];
02817
02818 }
02819 }
02820 s_wsfe(&io___205);
02821 do_fio(&c__1, c3, (ftnlen)3);
02822 do_fio(&c__1, (char *)&nbval[i__ - 1], (ftnlen)sizeof(integer));
02823 do_fio(&c__1, (char *)&nbmin[i__ - 1], (ftnlen)sizeof(integer));
02824 do_fio(&c__1, (char *)&nxval[i__ - 1], (ftnlen)sizeof(integer));
02825 e_wsfe();
02826 if (tstchk) {
02827 schkst_(&nn, nval, &maxtyp, dotype, iseed, &thresh, &c__6, a,
02828 &c__132, &a[17424], d__, &d__[132], &d__[264], &d__[
02829 396], &d__[528], &d__[660], &d__[792], &d__[924], &
02830 d__[1056], &d__[1188], &d__[1320], &a[34848], &c__132,
02831 &a[52272], &a[69696], &d__[1452], &a[87120], work, &
02832 c__87781, iwork, &c__89760, result, &info);
02833 if (info != 0) {
02834 s_wsfe(&io___206);
02835 do_fio(&c__1, "SCHKST", (ftnlen)6);
02836 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
02837 e_wsfe();
02838 }
02839 }
02840 if (tstdrv) {
02841 sdrvst_(&nn, nval, &c__18, dotype, iseed, &thresh, &c__6, a, &
02842 c__132, &d__[264], &d__[396], &d__[528], &d__[660], &
02843 d__[924], &d__[1056], &d__[1188], &d__[1320], &a[
02844 17424], &c__132, &a[34848], &d__[1452], &a[52272],
02845 work, &c__87781, iwork, &c__89760, result, &info);
02846 if (info != 0) {
02847 s_wsfe(&io___207);
02848 do_fio(&c__1, "SDRVST", (ftnlen)6);
02849 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
02850 e_wsfe();
02851 }
02852 }
02853
02854 }
02855
02856 } else if (lsamen_(&c__3, c3, "SSG")) {
02857
02858
02859
02860
02861
02862
02863
02864
02865
02866 maxtyp = 21;
02867 ntypes = min(maxtyp,ntypes);
02868 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
02869 xlaenv_(&c__9, &c__25);
02870 i__1 = nparms;
02871 for (i__ = 1; i__ <= i__1; ++i__) {
02872 xlaenv_(&c__1, &nbval[i__ - 1]);
02873 xlaenv_(&c__2, &nbmin[i__ - 1]);
02874 xlaenv_(&c__3, &nxval[i__ - 1]);
02875
02876 if (newsd == 0) {
02877 for (k = 1; k <= 4; ++k) {
02878 iseed[k - 1] = ioldsd[k - 1];
02879
02880 }
02881 }
02882 s_wsfe(&io___208);
02883 do_fio(&c__1, c3, (ftnlen)3);
02884 do_fio(&c__1, (char *)&nbval[i__ - 1], (ftnlen)sizeof(integer));
02885 do_fio(&c__1, (char *)&nbmin[i__ - 1], (ftnlen)sizeof(integer));
02886 do_fio(&c__1, (char *)&nxval[i__ - 1], (ftnlen)sizeof(integer));
02887 e_wsfe();
02888 if (tstchk) {
02889 sdrvsg_(&nn, nval, &maxtyp, dotype, iseed, &thresh, &c__6, a,
02890 &c__132, &a[17424], &c__132, &d__[264], &a[34848], &
02891 c__132, &a[52272], &a[69696], &a[87120], &a[104544],
02892 work, &c__87781, iwork, &c__89760, result, &info);
02893 if (info != 0) {
02894 s_wsfe(&io___209);
02895 do_fio(&c__1, "SDRVSG", (ftnlen)6);
02896 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
02897 e_wsfe();
02898 }
02899 }
02900
02901 }
02902
02903 } else if (lsamen_(&c__3, c3, "SBD") || lsamen_(&
02904 c__3, c3, "SVD")) {
02905
02906
02907
02908
02909
02910
02911
02912
02913
02914
02915 maxtyp = 16;
02916 ntypes = min(maxtyp,ntypes);
02917 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
02918 xlaenv_(&c__1, &c__1);
02919 xlaenv_(&c__9, &c__25);
02920
02921
02922
02923 if (tsterr && tstchk) {
02924 serrbd_("SBD", &c__6);
02925 }
02926 if (tsterr && tstdrv) {
02927 serred_("SBD", &c__6);
02928 }
02929
02930 i__1 = nparms;
02931 for (i__ = 1; i__ <= i__1; ++i__) {
02932 nrhs = nsval[i__ - 1];
02933 xlaenv_(&c__1, &nbval[i__ - 1]);
02934 xlaenv_(&c__2, &nbmin[i__ - 1]);
02935 xlaenv_(&c__3, &nxval[i__ - 1]);
02936 if (newsd == 0) {
02937 for (k = 1; k <= 4; ++k) {
02938 iseed[k - 1] = ioldsd[k - 1];
02939
02940 }
02941 }
02942 s_wsfe(&io___211);
02943 do_fio(&c__1, c3, (ftnlen)3);
02944 do_fio(&c__1, (char *)&nbval[i__ - 1], (ftnlen)sizeof(integer));
02945 do_fio(&c__1, (char *)&nbmin[i__ - 1], (ftnlen)sizeof(integer));
02946 do_fio(&c__1, (char *)&nxval[i__ - 1], (ftnlen)sizeof(integer));
02947 do_fio(&c__1, (char *)&nrhs, (ftnlen)sizeof(integer));
02948 e_wsfe();
02949 if (tstchk) {
02950 schkbd_(&nn, mval, nval, &maxtyp, dotype, &nrhs, iseed, &
02951 thresh, a, &c__132, d__, &d__[132], &d__[264], &d__[
02952 396], &a[17424], &c__132, &a[34848], &a[52272], &a[
02953 69696], &c__132, &a[87120], &c__132, &a[104544], &a[
02954 121968], work, &c__87781, iwork, &c__6, &info);
02955 if (info != 0) {
02956 s_wsfe(&io___212);
02957 do_fio(&c__1, "SCHKBD", (ftnlen)6);
02958 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
02959 e_wsfe();
02960 }
02961 }
02962 if (tstdrv) {
02963 sdrvbd_(&nn, mval, nval, &maxtyp, dotype, iseed, &thresh, a, &
02964 c__132, &a[17424], &c__132, &a[34848], &c__132, &a[
02965 52272], &a[69696], &a[87120], d__, &d__[132], &d__[
02966 264], work, &c__87781, iwork, &c__6, &info);
02967 }
02968
02969 }
02970
02971 } else if (lsamen_(&c__3, c3, "SEV")) {
02972
02973
02974
02975
02976
02977
02978 maxtyp = 21;
02979 ntypes = min(maxtyp,ntypes);
02980 if (ntypes <= 0) {
02981 s_wsfe(&io___213);
02982 do_fio(&c__1, c3, (ftnlen)3);
02983 e_wsfe();
02984 } else {
02985 if (tsterr) {
02986 serred_(c3, &c__6);
02987 }
02988 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
02989 sdrvev_(&nn, nval, &ntypes, dotype, iseed, &thresh, &c__6, a, &
02990 c__132, &a[17424], d__, &d__[132], &d__[264], &d__[396], &
02991 a[34848], &c__132, &a[52272], &c__132, &a[69696], &c__132,
02992 result, work, &c__87781, iwork, &info);
02993 if (info != 0) {
02994 s_wsfe(&io___214);
02995 do_fio(&c__1, "SGEEV", (ftnlen)5);
02996 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
02997 e_wsfe();
02998 }
02999 }
03000 s_wsfe(&io___215);
03001 e_wsfe();
03002 goto L10;
03003
03004 } else if (lsamen_(&c__3, c3, "SES")) {
03005
03006
03007
03008
03009
03010
03011 maxtyp = 21;
03012 ntypes = min(maxtyp,ntypes);
03013 if (ntypes <= 0) {
03014 s_wsfe(&io___216);
03015 do_fio(&c__1, c3, (ftnlen)3);
03016 e_wsfe();
03017 } else {
03018 if (tsterr) {
03019 serred_(c3, &c__6);
03020 }
03021 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
03022 sdrves_(&nn, nval, &ntypes, dotype, iseed, &thresh, &c__6, a, &
03023 c__132, &a[17424], &a[34848], d__, &d__[132], &d__[264], &
03024 d__[396], &a[52272], &c__132, result, work, &c__87781,
03025 iwork, logwrk, &info);
03026 if (info != 0) {
03027 s_wsfe(&io___217);
03028 do_fio(&c__1, "SGEES", (ftnlen)5);
03029 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03030 e_wsfe();
03031 }
03032 }
03033 s_wsfe(&io___218);
03034 e_wsfe();
03035 goto L10;
03036
03037 } else if (lsamen_(&c__3, c3, "SVX")) {
03038
03039
03040
03041
03042
03043
03044 maxtyp = 21;
03045 ntypes = min(maxtyp,ntypes);
03046 if (ntypes < 0) {
03047 s_wsfe(&io___219);
03048 do_fio(&c__1, c3, (ftnlen)3);
03049 e_wsfe();
03050 } else {
03051 if (tsterr) {
03052 serred_(c3, &c__6);
03053 }
03054 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
03055 sdrvvx_(&nn, nval, &ntypes, dotype, iseed, &thresh, &c__5, &c__6,
03056 a, &c__132, &a[17424], d__, &d__[132], &d__[264], &d__[
03057 396], &a[34848], &c__132, &a[52272], &c__132, &a[69696], &
03058 c__132, &d__[528], &d__[660], &d__[792], &d__[924], &d__[
03059 1056], &d__[1188], &d__[1320], &d__[1452], result, work, &
03060 c__87781, iwork, &info);
03061 if (info != 0) {
03062 s_wsfe(&io___220);
03063 do_fio(&c__1, "SGEEVX", (ftnlen)6);
03064 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03065 e_wsfe();
03066 }
03067 }
03068 s_wsfe(&io___221);
03069 e_wsfe();
03070 goto L10;
03071
03072 } else if (lsamen_(&c__3, c3, "SSX")) {
03073
03074
03075
03076
03077
03078
03079 maxtyp = 21;
03080 ntypes = min(maxtyp,ntypes);
03081 if (ntypes < 0) {
03082 s_wsfe(&io___222);
03083 do_fio(&c__1, c3, (ftnlen)3);
03084 e_wsfe();
03085 } else {
03086 if (tsterr) {
03087 serred_(c3, &c__6);
03088 }
03089 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
03090 sdrvsx_(&nn, nval, &ntypes, dotype, iseed, &thresh, &c__5, &c__6,
03091 a, &c__132, &a[17424], &a[34848], d__, &d__[132], &d__[
03092 264], &d__[396], &d__[528], &d__[660], &a[52272], &c__132,
03093 &a[69696], result, work, &c__87781, iwork, logwrk, &info)
03094 ;
03095 if (info != 0) {
03096 s_wsfe(&io___223);
03097 do_fio(&c__1, "SGEESX", (ftnlen)6);
03098 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03099 e_wsfe();
03100 }
03101 }
03102 s_wsfe(&io___224);
03103 e_wsfe();
03104 goto L10;
03105
03106 } else if (lsamen_(&c__3, c3, "SGG")) {
03107
03108
03109
03110
03111
03112
03113
03114
03115
03116
03117
03118 maxtyp = 26;
03119 ntypes = min(maxtyp,ntypes);
03120 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
03121 if (tstchk && tsterr) {
03122 serrgg_(c3, &c__6);
03123 }
03124 i__1 = nparms;
03125 for (i__ = 1; i__ <= i__1; ++i__) {
03126 xlaenv_(&c__1, &nbval[i__ - 1]);
03127 xlaenv_(&c__2, &nbmin[i__ - 1]);
03128 xlaenv_(&c__4, &nsval[i__ - 1]);
03129 xlaenv_(&c__8, &mxbval[i__ - 1]);
03130 xlaenv_(&c__5, &nbcol[i__ - 1]);
03131
03132 if (newsd == 0) {
03133 for (k = 1; k <= 4; ++k) {
03134 iseed[k - 1] = ioldsd[k - 1];
03135
03136 }
03137 }
03138 s_wsfe(&io___225);
03139 do_fio(&c__1, c3, (ftnlen)3);
03140 do_fio(&c__1, (char *)&nbval[i__ - 1], (ftnlen)sizeof(integer));
03141 do_fio(&c__1, (char *)&nbmin[i__ - 1], (ftnlen)sizeof(integer));
03142 do_fio(&c__1, (char *)&nsval[i__ - 1], (ftnlen)sizeof(integer));
03143 do_fio(&c__1, (char *)&mxbval[i__ - 1], (ftnlen)sizeof(integer));
03144 do_fio(&c__1, (char *)&nbcol[i__ - 1], (ftnlen)sizeof(integer));
03145 e_wsfe();
03146 tstdif = FALSE_;
03147 thrshn = 10.f;
03148 if (tstchk) {
03149 schkgg_(&nn, nval, &maxtyp, dotype, iseed, &thresh, &tstdif, &
03150 thrshn, &c__6, a, &c__132, &a[17424], &a[34848], &a[
03151 52272], &a[69696], &a[87120], &a[104544], &a[121968],
03152 &a[139392], &c__132, &a[156816], &a[174240], &a[
03153 191664], d__, &d__[132], &d__[264], &d__[396], &d__[
03154 528], &d__[660], &a[209088], &a[226512], work, &
03155 c__87781, logwrk, result, &info);
03156 if (info != 0) {
03157 s_wsfe(&io___228);
03158 do_fio(&c__1, "SCHKGG", (ftnlen)6);
03159 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03160 e_wsfe();
03161 }
03162 }
03163 xlaenv_(&c__1, &c__1);
03164 if (tstdrv) {
03165 sdrvgg_(&nn, nval, &maxtyp, dotype, iseed, &thresh, &thrshn, &
03166 c__6, a, &c__132, &a[17424], &a[34848], &a[52272], &a[
03167 69696], &a[87120], &a[104544], &c__132, &a[121968],
03168 d__, &d__[132], &d__[264], &d__[396], &d__[528], &d__[
03169 660], &a[209088], &a[226512], work, &c__87781, result,
03170 &info);
03171 if (info != 0) {
03172 s_wsfe(&io___229);
03173 do_fio(&c__1, "SDRVGG", (ftnlen)6);
03174 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03175 e_wsfe();
03176 }
03177 }
03178
03179 }
03180
03181 } else if (lsamen_(&c__3, c3, "SGS")) {
03182
03183
03184
03185
03186
03187
03188 maxtyp = 26;
03189 ntypes = min(maxtyp,ntypes);
03190 if (ntypes <= 0) {
03191 s_wsfe(&io___230);
03192 do_fio(&c__1, c3, (ftnlen)3);
03193 e_wsfe();
03194 } else {
03195 if (tsterr) {
03196 serrgg_(c3, &c__6);
03197 }
03198 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
03199 sdrges_(&nn, nval, &maxtyp, dotype, iseed, &thresh, &c__6, a, &
03200 c__132, &a[17424], &a[34848], &a[52272], &a[104544], &
03201 c__132, &a[121968], d__, &d__[132], &d__[264], work, &
03202 c__87781, result, logwrk, &info);
03203
03204 if (info != 0) {
03205 s_wsfe(&io___231);
03206 do_fio(&c__1, "SDRGES", (ftnlen)6);
03207 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03208 e_wsfe();
03209 }
03210 }
03211 s_wsfe(&io___232);
03212 e_wsfe();
03213 goto L10;
03214
03215 } else if (sgx) {
03216
03217
03218
03219
03220
03221
03222 maxtyp = 5;
03223 ntypes = maxtyp;
03224 if (nn < 0) {
03225 s_wsfe(&io___233);
03226 do_fio(&c__1, c3, (ftnlen)3);
03227 e_wsfe();
03228 } else {
03229 if (tsterr) {
03230 serrgg_(c3, &c__6);
03231 }
03232 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
03233 xlaenv_(&c__5, &c__2);
03234 sdrgsx_(&nn, &c__20, &thresh, &c__5, &c__6, a, &c__132, &a[17424],
03235 &a[34848], &a[52272], &a[69696], &a[87120], d__, &d__[
03236 132], &d__[264], c__, &c__400, &a[191664], work, &
03237 c__87781, iwork, &c__89760, logwrk, &info);
03238 if (info != 0) {
03239 s_wsfe(&io___235);
03240 do_fio(&c__1, "SDRGSX", (ftnlen)6);
03241 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03242 e_wsfe();
03243 }
03244 }
03245 s_wsfe(&io___236);
03246 e_wsfe();
03247 goto L10;
03248
03249 } else if (lsamen_(&c__3, c3, "SGV")) {
03250
03251
03252
03253
03254
03255
03256 maxtyp = 26;
03257 ntypes = min(maxtyp,ntypes);
03258 if (ntypes <= 0) {
03259 s_wsfe(&io___237);
03260 do_fio(&c__1, c3, (ftnlen)3);
03261 e_wsfe();
03262 } else {
03263 if (tsterr) {
03264 serrgg_(c3, &c__6);
03265 }
03266 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
03267 sdrgev_(&nn, nval, &maxtyp, dotype, iseed, &thresh, &c__6, a, &
03268 c__132, &a[17424], &a[34848], &a[52272], &a[104544], &
03269 c__132, &a[121968], &a[139392], &c__132, d__, &d__[132], &
03270 d__[264], &d__[396], &d__[528], &d__[660], work, &
03271 c__87781, result, &info);
03272 if (info != 0) {
03273 s_wsfe(&io___238);
03274 do_fio(&c__1, "SDRGEV", (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 (sxv) {
03284
03285
03286
03287
03288
03289
03290 maxtyp = 2;
03291 ntypes = maxtyp;
03292 if (nn < 0) {
03293 s_wsfe(&io___240);
03294 do_fio(&c__1, c3, (ftnlen)3);
03295 e_wsfe();
03296 } else {
03297 if (tsterr) {
03298 serrgg_(c3, &c__6);
03299 }
03300 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
03301 sdrgvx_(&nn, &thresh, &c__5, &c__6, a, &c__132, &a[17424], &a[
03302 34848], &a[52272], d__, &d__[132], &d__[264], &a[69696], &
03303 a[87120], iwork, &iwork[1], &d__[396], &d__[528], &d__[
03304 660], &d__[792], &d__[924], &d__[1056], work, &c__87781, &
03305 iwork[2], &c__89758, result, logwrk, &info);
03306
03307 if (info != 0) {
03308 s_wsfe(&io___241);
03309 do_fio(&c__1, "SDRGVX", (ftnlen)6);
03310 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03311 e_wsfe();
03312 }
03313 }
03314 s_wsfe(&io___242);
03315 e_wsfe();
03316 goto L10;
03317
03318 } else if (lsamen_(&c__3, c3, "SSB")) {
03319
03320
03321
03322
03323
03324 maxtyp = 15;
03325 ntypes = min(maxtyp,ntypes);
03326 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
03327 if (tsterr) {
03328 serrst_("SSB", &c__6);
03329 }
03330 schksb_(&nn, nval, &nk, kval, &maxtyp, dotype, iseed, &thresh, &c__6,
03331 a, &c__132, d__, &d__[132], &a[17424], &c__132, work, &
03332 c__87781, result, &info);
03333 if (info != 0) {
03334 s_wsfe(&io___243);
03335 do_fio(&c__1, "SCHKSB", (ftnlen)6);
03336 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03337 e_wsfe();
03338 }
03339
03340 } else if (lsamen_(&c__3, c3, "SBB")) {
03341
03342
03343
03344
03345
03346 maxtyp = 15;
03347 ntypes = min(maxtyp,ntypes);
03348 alareq_(c3, &ntypes, dotype, &maxtyp, &c__5, &c__6);
03349 i__1 = nparms;
03350 for (i__ = 1; i__ <= i__1; ++i__) {
03351 nrhs = nsval[i__ - 1];
03352
03353 if (newsd == 0) {
03354 for (k = 1; k <= 4; ++k) {
03355 iseed[k - 1] = ioldsd[k - 1];
03356
03357 }
03358 }
03359 s_wsfe(&io___244);
03360 do_fio(&c__1, c3, (ftnlen)3);
03361 do_fio(&c__1, (char *)&nrhs, (ftnlen)sizeof(integer));
03362 e_wsfe();
03363 schkbb_(&nn, mval, nval, &nk, kval, &maxtyp, dotype, &nrhs, iseed,
03364 &thresh, &c__6, a, &c__132, &a[17424], &c__264, d__, &
03365 d__[132], &a[52272], &c__132, &a[69696], &c__132, &a[
03366 87120], &c__132, &a[104544], work, &c__87781, result, &
03367 info);
03368 if (info != 0) {
03369 s_wsfe(&io___245);
03370 do_fio(&c__1, "SCHKBB", (ftnlen)6);
03371 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03372 e_wsfe();
03373 }
03374
03375 }
03376
03377 } else if (lsamen_(&c__3, c3, "GLM")) {
03378
03379
03380
03381
03382
03383 xlaenv_(&c__1, &c__1);
03384 if (tsterr) {
03385 serrgg_("GLM", &c__6);
03386 }
03387 sckglm_(&nn, mval, pval, nval, &ntypes, iseed, &thresh, &c__132, a, &
03388 a[17424], b, &b[17424], x, work, d__, &c__5, &c__6, &info);
03389 if (info != 0) {
03390 s_wsfe(&io___248);
03391 do_fio(&c__1, "SCKGLM", (ftnlen)6);
03392 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03393 e_wsfe();
03394 }
03395
03396 } else if (lsamen_(&c__3, c3, "GQR")) {
03397
03398
03399
03400
03401
03402 xlaenv_(&c__1, &c__1);
03403 if (tsterr) {
03404 serrgg_("GQR", &c__6);
03405 }
03406 sckgqr_(&nn, mval, &nn, pval, &nn, nval, &ntypes, iseed, &thresh, &
03407 c__132, a, &a[17424], &a[34848], &a[52272], taua, b, &b[17424]
03408 , &b[34848], &b[52272], &b[69696], taub, work, d__, &c__5, &
03409 c__6, &info);
03410 if (info != 0) {
03411 s_wsfe(&io___251);
03412 do_fio(&c__1, "SCKGQR", (ftnlen)6);
03413 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03414 e_wsfe();
03415 }
03416
03417 } else if (lsamen_(&c__3, c3, "GSV")) {
03418
03419
03420
03421
03422
03423 if (tsterr) {
03424 serrgg_("GSV", &c__6);
03425 }
03426 sckgsv_(&nn, mval, pval, nval, &ntypes, iseed, &thresh, &c__132, a, &
03427 a[17424], b, &b[17424], &a[34848], &b[34848], &a[52272], taua,
03428 taub, &b[52272], iwork, work, d__, &c__5, &c__6, &info);
03429 if (info != 0) {
03430 s_wsfe(&io___252);
03431 do_fio(&c__1, "SCKGSV", (ftnlen)6);
03432 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03433 e_wsfe();
03434 }
03435
03436 } else if (lsamen_(&c__3, c3, "LSE")) {
03437
03438
03439
03440
03441
03442 xlaenv_(&c__1, &c__1);
03443 if (tsterr) {
03444 serrgg_("LSE", &c__6);
03445 }
03446 scklse_(&nn, mval, pval, nval, &ntypes, iseed, &thresh, &c__132, a, &
03447 a[17424], b, &b[17424], x, work, d__, &c__5, &c__6, &info);
03448 if (info != 0) {
03449 s_wsfe(&io___253);
03450 do_fio(&c__1, "SCKLSE", (ftnlen)6);
03451 do_fio(&c__1, (char *)&info, (ftnlen)sizeof(integer));
03452 e_wsfe();
03453 }
03454
03455 } else {
03456 s_wsle(&io___254);
03457 e_wsle();
03458 s_wsle(&io___255);
03459 e_wsle();
03460 s_wsfe(&io___256);
03461 do_fio(&c__1, c3, (ftnlen)3);
03462 e_wsfe();
03463 }
03464 if (! (sgx || sxv)) {
03465 goto L190;
03466 }
03467 L380:
03468 s_wsfe(&io___257);
03469 e_wsfe();
03470 s2 = second_();
03471 s_wsfe(&io___259);
03472 r__1 = s2 - s1;
03473 do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
03474 e_wsfe();
03475
03476
03477
03478
03479
03480 return 0;
03481 }
03482
03483 int schkee_ () { MAIN__ (); return 0; }