00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013 #include "f2c.h"
00014 #include "blaswrap.h"
00015
00016
00017
00018 struct {
00019 integer infot, nout;
00020 logical ok, lerr;
00021 } infoc_;
00022
00023 #define infoc_1 infoc_
00024
00025 struct {
00026 char srnamt[32];
00027 } srnamc_;
00028
00029 #define srnamc_1 srnamc_
00030
00031
00032
00033 static integer c__2 = 2;
00034 static integer c__0 = 0;
00035 static integer c__1 = 1;
00036 static integer c_n1 = -1;
00037 static integer c__18 = 18;
00038 static integer c__3 = 3;
00039 static integer c__32 = 32;
00040 static logical c_true = TRUE_;
00041 static logical c_false = FALSE_;
00042 static integer c__20 = 20;
00043
00044 int derrgg_(char *path, integer *nunit)
00045 {
00046
00047 static char fmt_9999[] = "(1x,a3,\002 routines passed the tests of the e"
00048 "rror exits (\002,i3,\002 tests done)\002)";
00049 static char fmt_9998[] = "(\002 *** \002,a3,\002 routines failed the tes"
00050 "ts of the error \002,\002exits ***\002)";
00051
00052
00053 integer s_wsle(cilist *), e_wsle(void);
00054 int s_copy(char *, char *, ftnlen, ftnlen);
00055 integer s_wsfe(cilist *), do_fio(integer *, char *, ftnlen), e_wsfe(void);
00056
00057
00058 doublereal a[9] , b[9] ;
00059 integer i__, j, m;
00060 doublereal q[9] , u[9] , v[9]
00061 , w[18], z__[9] ;
00062 char c2[2];
00063 doublereal r1[3], r2[3], r3[3];
00064 logical bw[3];
00065 doublereal ls[3];
00066 integer iw[3], nt;
00067 doublereal rs[3], dif, rce[2];
00068 logical sel[3];
00069 doublereal tau[3], rcv[2];
00070 integer info, sdim;
00071 doublereal anrm, bnrm, tola, tolb;
00072 integer ifst, ilst;
00073 doublereal scale;
00074 extern int dgges_(char *, char *, char *, L_fp, integer *
00075 , doublereal *, integer *, doublereal *, integer *, integer *,
00076 doublereal *, doublereal *, doublereal *, doublereal *, integer *,
00077 doublereal *, integer *, doublereal *, integer *, logical *,
00078 integer *), dggev_(char *, char *,
00079 integer *, doublereal *, integer *, doublereal *, integer *,
00080 doublereal *, doublereal *, doublereal *, doublereal *, integer *,
00081 doublereal *, integer *, doublereal *, integer *, integer *), dgghrd_(char *, char *, integer *, integer *,
00082 integer *, doublereal *, integer *, doublereal *, integer *,
00083 doublereal *, integer *, doublereal *, integer *, integer *), dggglm_(integer *, integer *, integer *,
00084 doublereal *, integer *, doublereal *, integer *, doublereal *,
00085 doublereal *, doublereal *, doublereal *, integer *, integer *),
00086 dgglse_(integer *, integer *, integer *, doublereal *, integer *,
00087 doublereal *, integer *, doublereal *, doublereal *, doublereal *,
00088 doublereal *, integer *, integer *), dggqrf_(integer *, integer *
00089 , integer *, doublereal *, integer *, doublereal *, doublereal *,
00090 integer *, doublereal *, doublereal *, integer *, integer *),
00091 dggrqf_(integer *, integer *, integer *, doublereal *, integer *,
00092 doublereal *, doublereal *, integer *, doublereal *, doublereal *,
00093 integer *, integer *);
00094 integer ncycle;
00095 extern logical dlctes_(), lsamen_(integer *, char *, char *);
00096 extern int dggsvd_(char *, char *, char *, integer *,
00097 integer *, integer *, integer *, integer *, doublereal *, integer
00098 *, doublereal *, integer *, doublereal *, doublereal *,
00099 doublereal *, integer *, doublereal *, integer *, doublereal *,
00100 integer *, doublereal *, integer *, integer *), dggesx_(char *, char *, char *, L_fp, char *, integer *,
00101 doublereal *, integer *, doublereal *, integer *, integer *,
00102 doublereal *, doublereal *, doublereal *, doublereal *, integer *,
00103 doublereal *, integer *, doublereal *, doublereal *, doublereal *
00104 , integer *, integer *, integer *, logical *, integer *), dhgeqz_(char *, char *, char *, integer *
00105 , integer *, integer *, doublereal *, integer *, doublereal *,
00106 integer *, doublereal *, doublereal *, doublereal *, doublereal *,
00107 integer *, doublereal *, integer *, doublereal *, integer *,
00108 integer *), dtgevc_(char *, char *,
00109 logical *, integer *, doublereal *, integer *, doublereal *,
00110 integer *, doublereal *, integer *, doublereal *, integer *,
00111 integer *, integer *, doublereal *, integer *),
00112 chkxer_(char *, integer *, integer *, logical *, logical *), dggevx_(char *, char *, char *, char *, integer *,
00113 doublereal *, integer *, doublereal *, integer *, doublereal *,
00114 doublereal *, doublereal *, doublereal *, integer *, doublereal *,
00115 integer *, integer *, integer *, doublereal *, doublereal *,
00116 doublereal *, doublereal *, doublereal *, doublereal *,
00117 doublereal *, integer *, integer *, logical *, integer *), dtgexc_(logical *, logical *, integer *,
00118 doublereal *, integer *, doublereal *, integer *, doublereal *,
00119 integer *, doublereal *, integer *, integer *, integer *,
00120 doublereal *, integer *, integer *), dtgsen_(integer *, logical *,
00121 logical *, logical *, integer *, doublereal *, integer *,
00122 doublereal *, integer *, doublereal *, doublereal *, doublereal *,
00123 doublereal *, integer *, doublereal *, integer *, integer *,
00124 doublereal *, doublereal *, doublereal *, doublereal *, integer *,
00125 integer *, integer *, integer *), dtgsja_(char *, char *, char *,
00126 integer *, integer *, integer *, integer *, integer *,
00127 doublereal *, integer *, doublereal *, integer *, doublereal *,
00128 doublereal *, doublereal *, doublereal *, doublereal *, integer *,
00129 doublereal *, integer *, doublereal *, integer *, doublereal *,
00130 integer *, integer *), dtgsna_(char *,
00131 char *, logical *, integer *, doublereal *, integer *, doublereal
00132 *, integer *, doublereal *, integer *, doublereal *, integer *,
00133 doublereal *, doublereal *, integer *, integer *, doublereal *,
00134 integer *, integer *, integer *), dggsvp_(char *,
00135 char *, char *, integer *, integer *, integer *, doublereal *,
00136 integer *, doublereal *, integer *, doublereal *, doublereal *,
00137 integer *, integer *, doublereal *, integer *, doublereal *,
00138 integer *, doublereal *, integer *, integer *, doublereal *,
00139 doublereal *, integer *);
00140 extern logical dlctsx_();
00141 integer dummyk, dummyl;
00142 extern int dtgsyl_(char *, integer *, integer *, integer
00143 *, doublereal *, integer *, doublereal *, integer *, doublereal *,
00144 integer *, doublereal *, integer *, doublereal *, integer *,
00145 doublereal *, integer *, doublereal *, doublereal *, doublereal *,
00146 integer *, integer *, integer *);
00147
00148
00149 static cilist io___1 = { 0, 0, 0, 0, 0 };
00150 static cilist io___38 = { 0, 0, 0, fmt_9999, 0 };
00151 static cilist io___39 = { 0, 0, 0, fmt_9998, 0 };
00152
00153
00154
00155
00156
00157
00158
00159
00160
00161
00162
00163
00164
00165
00166
00167
00168
00169
00170
00171
00172
00173
00174
00175
00176
00177
00178
00179
00180
00181
00182
00183
00184
00185
00186
00187
00188
00189
00190
00191
00192
00193
00194
00195
00196 infoc_1.nout = *nunit;
00197 io___1.ciunit = infoc_1.nout;
00198 s_wsle(&io___1);
00199 e_wsle();
00200 s_copy(c2, path + 1, (ftnlen)2, (ftnlen)2);
00201
00202
00203
00204 for (j = 1; j <= 3; ++j) {
00205 sel[j - 1] = TRUE_;
00206 for (i__ = 1; i__ <= 3; ++i__) {
00207 a[i__ + j * 3 - 4] = 0.;
00208 b[i__ + j * 3 - 4] = 0.;
00209
00210 }
00211
00212 }
00213 for (i__ = 1; i__ <= 3; ++i__) {
00214 a[i__ + i__ * 3 - 4] = 1.;
00215 b[i__ + i__ * 3 - 4] = 1.;
00216
00217 }
00218 infoc_1.ok = TRUE_;
00219 tola = 1.;
00220 tolb = 1.;
00221 ifst = 1;
00222 ilst = 1;
00223 nt = 0;
00224
00225
00226
00227 if (lsamen_(&c__2, c2, "GG")) {
00228
00229
00230
00231 s_copy(srnamc_1.srnamt, "DGGHRD", (ftnlen)32, (ftnlen)6);
00232 infoc_1.infot = 1;
00233 dgghrd_("/", "N", &c__0, &c__1, &c__0, a, &c__1, b, &c__1, q, &c__1,
00234 z__, &c__1, &info);
00235 chkxer_("DGGHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00236 infoc_1.ok);
00237 infoc_1.infot = 2;
00238 dgghrd_("N", "/", &c__0, &c__1, &c__0, a, &c__1, b, &c__1, q, &c__1,
00239 z__, &c__1, &info);
00240 chkxer_("DGGHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00241 infoc_1.ok);
00242 infoc_1.infot = 3;
00243 dgghrd_("N", "N", &c_n1, &c__0, &c__0, a, &c__1, b, &c__1, q, &c__1,
00244 z__, &c__1, &info);
00245 chkxer_("DGGHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00246 infoc_1.ok);
00247 infoc_1.infot = 4;
00248 dgghrd_("N", "N", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, q, &c__1,
00249 z__, &c__1, &info);
00250 chkxer_("DGGHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00251 infoc_1.ok);
00252 infoc_1.infot = 5;
00253 dgghrd_("N", "N", &c__0, &c__1, &c__1, a, &c__1, b, &c__1, q, &c__1,
00254 z__, &c__1, &info);
00255 chkxer_("DGGHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00256 infoc_1.ok);
00257 infoc_1.infot = 7;
00258 dgghrd_("N", "N", &c__2, &c__1, &c__1, a, &c__1, b, &c__2, q, &c__1,
00259 z__, &c__1, &info);
00260 chkxer_("DGGHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00261 infoc_1.ok);
00262 infoc_1.infot = 9;
00263 dgghrd_("N", "N", &c__2, &c__1, &c__1, a, &c__2, b, &c__1, q, &c__1,
00264 z__, &c__1, &info);
00265 chkxer_("DGGHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00266 infoc_1.ok);
00267 infoc_1.infot = 11;
00268 dgghrd_("V", "N", &c__2, &c__1, &c__1, a, &c__2, b, &c__2, q, &c__1,
00269 z__, &c__1, &info);
00270 chkxer_("DGGHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00271 infoc_1.ok);
00272 infoc_1.infot = 13;
00273 dgghrd_("N", "V", &c__2, &c__1, &c__1, a, &c__2, b, &c__2, q, &c__1,
00274 z__, &c__1, &info);
00275 chkxer_("DGGHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00276 infoc_1.ok);
00277 nt += 9;
00278
00279
00280
00281 s_copy(srnamc_1.srnamt, "DHGEQZ", (ftnlen)32, (ftnlen)6);
00282 infoc_1.infot = 1;
00283 dhgeqz_("/", "N", "N", &c__0, &c__1, &c__0, a, &c__1, b, &c__1, r1,
00284 r2, r3, q, &c__1, z__, &c__1, w, &c__18, &info);
00285 chkxer_("DHGEQZ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00286 infoc_1.ok);
00287 infoc_1.infot = 2;
00288 dhgeqz_("E", "/", "N", &c__0, &c__1, &c__0, a, &c__1, b, &c__1, r1,
00289 r2, r3, q, &c__1, z__, &c__1, w, &c__18, &info);
00290 chkxer_("DHGEQZ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00291 infoc_1.ok);
00292 infoc_1.infot = 3;
00293 dhgeqz_("E", "N", "/", &c__0, &c__1, &c__0, a, &c__1, b, &c__1, r1,
00294 r2, r3, q, &c__1, z__, &c__1, w, &c__18, &info);
00295 chkxer_("DHGEQZ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00296 infoc_1.ok);
00297 infoc_1.infot = 4;
00298 dhgeqz_("E", "N", "N", &c_n1, &c__0, &c__0, a, &c__1, b, &c__1, r1,
00299 r2, r3, q, &c__1, z__, &c__1, w, &c__18, &info);
00300 chkxer_("DHGEQZ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00301 infoc_1.ok);
00302 infoc_1.infot = 5;
00303 dhgeqz_("E", "N", "N", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, r1,
00304 r2, r3, q, &c__1, z__, &c__1, w, &c__18, &info);
00305 chkxer_("DHGEQZ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00306 infoc_1.ok);
00307 infoc_1.infot = 6;
00308 dhgeqz_("E", "N", "N", &c__0, &c__1, &c__1, a, &c__1, b, &c__1, r1,
00309 r2, r3, q, &c__1, z__, &c__1, w, &c__18, &info);
00310 chkxer_("DHGEQZ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00311 infoc_1.ok);
00312 infoc_1.infot = 8;
00313 dhgeqz_("E", "N", "N", &c__2, &c__1, &c__1, a, &c__1, b, &c__2, r1,
00314 r2, r3, q, &c__1, z__, &c__1, w, &c__18, &info);
00315 chkxer_("DHGEQZ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00316 infoc_1.ok);
00317 infoc_1.infot = 10;
00318 dhgeqz_("E", "N", "N", &c__2, &c__1, &c__1, a, &c__2, b, &c__1, r1,
00319 r2, r3, q, &c__1, z__, &c__1, w, &c__18, &info);
00320 chkxer_("DHGEQZ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00321 infoc_1.ok);
00322 infoc_1.infot = 15;
00323 dhgeqz_("E", "V", "N", &c__2, &c__1, &c__1, a, &c__2, b, &c__2, r1,
00324 r2, r3, q, &c__1, z__, &c__1, w, &c__18, &info);
00325 chkxer_("DHGEQZ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00326 infoc_1.ok);
00327 infoc_1.infot = 17;
00328 dhgeqz_("E", "N", "V", &c__2, &c__1, &c__1, a, &c__2, b, &c__2, r1,
00329 r2, r3, q, &c__1, z__, &c__1, w, &c__18, &info);
00330 chkxer_("DHGEQZ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00331 infoc_1.ok);
00332 nt += 10;
00333
00334
00335
00336 s_copy(srnamc_1.srnamt, "DTGEVC", (ftnlen)32, (ftnlen)6);
00337 infoc_1.infot = 1;
00338 dtgevc_("/", "A", sel, &c__0, a, &c__1, b, &c__1, q, &c__1, z__, &
00339 c__1, &c__0, &m, w, &info);
00340 chkxer_("DTGEVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00341 infoc_1.ok);
00342 infoc_1.infot = 2;
00343 dtgevc_("R", "/", sel, &c__0, a, &c__1, b, &c__1, q, &c__1, z__, &
00344 c__1, &c__0, &m, w, &info);
00345 chkxer_("DTGEVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00346 infoc_1.ok);
00347 infoc_1.infot = 4;
00348 dtgevc_("R", "A", sel, &c_n1, a, &c__1, b, &c__1, q, &c__1, z__, &
00349 c__1, &c__0, &m, w, &info);
00350 chkxer_("DTGEVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00351 infoc_1.ok);
00352 infoc_1.infot = 6;
00353 dtgevc_("R", "A", sel, &c__2, a, &c__1, b, &c__2, q, &c__1, z__, &
00354 c__2, &c__0, &m, w, &info);
00355 chkxer_("DTGEVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00356 infoc_1.ok);
00357 infoc_1.infot = 8;
00358 dtgevc_("R", "A", sel, &c__2, a, &c__2, b, &c__1, q, &c__1, z__, &
00359 c__2, &c__0, &m, w, &info);
00360 chkxer_("DTGEVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00361 infoc_1.ok);
00362 infoc_1.infot = 10;
00363 dtgevc_("L", "A", sel, &c__2, a, &c__2, b, &c__2, q, &c__1, z__, &
00364 c__1, &c__0, &m, w, &info);
00365 chkxer_("DTGEVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00366 infoc_1.ok);
00367 infoc_1.infot = 12;
00368 dtgevc_("R", "A", sel, &c__2, a, &c__2, b, &c__2, q, &c__1, z__, &
00369 c__1, &c__0, &m, w, &info);
00370 chkxer_("DTGEVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00371 infoc_1.ok);
00372 infoc_1.infot = 13;
00373 dtgevc_("R", "A", sel, &c__2, a, &c__2, b, &c__2, q, &c__1, z__, &
00374 c__2, &c__1, &m, w, &info);
00375 chkxer_("DTGEVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00376 infoc_1.ok);
00377 nt += 8;
00378
00379
00380
00381 } else if (lsamen_(&c__3, path, "GSV")) {
00382
00383
00384
00385 s_copy(srnamc_1.srnamt, "DGGSVD", (ftnlen)32, (ftnlen)6);
00386 infoc_1.infot = 1;
00387 dggsvd_("/", "N", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
00388 c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
00389 info);
00390 chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00391 infoc_1.ok);
00392 infoc_1.infot = 2;
00393 dggsvd_("N", "/", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
00394 c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
00395 info);
00396 chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00397 infoc_1.ok);
00398 infoc_1.infot = 3;
00399 dggsvd_("N", "N", "/", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
00400 c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
00401 info);
00402 chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00403 infoc_1.ok);
00404 infoc_1.infot = 4;
00405 dggsvd_("N", "N", "N", &c_n1, &c__0, &c__0, &dummyk, &dummyl, a, &
00406 c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
00407 info);
00408 chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00409 infoc_1.ok);
00410 infoc_1.infot = 5;
00411 dggsvd_("N", "N", "N", &c__0, &c_n1, &c__0, &dummyk, &dummyl, a, &
00412 c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
00413 info);
00414 chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00415 infoc_1.ok);
00416 infoc_1.infot = 6;
00417 dggsvd_("N", "N", "N", &c__0, &c__0, &c_n1, &dummyk, &dummyl, a, &
00418 c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
00419 info);
00420 chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00421 infoc_1.ok);
00422 infoc_1.infot = 10;
00423 dggsvd_("N", "N", "N", &c__2, &c__1, &c__1, &dummyk, &dummyl, a, &
00424 c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
00425 info);
00426 chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00427 infoc_1.ok);
00428 infoc_1.infot = 12;
00429 dggsvd_("N", "N", "N", &c__1, &c__1, &c__2, &dummyk, &dummyl, a, &
00430 c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
00431 info);
00432 chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00433 infoc_1.ok);
00434 infoc_1.infot = 16;
00435 dggsvd_("U", "N", "N", &c__2, &c__2, &c__2, &dummyk, &dummyl, a, &
00436 c__2, b, &c__2, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
00437 info);
00438 chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00439 infoc_1.ok);
00440 infoc_1.infot = 18;
00441 dggsvd_("N", "V", "N", &c__1, &c__1, &c__2, &dummyk, &dummyl, a, &
00442 c__1, b, &c__2, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
00443 info);
00444 chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00445 infoc_1.ok);
00446 infoc_1.infot = 20;
00447 dggsvd_("N", "N", "Q", &c__1, &c__2, &c__1, &dummyk, &dummyl, a, &
00448 c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
00449 info);
00450 chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00451 infoc_1.ok);
00452 nt += 11;
00453
00454
00455
00456 s_copy(srnamc_1.srnamt, "DGGSVP", (ftnlen)32, (ftnlen)6);
00457 infoc_1.infot = 1;
00458 dggsvp_("/", "N", "N", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, &tola,
00459 &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
00460 tau, w, &info);
00461 chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00462 infoc_1.ok);
00463 infoc_1.infot = 2;
00464 dggsvp_("N", "/", "N", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, &tola,
00465 &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
00466 tau, w, &info);
00467 chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00468 infoc_1.ok);
00469 infoc_1.infot = 3;
00470 dggsvp_("N", "N", "/", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, &tola,
00471 &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
00472 tau, w, &info);
00473 chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00474 infoc_1.ok);
00475 infoc_1.infot = 4;
00476 dggsvp_("N", "N", "N", &c_n1, &c__0, &c__0, a, &c__1, b, &c__1, &tola,
00477 &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
00478 tau, w, &info);
00479 chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00480 infoc_1.ok);
00481 infoc_1.infot = 5;
00482 dggsvp_("N", "N", "N", &c__0, &c_n1, &c__0, a, &c__1, b, &c__1, &tola,
00483 &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
00484 tau, w, &info);
00485 chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00486 infoc_1.ok);
00487 infoc_1.infot = 6;
00488 dggsvp_("N", "N", "N", &c__0, &c__0, &c_n1, a, &c__1, b, &c__1, &tola,
00489 &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
00490 tau, w, &info);
00491 chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00492 infoc_1.ok);
00493 infoc_1.infot = 8;
00494 dggsvp_("N", "N", "N", &c__2, &c__1, &c__1, a, &c__1, b, &c__1, &tola,
00495 &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
00496 tau, w, &info);
00497 chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00498 infoc_1.ok);
00499 infoc_1.infot = 10;
00500 dggsvp_("N", "N", "N", &c__1, &c__2, &c__1, a, &c__1, b, &c__1, &tola,
00501 &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
00502 tau, w, &info);
00503 chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00504 infoc_1.ok);
00505 infoc_1.infot = 16;
00506 dggsvp_("U", "N", "N", &c__2, &c__2, &c__2, a, &c__2, b, &c__2, &tola,
00507 &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
00508 tau, w, &info);
00509 chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00510 infoc_1.ok);
00511 infoc_1.infot = 18;
00512 dggsvp_("N", "V", "N", &c__1, &c__2, &c__1, a, &c__1, b, &c__2, &tola,
00513 &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
00514 tau, w, &info);
00515 chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00516 infoc_1.ok);
00517 infoc_1.infot = 20;
00518 dggsvp_("N", "N", "Q", &c__1, &c__1, &c__2, a, &c__1, b, &c__1, &tola,
00519 &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
00520 tau, w, &info);
00521 chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00522 infoc_1.ok);
00523 nt += 11;
00524
00525
00526
00527 s_copy(srnamc_1.srnamt, "DTGSJA", (ftnlen)32, (ftnlen)6);
00528 infoc_1.infot = 1;
00529 dtgsja_("/", "N", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
00530 c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
00531 c__1, w, &ncycle, &info);
00532 chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00533 infoc_1.ok);
00534 infoc_1.infot = 2;
00535 dtgsja_("N", "/", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
00536 c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
00537 c__1, w, &ncycle, &info);
00538 chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00539 infoc_1.ok);
00540 infoc_1.infot = 3;
00541 dtgsja_("N", "N", "/", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
00542 c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
00543 c__1, w, &ncycle, &info);
00544 chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00545 infoc_1.ok);
00546 infoc_1.infot = 4;
00547 dtgsja_("N", "N", "N", &c_n1, &c__0, &c__0, &dummyk, &dummyl, a, &
00548 c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
00549 c__1, w, &ncycle, &info);
00550 chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00551 infoc_1.ok);
00552 infoc_1.infot = 5;
00553 dtgsja_("N", "N", "N", &c__0, &c_n1, &c__0, &dummyk, &dummyl, a, &
00554 c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
00555 c__1, w, &ncycle, &info);
00556 chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00557 infoc_1.ok);
00558 infoc_1.infot = 6;
00559 dtgsja_("N", "N", "N", &c__0, &c__0, &c_n1, &dummyk, &dummyl, a, &
00560 c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
00561 c__1, w, &ncycle, &info);
00562 chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00563 infoc_1.ok);
00564 infoc_1.infot = 10;
00565 dtgsja_("N", "N", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
00566 c__0, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
00567 c__1, w, &ncycle, &info);
00568 chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00569 infoc_1.ok);
00570 infoc_1.infot = 12;
00571 dtgsja_("N", "N", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
00572 c__1, b, &c__0, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
00573 c__1, w, &ncycle, &info);
00574 chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00575 infoc_1.ok);
00576 infoc_1.infot = 18;
00577 dtgsja_("U", "N", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
00578 c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__0, v, &c__1, q, &
00579 c__1, w, &ncycle, &info);
00580 chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00581 infoc_1.ok);
00582 infoc_1.infot = 20;
00583 dtgsja_("N", "V", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
00584 c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__0, q, &
00585 c__1, w, &ncycle, &info);
00586 chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00587 infoc_1.ok);
00588 infoc_1.infot = 22;
00589 dtgsja_("N", "N", "Q", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
00590 c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
00591 c__0, w, &ncycle, &info);
00592 chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00593 infoc_1.ok);
00594 nt += 11;
00595
00596
00597
00598 } else if (lsamen_(&c__3, path, "GLM")) {
00599
00600
00601
00602 s_copy(srnamc_1.srnamt, "DGGGLM", (ftnlen)32, (ftnlen)6);
00603 infoc_1.infot = 1;
00604 dggglm_(&c_n1, &c__0, &c__0, a, &c__1, b, &c__1, r1, r2, r3, w, &
00605 c__18, &info);
00606 chkxer_("DGGGLM", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00607 infoc_1.ok);
00608 infoc_1.infot = 2;
00609 dggglm_(&c__0, &c_n1, &c__0, a, &c__1, b, &c__1, r1, r2, r3, w, &
00610 c__18, &info);
00611 chkxer_("DGGGLM", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00612 infoc_1.ok);
00613 infoc_1.infot = 2;
00614 dggglm_(&c__0, &c__1, &c__0, a, &c__1, b, &c__1, r1, r2, r3, w, &
00615 c__18, &info);
00616 chkxer_("DGGGLM", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00617 infoc_1.ok);
00618 infoc_1.infot = 3;
00619 dggglm_(&c__0, &c__0, &c_n1, a, &c__1, b, &c__1, r1, r2, r3, w, &
00620 c__18, &info);
00621 chkxer_("DGGGLM", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00622 infoc_1.ok);
00623 infoc_1.infot = 3;
00624 dggglm_(&c__1, &c__0, &c__0, a, &c__1, b, &c__1, r1, r2, r3, w, &
00625 c__18, &info);
00626 chkxer_("DGGGLM", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00627 infoc_1.ok);
00628 infoc_1.infot = 5;
00629 dggglm_(&c__0, &c__0, &c__0, a, &c__0, b, &c__1, r1, r2, r3, w, &
00630 c__18, &info);
00631 chkxer_("DGGGLM", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00632 infoc_1.ok);
00633 infoc_1.infot = 7;
00634 dggglm_(&c__0, &c__0, &c__0, a, &c__1, b, &c__0, r1, r2, r3, w, &
00635 c__18, &info);
00636 chkxer_("DGGGLM", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00637 infoc_1.ok);
00638 infoc_1.infot = 12;
00639 dggglm_(&c__1, &c__1, &c__1, a, &c__1, b, &c__1, r1, r2, r3, w, &c__1,
00640 &info);
00641 chkxer_("DGGGLM", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00642 infoc_1.ok);
00643 nt += 8;
00644
00645
00646
00647 } else if (lsamen_(&c__3, path, "LSE")) {
00648
00649
00650
00651 s_copy(srnamc_1.srnamt, "DGGLSE", (ftnlen)32, (ftnlen)6);
00652 infoc_1.infot = 1;
00653 dgglse_(&c_n1, &c__0, &c__0, a, &c__1, b, &c__1, r1, r2, r3, w, &
00654 c__18, &info);
00655 chkxer_("DGGLSE", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00656 infoc_1.ok);
00657 infoc_1.infot = 2;
00658 dgglse_(&c__0, &c_n1, &c__0, a, &c__1, b, &c__1, r1, r2, r3, w, &
00659 c__18, &info);
00660 chkxer_("DGGLSE", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00661 infoc_1.ok);
00662 infoc_1.infot = 3;
00663 dgglse_(&c__0, &c__0, &c_n1, a, &c__1, b, &c__1, r1, r2, r3, w, &
00664 c__18, &info);
00665 chkxer_("DGGLSE", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00666 infoc_1.ok);
00667 infoc_1.infot = 3;
00668 dgglse_(&c__0, &c__0, &c__1, a, &c__1, b, &c__1, r1, r2, r3, w, &
00669 c__18, &info);
00670 chkxer_("DGGLSE", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00671 infoc_1.ok);
00672 infoc_1.infot = 3;
00673 dgglse_(&c__0, &c__1, &c__0, a, &c__1, b, &c__1, r1, r2, r3, w, &
00674 c__18, &info);
00675 chkxer_("DGGLSE", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00676 infoc_1.ok);
00677 infoc_1.infot = 5;
00678 dgglse_(&c__0, &c__0, &c__0, a, &c__0, b, &c__1, r1, r2, r3, w, &
00679 c__18, &info);
00680 chkxer_("DGGLSE", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00681 infoc_1.ok);
00682 infoc_1.infot = 7;
00683 dgglse_(&c__0, &c__0, &c__0, a, &c__1, b, &c__0, r1, r2, r3, w, &
00684 c__18, &info);
00685 chkxer_("DGGLSE", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00686 infoc_1.ok);
00687 infoc_1.infot = 12;
00688 dgglse_(&c__1, &c__1, &c__1, a, &c__1, b, &c__1, r1, r2, r3, w, &c__1,
00689 &info);
00690 chkxer_("DGGLSE", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00691 infoc_1.ok);
00692 nt += 8;
00693
00694
00695
00696 } else if (lsamen_(&c__3, path, "GQR")) {
00697
00698
00699
00700 s_copy(srnamc_1.srnamt, "DGGQRF", (ftnlen)32, (ftnlen)6);
00701 infoc_1.infot = 1;
00702 dggqrf_(&c_n1, &c__0, &c__0, a, &c__1, r1, b, &c__1, r2, w, &c__18, &
00703 info);
00704 chkxer_("DGGQRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00705 infoc_1.ok);
00706 infoc_1.infot = 2;
00707 dggqrf_(&c__0, &c_n1, &c__0, a, &c__1, r1, b, &c__1, r2, w, &c__18, &
00708 info);
00709 chkxer_("DGGQRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00710 infoc_1.ok);
00711 infoc_1.infot = 3;
00712 dggqrf_(&c__0, &c__0, &c_n1, a, &c__1, r1, b, &c__1, r2, w, &c__18, &
00713 info);
00714 chkxer_("DGGQRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00715 infoc_1.ok);
00716 infoc_1.infot = 5;
00717 dggqrf_(&c__0, &c__0, &c__0, a, &c__0, r1, b, &c__1, r2, w, &c__18, &
00718 info);
00719 chkxer_("DGGQRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00720 infoc_1.ok);
00721 infoc_1.infot = 8;
00722 dggqrf_(&c__0, &c__0, &c__0, a, &c__1, r1, b, &c__0, r2, w, &c__18, &
00723 info);
00724 chkxer_("DGGQRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00725 infoc_1.ok);
00726 infoc_1.infot = 11;
00727 dggqrf_(&c__1, &c__1, &c__2, a, &c__1, r1, b, &c__1, r2, w, &c__1, &
00728 info);
00729 chkxer_("DGGQRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00730 infoc_1.ok);
00731 nt += 6;
00732
00733
00734
00735 s_copy(srnamc_1.srnamt, "DGGRQF", (ftnlen)32, (ftnlen)6);
00736 infoc_1.infot = 1;
00737 dggrqf_(&c_n1, &c__0, &c__0, a, &c__1, r1, b, &c__1, r2, w, &c__18, &
00738 info);
00739 chkxer_("DGGRQF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00740 infoc_1.ok);
00741 infoc_1.infot = 2;
00742 dggrqf_(&c__0, &c_n1, &c__0, a, &c__1, r1, b, &c__1, r2, w, &c__18, &
00743 info);
00744 chkxer_("DGGRQF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00745 infoc_1.ok);
00746 infoc_1.infot = 3;
00747 dggrqf_(&c__0, &c__0, &c_n1, a, &c__1, r1, b, &c__1, r2, w, &c__18, &
00748 info);
00749 chkxer_("DGGRQF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00750 infoc_1.ok);
00751 infoc_1.infot = 5;
00752 dggrqf_(&c__0, &c__0, &c__0, a, &c__0, r1, b, &c__1, r2, w, &c__18, &
00753 info);
00754 chkxer_("DGGRQF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00755 infoc_1.ok);
00756 infoc_1.infot = 8;
00757 dggrqf_(&c__0, &c__0, &c__0, a, &c__1, r1, b, &c__0, r2, w, &c__18, &
00758 info);
00759 chkxer_("DGGRQF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00760 infoc_1.ok);
00761 infoc_1.infot = 11;
00762 dggrqf_(&c__1, &c__1, &c__2, a, &c__1, r1, b, &c__1, r2, w, &c__1, &
00763 info);
00764 chkxer_("DGGRQF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00765 infoc_1.ok);
00766 nt += 6;
00767
00768
00769
00770 } else if (lsamen_(&c__3, path, "DGS") || lsamen_(&
00771 c__3, path, "DGV") || lsamen_(&c__3, path,
00772 "DGX") || lsamen_(&c__3, path, "DXV")) {
00773
00774
00775
00776 s_copy(srnamc_1.srnamt, "DGGES ", (ftnlen)32, (ftnlen)6);
00777 infoc_1.infot = 1;
00778 dgges_("/", "N", "S", (L_fp)dlctes_, &c__1, a, &c__1, b, &c__1, &sdim,
00779 r1, r2, r3, q, &c__1, u, &c__1, w, &c__1, bw, &info);
00780 chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00781 infoc_1.ok);
00782 infoc_1.infot = 2;
00783 dgges_("N", "/", "S", (L_fp)dlctes_, &c__1, a, &c__1, b, &c__1, &sdim,
00784 r1, r2, r3, q, &c__1, u, &c__1, w, &c__1, bw, &info);
00785 chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00786 infoc_1.ok);
00787 infoc_1.infot = 3;
00788 dgges_("N", "V", "/", (L_fp)dlctes_, &c__1, a, &c__1, b, &c__1, &sdim,
00789 r1, r2, r3, q, &c__1, u, &c__1, w, &c__1, bw, &info);
00790 chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00791 infoc_1.ok);
00792 infoc_1.infot = 5;
00793 dgges_("N", "V", "S", (L_fp)dlctes_, &c_n1, a, &c__1, b, &c__1, &sdim,
00794 r1, r2, r3, q, &c__1, u, &c__1, w, &c__1, bw, &info);
00795 chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00796 infoc_1.ok);
00797 infoc_1.infot = 7;
00798 dgges_("N", "V", "S", (L_fp)dlctes_, &c__1, a, &c__0, b, &c__1, &sdim,
00799 r1, r2, r3, q, &c__1, u, &c__1, w, &c__1, bw, &info);
00800 chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00801 infoc_1.ok);
00802 infoc_1.infot = 9;
00803 dgges_("N", "V", "S", (L_fp)dlctes_, &c__1, a, &c__1, b, &c__0, &sdim,
00804 r1, r2, r3, q, &c__1, u, &c__1, w, &c__1, bw, &info);
00805 chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00806 infoc_1.ok);
00807 infoc_1.infot = 15;
00808 dgges_("N", "V", "S", (L_fp)dlctes_, &c__1, a, &c__1, b, &c__1, &sdim,
00809 r1, r2, r3, q, &c__0, u, &c__1, w, &c__1, bw, &info);
00810 chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00811 infoc_1.ok);
00812 infoc_1.infot = 15;
00813 dgges_("V", "V", "S", (L_fp)dlctes_, &c__2, a, &c__2, b, &c__2, &sdim,
00814 r1, r2, r3, q, &c__1, u, &c__2, w, &c__1, bw, &info);
00815 chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00816 infoc_1.ok);
00817 infoc_1.infot = 17;
00818 dgges_("N", "V", "S", (L_fp)dlctes_, &c__1, a, &c__1, b, &c__1, &sdim,
00819 r1, r2, r3, q, &c__1, u, &c__0, w, &c__1, bw, &info);
00820 chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00821 infoc_1.ok);
00822 infoc_1.infot = 17;
00823 dgges_("V", "V", "S", (L_fp)dlctes_, &c__2, a, &c__2, b, &c__2, &sdim,
00824 r1, r2, r3, q, &c__2, u, &c__1, w, &c__1, bw, &info);
00825 chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00826 infoc_1.ok);
00827 infoc_1.infot = 19;
00828 dgges_("V", "V", "S", (L_fp)dlctes_, &c__2, a, &c__2, b, &c__2, &sdim,
00829 r1, r2, r3, q, &c__2, u, &c__2, w, &c__1, bw, &info);
00830 chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00831 infoc_1.ok);
00832 nt += 11;
00833
00834
00835
00836 s_copy(srnamc_1.srnamt, "DGGESX", (ftnlen)32, (ftnlen)6);
00837 infoc_1.infot = 1;
00838 dggesx_("/", "N", "S", (L_fp)dlctsx_, "N", &c__1, a, &c__1, b, &c__1,
00839 &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
00840 &c__1, bw, &info)
00841 ;
00842 chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00843 infoc_1.ok);
00844 infoc_1.infot = 2;
00845 dggesx_("N", "/", "S", (L_fp)dlctsx_, "N", &c__1, a, &c__1, b, &c__1,
00846 &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
00847 &c__1, bw, &info)
00848 ;
00849 chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00850 infoc_1.ok);
00851 infoc_1.infot = 3;
00852 dggesx_("V", "V", "/", (L_fp)dlctsx_, "N", &c__1, a, &c__1, b, &c__1,
00853 &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
00854 &c__1, bw, &info)
00855 ;
00856 chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00857 infoc_1.ok);
00858 infoc_1.infot = 5;
00859 dggesx_("V", "V", "S", (L_fp)dlctsx_, "/", &c__1, a, &c__1, b, &c__1,
00860 &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
00861 &c__1, bw, &info)
00862 ;
00863 chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00864 infoc_1.ok);
00865 infoc_1.infot = 6;
00866 dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c_n1, a, &c__1, b, &c__1,
00867 &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
00868 &c__1, bw, &info)
00869 ;
00870 chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00871 infoc_1.ok);
00872 infoc_1.infot = 8;
00873 dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__1, a, &c__0, b, &c__1,
00874 &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
00875 &c__1, bw, &info)
00876 ;
00877 chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00878 infoc_1.ok);
00879 infoc_1.infot = 10;
00880 dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__1, a, &c__1, b, &c__0,
00881 &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
00882 &c__1, bw, &info)
00883 ;
00884 chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00885 infoc_1.ok);
00886 infoc_1.infot = 16;
00887 dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__1, a, &c__1, b, &c__1,
00888 &sdim, r1, r2, r3, q, &c__0, u, &c__1, rce, rcv, w, &c__1, iw,
00889 &c__1, bw, &info)
00890 ;
00891 chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00892 infoc_1.ok);
00893 infoc_1.infot = 16;
00894 dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__2, a, &c__2, b, &c__2,
00895 &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
00896 &c__1, bw, &info)
00897 ;
00898 chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00899 infoc_1.ok);
00900 infoc_1.infot = 18;
00901 dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__1, a, &c__1, b, &c__1,
00902 &sdim, r1, r2, r3, q, &c__1, u, &c__0, rce, rcv, w, &c__1, iw,
00903 &c__1, bw, &info)
00904 ;
00905 chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00906 infoc_1.ok);
00907 infoc_1.infot = 18;
00908 dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__2, a, &c__2, b, &c__2,
00909 &sdim, r1, r2, r3, q, &c__2, u, &c__1, rce, rcv, w, &c__1, iw,
00910 &c__1, bw, &info)
00911 ;
00912 chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00913 infoc_1.ok);
00914 infoc_1.infot = 22;
00915 dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__2, a, &c__2, b, &c__2,
00916 &sdim, r1, r2, r3, q, &c__2, u, &c__2, rce, rcv, w, &c__1, iw,
00917 &c__1, bw, &info)
00918 ;
00919 chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00920 infoc_1.ok);
00921 infoc_1.infot = 24;
00922 dggesx_("V", "V", "S", (L_fp)dlctsx_, "V", &c__1, a, &c__1, b, &c__1,
00923 &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__32,
00924 iw, &c__0, bw, &info);
00925 chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00926 infoc_1.ok);
00927 nt += 13;
00928
00929
00930
00931 s_copy(srnamc_1.srnamt, "DGGEV ", (ftnlen)32, (ftnlen)6);
00932 infoc_1.infot = 1;
00933 dggev_("/", "N", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q, &c__1, u, &
00934 c__1, w, &c__1, &info);
00935 chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00936 infoc_1.ok);
00937 infoc_1.infot = 2;
00938 dggev_("N", "/", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q, &c__1, u, &
00939 c__1, w, &c__1, &info);
00940 chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00941 infoc_1.ok);
00942 infoc_1.infot = 3;
00943 dggev_("V", "V", &c_n1, a, &c__1, b, &c__1, r1, r2, r3, q, &c__1, u, &
00944 c__1, w, &c__1, &info);
00945 chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00946 infoc_1.ok);
00947 infoc_1.infot = 5;
00948 dggev_("V", "V", &c__1, a, &c__0, b, &c__1, r1, r2, r3, q, &c__1, u, &
00949 c__1, w, &c__1, &info);
00950 chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00951 infoc_1.ok);
00952 infoc_1.infot = 7;
00953 dggev_("V", "V", &c__1, a, &c__1, b, &c__0, r1, r2, r3, q, &c__1, u, &
00954 c__1, w, &c__1, &info);
00955 chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00956 infoc_1.ok);
00957 infoc_1.infot = 12;
00958 dggev_("N", "V", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q, &c__0, u, &
00959 c__1, w, &c__1, &info);
00960 chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00961 infoc_1.ok);
00962 infoc_1.infot = 12;
00963 dggev_("V", "V", &c__2, a, &c__2, b, &c__2, r1, r2, r3, q, &c__1, u, &
00964 c__2, w, &c__1, &info);
00965 chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00966 infoc_1.ok);
00967 infoc_1.infot = 14;
00968 dggev_("V", "N", &c__2, a, &c__2, b, &c__2, r1, r2, r3, q, &c__2, u, &
00969 c__0, w, &c__1, &info);
00970 chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00971 infoc_1.ok);
00972 infoc_1.infot = 14;
00973 dggev_("V", "V", &c__2, a, &c__2, b, &c__2, r1, r2, r3, q, &c__2, u, &
00974 c__1, w, &c__1, &info);
00975 chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00976 infoc_1.ok);
00977 infoc_1.infot = 16;
00978 dggev_("V", "V", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q, &c__1, u, &
00979 c__1, w, &c__1, &info);
00980 chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00981 infoc_1.ok);
00982 nt += 10;
00983
00984
00985
00986 s_copy(srnamc_1.srnamt, "DGGEVX", (ftnlen)32, (ftnlen)6);
00987 infoc_1.infot = 1;
00988 dggevx_("/", "N", "N", "N", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q,
00989 &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
00990 w, &c__1, iw, bw, &info);
00991 chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00992 infoc_1.ok);
00993 infoc_1.infot = 2;
00994 dggevx_("N", "/", "N", "N", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q,
00995 &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
00996 w, &c__1, iw, bw, &info);
00997 chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
00998 infoc_1.ok);
00999 infoc_1.infot = 3;
01000 dggevx_("N", "N", "/", "N", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q,
01001 &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
01002 w, &c__1, iw, bw, &info);
01003 chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01004 infoc_1.ok);
01005 infoc_1.infot = 4;
01006 dggevx_("N", "N", "N", "/", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q,
01007 &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
01008 w, &c__1, iw, bw, &info);
01009 chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01010 infoc_1.ok);
01011 infoc_1.infot = 5;
01012 dggevx_("N", "N", "N", "N", &c_n1, a, &c__1, b, &c__1, r1, r2, r3, q,
01013 &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
01014 w, &c__1, iw, bw, &info);
01015 chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01016 infoc_1.ok);
01017 infoc_1.infot = 7;
01018 dggevx_("N", "N", "N", "N", &c__1, a, &c__0, b, &c__1, r1, r2, r3, q,
01019 &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
01020 w, &c__1, iw, bw, &info);
01021 chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01022 infoc_1.ok);
01023 infoc_1.infot = 9;
01024 dggevx_("N", "N", "N", "N", &c__1, a, &c__1, b, &c__0, r1, r2, r3, q,
01025 &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
01026 w, &c__1, iw, bw, &info);
01027 chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01028 infoc_1.ok);
01029 infoc_1.infot = 14;
01030 dggevx_("N", "N", "N", "N", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q,
01031 &c__0, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
01032 w, &c__1, iw, bw, &info);
01033 chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01034 infoc_1.ok);
01035 infoc_1.infot = 14;
01036 dggevx_("N", "V", "N", "N", &c__2, a, &c__2, b, &c__2, r1, r2, r3, q,
01037 &c__1, u, &c__2, &c__1, &c__2, ls, rs, &anrm, &bnrm, rce, rcv,
01038 w, &c__1, iw, bw, &info);
01039 chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01040 infoc_1.ok);
01041 infoc_1.infot = 16;
01042 dggevx_("N", "N", "N", "N", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q,
01043 &c__1, u, &c__0, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
01044 w, &c__1, iw, bw, &info);
01045 chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01046 infoc_1.ok);
01047 infoc_1.infot = 16;
01048 dggevx_("N", "N", "V", "N", &c__2, a, &c__2, b, &c__2, r1, r2, r3, q,
01049 &c__2, u, &c__1, &c__1, &c__2, ls, rs, &anrm, &bnrm, rce, rcv,
01050 w, &c__1, iw, bw, &info);
01051 chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01052 infoc_1.ok);
01053 infoc_1.infot = 26;
01054 dggevx_("N", "N", "V", "N", &c__2, a, &c__2, b, &c__2, r1, r2, r3, q,
01055 &c__2, u, &c__2, &c__1, &c__2, ls, rs, &anrm, &bnrm, rce, rcv,
01056 w, &c__1, iw, bw, &info);
01057 chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01058 infoc_1.ok);
01059 nt += 12;
01060
01061
01062
01063 s_copy(srnamc_1.srnamt, "DTGEXC", (ftnlen)32, (ftnlen)6);
01064 infoc_1.infot = 3;
01065 dtgexc_(&c_true, &c_true, &c_n1, a, &c__1, b, &c__1, q, &c__1, z__, &
01066 c__1, &ifst, &ilst, w, &c__1, &info);
01067 chkxer_("DTGEXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01068 infoc_1.ok);
01069 infoc_1.infot = 5;
01070 dtgexc_(&c_true, &c_true, &c__1, a, &c__0, b, &c__1, q, &c__1, z__, &
01071 c__1, &ifst, &ilst, w, &c__1, &info);
01072 chkxer_("DTGEXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01073 infoc_1.ok);
01074 infoc_1.infot = 7;
01075 dtgexc_(&c_true, &c_true, &c__1, a, &c__1, b, &c__0, q, &c__1, z__, &
01076 c__1, &ifst, &ilst, w, &c__1, &info);
01077 chkxer_("DTGEXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01078 infoc_1.ok);
01079 infoc_1.infot = 9;
01080 dtgexc_(&c_false, &c_true, &c__1, a, &c__1, b, &c__1, q, &c__0, z__, &
01081 c__1, &ifst, &ilst, w, &c__1, &info);
01082 chkxer_("DTGEXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01083 infoc_1.ok);
01084 infoc_1.infot = 9;
01085 dtgexc_(&c_true, &c_true, &c__1, a, &c__1, b, &c__1, q, &c__0, z__, &
01086 c__1, &ifst, &ilst, w, &c__1, &info);
01087 chkxer_("DTGEXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01088 infoc_1.ok);
01089 infoc_1.infot = 11;
01090 dtgexc_(&c_true, &c_false, &c__1, a, &c__1, b, &c__1, q, &c__1, z__, &
01091 c__0, &ifst, &ilst, w, &c__1, &info);
01092 chkxer_("DTGEXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01093 infoc_1.ok);
01094 infoc_1.infot = 11;
01095 dtgexc_(&c_true, &c_true, &c__1, a, &c__1, b, &c__1, q, &c__1, z__, &
01096 c__0, &ifst, &ilst, w, &c__1, &info);
01097 chkxer_("DTGEXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01098 infoc_1.ok);
01099 infoc_1.infot = 15;
01100 dtgexc_(&c_true, &c_true, &c__1, a, &c__1, b, &c__1, q, &c__1, z__, &
01101 c__1, &ifst, &ilst, w, &c__0, &info);
01102 chkxer_("DTGEXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01103 infoc_1.ok);
01104 nt += 8;
01105
01106
01107
01108 s_copy(srnamc_1.srnamt, "DTGSEN", (ftnlen)32, (ftnlen)6);
01109 infoc_1.infot = 1;
01110 dtgsen_(&c_n1, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
01111 r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
01112 iw, &c__1, &info);
01113 chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01114 infoc_1.ok);
01115 infoc_1.infot = 5;
01116 dtgsen_(&c__1, &c_true, &c_true, sel, &c_n1, a, &c__1, b, &c__1, r1,
01117 r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
01118 iw, &c__1, &info);
01119 chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01120 infoc_1.ok);
01121 infoc_1.infot = 7;
01122 dtgsen_(&c__1, &c_true, &c_true, sel, &c__1, a, &c__0, b, &c__1, r1,
01123 r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
01124 iw, &c__1, &info);
01125 chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01126 infoc_1.ok);
01127 infoc_1.infot = 9;
01128 dtgsen_(&c__1, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__0, r1,
01129 r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
01130 iw, &c__1, &info);
01131 chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01132 infoc_1.ok);
01133 infoc_1.infot = 14;
01134 dtgsen_(&c__1, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
01135 r2, r3, q, &c__0, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
01136 iw, &c__1, &info);
01137 chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01138 infoc_1.ok);
01139 infoc_1.infot = 16;
01140 dtgsen_(&c__1, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
01141 r2, r3, q, &c__1, z__, &c__0, &m, &tola, &tolb, rcv, w, &c__1,
01142 iw, &c__1, &info);
01143 chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01144 infoc_1.ok);
01145 infoc_1.infot = 22;
01146 dtgsen_(&c__0, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
01147 r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
01148 iw, &c__1, &info);
01149 chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01150 infoc_1.ok);
01151 infoc_1.infot = 22;
01152 dtgsen_(&c__1, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
01153 r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
01154 iw, &c__1, &info);
01155 chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01156 infoc_1.ok);
01157 infoc_1.infot = 22;
01158 dtgsen_(&c__2, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
01159 r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
01160 iw, &c__1, &info);
01161 chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01162 infoc_1.ok);
01163 infoc_1.infot = 24;
01164 dtgsen_(&c__0, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
01165 r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &
01166 c__20, iw, &c__0, &info);
01167 chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01168 infoc_1.ok);
01169 infoc_1.infot = 24;
01170 dtgsen_(&c__1, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
01171 r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &
01172 c__20, iw, &c__0, &info);
01173 chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01174 infoc_1.ok);
01175 infoc_1.infot = 24;
01176 dtgsen_(&c__2, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
01177 r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &
01178 c__20, iw, &c__1, &info);
01179 chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01180 infoc_1.ok);
01181 nt += 12;
01182
01183
01184
01185 s_copy(srnamc_1.srnamt, "DTGSNA", (ftnlen)32, (ftnlen)6);
01186 infoc_1.infot = 1;
01187 dtgsna_("/", "A", sel, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &c__1,
01188 r1, r2, &c__1, &m, w, &c__1, iw, &info);
01189 chkxer_("DTGSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01190 infoc_1.ok);
01191 infoc_1.infot = 2;
01192 dtgsna_("B", "/", sel, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &c__1,
01193 r1, r2, &c__1, &m, w, &c__1, iw, &info);
01194 chkxer_("DTGSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01195 infoc_1.ok);
01196 infoc_1.infot = 4;
01197 dtgsna_("B", "A", sel, &c_n1, a, &c__1, b, &c__1, q, &c__1, u, &c__1,
01198 r1, r2, &c__1, &m, w, &c__1, iw, &info);
01199 chkxer_("DTGSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01200 infoc_1.ok);
01201 infoc_1.infot = 6;
01202 dtgsna_("B", "A", sel, &c__1, a, &c__0, b, &c__1, q, &c__1, u, &c__1,
01203 r1, r2, &c__1, &m, w, &c__1, iw, &info);
01204 chkxer_("DTGSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01205 infoc_1.ok);
01206 infoc_1.infot = 8;
01207 dtgsna_("B", "A", sel, &c__1, a, &c__1, b, &c__0, q, &c__1, u, &c__1,
01208 r1, r2, &c__1, &m, w, &c__1, iw, &info);
01209 chkxer_("DTGSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01210 infoc_1.ok);
01211 infoc_1.infot = 10;
01212 dtgsna_("E", "A", sel, &c__1, a, &c__1, b, &c__1, q, &c__0, u, &c__1,
01213 r1, r2, &c__1, &m, w, &c__1, iw, &info);
01214 chkxer_("DTGSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01215 infoc_1.ok);
01216 infoc_1.infot = 12;
01217 dtgsna_("E", "A", sel, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &c__0,
01218 r1, r2, &c__1, &m, w, &c__1, iw, &info);
01219 chkxer_("DTGSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01220 infoc_1.ok);
01221 infoc_1.infot = 15;
01222 dtgsna_("E", "A", sel, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &c__1,
01223 r1, r2, &c__0, &m, w, &c__1, iw, &info);
01224 chkxer_("DTGSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01225 infoc_1.ok);
01226 infoc_1.infot = 18;
01227 dtgsna_("E", "A", sel, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &c__1,
01228 r1, r2, &c__1, &m, w, &c__0, iw, &info);
01229 chkxer_("DTGSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01230 infoc_1.ok);
01231 nt += 9;
01232
01233
01234
01235 s_copy(srnamc_1.srnamt, "DTGSYL", (ftnlen)32, (ftnlen)6);
01236 infoc_1.infot = 1;
01237 dtgsyl_("/", &c__0, &c__1, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &
01238 c__1, v, &c__1, z__, &c__1, &scale, &dif, w, &c__1, iw, &info);
01239 chkxer_("DTGSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01240 infoc_1.ok);
01241 infoc_1.infot = 2;
01242 dtgsyl_("N", &c_n1, &c__1, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &
01243 c__1, v, &c__1, z__, &c__1, &scale, &dif, w, &c__1, iw, &info);
01244 chkxer_("DTGSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01245 infoc_1.ok);
01246 infoc_1.infot = 3;
01247 dtgsyl_("N", &c__0, &c__0, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &
01248 c__1, v, &c__1, z__, &c__1, &scale, &dif, w, &c__1, iw, &info);
01249 chkxer_("DTGSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01250 infoc_1.ok);
01251 infoc_1.infot = 4;
01252 dtgsyl_("N", &c__0, &c__1, &c__0, a, &c__1, b, &c__1, q, &c__1, u, &
01253 c__1, v, &c__1, z__, &c__1, &scale, &dif, w, &c__1, iw, &info);
01254 chkxer_("DTGSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01255 infoc_1.ok);
01256 infoc_1.infot = 6;
01257 dtgsyl_("N", &c__0, &c__1, &c__1, a, &c__0, b, &c__1, q, &c__1, u, &
01258 c__1, v, &c__1, z__, &c__1, &scale, &dif, w, &c__1, iw, &info);
01259 chkxer_("DTGSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01260 infoc_1.ok);
01261 infoc_1.infot = 8;
01262 dtgsyl_("N", &c__0, &c__1, &c__1, a, &c__1, b, &c__0, q, &c__1, u, &
01263 c__1, v, &c__1, z__, &c__1, &scale, &dif, w, &c__1, iw, &info);
01264 chkxer_("DTGSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01265 infoc_1.ok);
01266 infoc_1.infot = 10;
01267 dtgsyl_("N", &c__0, &c__1, &c__1, a, &c__1, b, &c__1, q, &c__0, u, &
01268 c__1, v, &c__1, z__, &c__1, &scale, &dif, w, &c__1, iw, &info);
01269 chkxer_("DTGSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01270 infoc_1.ok);
01271 infoc_1.infot = 12;
01272 dtgsyl_("N", &c__0, &c__1, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &
01273 c__0, v, &c__1, z__, &c__1, &scale, &dif, w, &c__1, iw, &info);
01274 chkxer_("DTGSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01275 infoc_1.ok);
01276 infoc_1.infot = 14;
01277 dtgsyl_("N", &c__0, &c__1, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &
01278 c__1, v, &c__0, z__, &c__1, &scale, &dif, w, &c__1, iw, &info);
01279 chkxer_("DTGSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01280 infoc_1.ok);
01281 infoc_1.infot = 16;
01282 dtgsyl_("N", &c__0, &c__1, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &
01283 c__1, v, &c__1, z__, &c__0, &scale, &dif, w, &c__1, iw, &info);
01284 chkxer_("DTGSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01285 infoc_1.ok);
01286 infoc_1.infot = 20;
01287 dtgsyl_("N", &c__1, &c__1, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &
01288 c__1, v, &c__1, z__, &c__1, &scale, &dif, w, &c__1, iw, &info);
01289 chkxer_("DTGSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01290 infoc_1.ok);
01291 infoc_1.infot = 20;
01292 dtgsyl_("N", &c__2, &c__1, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &
01293 c__1, v, &c__1, z__, &c__1, &scale, &dif, w, &c__1, iw, &info);
01294 chkxer_("DTGSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
01295 infoc_1.ok);
01296 nt += 12;
01297 }
01298
01299
01300
01301 if (infoc_1.ok) {
01302 io___38.ciunit = infoc_1.nout;
01303 s_wsfe(&io___38);
01304 do_fio(&c__1, path, (ftnlen)3);
01305 do_fio(&c__1, (char *)&nt, (ftnlen)sizeof(integer));
01306 e_wsfe();
01307 } else {
01308 io___39.ciunit = infoc_1.nout;
01309 s_wsfe(&io___39);
01310 do_fio(&c__1, path, (ftnlen)3);
01311 e_wsfe();
01312 }
01313
01314
01315 return 0;
01316
01317
01318
01319 }