PageRenderTime 56ms CodeModel.GetById 12ms RepoModel.GetById 0ms app.codeStats 0ms

/contrib/clapack-3.2.1-CMAKE/TESTING/EIG/derrgg.c

http://github.com/openmeeg/openmeeg
C | 1319 lines | 1160 code | 80 blank | 79 comment | 19 complexity | d630fd3674d9913f628aae541ed27c50 MD5 | raw file
Possible License(s): BSD-3-Clause

Large files files are truncated, but you can click here to view the full file

  1. /* derrgg.f -- translated by f2c (version 20061008).
  2. You must link the resulting object file with libf2c:
  3. on Microsoft Windows system, link with libf2c.lib;
  4. on Linux or Unix systems, link with .../path/to/libf2c.a -lm
  5. or, if you install libf2c.a in a standard place, with -lf2c -lm
  6. -- in that order, at the end of the command line, as in
  7. cc *.o -lf2c -lm
  8. Source for libf2c is in /netlib/f2c/libf2c.zip, e.g.,
  9. http://www.netlib.org/f2c/libf2c.zip
  10. */
  11. #include "f2c.h"
  12. #include "blaswrap.h"
  13. /* Common Block Declarations */
  14. struct {
  15. integer infot, nout;
  16. logical ok, lerr;
  17. } infoc_;
  18. #define infoc_1 infoc_
  19. struct {
  20. char srnamt[32];
  21. } srnamc_;
  22. #define srnamc_1 srnamc_
  23. /* Table of constant values */
  24. static integer c__2 = 2;
  25. static integer c__0 = 0;
  26. static integer c__1 = 1;
  27. static integer c_n1 = -1;
  28. static integer c__18 = 18;
  29. static integer c__3 = 3;
  30. static integer c__32 = 32;
  31. static logical c_true = TRUE_;
  32. static logical c_false = FALSE_;
  33. static integer c__20 = 20;
  34. /* Subroutine */ int derrgg_(char *path, integer *nunit)
  35. {
  36. /* Format strings */
  37. static char fmt_9999[] = "(1x,a3,\002 routines passed the tests of the e"
  38. "rror exits (\002,i3,\002 tests done)\002)";
  39. static char fmt_9998[] = "(\002 *** \002,a3,\002 routines failed the tes"
  40. "ts of the error \002,\002exits ***\002)";
  41. /* Builtin functions */
  42. integer s_wsle(cilist *), e_wsle(void);
  43. /* Subroutine */ int s_copy(char *, char *, ftnlen, ftnlen);
  44. integer s_wsfe(cilist *), do_fio(integer *, char *, ftnlen), e_wsfe(void);
  45. /* Local variables */
  46. doublereal a[9] /* was [3][3] */, b[9] /* was [3][3] */;
  47. integer i__, j, m;
  48. doublereal q[9] /* was [3][3] */, u[9] /* was [3][3] */, v[9] /*
  49. was [3][3] */, w[18], z__[9] /* was [3][3] */;
  50. char c2[2];
  51. doublereal r1[3], r2[3], r3[3];
  52. logical bw[3];
  53. doublereal ls[3];
  54. integer iw[3], nt;
  55. doublereal rs[3], dif, rce[2];
  56. logical sel[3];
  57. doublereal tau[3], rcv[2];
  58. integer info, sdim;
  59. doublereal anrm, bnrm, tola, tolb;
  60. integer ifst, ilst;
  61. doublereal scale;
  62. extern /* Subroutine */ int dgges_(char *, char *, char *, L_fp, integer *
  63. , doublereal *, integer *, doublereal *, integer *, integer *,
  64. doublereal *, doublereal *, doublereal *, doublereal *, integer *,
  65. doublereal *, integer *, doublereal *, integer *, logical *,
  66. integer *), dggev_(char *, char *,
  67. integer *, doublereal *, integer *, doublereal *, integer *,
  68. doublereal *, doublereal *, doublereal *, doublereal *, integer *,
  69. doublereal *, integer *, doublereal *, integer *, integer *), dgghrd_(char *, char *, integer *, integer *,
  70. integer *, doublereal *, integer *, doublereal *, integer *,
  71. doublereal *, integer *, doublereal *, integer *, integer *), dggglm_(integer *, integer *, integer *,
  72. doublereal *, integer *, doublereal *, integer *, doublereal *,
  73. doublereal *, doublereal *, doublereal *, integer *, integer *),
  74. dgglse_(integer *, integer *, integer *, doublereal *, integer *,
  75. doublereal *, integer *, doublereal *, doublereal *, doublereal *,
  76. doublereal *, integer *, integer *), dggqrf_(integer *, integer *
  77. , integer *, doublereal *, integer *, doublereal *, doublereal *,
  78. integer *, doublereal *, doublereal *, integer *, integer *),
  79. dggrqf_(integer *, integer *, integer *, doublereal *, integer *,
  80. doublereal *, doublereal *, integer *, doublereal *, doublereal *,
  81. integer *, integer *);
  82. integer ncycle;
  83. extern logical dlctes_(), lsamen_(integer *, char *, char *);
  84. extern /* Subroutine */ int dggsvd_(char *, char *, char *, integer *,
  85. integer *, integer *, integer *, integer *, doublereal *, integer
  86. *, doublereal *, integer *, doublereal *, doublereal *,
  87. doublereal *, integer *, doublereal *, integer *, doublereal *,
  88. integer *, doublereal *, integer *, integer *), dggesx_(char *, char *, char *, L_fp, char *, integer *,
  89. doublereal *, integer *, doublereal *, integer *, integer *,
  90. doublereal *, doublereal *, doublereal *, doublereal *, integer *,
  91. doublereal *, integer *, doublereal *, doublereal *, doublereal *
  92. , integer *, integer *, integer *, logical *, integer *), dhgeqz_(char *, char *, char *, integer *
  93. , integer *, integer *, doublereal *, integer *, doublereal *,
  94. integer *, doublereal *, doublereal *, doublereal *, doublereal *,
  95. integer *, doublereal *, integer *, doublereal *, integer *,
  96. integer *), dtgevc_(char *, char *,
  97. logical *, integer *, doublereal *, integer *, doublereal *,
  98. integer *, doublereal *, integer *, doublereal *, integer *,
  99. integer *, integer *, doublereal *, integer *),
  100. chkxer_(char *, integer *, integer *, logical *, logical *), dggevx_(char *, char *, char *, char *, integer *,
  101. doublereal *, integer *, doublereal *, integer *, doublereal *,
  102. doublereal *, doublereal *, doublereal *, integer *, doublereal *,
  103. integer *, integer *, integer *, doublereal *, doublereal *,
  104. doublereal *, doublereal *, doublereal *, doublereal *,
  105. doublereal *, integer *, integer *, logical *, integer *), dtgexc_(logical *, logical *, integer *,
  106. doublereal *, integer *, doublereal *, integer *, doublereal *,
  107. integer *, doublereal *, integer *, integer *, integer *,
  108. doublereal *, integer *, integer *), dtgsen_(integer *, logical *,
  109. logical *, logical *, integer *, doublereal *, integer *,
  110. doublereal *, integer *, doublereal *, doublereal *, doublereal *,
  111. doublereal *, integer *, doublereal *, integer *, integer *,
  112. doublereal *, doublereal *, doublereal *, doublereal *, integer *,
  113. integer *, integer *, integer *), dtgsja_(char *, char *, char *,
  114. integer *, integer *, integer *, integer *, integer *,
  115. doublereal *, integer *, doublereal *, integer *, doublereal *,
  116. doublereal *, doublereal *, doublereal *, doublereal *, integer *,
  117. doublereal *, integer *, doublereal *, integer *, doublereal *,
  118. integer *, integer *), dtgsna_(char *,
  119. char *, logical *, integer *, doublereal *, integer *, doublereal
  120. *, integer *, doublereal *, integer *, doublereal *, integer *,
  121. doublereal *, doublereal *, integer *, integer *, doublereal *,
  122. integer *, integer *, integer *), dggsvp_(char *,
  123. char *, char *, integer *, integer *, integer *, doublereal *,
  124. integer *, doublereal *, integer *, doublereal *, doublereal *,
  125. integer *, integer *, doublereal *, integer *, doublereal *,
  126. integer *, doublereal *, integer *, integer *, doublereal *,
  127. doublereal *, integer *);
  128. extern logical dlctsx_();
  129. integer dummyk, dummyl;
  130. extern /* Subroutine */ int dtgsyl_(char *, integer *, integer *, integer
  131. *, doublereal *, integer *, doublereal *, integer *, doublereal *,
  132. integer *, doublereal *, integer *, doublereal *, integer *,
  133. doublereal *, integer *, doublereal *, doublereal *, doublereal *,
  134. integer *, integer *, integer *);
  135. /* Fortran I/O blocks */
  136. static cilist io___1 = { 0, 0, 0, 0, 0 };
  137. static cilist io___38 = { 0, 0, 0, fmt_9999, 0 };
  138. static cilist io___39 = { 0, 0, 0, fmt_9998, 0 };
  139. /* -- LAPACK test routine (version 3.1) -- */
  140. /* Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. */
  141. /* November 2006 */
  142. /* .. Scalar Arguments .. */
  143. /* .. */
  144. /* Purpose */
  145. /* ======= */
  146. /* DERRGG tests the error exits for DGGES, DGGESX, DGGEV, DGGEVX, */
  147. /* DGGGLM, DGGHRD, DGGLSE, DGGQRF, DGGRQF, DGGSVD, DGGSVP, DHGEQZ, */
  148. /* DTGEVC, DTGEXC, DTGSEN, DTGSJA, DTGSNA, and DTGSYL. */
  149. /* Arguments */
  150. /* ========= */
  151. /* PATH (input) CHARACTER*3 */
  152. /* The LAPACK path name for the routines to be tested. */
  153. /* NUNIT (input) INTEGER */
  154. /* The unit number for output. */
  155. /* ===================================================================== */
  156. /* .. Parameters .. */
  157. /* .. */
  158. /* .. Local Scalars .. */
  159. /* .. */
  160. /* .. Local Arrays .. */
  161. /* .. */
  162. /* .. External Functions .. */
  163. /* .. */
  164. /* .. External Subroutines .. */
  165. /* .. */
  166. /* .. Scalars in Common .. */
  167. /* .. */
  168. /* .. Common blocks .. */
  169. /* .. */
  170. /* .. Executable Statements .. */
  171. infoc_1.nout = *nunit;
  172. io___1.ciunit = infoc_1.nout;
  173. s_wsle(&io___1);
  174. e_wsle();
  175. s_copy(c2, path + 1, (ftnlen)2, (ftnlen)2);
  176. /* Set the variables to innocuous values. */
  177. for (j = 1; j <= 3; ++j) {
  178. sel[j - 1] = TRUE_;
  179. for (i__ = 1; i__ <= 3; ++i__) {
  180. a[i__ + j * 3 - 4] = 0.;
  181. b[i__ + j * 3 - 4] = 0.;
  182. /* L10: */
  183. }
  184. /* L20: */
  185. }
  186. for (i__ = 1; i__ <= 3; ++i__) {
  187. a[i__ + i__ * 3 - 4] = 1.;
  188. b[i__ + i__ * 3 - 4] = 1.;
  189. /* L30: */
  190. }
  191. infoc_1.ok = TRUE_;
  192. tola = 1.;
  193. tolb = 1.;
  194. ifst = 1;
  195. ilst = 1;
  196. nt = 0;
  197. /* Test error exits for the GG path. */
  198. if (lsamen_(&c__2, c2, "GG")) {
  199. /* DGGHRD */
  200. s_copy(srnamc_1.srnamt, "DGGHRD", (ftnlen)32, (ftnlen)6);
  201. infoc_1.infot = 1;
  202. dgghrd_("/", "N", &c__0, &c__1, &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 = 2;
  207. dgghrd_("N", "/", &c__0, &c__1, &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 = 3;
  212. dgghrd_("N", "N", &c_n1, &c__0, &c__0, 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 = 4;
  217. dgghrd_("N", "N", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, 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 = 5;
  222. dgghrd_("N", "N", &c__0, &c__1, &c__1, a, &c__1, 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 = 7;
  227. dgghrd_("N", "N", &c__2, &c__1, &c__1, a, &c__1, 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 = 9;
  232. dgghrd_("N", "N", &c__2, &c__1, &c__1, a, &c__2, b, &c__1, 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. infoc_1.infot = 11;
  237. dgghrd_("V", "N", &c__2, &c__1, &c__1, a, &c__2, b, &c__2, q, &c__1,
  238. z__, &c__1, &info);
  239. chkxer_("DGGHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  240. infoc_1.ok);
  241. infoc_1.infot = 13;
  242. dgghrd_("N", "V", &c__2, &c__1, &c__1, a, &c__2, b, &c__2, q, &c__1,
  243. z__, &c__1, &info);
  244. chkxer_("DGGHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  245. infoc_1.ok);
  246. nt += 9;
  247. /* DHGEQZ */
  248. s_copy(srnamc_1.srnamt, "DHGEQZ", (ftnlen)32, (ftnlen)6);
  249. infoc_1.infot = 1;
  250. dhgeqz_("/", "N", "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 = 2;
  255. dhgeqz_("E", "/", "N", &c__0, &c__1, &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 = 3;
  260. dhgeqz_("E", "N", "/", &c__0, &c__1, &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 = 4;
  265. dhgeqz_("E", "N", "N", &c_n1, &c__0, &c__0, 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 = 5;
  270. dhgeqz_("E", "N", "N", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, 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 = 6;
  275. dhgeqz_("E", "N", "N", &c__0, &c__1, &c__1, a, &c__1, 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 = 8;
  280. dhgeqz_("E", "N", "N", &c__2, &c__1, &c__1, a, &c__1, 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 = 10;
  285. dhgeqz_("E", "N", "N", &c__2, &c__1, &c__1, a, &c__2, b, &c__1, 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. infoc_1.infot = 15;
  290. dhgeqz_("E", "V", "N", &c__2, &c__1, &c__1, a, &c__2, b, &c__2, r1,
  291. r2, r3, q, &c__1, z__, &c__1, w, &c__18, &info);
  292. chkxer_("DHGEQZ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  293. infoc_1.ok);
  294. infoc_1.infot = 17;
  295. dhgeqz_("E", "N", "V", &c__2, &c__1, &c__1, a, &c__2, b, &c__2, r1,
  296. r2, r3, q, &c__1, z__, &c__1, w, &c__18, &info);
  297. chkxer_("DHGEQZ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  298. infoc_1.ok);
  299. nt += 10;
  300. /* DTGEVC */
  301. s_copy(srnamc_1.srnamt, "DTGEVC", (ftnlen)32, (ftnlen)6);
  302. infoc_1.infot = 1;
  303. dtgevc_("/", "A", sel, &c__0, 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 = 2;
  308. dtgevc_("R", "/", sel, &c__0, a, &c__1, b, &c__1, q, &c__1, z__, &
  309. c__1, &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 = 4;
  313. dtgevc_("R", "A", sel, &c_n1, a, &c__1, b, &c__1, q, &c__1, z__, &
  314. c__1, &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 = 6;
  318. dtgevc_("R", "A", sel, &c__2, a, &c__1, b, &c__2, q, &c__1, z__, &
  319. c__2, &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 = 8;
  323. dtgevc_("R", "A", sel, &c__2, a, &c__2, b, &c__1, q, &c__1, z__, &
  324. c__2, &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 = 10;
  328. dtgevc_("L", "A", sel, &c__2, a, &c__2, b, &c__2, q, &c__1, z__, &
  329. c__1, &c__0, &m, w, &info);
  330. chkxer_("DTGEVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  331. infoc_1.ok);
  332. infoc_1.infot = 12;
  333. dtgevc_("R", "A", sel, &c__2, a, &c__2, b, &c__2, q, &c__1, z__, &
  334. c__1, &c__0, &m, w, &info);
  335. chkxer_("DTGEVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  336. infoc_1.ok);
  337. infoc_1.infot = 13;
  338. dtgevc_("R", "A", sel, &c__2, a, &c__2, b, &c__2, q, &c__1, z__, &
  339. c__2, &c__1, &m, w, &info);
  340. chkxer_("DTGEVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  341. infoc_1.ok);
  342. nt += 8;
  343. /* Test error exits for the GSV path. */
  344. } else if (lsamen_(&c__3, path, "GSV")) {
  345. /* DGGSVD */
  346. s_copy(srnamc_1.srnamt, "DGGSVD", (ftnlen)32, (ftnlen)6);
  347. infoc_1.infot = 1;
  348. dggsvd_("/", "N", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  349. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  350. info);
  351. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  352. infoc_1.ok);
  353. infoc_1.infot = 2;
  354. dggsvd_("N", "/", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  355. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  356. info);
  357. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  358. infoc_1.ok);
  359. infoc_1.infot = 3;
  360. dggsvd_("N", "N", "/", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  361. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  362. info);
  363. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  364. infoc_1.ok);
  365. infoc_1.infot = 4;
  366. dggsvd_("N", "N", "N", &c_n1, &c__0, &c__0, &dummyk, &dummyl, a, &
  367. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  368. info);
  369. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  370. infoc_1.ok);
  371. infoc_1.infot = 5;
  372. dggsvd_("N", "N", "N", &c__0, &c_n1, &c__0, &dummyk, &dummyl, a, &
  373. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  374. info);
  375. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  376. infoc_1.ok);
  377. infoc_1.infot = 6;
  378. dggsvd_("N", "N", "N", &c__0, &c__0, &c_n1, &dummyk, &dummyl, a, &
  379. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  380. info);
  381. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  382. infoc_1.ok);
  383. infoc_1.infot = 10;
  384. dggsvd_("N", "N", "N", &c__2, &c__1, &c__1, &dummyk, &dummyl, a, &
  385. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  386. info);
  387. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  388. infoc_1.ok);
  389. infoc_1.infot = 12;
  390. dggsvd_("N", "N", "N", &c__1, &c__1, &c__2, &dummyk, &dummyl, a, &
  391. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  392. info);
  393. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  394. infoc_1.ok);
  395. infoc_1.infot = 16;
  396. dggsvd_("U", "N", "N", &c__2, &c__2, &c__2, &dummyk, &dummyl, a, &
  397. c__2, b, &c__2, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  398. info);
  399. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  400. infoc_1.ok);
  401. infoc_1.infot = 18;
  402. dggsvd_("N", "V", "N", &c__1, &c__1, &c__2, &dummyk, &dummyl, a, &
  403. c__1, b, &c__2, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  404. info);
  405. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  406. infoc_1.ok);
  407. infoc_1.infot = 20;
  408. dggsvd_("N", "N", "Q", &c__1, &c__2, &c__1, &dummyk, &dummyl, a, &
  409. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  410. info);
  411. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  412. infoc_1.ok);
  413. nt += 11;
  414. /* DGGSVP */
  415. s_copy(srnamc_1.srnamt, "DGGSVP", (ftnlen)32, (ftnlen)6);
  416. infoc_1.infot = 1;
  417. dggsvp_("/", "N", "N", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, &tola,
  418. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  419. tau, w, &info);
  420. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  421. infoc_1.ok);
  422. infoc_1.infot = 2;
  423. dggsvp_("N", "/", "N", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, &tola,
  424. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  425. tau, w, &info);
  426. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  427. infoc_1.ok);
  428. infoc_1.infot = 3;
  429. dggsvp_("N", "N", "/", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, &tola,
  430. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  431. tau, w, &info);
  432. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  433. infoc_1.ok);
  434. infoc_1.infot = 4;
  435. dggsvp_("N", "N", "N", &c_n1, &c__0, &c__0, a, &c__1, b, &c__1, &tola,
  436. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  437. tau, w, &info);
  438. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  439. infoc_1.ok);
  440. infoc_1.infot = 5;
  441. dggsvp_("N", "N", "N", &c__0, &c_n1, &c__0, a, &c__1, b, &c__1, &tola,
  442. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  443. tau, w, &info);
  444. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  445. infoc_1.ok);
  446. infoc_1.infot = 6;
  447. dggsvp_("N", "N", "N", &c__0, &c__0, &c_n1, a, &c__1, b, &c__1, &tola,
  448. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  449. tau, w, &info);
  450. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  451. infoc_1.ok);
  452. infoc_1.infot = 8;
  453. dggsvp_("N", "N", "N", &c__2, &c__1, &c__1, a, &c__1, b, &c__1, &tola,
  454. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  455. tau, w, &info);
  456. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  457. infoc_1.ok);
  458. infoc_1.infot = 10;
  459. dggsvp_("N", "N", "N", &c__1, &c__2, &c__1, a, &c__1, b, &c__1, &tola,
  460. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  461. tau, w, &info);
  462. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  463. infoc_1.ok);
  464. infoc_1.infot = 16;
  465. dggsvp_("U", "N", "N", &c__2, &c__2, &c__2, a, &c__2, b, &c__2, &tola,
  466. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  467. tau, w, &info);
  468. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  469. infoc_1.ok);
  470. infoc_1.infot = 18;
  471. dggsvp_("N", "V", "N", &c__1, &c__2, &c__1, a, &c__1, b, &c__2, &tola,
  472. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  473. tau, w, &info);
  474. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  475. infoc_1.ok);
  476. infoc_1.infot = 20;
  477. dggsvp_("N", "N", "Q", &c__1, &c__1, &c__2, a, &c__1, b, &c__1, &tola,
  478. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  479. tau, w, &info);
  480. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  481. infoc_1.ok);
  482. nt += 11;
  483. /* DTGSJA */
  484. s_copy(srnamc_1.srnamt, "DTGSJA", (ftnlen)32, (ftnlen)6);
  485. infoc_1.infot = 1;
  486. dtgsja_("/", "N", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  487. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  488. c__1, w, &ncycle, &info);
  489. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  490. infoc_1.ok);
  491. infoc_1.infot = 2;
  492. dtgsja_("N", "/", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  493. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  494. c__1, w, &ncycle, &info);
  495. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  496. infoc_1.ok);
  497. infoc_1.infot = 3;
  498. dtgsja_("N", "N", "/", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  499. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  500. c__1, w, &ncycle, &info);
  501. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  502. infoc_1.ok);
  503. infoc_1.infot = 4;
  504. dtgsja_("N", "N", "N", &c_n1, &c__0, &c__0, &dummyk, &dummyl, a, &
  505. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  506. c__1, w, &ncycle, &info);
  507. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  508. infoc_1.ok);
  509. infoc_1.infot = 5;
  510. dtgsja_("N", "N", "N", &c__0, &c_n1, &c__0, &dummyk, &dummyl, a, &
  511. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  512. c__1, w, &ncycle, &info);
  513. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  514. infoc_1.ok);
  515. infoc_1.infot = 6;
  516. dtgsja_("N", "N", "N", &c__0, &c__0, &c_n1, &dummyk, &dummyl, a, &
  517. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  518. c__1, w, &ncycle, &info);
  519. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  520. infoc_1.ok);
  521. infoc_1.infot = 10;
  522. dtgsja_("N", "N", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  523. c__0, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  524. c__1, w, &ncycle, &info);
  525. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  526. infoc_1.ok);
  527. infoc_1.infot = 12;
  528. dtgsja_("N", "N", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  529. c__1, b, &c__0, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  530. c__1, w, &ncycle, &info);
  531. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  532. infoc_1.ok);
  533. infoc_1.infot = 18;
  534. dtgsja_("U", "N", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  535. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__0, v, &c__1, q, &
  536. c__1, w, &ncycle, &info);
  537. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  538. infoc_1.ok);
  539. infoc_1.infot = 20;
  540. dtgsja_("N", "V", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  541. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__0, q, &
  542. c__1, w, &ncycle, &info);
  543. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  544. infoc_1.ok);
  545. infoc_1.infot = 22;
  546. dtgsja_("N", "N", "Q", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  547. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  548. c__0, w, &ncycle, &info);
  549. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  550. infoc_1.ok);
  551. nt += 11;
  552. /* Test error exits for the GLM path. */
  553. } else if (lsamen_(&c__3, path, "GLM")) {
  554. /* DGGGLM */
  555. s_copy(srnamc_1.srnamt, "DGGGLM", (ftnlen)32, (ftnlen)6);
  556. infoc_1.infot = 1;
  557. dggglm_(&c_n1, &c__0, &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 = 2;
  562. dggglm_(&c__0, &c_n1, &c__0, 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 = 2;
  567. dggglm_(&c__0, &c__1, &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 = 3;
  572. dggglm_(&c__0, &c__0, &c_n1, a, &c__1, 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 = 3;
  577. dggglm_(&c__1, &c__0, &c__0, a, &c__1, b, &c__1, 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 = 5;
  582. dggglm_(&c__0, &c__0, &c__0, a, &c__0, b, &c__1, r1, r2, r3, w, &
  583. c__18, &info);
  584. chkxer_("DGGGLM", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  585. infoc_1.ok);
  586. infoc_1.infot = 7;
  587. dggglm_(&c__0, &c__0, &c__0, a, &c__1, b, &c__0, r1, r2, r3, w, &
  588. c__18, &info);
  589. chkxer_("DGGGLM", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  590. infoc_1.ok);
  591. infoc_1.infot = 12;
  592. dggglm_(&c__1, &c__1, &c__1, a, &c__1, b, &c__1, r1, r2, r3, w, &c__1,
  593. &info);
  594. chkxer_("DGGGLM", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  595. infoc_1.ok);
  596. nt += 8;
  597. /* Test error exits for the LSE path. */
  598. } else if (lsamen_(&c__3, path, "LSE")) {
  599. /* DGGLSE */
  600. s_copy(srnamc_1.srnamt, "DGGLSE", (ftnlen)32, (ftnlen)6);
  601. infoc_1.infot = 1;
  602. dgglse_(&c_n1, &c__0, &c__0, 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 = 2;
  607. dgglse_(&c__0, &c_n1, &c__0, 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__0, &c_n1, 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 = 3;
  617. dgglse_(&c__0, &c__0, &c__1, a, &c__1, 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 = 3;
  622. dgglse_(&c__0, &c__1, &c__0, a, &c__1, b, &c__1, 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 = 5;
  627. dgglse_(&c__0, &c__0, &c__0, a, &c__0, b, &c__1, r1, r2, r3, w, &
  628. c__18, &info);
  629. chkxer_("DGGLSE", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  630. infoc_1.ok);
  631. infoc_1.infot = 7;
  632. dgglse_(&c__0, &c__0, &c__0, a, &c__1, b, &c__0, r1, r2, r3, w, &
  633. c__18, &info);
  634. chkxer_("DGGLSE", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  635. infoc_1.ok);
  636. infoc_1.infot = 12;
  637. dgglse_(&c__1, &c__1, &c__1, a, &c__1, b, &c__1, r1, r2, r3, w, &c__1,
  638. &info);
  639. chkxer_("DGGLSE", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  640. infoc_1.ok);
  641. nt += 8;
  642. /* Test error exits for the GQR path. */
  643. } else if (lsamen_(&c__3, path, "GQR")) {
  644. /* DGGQRF */
  645. s_copy(srnamc_1.srnamt, "DGGQRF", (ftnlen)32, (ftnlen)6);
  646. infoc_1.infot = 1;
  647. dggqrf_(&c_n1, &c__0, &c__0, 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 = 2;
  652. dggqrf_(&c__0, &c_n1, &c__0, a, &c__1, 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 = 3;
  657. dggqrf_(&c__0, &c__0, &c_n1, a, &c__1, r1, b, &c__1, 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 = 5;
  662. dggqrf_(&c__0, &c__0, &c__0, a, &c__0, r1, b, &c__1, r2, w, &c__18, &
  663. info);
  664. chkxer_("DGGQRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  665. infoc_1.ok);
  666. infoc_1.infot = 8;
  667. dggqrf_(&c__0, &c__0, &c__0, a, &c__1, r1, b, &c__0, r2, w, &c__18, &
  668. info);
  669. chkxer_("DGGQRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  670. infoc_1.ok);
  671. infoc_1.infot = 11;
  672. dggqrf_(&c__1, &c__1, &c__2, a, &c__1, r1, b, &c__1, r2, w, &c__1, &
  673. info);
  674. chkxer_("DGGQRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  675. infoc_1.ok);
  676. nt += 6;
  677. /* DGGRQF */
  678. s_copy(srnamc_1.srnamt, "DGGRQF", (ftnlen)32, (ftnlen)6);
  679. infoc_1.infot = 1;
  680. dggrqf_(&c_n1, &c__0, &c__0, 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 = 2;
  685. dggrqf_(&c__0, &c_n1, &c__0, a, &c__1, 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 = 3;
  690. dggrqf_(&c__0, &c__0, &c_n1, a, &c__1, r1, b, &c__1, 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 = 5;
  695. dggrqf_(&c__0, &c__0, &c__0, a, &c__0, r1, b, &c__1, r2, w, &c__18, &
  696. info);
  697. chkxer_("DGGRQF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  698. infoc_1.ok);
  699. infoc_1.infot = 8;
  700. dggrqf_(&c__0, &c__0, &c__0, a, &c__1, r1, b, &c__0, r2, w, &c__18, &
  701. info);
  702. chkxer_("DGGRQF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  703. infoc_1.ok);
  704. infoc_1.infot = 11;
  705. dggrqf_(&c__1, &c__1, &c__2, a, &c__1, r1, b, &c__1, r2, w, &c__1, &
  706. info);
  707. chkxer_("DGGRQF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  708. infoc_1.ok);
  709. nt += 6;
  710. /* Test error exits for the DGS, DGV, DGX, and DXV paths. */
  711. } else if (lsamen_(&c__3, path, "DGS") || lsamen_(&
  712. c__3, path, "DGV") || lsamen_(&c__3, path,
  713. "DGX") || lsamen_(&c__3, path, "DXV")) {
  714. /* DGGES */
  715. s_copy(srnamc_1.srnamt, "DGGES ", (ftnlen)32, (ftnlen)6);
  716. infoc_1.infot = 1;
  717. dgges_("/", "N", "S", (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 = 2;
  722. dgges_("N", "/", "S", (L_fp)dlctes_, &c__1, 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 = 3;
  727. dgges_("N", "V", "/", (L_fp)dlctes_, &c__1, a, &c__1, 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 = 5;
  732. dgges_("N", "V", "S", (L_fp)dlctes_, &c_n1, a, &c__1, b, &c__1, &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 = 7;
  737. dgges_("N", "V", "S", (L_fp)dlctes_, &c__1, a, &c__0, b, &c__1, &sdim,
  738. r1, r2, r3, q, &c__1, 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 = 9;
  742. dgges_("N", "V", "S", (L_fp)dlctes_, &c__1, a, &c__1, b, &c__0, &sdim,
  743. r1, r2, r3, q, &c__1, u, &c__1, 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 = 15;
  747. dgges_("N", "V", "S", (L_fp)dlctes_, &c__1, a, &c__1, b, &c__1, &sdim,
  748. r1, r2, r3, q, &c__0, u, &c__1, 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 = 15;
  752. dgges_("V", "V", "S", (L_fp)dlctes_, &c__2, a, &c__2, b, &c__2, &sdim,
  753. r1, r2, r3, q, &c__1, u, &c__2, 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 = 17;
  757. dgges_("N", "V", "S", (L_fp)dlctes_, &c__1, a, &c__1, b, &c__1, &sdim,
  758. r1, r2, r3, q, &c__1, u, &c__0, w, &c__1, bw, &info);
  759. chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  760. infoc_1.ok);
  761. infoc_1.infot = 17;
  762. dgges_("V", "V", "S", (L_fp)dlctes_, &c__2, a, &c__2, b, &c__2, &sdim,
  763. r1, r2, r3, q, &c__2, u, &c__1, w, &c__1, bw, &info);
  764. chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  765. infoc_1.ok);
  766. infoc_1.infot = 19;
  767. dgges_("V", "V", "S", (L_fp)dlctes_, &c__2, a, &c__2, b, &c__2, &sdim,
  768. r1, r2, r3, q, &c__2, u, &c__2, w, &c__1, bw, &info);
  769. chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  770. infoc_1.ok);
  771. nt += 11;
  772. /* DGGESX */
  773. s_copy(srnamc_1.srnamt, "DGGESX", (ftnlen)32, (ftnlen)6);
  774. infoc_1.infot = 1;
  775. dggesx_("/", "N", "S", (L_fp)dlctsx_, "N", &c__1, a, &c__1, b, &c__1,
  776. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
  777. &c__1, bw, &info)
  778. ;
  779. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  780. infoc_1.ok);
  781. infoc_1.infot = 2;
  782. dggesx_("N", "/", "S", (L_fp)dlctsx_, "N", &c__1, a, &c__1, b, &c__1,
  783. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
  784. &c__1, bw, &info)
  785. ;
  786. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  787. infoc_1.ok);
  788. infoc_1.infot = 3;
  789. dggesx_("V", "V", "/", (L_fp)dlctsx_, "N", &c__1, a, &c__1, b, &c__1,
  790. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
  791. &c__1, bw, &info)
  792. ;
  793. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  794. infoc_1.ok);
  795. infoc_1.infot = 5;
  796. dggesx_("V", "V", "S", (L_fp)dlctsx_, "/", &c__1, a, &c__1, b, &c__1,
  797. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
  798. &c__1, bw, &info)
  799. ;
  800. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  801. infoc_1.ok);
  802. infoc_1.infot = 6;
  803. dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c_n1, a, &c__1, b, &c__1,
  804. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
  805. &c__1, bw, &info)
  806. ;
  807. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  808. infoc_1.ok);
  809. infoc_1.infot = 8;
  810. dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__1, a, &c__0, b, &c__1,
  811. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
  812. &c__1, bw, &info)
  813. ;
  814. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  815. infoc_1.ok);
  816. infoc_1.infot = 10;
  817. dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__1, a, &c__1, b, &c__0,
  818. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
  819. &c__1, bw, &info)
  820. ;
  821. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  822. infoc_1.ok);
  823. infoc_1.infot = 16;
  824. dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__1, a, &c__1, b, &c__1,
  825. &sdim, r1, r2, r3, q, &c__0, u, &c__1, rce, rcv, w, &c__1, iw,
  826. &c__1, bw, &info)
  827. ;
  828. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  829. infoc_1.ok);
  830. infoc_1.infot = 16;
  831. dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__2, a, &c__2, b, &c__2,
  832. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
  833. &c__1, bw, &info)
  834. ;
  835. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  836. infoc_1.ok);
  837. infoc_1.infot = 18;
  838. dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__1, a, &c__1, b, &c__1,
  839. &sdim, r1, r2, r3, q, &c__1, u, &c__0, rce, rcv, w, &c__1, iw,
  840. &c__1, bw, &info)
  841. ;
  842. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  843. infoc_1.ok);
  844. infoc_1.infot = 18;
  845. dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__2, a, &c__2, b, &c__2,
  846. &sdim, r1, r2, r3, q, &c__2, u, &c__1, rce, rcv, w, &c__1, iw,
  847. &c__1, bw, &info)
  848. ;
  849. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  850. infoc_1.ok);
  851. infoc_1.infot = 22;
  852. dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__2, a, &c__2, b, &c__2,
  853. &sdim, r1, r2, r3, q, &c__2, u, &c__2, rce, rcv, w, &c__1, iw,
  854. &c__1, bw, &info)
  855. ;
  856. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  857. infoc_1.ok);
  858. infoc_1.infot = 24;
  859. dggesx_("V", "V", "S", (L_fp)dlctsx_, "V", &c__1, a, &c__1, b, &c__1,
  860. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__32,
  861. iw, &c__0, bw, &info);
  862. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  863. infoc_1.ok);
  864. nt += 13;
  865. /* DGGEV */
  866. s_copy(srnamc_1.srnamt, "DGGEV ", (ftnlen)32, (ftnlen)6);
  867. infoc_1.infot = 1;
  868. dggev_("/", "N", &c__1, 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 = 2;
  873. dggev_("N", "/", &c__1, a, &c__1, 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 = 3;
  878. dggev_("V", "V", &c_n1, a, &c__1, b, &c__1, 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 = 5;
  883. dggev_("V", "V", &c__1, a, &c__0, b, &c__1, r1, r2, r3, q, &c__1, 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 = 7;
  888. dggev_("V", "V", &c__1, a, &c__1, b, &c__0, r1, r2, r3, q, &c__1, u, &
  889. c__1, w, &c__1, &info);
  890. chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  891. infoc_1.ok);
  892. infoc_1.infot = 12;
  893. dggev_("N", "V", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q, &c__0, u, &
  894. c__1, w, &c__1, &info);
  895. chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  896. infoc_1.ok);
  897. infoc_1.infot = 12;
  898. dggev_("V", "V", &c__2, a, &c__2, b, &c__2, r1, r2, r3, q, &c__1, u, &
  899. c__2, w, &c__1, &info);
  900. chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  901. infoc_1.ok);
  902. infoc_1.infot = 14;
  903. dggev_("V", "N", &c__2, a, &c__2, b, &c__2, r1, r2, r3, q, &c__2, u, &
  904. c__0, w, &c__1, &info);
  905. chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  906. infoc_1.ok);
  907. infoc_1.infot = 14;
  908. dggev_("V", "V", &c__2, a, &c__2, b, &c__2, r1, r2, r3, q, &c__2, u, &
  909. c__1, w, &c__1, &info);
  910. chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  911. infoc_1.ok);
  912. infoc_1.infot = 16;
  913. dggev_("V", "V", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q, &c__1, u, &
  914. c__1, w, &c__1, &info);
  915. chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  916. infoc_1.ok);
  917. nt += 10;
  918. /* DGGEVX */
  919. s_copy(srnamc_1.srnamt, "DGGEVX", (ftnlen)32, (ftnlen)6);
  920. infoc_1.infot = 1;
  921. dggevx_("/", "N", "N", "N", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q,
  922. &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  923. w, &c__1, iw, bw, &info);
  924. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  925. infoc_1.ok);
  926. infoc_1.infot = 2;
  927. dggevx_("N", "/", "N", "N", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q,
  928. &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  929. w, &c__1, iw, bw, &info);
  930. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  931. infoc_1.ok);
  932. infoc_1.infot = 3;
  933. dggevx_("N", "N", "/", "N", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q,
  934. &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  935. w, &c__1, iw, bw, &info);
  936. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  937. infoc_1.ok);
  938. infoc_1.infot = 4;
  939. dggevx_("N", "N", "N", "/", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q,
  940. &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  941. w, &c__1, iw, bw, &info);
  942. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  943. infoc_1.ok);
  944. infoc_1.infot = 5;
  945. dggevx_("N", "N", "N", "N", &c_n1, a, &c__1, b, &c__1, r1, r2, r3, q,
  946. &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  947. w, &c__1, iw, bw, &info);
  948. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  949. infoc_1.ok);
  950. infoc_1.infot = 7;
  951. dggevx_("N", "N", "N", "N", &c__1, a, &c__0, b, &c__1, r1, r2, r3, q,
  952. &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  953. w, &c__1, iw, bw, &info);
  954. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  955. infoc_1.ok);
  956. infoc_1.infot = 9;
  957. dggevx_("N", "N", "N", "N", &c__1, a, &c__1, b, &c__0, r1, r2, r3, q,
  958. &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  959. w, &c__1, iw, bw, &info);
  960. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  961. infoc_1.ok);
  962. infoc_1.infot = 14;
  963. dggevx_("N", "N", "N", "N", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q,
  964. &c__0, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  965. w, &c__1, iw, bw, &info);
  966. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  967. infoc_1.ok);
  968. infoc_1.infot = 14;
  969. dggevx_("N", "V", "N", "N", &c__2, a, &c__2, b, &c__2, r1, r2, r3, q,
  970. &c__1, u, &c__2, &c__1, &c__2, ls, rs, &anrm, &bnrm, rce, rcv,
  971. w, &c__1, iw, bw, &info);
  972. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  973. infoc_1.ok);
  974. infoc_1.infot = 16;
  975. dggevx_("N", "N", "N", "N", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q,
  976. &c__1, u, &c__0, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  977. w, &c__1, iw, bw, &info);
  978. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  979. infoc_1.ok);
  980. infoc_1.infot = 16;
  981. dggevx_("N", "N", "V", "N", &c__2, a, &c__2, b, &c__2, r1, r2, r3, q,
  982. &c__2, u, &c__1, &c__1, &c__2, ls, rs, &anrm, &bnrm, rce, rcv,
  983. w, &c__1, iw, bw, &info);
  984. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  985. infoc_1.ok);
  986. infoc_1.infot = 26;
  987. dggevx_("N", "N", "V", "N", &c__2, a, &c__2, b, &c__2, r1, r2, r3, q,
  988. &c__2, u, &c__2, &c__1, &c__2, ls, rs, &anrm, &bnrm, rce, rcv,
  989. w, &c__1, iw, bw, &info);
  990. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  991. infoc_1.ok);
  992. nt += 12;
  993. /* DTGEXC */
  994. s_copy(srnamc_1.srnamt, "DTGEXC", (ftnlen)32, (ftnlen)6);
  995. infoc_1.infot = 3;
  996. dtgexc_(&c_true, &c_true, &c_n1, a, &c__1, b, &c__1, 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 = 5;
  1001. dtgexc_(&c_true, &c_true, &c__1, a, &c__0, b, &c__1, q, &c__1, 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 = 7;
  1006. dtgexc_(&c_true, &c_true, &c__1, a, &c__1, b, &c__0, q, &c__1, 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 = 9;
  1011. dtgexc_(&c_false, &c_true, &c__1, a, &c__1, b, &c__1, q, &c__0, z__, &
  1012. c__1, &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 = 9;
  1016. dtgexc_(&c_true, &c_true, &c__1, a, &c__1, b, &c__1, q, &c__0, z__, &
  1017. c__1, &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 = 11;
  1021. dtgexc_(&c_true, &c_false, &c__1, a, &c__1, b, &c__1, q, &c__1, z__, &
  1022. c__0, &ifst, &ilst, w, &c__1, &info);
  1023. chkxer_("DTGEXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1024. infoc_1.ok);
  1025. infoc_1.infot = 11;
  1026. dtgexc_(&c_true, &c_true, &c__1, a, &c__1, b, &c__1, q, &c__1, z__, &
  1027. c__0, &ifst, &ilst, w, &c__1, &info);
  1028. chkxer_("DTGEXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1029. infoc_1.ok);
  1030. infoc_1.infot = 15;
  1031. dtgexc_(&c_true, &c_true, &c__1, a, &c__1, b, &c__1, q, &c__1, z__, &
  1032. c__1, &ifst, &ilst, w, &c__0, &info);
  1033. chkxer_("DTGEXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1034. infoc_1.ok);
  1035. nt += 8;
  1036. /* DTGSEN */
  1037. s_copy(srnamc_1.srnamt, "DTGSEN", (ftnlen)32, (ftnlen)6);
  1038. infoc_1.infot = 1;
  1039. dtgsen_(&c_n1, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1040. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
  1041. iw, &c__1, &info);
  1042. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1043. infoc_1.ok);
  1044. infoc_1.infot = 5;
  1045. dtgsen_(&c__1, &c_true, &c_true, sel, &c_n1, a, &c__1, b, &c__1, r1,
  1046. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
  1047. iw, &c__1, &info);
  1048. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1049. infoc_1.ok);
  1050. infoc_1.infot = 7;
  1051. dtgsen_(&c__1, &c_true, &c_true, sel, &c__1, a, &c__0, b, &c__1, r1,
  1052. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
  1053. iw, &c__1, &info);
  1054. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1055. infoc_1.ok);
  1056. infoc_1.infot = 9;
  1057. dtgsen_(&c__1, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__0, r1,
  1058. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
  1059. iw, &c__1, &info);
  1060. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1061. infoc_1.ok);
  1062. infoc_1.infot = 14;
  1063. dtgsen_(&c__1, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1064. r2, r3, q, &c__0, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
  1065. iw, &c__1, &info);
  1066. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1067. infoc_1.ok);
  1068. infoc_1.infot = 16;
  1069. dtgsen_(&c__1, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1070. r2, r3, q, &c__1, z__, &c__0, &m, &tola, &tolb, rcv, w, &c__1,
  1071. iw, &c__1, &info);
  1072. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1073. infoc_1.ok);
  1074. infoc_1.infot = 22;
  1075. dtgsen_(&c__0, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1076. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
  1077. iw, &c__1, &info);
  1078. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1079. infoc_1.ok);
  1080. infoc_1.infot = 22;
  1081. dtgsen_(&c__1, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1082. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
  1083. iw, &c__1, &info);
  1084. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1085. infoc_1.ok);
  1086. infoc_1.infot = 22;
  1087. dtgsen_(&c__2, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1088. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
  1089. iw, &c__1, &info);
  1090. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1091. infoc_1.ok);
  1092. infoc_1.infot = 24;
  1093. dtgsen_(&c__0, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1094. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &
  1095. c__20, iw, &c__0, &info);
  1096. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1097. infoc_1.ok);
  1098. infoc_1.infot = 24;
  1099. dtgsen_(&c__1, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1100. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &
  1101. c__20, iw, &c__0, &info);
  1102. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1103. infoc_1.ok);
  1104. infoc_1.infot = 24;
  1105. dtgsen_(&c__2, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1106. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &
  1107. c__20, iw, &c__1, &info);
  1108. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1109. infoc_1.ok);
  1110. nt += 12;
  1111. /* DTGSNA */
  1112. s_copy(srnamc_1.srnamt, "DTGSNA", (ftnlen)32, (ftnlen)6);
  1113. infoc_1.infot = 1;
  1114. dtgsna_("/", "A", sel, &c__1, 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 = 2;
  1119. dtgsna_("B", "/", sel, &c__1, a, &c__1, 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 = 4;
  1124. dtgsna_("B", "A", sel, &c_n1, a, &c__1, b, &c__1, 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 = 6;
  1129. dtgsna_("B", "A", sel, &c__1, a, &c__0, b, &c__1, q, &c__1, 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 = 8;
  1134. dtgsna_("B", "A", sel, &c__1, a, &c__1, b, &c__0, q, &c__1, u, &c__1,
  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 = 10;
  1139. dtgsna_("E", "A", sel, &c__1, a, &c__1, b, &c__1, q, &c__0, u, &c__1,
  1140. r1, r2, &c__1, &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 = 12;
  1144. dtgsna_("E", "A", sel, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &c__0,
  1145. r1, r2, &c__1, &m, w, &c__1, iw, &info);
  1146. chkxer_("DTGSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1147. infoc_1.ok);
  1148. infoc_1.infot = 15;
  1149. dtgsna_("E", "A", sel, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &c__1,
  1150. r1, r2, &c__0, &m, w, &c__1, iw, &info);
  1151. chkxer_("DTGSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1152. infoc_1.ok);
  1153. infoc_1.infot = 18;
  1154. dtgsna_("E", "A", sel, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &c__1,
  1155. r1, r2, &c__1, &m, w, &c__0, iw, &info);
  1156. chkxer_("DTGSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1157. infoc_1.ok);
  1158. nt += 9;
  1159. /* DTGSYL */
  1160. s_copy(srnamc_1.srnamt, "DTGSYL", (ftnlen)32, (ftnlen)6);
  1161. infoc_1.infot = 1;
  1162. dtgsyl_("/", &c__0, &c__1, &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 = 2;
  1167. dtgsyl_("N", &c_n1, &c__1, &c__1, 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 = 3;
  1172. dtgsyl_("N", &c__0, &c__0, &c__1, a, &c__1, 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 = 4;
  1177. dtgsyl_("N", &c__0, &c__1, &c__0, a, &c__1, b, &c__1, 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 = 6;
  1182. dtgsyl_("N", &c__0, &c__1, &c__1, a, &c__0, b, &c__1, q, &c__1, 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 = 8;
  1187. dtgsyl_("N", &c__0, &c__1, &c__1, a, &c__1, b, &c__0, q, &c__1, u, &
  1188. c__1, 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 = 10;
  1192. dtgsyl_("N", &c__0, &c__1, &c__1, a, &c__1, b, &c__1, q, &c__0, u, &
  1193. c__1, v, &c__1, z__, &c__1, &scale, &dif, w,

Large files files are truncated, but you can click here to view the full file