/nvxs/CLAPACK/TESTING/EIG/derrgg.c

https://github.com/nagadomi/animeface-2009 · C · 1307 lines · 1160 code · 79 blank · 68 comment · 19 complexity · e1599ac521d9b2895ae522fc699deebc MD5 · raw file

Large files are truncated click here to view the full file

  1. #include "f2c.h"
  2. #include "blaswrap.h"
  3. /* Common Block Declarations */
  4. struct {
  5. integer infot, nout;
  6. logical ok, lerr;
  7. } infoc_;
  8. #define infoc_1 infoc_
  9. struct {
  10. char srnamt[6];
  11. } srnamc_;
  12. #define srnamc_1 srnamc_
  13. /* Table of constant values */
  14. static integer c__2 = 2;
  15. static integer c__0 = 0;
  16. static integer c__1 = 1;
  17. static integer c_n1 = -1;
  18. static integer c__18 = 18;
  19. static integer c__3 = 3;
  20. static integer c__32 = 32;
  21. static logical c_true = TRUE_;
  22. static logical c_false = FALSE_;
  23. static integer c__20 = 20;
  24. /* Subroutine */ int derrgg_(char *path, integer *nunit)
  25. {
  26. /* Format strings */
  27. static char fmt_9999[] = "(1x,a3,\002 routines passed the tests of the e"
  28. "rror exits (\002,i3,\002 tests done)\002)";
  29. static char fmt_9998[] = "(\002 *** \002,a3,\002 routines failed the tes"
  30. "ts of the error \002,\002exits ***\002)";
  31. /* Builtin functions */
  32. integer s_wsle(cilist *), e_wsle(void);
  33. /* Subroutine */ int s_copy(char *, char *, ftnlen, ftnlen);
  34. integer s_wsfe(cilist *), do_fio(integer *, char *, ftnlen), e_wsfe(void);
  35. /* Local variables */
  36. doublereal a[9] /* was [3][3] */, b[9] /* was [3][3] */;
  37. integer i__, j, m;
  38. doublereal q[9] /* was [3][3] */, u[9] /* was [3][3] */, v[9] /*
  39. was [3][3] */, w[18], z__[9] /* was [3][3] */;
  40. char c2[2];
  41. doublereal r1[3], r2[3], r3[3];
  42. logical bw[3];
  43. doublereal ls[3];
  44. integer iw[3], nt;
  45. doublereal rs[3], dif, rce[2];
  46. logical sel[3];
  47. doublereal tau[3], rcv[2];
  48. integer info, sdim;
  49. doublereal anrm, bnrm, tola, tolb;
  50. integer ifst, ilst;
  51. doublereal scale;
  52. extern /* Subroutine */ int dgges_(char *, char *, char *, L_fp, integer *
  53. , doublereal *, integer *, doublereal *, integer *, integer *,
  54. doublereal *, doublereal *, doublereal *, doublereal *, integer *,
  55. doublereal *, integer *, doublereal *, integer *, logical *,
  56. integer *), dggev_(char *, char *,
  57. integer *, doublereal *, integer *, doublereal *, integer *,
  58. doublereal *, doublereal *, doublereal *, doublereal *, integer *,
  59. doublereal *, integer *, doublereal *, integer *, integer *), dgghrd_(char *, char *, integer *, integer *,
  60. integer *, doublereal *, integer *, doublereal *, integer *,
  61. doublereal *, integer *, doublereal *, integer *, integer *), dggglm_(integer *, integer *, integer *,
  62. doublereal *, integer *, doublereal *, integer *, doublereal *,
  63. doublereal *, doublereal *, doublereal *, integer *, integer *),
  64. dgglse_(integer *, integer *, integer *, doublereal *, integer *,
  65. doublereal *, integer *, doublereal *, doublereal *, doublereal *,
  66. doublereal *, integer *, integer *), dggqrf_(integer *, integer *
  67. , integer *, doublereal *, integer *, doublereal *, doublereal *,
  68. integer *, doublereal *, doublereal *, integer *, integer *),
  69. dggrqf_(integer *, integer *, integer *, doublereal *, integer *,
  70. doublereal *, doublereal *, integer *, doublereal *, doublereal *,
  71. integer *, integer *);
  72. integer ncycle;
  73. extern logical dlctes_(), lsamen_(integer *, char *, char *);
  74. extern /* Subroutine */ int dggsvd_(char *, char *, char *, integer *,
  75. integer *, integer *, integer *, integer *, doublereal *, integer
  76. *, doublereal *, integer *, doublereal *, doublereal *,
  77. doublereal *, integer *, doublereal *, integer *, doublereal *,
  78. integer *, doublereal *, integer *, integer *), dggesx_(char *, char *, char *, L_fp, char *, integer *,
  79. doublereal *, integer *, doublereal *, integer *, integer *,
  80. doublereal *, doublereal *, doublereal *, doublereal *, integer *,
  81. doublereal *, integer *, doublereal *, doublereal *, doublereal *
  82. , integer *, integer *, integer *, logical *, integer *), dhgeqz_(char *, char *, char *, integer *
  83. , integer *, integer *, doublereal *, integer *, doublereal *,
  84. integer *, doublereal *, doublereal *, doublereal *, doublereal *,
  85. integer *, doublereal *, integer *, doublereal *, integer *,
  86. integer *), dtgevc_(char *, char *,
  87. logical *, integer *, doublereal *, integer *, doublereal *,
  88. integer *, doublereal *, integer *, doublereal *, integer *,
  89. integer *, integer *, doublereal *, integer *),
  90. chkxer_(char *, integer *, integer *, logical *, logical *), dggevx_(char *, char *, char *, char *, integer *,
  91. doublereal *, integer *, doublereal *, integer *, doublereal *,
  92. doublereal *, doublereal *, doublereal *, integer *, doublereal *,
  93. integer *, integer *, integer *, doublereal *, doublereal *,
  94. doublereal *, doublereal *, doublereal *, doublereal *,
  95. doublereal *, integer *, integer *, logical *, integer *), dtgexc_(logical *, logical *, integer *,
  96. doublereal *, integer *, doublereal *, integer *, doublereal *,
  97. integer *, doublereal *, integer *, integer *, integer *,
  98. doublereal *, integer *, integer *), dtgsen_(integer *, logical *,
  99. logical *, logical *, integer *, doublereal *, integer *,
  100. doublereal *, integer *, doublereal *, doublereal *, doublereal *,
  101. doublereal *, integer *, doublereal *, integer *, integer *,
  102. doublereal *, doublereal *, doublereal *, doublereal *, integer *,
  103. integer *, integer *, integer *), dtgsja_(char *, char *, char *,
  104. integer *, integer *, integer *, integer *, integer *,
  105. doublereal *, integer *, doublereal *, integer *, doublereal *,
  106. doublereal *, doublereal *, doublereal *, doublereal *, integer *,
  107. doublereal *, integer *, doublereal *, integer *, doublereal *,
  108. integer *, integer *), dtgsna_(char *,
  109. char *, logical *, integer *, doublereal *, integer *, doublereal
  110. *, integer *, doublereal *, integer *, doublereal *, integer *,
  111. doublereal *, doublereal *, integer *, integer *, doublereal *,
  112. integer *, integer *, integer *), dggsvp_(char *,
  113. char *, char *, integer *, integer *, integer *, doublereal *,
  114. integer *, doublereal *, integer *, doublereal *, doublereal *,
  115. integer *, integer *, doublereal *, integer *, doublereal *,
  116. integer *, doublereal *, integer *, integer *, doublereal *,
  117. doublereal *, integer *);
  118. extern logical dlctsx_();
  119. integer dummyk, dummyl;
  120. extern /* Subroutine */ int dtgsyl_(char *, integer *, integer *, integer
  121. *, doublereal *, integer *, doublereal *, integer *, doublereal *,
  122. integer *, doublereal *, integer *, doublereal *, integer *,
  123. doublereal *, integer *, doublereal *, doublereal *, doublereal *,
  124. integer *, integer *, integer *);
  125. /* Fortran I/O blocks */
  126. static cilist io___1 = { 0, 0, 0, 0, 0 };
  127. static cilist io___38 = { 0, 0, 0, fmt_9999, 0 };
  128. static cilist io___39 = { 0, 0, 0, fmt_9998, 0 };
  129. /* -- LAPACK test routine (version 3.1) -- */
  130. /* Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. */
  131. /* November 2006 */
  132. /* .. Scalar Arguments .. */
  133. /* .. */
  134. /* Purpose */
  135. /* ======= */
  136. /* DERRGG tests the error exits for DGGES, DGGESX, DGGEV, DGGEVX, */
  137. /* DGGGLM, DGGHRD, DGGLSE, DGGQRF, DGGRQF, DGGSVD, DGGSVP, DHGEQZ, */
  138. /* DTGEVC, DTGEXC, DTGSEN, DTGSJA, DTGSNA, and DTGSYL. */
  139. /* Arguments */
  140. /* ========= */
  141. /* PATH (input) CHARACTER*3 */
  142. /* The LAPACK path name for the routines to be tested. */
  143. /* NUNIT (input) INTEGER */
  144. /* The unit number for output. */
  145. /* ===================================================================== */
  146. /* .. Parameters .. */
  147. /* .. */
  148. /* .. Local Scalars .. */
  149. /* .. */
  150. /* .. Local Arrays .. */
  151. /* .. */
  152. /* .. External Functions .. */
  153. /* .. */
  154. /* .. External Subroutines .. */
  155. /* .. */
  156. /* .. Scalars in Common .. */
  157. /* .. */
  158. /* .. Common blocks .. */
  159. /* .. */
  160. /* .. Executable Statements .. */
  161. infoc_1.nout = *nunit;
  162. io___1.ciunit = infoc_1.nout;
  163. s_wsle(&io___1);
  164. e_wsle();
  165. s_copy(c2, path + 1, (ftnlen)2, (ftnlen)2);
  166. /* Set the variables to innocuous values. */
  167. for (j = 1; j <= 3; ++j) {
  168. sel[j - 1] = TRUE_;
  169. for (i__ = 1; i__ <= 3; ++i__) {
  170. a[i__ + j * 3 - 4] = 0.;
  171. b[i__ + j * 3 - 4] = 0.;
  172. /* L10: */
  173. }
  174. /* L20: */
  175. }
  176. for (i__ = 1; i__ <= 3; ++i__) {
  177. a[i__ + i__ * 3 - 4] = 1.;
  178. b[i__ + i__ * 3 - 4] = 1.;
  179. /* L30: */
  180. }
  181. infoc_1.ok = TRUE_;
  182. tola = 1.;
  183. tolb = 1.;
  184. ifst = 1;
  185. ilst = 1;
  186. nt = 0;
  187. /* Test error exits for the GG path. */
  188. if (lsamen_(&c__2, c2, "GG")) {
  189. /* DGGHRD */
  190. s_copy(srnamc_1.srnamt, "DGGHRD", (ftnlen)6, (ftnlen)6);
  191. infoc_1.infot = 1;
  192. dgghrd_("/", "N", &c__0, &c__1, &c__0, a, &c__1, b, &c__1, q, &c__1,
  193. z__, &c__1, &info);
  194. chkxer_("DGGHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  195. infoc_1.ok);
  196. infoc_1.infot = 2;
  197. dgghrd_("N", "/", &c__0, &c__1, &c__0, a, &c__1, b, &c__1, q, &c__1,
  198. z__, &c__1, &info);
  199. chkxer_("DGGHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  200. infoc_1.ok);
  201. infoc_1.infot = 3;
  202. dgghrd_("N", "N", &c_n1, &c__0, &c__0, a, &c__1, b, &c__1, q, &c__1,
  203. z__, &c__1, &info);
  204. chkxer_("DGGHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  205. infoc_1.ok);
  206. infoc_1.infot = 4;
  207. dgghrd_("N", "N", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, q, &c__1,
  208. z__, &c__1, &info);
  209. chkxer_("DGGHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  210. infoc_1.ok);
  211. infoc_1.infot = 5;
  212. dgghrd_("N", "N", &c__0, &c__1, &c__1, a, &c__1, b, &c__1, q, &c__1,
  213. z__, &c__1, &info);
  214. chkxer_("DGGHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  215. infoc_1.ok);
  216. infoc_1.infot = 7;
  217. dgghrd_("N", "N", &c__2, &c__1, &c__1, a, &c__1, b, &c__2, q, &c__1,
  218. z__, &c__1, &info);
  219. chkxer_("DGGHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  220. infoc_1.ok);
  221. infoc_1.infot = 9;
  222. dgghrd_("N", "N", &c__2, &c__1, &c__1, a, &c__2, b, &c__1, q, &c__1,
  223. z__, &c__1, &info);
  224. chkxer_("DGGHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  225. infoc_1.ok);
  226. infoc_1.infot = 11;
  227. dgghrd_("V", "N", &c__2, &c__1, &c__1, a, &c__2, b, &c__2, q, &c__1,
  228. z__, &c__1, &info);
  229. chkxer_("DGGHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  230. infoc_1.ok);
  231. infoc_1.infot = 13;
  232. dgghrd_("N", "V", &c__2, &c__1, &c__1, a, &c__2, b, &c__2, q, &c__1,
  233. z__, &c__1, &info);
  234. chkxer_("DGGHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  235. infoc_1.ok);
  236. nt += 9;
  237. /* DHGEQZ */
  238. s_copy(srnamc_1.srnamt, "DHGEQZ", (ftnlen)6, (ftnlen)6);
  239. infoc_1.infot = 1;
  240. dhgeqz_("/", "N", "N", &c__0, &c__1, &c__0, a, &c__1, b, &c__1, r1,
  241. r2, r3, q, &c__1, z__, &c__1, w, &c__18, &info);
  242. chkxer_("DHGEQZ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  243. infoc_1.ok);
  244. infoc_1.infot = 2;
  245. dhgeqz_("E", "/", "N", &c__0, &c__1, &c__0, a, &c__1, b, &c__1, r1,
  246. r2, r3, q, &c__1, z__, &c__1, w, &c__18, &info);
  247. chkxer_("DHGEQZ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  248. infoc_1.ok);
  249. infoc_1.infot = 3;
  250. dhgeqz_("E", "N", "/", &c__0, &c__1, &c__0, a, &c__1, b, &c__1, r1,
  251. r2, r3, q, &c__1, z__, &c__1, w, &c__18, &info);
  252. chkxer_("DHGEQZ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  253. infoc_1.ok);
  254. infoc_1.infot = 4;
  255. dhgeqz_("E", "N", "N", &c_n1, &c__0, &c__0, a, &c__1, b, &c__1, r1,
  256. r2, r3, q, &c__1, z__, &c__1, w, &c__18, &info);
  257. chkxer_("DHGEQZ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  258. infoc_1.ok);
  259. infoc_1.infot = 5;
  260. dhgeqz_("E", "N", "N", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, r1,
  261. r2, r3, q, &c__1, z__, &c__1, w, &c__18, &info);
  262. chkxer_("DHGEQZ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  263. infoc_1.ok);
  264. infoc_1.infot = 6;
  265. dhgeqz_("E", "N", "N", &c__0, &c__1, &c__1, a, &c__1, b, &c__1, r1,
  266. r2, r3, q, &c__1, z__, &c__1, w, &c__18, &info);
  267. chkxer_("DHGEQZ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  268. infoc_1.ok);
  269. infoc_1.infot = 8;
  270. dhgeqz_("E", "N", "N", &c__2, &c__1, &c__1, a, &c__1, b, &c__2, r1,
  271. r2, r3, q, &c__1, z__, &c__1, w, &c__18, &info);
  272. chkxer_("DHGEQZ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  273. infoc_1.ok);
  274. infoc_1.infot = 10;
  275. dhgeqz_("E", "N", "N", &c__2, &c__1, &c__1, a, &c__2, b, &c__1, r1,
  276. r2, r3, q, &c__1, z__, &c__1, w, &c__18, &info);
  277. chkxer_("DHGEQZ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  278. infoc_1.ok);
  279. infoc_1.infot = 15;
  280. dhgeqz_("E", "V", "N", &c__2, &c__1, &c__1, a, &c__2, b, &c__2, r1,
  281. r2, r3, q, &c__1, z__, &c__1, w, &c__18, &info);
  282. chkxer_("DHGEQZ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  283. infoc_1.ok);
  284. infoc_1.infot = 17;
  285. dhgeqz_("E", "N", "V", &c__2, &c__1, &c__1, a, &c__2, b, &c__2, r1,
  286. r2, r3, q, &c__1, z__, &c__1, w, &c__18, &info);
  287. chkxer_("DHGEQZ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  288. infoc_1.ok);
  289. nt += 10;
  290. /* DTGEVC */
  291. s_copy(srnamc_1.srnamt, "DTGEVC", (ftnlen)6, (ftnlen)6);
  292. infoc_1.infot = 1;
  293. dtgevc_("/", "A", sel, &c__0, a, &c__1, b, &c__1, q, &c__1, z__, &
  294. c__1, &c__0, &m, w, &info);
  295. chkxer_("DTGEVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  296. infoc_1.ok);
  297. infoc_1.infot = 2;
  298. dtgevc_("R", "/", sel, &c__0, a, &c__1, b, &c__1, q, &c__1, z__, &
  299. c__1, &c__0, &m, w, &info);
  300. chkxer_("DTGEVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  301. infoc_1.ok);
  302. infoc_1.infot = 4;
  303. dtgevc_("R", "A", sel, &c_n1, a, &c__1, b, &c__1, q, &c__1, z__, &
  304. c__1, &c__0, &m, w, &info);
  305. chkxer_("DTGEVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  306. infoc_1.ok);
  307. infoc_1.infot = 6;
  308. dtgevc_("R", "A", sel, &c__2, a, &c__1, b, &c__2, q, &c__1, z__, &
  309. c__2, &c__0, &m, w, &info);
  310. chkxer_("DTGEVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  311. infoc_1.ok);
  312. infoc_1.infot = 8;
  313. dtgevc_("R", "A", sel, &c__2, a, &c__2, b, &c__1, q, &c__1, z__, &
  314. c__2, &c__0, &m, w, &info);
  315. chkxer_("DTGEVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  316. infoc_1.ok);
  317. infoc_1.infot = 10;
  318. dtgevc_("L", "A", sel, &c__2, a, &c__2, b, &c__2, q, &c__1, z__, &
  319. c__1, &c__0, &m, w, &info);
  320. chkxer_("DTGEVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  321. infoc_1.ok);
  322. infoc_1.infot = 12;
  323. dtgevc_("R", "A", sel, &c__2, a, &c__2, b, &c__2, q, &c__1, z__, &
  324. c__1, &c__0, &m, w, &info);
  325. chkxer_("DTGEVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  326. infoc_1.ok);
  327. infoc_1.infot = 13;
  328. dtgevc_("R", "A", sel, &c__2, a, &c__2, b, &c__2, q, &c__1, z__, &
  329. c__2, &c__1, &m, w, &info);
  330. chkxer_("DTGEVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  331. infoc_1.ok);
  332. nt += 8;
  333. /* Test error exits for the GSV path. */
  334. } else if (lsamen_(&c__3, path, "GSV")) {
  335. /* DGGSVD */
  336. s_copy(srnamc_1.srnamt, "DGGSVD", (ftnlen)6, (ftnlen)6);
  337. infoc_1.infot = 1;
  338. dggsvd_("/", "N", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  339. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  340. info);
  341. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  342. infoc_1.ok);
  343. infoc_1.infot = 2;
  344. dggsvd_("N", "/", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  345. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  346. info);
  347. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  348. infoc_1.ok);
  349. infoc_1.infot = 3;
  350. dggsvd_("N", "N", "/", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  351. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  352. info);
  353. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  354. infoc_1.ok);
  355. infoc_1.infot = 4;
  356. dggsvd_("N", "N", "N", &c_n1, &c__0, &c__0, &dummyk, &dummyl, a, &
  357. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  358. info);
  359. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  360. infoc_1.ok);
  361. infoc_1.infot = 5;
  362. dggsvd_("N", "N", "N", &c__0, &c_n1, &c__0, &dummyk, &dummyl, a, &
  363. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  364. info);
  365. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  366. infoc_1.ok);
  367. infoc_1.infot = 6;
  368. dggsvd_("N", "N", "N", &c__0, &c__0, &c_n1, &dummyk, &dummyl, a, &
  369. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  370. info);
  371. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  372. infoc_1.ok);
  373. infoc_1.infot = 10;
  374. dggsvd_("N", "N", "N", &c__2, &c__1, &c__1, &dummyk, &dummyl, a, &
  375. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  376. info);
  377. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  378. infoc_1.ok);
  379. infoc_1.infot = 12;
  380. dggsvd_("N", "N", "N", &c__1, &c__1, &c__2, &dummyk, &dummyl, a, &
  381. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  382. info);
  383. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  384. infoc_1.ok);
  385. infoc_1.infot = 16;
  386. dggsvd_("U", "N", "N", &c__2, &c__2, &c__2, &dummyk, &dummyl, a, &
  387. c__2, b, &c__2, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  388. info);
  389. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  390. infoc_1.ok);
  391. infoc_1.infot = 18;
  392. dggsvd_("N", "V", "N", &c__1, &c__1, &c__2, &dummyk, &dummyl, a, &
  393. c__1, b, &c__2, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  394. info);
  395. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  396. infoc_1.ok);
  397. infoc_1.infot = 20;
  398. dggsvd_("N", "N", "Q", &c__1, &c__2, &c__1, &dummyk, &dummyl, a, &
  399. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  400. info);
  401. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  402. infoc_1.ok);
  403. nt += 11;
  404. /* DGGSVP */
  405. s_copy(srnamc_1.srnamt, "DGGSVP", (ftnlen)6, (ftnlen)6);
  406. infoc_1.infot = 1;
  407. dggsvp_("/", "N", "N", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, &tola,
  408. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  409. tau, w, &info);
  410. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  411. infoc_1.ok);
  412. infoc_1.infot = 2;
  413. dggsvp_("N", "/", "N", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, &tola,
  414. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  415. tau, w, &info);
  416. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  417. infoc_1.ok);
  418. infoc_1.infot = 3;
  419. dggsvp_("N", "N", "/", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, &tola,
  420. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  421. tau, w, &info);
  422. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  423. infoc_1.ok);
  424. infoc_1.infot = 4;
  425. dggsvp_("N", "N", "N", &c_n1, &c__0, &c__0, a, &c__1, b, &c__1, &tola,
  426. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  427. tau, w, &info);
  428. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  429. infoc_1.ok);
  430. infoc_1.infot = 5;
  431. dggsvp_("N", "N", "N", &c__0, &c_n1, &c__0, a, &c__1, b, &c__1, &tola,
  432. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  433. tau, w, &info);
  434. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  435. infoc_1.ok);
  436. infoc_1.infot = 6;
  437. dggsvp_("N", "N", "N", &c__0, &c__0, &c_n1, a, &c__1, b, &c__1, &tola,
  438. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  439. tau, w, &info);
  440. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  441. infoc_1.ok);
  442. infoc_1.infot = 8;
  443. dggsvp_("N", "N", "N", &c__2, &c__1, &c__1, a, &c__1, b, &c__1, &tola,
  444. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  445. tau, w, &info);
  446. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  447. infoc_1.ok);
  448. infoc_1.infot = 10;
  449. dggsvp_("N", "N", "N", &c__1, &c__2, &c__1, a, &c__1, b, &c__1, &tola,
  450. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  451. tau, w, &info);
  452. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  453. infoc_1.ok);
  454. infoc_1.infot = 16;
  455. dggsvp_("U", "N", "N", &c__2, &c__2, &c__2, a, &c__2, b, &c__2, &tola,
  456. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  457. tau, w, &info);
  458. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  459. infoc_1.ok);
  460. infoc_1.infot = 18;
  461. dggsvp_("N", "V", "N", &c__1, &c__2, &c__1, a, &c__1, b, &c__2, &tola,
  462. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  463. tau, w, &info);
  464. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  465. infoc_1.ok);
  466. infoc_1.infot = 20;
  467. dggsvp_("N", "N", "Q", &c__1, &c__1, &c__2, a, &c__1, b, &c__1, &tola,
  468. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  469. tau, w, &info);
  470. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  471. infoc_1.ok);
  472. nt += 11;
  473. /* DTGSJA */
  474. s_copy(srnamc_1.srnamt, "DTGSJA", (ftnlen)6, (ftnlen)6);
  475. infoc_1.infot = 1;
  476. dtgsja_("/", "N", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  477. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  478. c__1, w, &ncycle, &info);
  479. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  480. infoc_1.ok);
  481. infoc_1.infot = 2;
  482. dtgsja_("N", "/", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  483. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  484. c__1, w, &ncycle, &info);
  485. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  486. infoc_1.ok);
  487. infoc_1.infot = 3;
  488. dtgsja_("N", "N", "/", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  489. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  490. c__1, w, &ncycle, &info);
  491. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  492. infoc_1.ok);
  493. infoc_1.infot = 4;
  494. dtgsja_("N", "N", "N", &c_n1, &c__0, &c__0, &dummyk, &dummyl, a, &
  495. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  496. c__1, w, &ncycle, &info);
  497. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  498. infoc_1.ok);
  499. infoc_1.infot = 5;
  500. dtgsja_("N", "N", "N", &c__0, &c_n1, &c__0, &dummyk, &dummyl, a, &
  501. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  502. c__1, w, &ncycle, &info);
  503. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  504. infoc_1.ok);
  505. infoc_1.infot = 6;
  506. dtgsja_("N", "N", "N", &c__0, &c__0, &c_n1, &dummyk, &dummyl, a, &
  507. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  508. c__1, w, &ncycle, &info);
  509. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  510. infoc_1.ok);
  511. infoc_1.infot = 10;
  512. dtgsja_("N", "N", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  513. c__0, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  514. c__1, w, &ncycle, &info);
  515. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  516. infoc_1.ok);
  517. infoc_1.infot = 12;
  518. dtgsja_("N", "N", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  519. c__1, b, &c__0, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  520. c__1, w, &ncycle, &info);
  521. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  522. infoc_1.ok);
  523. infoc_1.infot = 18;
  524. dtgsja_("U", "N", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  525. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__0, v, &c__1, q, &
  526. c__1, w, &ncycle, &info);
  527. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  528. infoc_1.ok);
  529. infoc_1.infot = 20;
  530. dtgsja_("N", "V", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  531. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__0, q, &
  532. c__1, w, &ncycle, &info);
  533. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  534. infoc_1.ok);
  535. infoc_1.infot = 22;
  536. dtgsja_("N", "N", "Q", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  537. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  538. c__0, w, &ncycle, &info);
  539. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  540. infoc_1.ok);
  541. nt += 11;
  542. /* Test error exits for the GLM path. */
  543. } else if (lsamen_(&c__3, path, "GLM")) {
  544. /* DGGGLM */
  545. s_copy(srnamc_1.srnamt, "DGGGLM", (ftnlen)6, (ftnlen)6);
  546. infoc_1.infot = 1;
  547. dggglm_(&c_n1, &c__0, &c__0, a, &c__1, b, &c__1, r1, r2, r3, w, &
  548. c__18, &info);
  549. chkxer_("DGGGLM", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  550. infoc_1.ok);
  551. infoc_1.infot = 2;
  552. dggglm_(&c__0, &c_n1, &c__0, a, &c__1, b, &c__1, r1, r2, r3, w, &
  553. c__18, &info);
  554. chkxer_("DGGGLM", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  555. infoc_1.ok);
  556. infoc_1.infot = 2;
  557. dggglm_(&c__0, &c__1, &c__0, a, &c__1, b, &c__1, r1, r2, r3, w, &
  558. c__18, &info);
  559. chkxer_("DGGGLM", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  560. infoc_1.ok);
  561. infoc_1.infot = 3;
  562. dggglm_(&c__0, &c__0, &c_n1, a, &c__1, b, &c__1, r1, r2, r3, w, &
  563. c__18, &info);
  564. chkxer_("DGGGLM", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  565. infoc_1.ok);
  566. infoc_1.infot = 3;
  567. dggglm_(&c__1, &c__0, &c__0, a, &c__1, b, &c__1, r1, r2, r3, w, &
  568. c__18, &info);
  569. chkxer_("DGGGLM", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  570. infoc_1.ok);
  571. infoc_1.infot = 5;
  572. dggglm_(&c__0, &c__0, &c__0, a, &c__0, b, &c__1, r1, r2, r3, w, &
  573. c__18, &info);
  574. chkxer_("DGGGLM", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  575. infoc_1.ok);
  576. infoc_1.infot = 7;
  577. dggglm_(&c__0, &c__0, &c__0, a, &c__1, b, &c__0, r1, r2, r3, w, &
  578. c__18, &info);
  579. chkxer_("DGGGLM", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  580. infoc_1.ok);
  581. infoc_1.infot = 12;
  582. dggglm_(&c__1, &c__1, &c__1, a, &c__1, b, &c__1, r1, r2, r3, w, &c__1,
  583. &info);
  584. chkxer_("DGGGLM", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  585. infoc_1.ok);
  586. nt += 8;
  587. /* Test error exits for the LSE path. */
  588. } else if (lsamen_(&c__3, path, "LSE")) {
  589. /* DGGLSE */
  590. s_copy(srnamc_1.srnamt, "DGGLSE", (ftnlen)6, (ftnlen)6);
  591. infoc_1.infot = 1;
  592. dgglse_(&c_n1, &c__0, &c__0, a, &c__1, b, &c__1, r1, r2, r3, w, &
  593. c__18, &info);
  594. chkxer_("DGGLSE", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  595. infoc_1.ok);
  596. infoc_1.infot = 2;
  597. dgglse_(&c__0, &c_n1, &c__0, a, &c__1, b, &c__1, r1, r2, r3, w, &
  598. c__18, &info);
  599. chkxer_("DGGLSE", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  600. infoc_1.ok);
  601. infoc_1.infot = 3;
  602. dgglse_(&c__0, &c__0, &c_n1, a, &c__1, b, &c__1, r1, r2, r3, w, &
  603. c__18, &info);
  604. chkxer_("DGGLSE", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  605. infoc_1.ok);
  606. infoc_1.infot = 3;
  607. dgglse_(&c__0, &c__0, &c__1, a, &c__1, b, &c__1, r1, r2, r3, w, &
  608. c__18, &info);
  609. chkxer_("DGGLSE", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  610. infoc_1.ok);
  611. infoc_1.infot = 3;
  612. dgglse_(&c__0, &c__1, &c__0, a, &c__1, b, &c__1, r1, r2, r3, w, &
  613. c__18, &info);
  614. chkxer_("DGGLSE", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  615. infoc_1.ok);
  616. infoc_1.infot = 5;
  617. dgglse_(&c__0, &c__0, &c__0, a, &c__0, b, &c__1, r1, r2, r3, w, &
  618. c__18, &info);
  619. chkxer_("DGGLSE", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  620. infoc_1.ok);
  621. infoc_1.infot = 7;
  622. dgglse_(&c__0, &c__0, &c__0, a, &c__1, b, &c__0, r1, r2, r3, w, &
  623. c__18, &info);
  624. chkxer_("DGGLSE", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  625. infoc_1.ok);
  626. infoc_1.infot = 12;
  627. dgglse_(&c__1, &c__1, &c__1, a, &c__1, b, &c__1, r1, r2, r3, w, &c__1,
  628. &info);
  629. chkxer_("DGGLSE", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  630. infoc_1.ok);
  631. nt += 8;
  632. /* Test error exits for the GQR path. */
  633. } else if (lsamen_(&c__3, path, "GQR")) {
  634. /* DGGQRF */
  635. s_copy(srnamc_1.srnamt, "DGGQRF", (ftnlen)6, (ftnlen)6);
  636. infoc_1.infot = 1;
  637. dggqrf_(&c_n1, &c__0, &c__0, a, &c__1, r1, b, &c__1, r2, w, &c__18, &
  638. info);
  639. chkxer_("DGGQRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  640. infoc_1.ok);
  641. infoc_1.infot = 2;
  642. dggqrf_(&c__0, &c_n1, &c__0, a, &c__1, r1, b, &c__1, r2, w, &c__18, &
  643. info);
  644. chkxer_("DGGQRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  645. infoc_1.ok);
  646. infoc_1.infot = 3;
  647. dggqrf_(&c__0, &c__0, &c_n1, a, &c__1, r1, b, &c__1, r2, w, &c__18, &
  648. info);
  649. chkxer_("DGGQRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  650. infoc_1.ok);
  651. infoc_1.infot = 5;
  652. dggqrf_(&c__0, &c__0, &c__0, a, &c__0, r1, b, &c__1, r2, w, &c__18, &
  653. info);
  654. chkxer_("DGGQRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  655. infoc_1.ok);
  656. infoc_1.infot = 8;
  657. dggqrf_(&c__0, &c__0, &c__0, a, &c__1, r1, b, &c__0, r2, w, &c__18, &
  658. info);
  659. chkxer_("DGGQRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  660. infoc_1.ok);
  661. infoc_1.infot = 11;
  662. dggqrf_(&c__1, &c__1, &c__2, a, &c__1, r1, b, &c__1, r2, w, &c__1, &
  663. info);
  664. chkxer_("DGGQRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  665. infoc_1.ok);
  666. nt += 6;
  667. /* DGGRQF */
  668. s_copy(srnamc_1.srnamt, "DGGRQF", (ftnlen)6, (ftnlen)6);
  669. infoc_1.infot = 1;
  670. dggrqf_(&c_n1, &c__0, &c__0, a, &c__1, r1, b, &c__1, r2, w, &c__18, &
  671. info);
  672. chkxer_("DGGRQF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  673. infoc_1.ok);
  674. infoc_1.infot = 2;
  675. dggrqf_(&c__0, &c_n1, &c__0, a, &c__1, r1, b, &c__1, r2, w, &c__18, &
  676. info);
  677. chkxer_("DGGRQF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  678. infoc_1.ok);
  679. infoc_1.infot = 3;
  680. dggrqf_(&c__0, &c__0, &c_n1, a, &c__1, r1, b, &c__1, r2, w, &c__18, &
  681. info);
  682. chkxer_("DGGRQF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  683. infoc_1.ok);
  684. infoc_1.infot = 5;
  685. dggrqf_(&c__0, &c__0, &c__0, a, &c__0, r1, b, &c__1, r2, w, &c__18, &
  686. info);
  687. chkxer_("DGGRQF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  688. infoc_1.ok);
  689. infoc_1.infot = 8;
  690. dggrqf_(&c__0, &c__0, &c__0, a, &c__1, r1, b, &c__0, r2, w, &c__18, &
  691. info);
  692. chkxer_("DGGRQF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  693. infoc_1.ok);
  694. infoc_1.infot = 11;
  695. dggrqf_(&c__1, &c__1, &c__2, a, &c__1, r1, b, &c__1, r2, w, &c__1, &
  696. info);
  697. chkxer_("DGGRQF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  698. infoc_1.ok);
  699. nt += 6;
  700. /* Test error exits for the DGS, DGV, DGX, and DXV paths. */
  701. } else if (lsamen_(&c__3, path, "DGS") || lsamen_(&
  702. c__3, path, "DGV") || lsamen_(&c__3, path,
  703. "DGX") || lsamen_(&c__3, path, "DXV")) {
  704. /* DGGES */
  705. s_copy(srnamc_1.srnamt, "DGGES ", (ftnlen)6, (ftnlen)6);
  706. infoc_1.infot = 1;
  707. dgges_("/", "N", "S", (L_fp)dlctes_, &c__1, a, &c__1, b, &c__1, &sdim,
  708. r1, r2, r3, q, &c__1, u, &c__1, w, &c__1, bw, &info);
  709. chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  710. infoc_1.ok);
  711. infoc_1.infot = 2;
  712. dgges_("N", "/", "S", (L_fp)dlctes_, &c__1, a, &c__1, b, &c__1, &sdim,
  713. r1, r2, r3, q, &c__1, u, &c__1, w, &c__1, bw, &info);
  714. chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  715. infoc_1.ok);
  716. infoc_1.infot = 3;
  717. dgges_("N", "V", "/", (L_fp)dlctes_, &c__1, a, &c__1, b, &c__1, &sdim,
  718. r1, r2, r3, q, &c__1, u, &c__1, w, &c__1, bw, &info);
  719. chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  720. infoc_1.ok);
  721. infoc_1.infot = 5;
  722. dgges_("N", "V", "S", (L_fp)dlctes_, &c_n1, a, &c__1, b, &c__1, &sdim,
  723. r1, r2, r3, q, &c__1, u, &c__1, w, &c__1, bw, &info);
  724. chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  725. infoc_1.ok);
  726. infoc_1.infot = 7;
  727. dgges_("N", "V", "S", (L_fp)dlctes_, &c__1, a, &c__0, b, &c__1, &sdim,
  728. r1, r2, r3, q, &c__1, u, &c__1, w, &c__1, bw, &info);
  729. chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  730. infoc_1.ok);
  731. infoc_1.infot = 9;
  732. dgges_("N", "V", "S", (L_fp)dlctes_, &c__1, a, &c__1, b, &c__0, &sdim,
  733. r1, r2, r3, q, &c__1, u, &c__1, w, &c__1, bw, &info);
  734. chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  735. infoc_1.ok);
  736. infoc_1.infot = 15;
  737. dgges_("N", "V", "S", (L_fp)dlctes_, &c__1, a, &c__1, b, &c__1, &sdim,
  738. r1, r2, r3, q, &c__0, u, &c__1, w, &c__1, bw, &info);
  739. chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  740. infoc_1.ok);
  741. infoc_1.infot = 15;
  742. dgges_("V", "V", "S", (L_fp)dlctes_, &c__2, a, &c__2, b, &c__2, &sdim,
  743. r1, r2, r3, q, &c__1, u, &c__2, w, &c__1, bw, &info);
  744. chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  745. infoc_1.ok);
  746. infoc_1.infot = 17;
  747. dgges_("N", "V", "S", (L_fp)dlctes_, &c__1, a, &c__1, b, &c__1, &sdim,
  748. r1, r2, r3, q, &c__1, u, &c__0, w, &c__1, bw, &info);
  749. chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  750. infoc_1.ok);
  751. infoc_1.infot = 17;
  752. dgges_("V", "V", "S", (L_fp)dlctes_, &c__2, a, &c__2, b, &c__2, &sdim,
  753. r1, r2, r3, q, &c__2, u, &c__1, w, &c__1, bw, &info);
  754. chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  755. infoc_1.ok);
  756. infoc_1.infot = 19;
  757. dgges_("V", "V", "S", (L_fp)dlctes_, &c__2, a, &c__2, b, &c__2, &sdim,
  758. r1, r2, r3, q, &c__2, u, &c__2, w, &c__1, bw, &info);
  759. chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  760. infoc_1.ok);
  761. nt += 11;
  762. /* DGGESX */
  763. s_copy(srnamc_1.srnamt, "DGGESX", (ftnlen)6, (ftnlen)6);
  764. infoc_1.infot = 1;
  765. dggesx_("/", "N", "S", (L_fp)dlctsx_, "N", &c__1, a, &c__1, b, &c__1,
  766. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
  767. &c__1, bw, &info)
  768. ;
  769. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  770. infoc_1.ok);
  771. infoc_1.infot = 2;
  772. dggesx_("N", "/", "S", (L_fp)dlctsx_, "N", &c__1, a, &c__1, b, &c__1,
  773. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
  774. &c__1, bw, &info)
  775. ;
  776. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  777. infoc_1.ok);
  778. infoc_1.infot = 3;
  779. dggesx_("V", "V", "/", (L_fp)dlctsx_, "N", &c__1, a, &c__1, b, &c__1,
  780. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
  781. &c__1, bw, &info)
  782. ;
  783. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  784. infoc_1.ok);
  785. infoc_1.infot = 5;
  786. dggesx_("V", "V", "S", (L_fp)dlctsx_, "/", &c__1, a, &c__1, b, &c__1,
  787. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
  788. &c__1, bw, &info)
  789. ;
  790. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  791. infoc_1.ok);
  792. infoc_1.infot = 6;
  793. dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c_n1, a, &c__1, b, &c__1,
  794. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
  795. &c__1, bw, &info)
  796. ;
  797. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  798. infoc_1.ok);
  799. infoc_1.infot = 8;
  800. dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__1, a, &c__0, b, &c__1,
  801. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
  802. &c__1, bw, &info)
  803. ;
  804. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  805. infoc_1.ok);
  806. infoc_1.infot = 10;
  807. dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__1, a, &c__1, b, &c__0,
  808. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
  809. &c__1, bw, &info)
  810. ;
  811. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  812. infoc_1.ok);
  813. infoc_1.infot = 16;
  814. dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__1, a, &c__1, b, &c__1,
  815. &sdim, r1, r2, r3, q, &c__0, u, &c__1, rce, rcv, w, &c__1, iw,
  816. &c__1, bw, &info)
  817. ;
  818. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  819. infoc_1.ok);
  820. infoc_1.infot = 16;
  821. dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__2, a, &c__2, b, &c__2,
  822. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
  823. &c__1, bw, &info)
  824. ;
  825. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  826. infoc_1.ok);
  827. infoc_1.infot = 18;
  828. dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__1, a, &c__1, b, &c__1,
  829. &sdim, r1, r2, r3, q, &c__1, u, &c__0, rce, rcv, w, &c__1, iw,
  830. &c__1, bw, &info)
  831. ;
  832. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  833. infoc_1.ok);
  834. infoc_1.infot = 18;
  835. dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__2, a, &c__2, b, &c__2,
  836. &sdim, r1, r2, r3, q, &c__2, u, &c__1, rce, rcv, w, &c__1, iw,
  837. &c__1, bw, &info)
  838. ;
  839. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  840. infoc_1.ok);
  841. infoc_1.infot = 22;
  842. dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__2, a, &c__2, b, &c__2,
  843. &sdim, r1, r2, r3, q, &c__2, u, &c__2, rce, rcv, w, &c__1, iw,
  844. &c__1, bw, &info)
  845. ;
  846. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  847. infoc_1.ok);
  848. infoc_1.infot = 24;
  849. dggesx_("V", "V", "S", (L_fp)dlctsx_, "V", &c__1, a, &c__1, b, &c__1,
  850. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__32,
  851. iw, &c__0, bw, &info);
  852. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  853. infoc_1.ok);
  854. nt += 13;
  855. /* DGGEV */
  856. s_copy(srnamc_1.srnamt, "DGGEV ", (ftnlen)6, (ftnlen)6);
  857. infoc_1.infot = 1;
  858. dggev_("/", "N", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q, &c__1, u, &
  859. c__1, w, &c__1, &info);
  860. chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  861. infoc_1.ok);
  862. infoc_1.infot = 2;
  863. dggev_("N", "/", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q, &c__1, u, &
  864. c__1, w, &c__1, &info);
  865. chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  866. infoc_1.ok);
  867. infoc_1.infot = 3;
  868. dggev_("V", "V", &c_n1, a, &c__1, b, &c__1, r1, r2, r3, q, &c__1, u, &
  869. c__1, w, &c__1, &info);
  870. chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  871. infoc_1.ok);
  872. infoc_1.infot = 5;
  873. dggev_("V", "V", &c__1, a, &c__0, b, &c__1, r1, r2, r3, q, &c__1, u, &
  874. c__1, w, &c__1, &info);
  875. chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  876. infoc_1.ok);
  877. infoc_1.infot = 7;
  878. dggev_("V", "V", &c__1, a, &c__1, b, &c__0, r1, r2, r3, q, &c__1, u, &
  879. c__1, w, &c__1, &info);
  880. chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  881. infoc_1.ok);
  882. infoc_1.infot = 12;
  883. dggev_("N", "V", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q, &c__0, u, &
  884. c__1, w, &c__1, &info);
  885. chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  886. infoc_1.ok);
  887. infoc_1.infot = 12;
  888. dggev_("V", "V", &c__2, a, &c__2, b, &c__2, r1, r2, r3, q, &c__1, u, &
  889. c__2, w, &c__1, &info);
  890. chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  891. infoc_1.ok);
  892. infoc_1.infot = 14;
  893. dggev_("V", "N", &c__2, a, &c__2, b, &c__2, r1, r2, r3, q, &c__2, u, &
  894. c__0, w, &c__1, &info);
  895. chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  896. infoc_1.ok);
  897. infoc_1.infot = 14;
  898. dggev_("V", "V", &c__2, a, &c__2, b, &c__2, r1, r2, r3, q, &c__2, u, &
  899. c__1, w, &c__1, &info);
  900. chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  901. infoc_1.ok);
  902. infoc_1.infot = 16;
  903. dggev_("V", "V", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q, &c__1, u, &
  904. c__1, w, &c__1, &info);
  905. chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  906. infoc_1.ok);
  907. nt += 10;
  908. /* DGGEVX */
  909. s_copy(srnamc_1.srnamt, "DGGEVX", (ftnlen)6, (ftnlen)6);
  910. infoc_1.infot = 1;
  911. dggevx_("/", "N", "N", "N", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q,
  912. &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  913. w, &c__1, iw, bw, &info);
  914. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  915. infoc_1.ok);
  916. infoc_1.infot = 2;
  917. dggevx_("N", "/", "N", "N", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q,
  918. &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  919. w, &c__1, iw, bw, &info);
  920. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  921. infoc_1.ok);
  922. infoc_1.infot = 3;
  923. dggevx_("N", "N", "/", "N", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q,
  924. &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  925. w, &c__1, iw, bw, &info);
  926. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  927. infoc_1.ok);
  928. infoc_1.infot = 4;
  929. dggevx_("N", "N", "N", "/", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q,
  930. &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  931. w, &c__1, iw, bw, &info);
  932. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  933. infoc_1.ok);
  934. infoc_1.infot = 5;
  935. dggevx_("N", "N", "N", "N", &c_n1, a, &c__1, b, &c__1, r1, r2, r3, q,
  936. &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  937. w, &c__1, iw, bw, &info);
  938. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  939. infoc_1.ok);
  940. infoc_1.infot = 7;
  941. dggevx_("N", "N", "N", "N", &c__1, a, &c__0, b, &c__1, r1, r2, r3, q,
  942. &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  943. w, &c__1, iw, bw, &info);
  944. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  945. infoc_1.ok);
  946. infoc_1.infot = 9;
  947. dggevx_("N", "N", "N", "N", &c__1, a, &c__1, b, &c__0, r1, r2, r3, q,
  948. &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  949. w, &c__1, iw, bw, &info);
  950. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  951. infoc_1.ok);
  952. infoc_1.infot = 14;
  953. dggevx_("N", "N", "N", "N", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q,
  954. &c__0, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  955. w, &c__1, iw, bw, &info);
  956. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  957. infoc_1.ok);
  958. infoc_1.infot = 14;
  959. dggevx_("N", "V", "N", "N", &c__2, a, &c__2, b, &c__2, r1, r2, r3, q,
  960. &c__1, u, &c__2, &c__1, &c__2, ls, rs, &anrm, &bnrm, rce, rcv,
  961. w, &c__1, iw, bw, &info);
  962. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  963. infoc_1.ok);
  964. infoc_1.infot = 16;
  965. dggevx_("N", "N", "N", "N", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q,
  966. &c__1, u, &c__0, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  967. w, &c__1, iw, bw, &info);
  968. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  969. infoc_1.ok);
  970. infoc_1.infot = 16;
  971. dggevx_("N", "N", "V", "N", &c__2, a, &c__2, b, &c__2, r1, r2, r3, q,
  972. &c__2, u, &c__1, &c__1, &c__2, ls, rs, &anrm, &bnrm, rce, rcv,
  973. w, &c__1, iw, bw, &info);
  974. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  975. infoc_1.ok);
  976. infoc_1.infot = 26;
  977. dggevx_("N", "N", "V", "N", &c__2, a, &c__2, b, &c__2, r1, r2, r3, q,
  978. &c__2, u, &c__2, &c__1, &c__2, ls, rs, &anrm, &bnrm, rce, rcv,
  979. w, &c__1, iw, bw, &info);
  980. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  981. infoc_1.ok);
  982. nt += 12;
  983. /* DTGEXC */
  984. s_copy(srnamc_1.srnamt, "DTGEXC", (ftnlen)6, (ftnlen)6);
  985. infoc_1.infot = 3;
  986. dtgexc_(&c_true, &c_true, &c_n1, a, &c__1, b, &c__1, q, &c__1, z__, &
  987. c__1, &ifst, &ilst, w, &c__1, &info);
  988. chkxer_("DTGEXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  989. infoc_1.ok);
  990. infoc_1.infot = 5;
  991. dtgexc_(&c_true, &c_true, &c__1, a, &c__0, b, &c__1, q, &c__1, z__, &
  992. c__1, &ifst, &ilst, w, &c__1, &info);
  993. chkxer_("DTGEXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  994. infoc_1.ok);
  995. infoc_1.infot = 7;
  996. dtgexc_(&c_true, &c_true, &c__1, a, &c__1, b, &c__0, q, &c__1, z__, &
  997. c__1, &ifst, &ilst, w, &c__1, &info);
  998. chkxer_("DTGEXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  999. infoc_1.ok);
  1000. infoc_1.infot = 9;
  1001. dtgexc_(&c_false, &c_true, &c__1, a, &c__1, b, &c__1, q, &c__0, z__, &
  1002. c__1, &ifst, &ilst, w, &c__1, &info);
  1003. chkxer_("DTGEXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1004. infoc_1.ok);
  1005. infoc_1.infot = 9;
  1006. dtgexc_(&c_true, &c_true, &c__1, a, &c__1, b, &c__1, q, &c__0, z__, &
  1007. c__1, &ifst, &ilst, w, &c__1, &info);
  1008. chkxer_("DTGEXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1009. infoc_1.ok);
  1010. infoc_1.infot = 11;
  1011. dtgexc_(&c_true, &c_false, &c__1, a, &c__1, b, &c__1, q, &c__1, z__, &
  1012. c__0, &ifst, &ilst, w, &c__1, &info);
  1013. chkxer_("DTGEXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1014. infoc_1.ok);
  1015. infoc_1.infot = 11;
  1016. dtgexc_(&c_true, &c_true, &c__1, a, &c__1, b, &c__1, q, &c__1, z__, &
  1017. c__0, &ifst, &ilst, w, &c__1, &info);
  1018. chkxer_("DTGEXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1019. infoc_1.ok);
  1020. infoc_1.infot = 15;
  1021. dtgexc_(&c_true, &c_true, &c__1, a, &c__1, b, &c__1, q, &c__1, z__, &
  1022. c__1, &ifst, &ilst, w, &c__0, &info);
  1023. chkxer_("DTGEXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1024. infoc_1.ok);
  1025. nt += 8;
  1026. /* DTGSEN */
  1027. s_copy(srnamc_1.srnamt, "DTGSEN", (ftnlen)6, (ftnlen)6);
  1028. infoc_1.infot = 1;
  1029. dtgsen_(&c_n1, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1030. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
  1031. iw, &c__1, &info);
  1032. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1033. infoc_1.ok);
  1034. infoc_1.infot = 5;
  1035. dtgsen_(&c__1, &c_true, &c_true, sel, &c_n1, a, &c__1, b, &c__1, r1,
  1036. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
  1037. iw, &c__1, &info);
  1038. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1039. infoc_1.ok);
  1040. infoc_1.infot = 7;
  1041. dtgsen_(&c__1, &c_true, &c_true, sel, &c__1, a, &c__0, b, &c__1, r1,
  1042. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
  1043. iw, &c__1, &info);
  1044. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1045. infoc_1.ok);
  1046. infoc_1.infot = 9;
  1047. dtgsen_(&c__1, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__0, r1,
  1048. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
  1049. iw, &c__1, &info);
  1050. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1051. infoc_1.ok);
  1052. infoc_1.infot = 14;
  1053. dtgsen_(&c__1, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1054. r2, r3, q, &c__0, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
  1055. iw, &c__1, &info);
  1056. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1057. infoc_1.ok);
  1058. infoc_1.infot = 16;
  1059. dtgsen_(&c__1, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1060. r2, r3, q, &c__1, z__, &c__0, &m, &tola, &tolb, rcv, w, &c__1,
  1061. iw, &c__1, &info);
  1062. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1063. infoc_1.ok);
  1064. infoc_1.infot = 22;
  1065. dtgsen_(&c__0, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1066. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
  1067. iw, &c__1, &info);
  1068. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1069. infoc_1.ok);
  1070. infoc_1.infot = 22;
  1071. dtgsen_(&c__1, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1072. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
  1073. iw, &c__1, &info);
  1074. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1075. infoc_1.ok);
  1076. infoc_1.infot = 22;
  1077. dtgsen_(&c__2, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1078. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
  1079. iw, &c__1, &info);
  1080. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1081. infoc_1.ok);
  1082. infoc_1.infot = 24;
  1083. dtgsen_(&c__0, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1084. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &
  1085. c__20, iw, &c__0, &info);
  1086. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1087. infoc_1.ok);
  1088. infoc_1.infot = 24;
  1089. dtgsen_(&c__1, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1090. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &
  1091. c__20, iw, &c__0, &info);
  1092. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1093. infoc_1.ok);
  1094. infoc_1.infot = 24;
  1095. dtgsen_(&c__2, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1096. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &
  1097. c__20, iw, &c__1, &info);
  1098. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1099. infoc_1.ok);
  1100. nt += 12;
  1101. /* DTGSNA */
  1102. s_copy(srnamc_1.srnamt, "DTGSNA", (ftnlen)6, (ftnlen)6);
  1103. infoc_1.infot = 1;
  1104. dtgsna_("/", "A", sel, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &c__1,
  1105. r1, r2, &c__1, &m, w, &c__1, iw, &info);
  1106. chkxer_("DTGSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1107. infoc_1.ok);
  1108. infoc_1.infot = 2;
  1109. dtgsna_("B", "/", sel, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &c__1,
  1110. r1, r2, &c__1, &m, w, &c__1, iw, &info);
  1111. chkxer_("DTGSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1112. infoc_1.ok);
  1113. infoc_1.infot = 4;
  1114. dtgsna_("B", "A", sel, &c_n1, a, &c__1, b, &c__1, q, &c__1, u, &c__1,
  1115. r1, r2, &c__1, &m, w, &c__1, iw, &info);
  1116. chkxer_("DTGSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1117. infoc_1.ok);
  1118. infoc_1.infot = 6;
  1119. dtgsna_("B", "A", sel, &c__1, a, &c__0, b, &c__1, q, &c__1, u, &c__1,
  1120. r1, r2, &c__1, &m, w, &c__1, iw, &info);
  1121. chkxer_("DTGSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1122. infoc_1.ok);
  1123. infoc_1.infot = 8;
  1124. dtgsna_("B", "A", sel, &c__1, a, &c__1, b, &c__0, q, &c__1, u, &c__1,
  1125. r1, r2, &c__1, &m, w, &c__1, iw, &info);
  1126. chkxer_("DTGSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1127. infoc_1.ok);
  1128. infoc_1.infot = 10;
  1129. dtgsna_("E", "A", sel, &c__1, a, &c__1, b, &c__1, q, &c__0, u, &c__1,
  1130. r1, r2, &c__1, &m, w, &c__1, iw, &info);
  1131. chkxer_("DTGSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1132. infoc_1.ok);
  1133. infoc_1.infot = 12;
  1134. dtgsna_("E", "A", sel, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &c__0,
  1135. r1, r2, &c__1, &m, w, &c__1, iw, &info);
  1136. chkxer_("DTGSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1137. infoc_1.ok);
  1138. infoc_1.infot = 15;
  1139. dtgsna_("E", "A", sel, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &c__1,
  1140. r1, r2, &c__0, &m, w, &c__1, iw, &info);
  1141. chkxer_("DTGSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1142. infoc_1.ok);
  1143. infoc_1.infot = 18;
  1144. dtgsna_("E", "A", sel, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &c__1,
  1145. r1, r2, &c__1, &m, w, &c__0, iw, &info);
  1146. chkxer_("DTGSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1147. infoc_1.ok);
  1148. nt += 9;
  1149. /* DTGSYL */
  1150. s_copy(srnamc_1.srnamt, "DTGSYL", (ftnlen)6, (ftnlen)6);
  1151. infoc_1.infot = 1;
  1152. dtgsyl_("/", &c__0, &c__1, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &
  1153. c__1, v, &c__1, z__, &c__1, &scale, &dif, w, &c__1, iw, &info);
  1154. chkxer_("DTGSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1155. infoc_1.ok);
  1156. infoc_1.infot = 2;
  1157. dtgsyl_("N", &c_n1, &c__1, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &
  1158. c__1, v, &c__1, z__, &c__1, &scale, &dif, w, &c__1, iw, &info);
  1159. chkxer_("DTGSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1160. infoc_1.ok);
  1161. infoc_1.infot = 3;
  1162. dtgsyl_("N", &c__0, &c__0, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &
  1163. c__1, v, &c__1, z__, &c__1, &scale, &dif, w, &c__1, iw, &info);
  1164. chkxer_("DTGSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1165. infoc_1.ok);
  1166. infoc_1.infot = 4;
  1167. dtgsyl_("N", &c__0, &c__1, &c__0, a, &c__1, b, &c__1, q, &c__1, u, &
  1168. c__1, v, &c__1, z__, &c__1, &scale, &dif, w, &c__1, iw, &info);
  1169. chkxer_("DTGSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1170. infoc_1.ok);
  1171. infoc_1.infot = 6;
  1172. dtgsyl_("N", &c__0, &c__1, &c__1, a, &c__0, b, &c__1, q, &c__1, u, &
  1173. c__1, v, &c__1, z__, &c__1, &scale, &dif, w, &c__1, iw, &info);
  1174. chkxer_("DTGSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1175. infoc_1.ok);
  1176. infoc_1.infot = 8;
  1177. dtgsyl_("N", &c__0, &c__1, &c__1, a, &c__1, b, &c__0, q, &c__1, u, &
  1178. c__1, v, &c__1, z__, &c__1, &scale, &dif, w, &c__1, iw, &info);
  1179. chkxer_("DTGSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1180. infoc_1.ok);
  1181. infoc_1.infot = 10;
  1182. dtgsyl_("N", &c__0, &c__1, &c__1, a, &c__1, b, &c__1, q, &c__0, u, &
  1183. c__1, v, &c__1, z__, &c__1, &scale, &dif, w, &c__1, iw, &info);
  1184. chkxer_("DTGSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1185. infoc_1.ok);
  1186. infoc_1.infot = 12;
  1187. dtgsyl_("N", &c__0, &c__1, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &
  1188. c__0, v, &c__1, z__, &c__1, &scale, &dif, w, &c__1, iw, &info);
  1189. chkxer_("DTGSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1190. infoc_1.ok);
  1191. infoc_1.infot = 14;
  1192. dtgsyl_("N", &c__0, &c__1, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &
  1193. c__1, v, &c__0, z__, &c__1, &scale, &dif, w, &c__