PageRenderTime 31ms CodeModel.GetById 21ms RepoModel.GetById 0ms app.codeStats 1ms

/thirdPartyCode/CLAPACK/TESTING/EIG/derrgg.c

https://bitbucket.org/rutad/libs
C | 1310 lines | 1169 code | 78 blank | 63 comment | 19 complexity | 75646729505c8bcd8133c53f0323da2b MD5 | raw file
Possible License(s): BSD-3-Clause, MPL-2.0-no-copyleft-exception, GPL-3.0

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

  1. #include "blaswrap.h"
  2. /* -- translated by f2c (version 19990503).
  3. You must link the resulting object file with the libraries:
  4. -lf2c -lm (in that order)
  5. */
  6. #include "f2c.h"
  7. /* Common Block Declarations */
  8. struct {
  9. integer infot, nout;
  10. logical ok, lerr;
  11. } infoc_;
  12. #define infoc_1 infoc_
  13. struct {
  14. char srnamt[6];
  15. } srnamc_;
  16. #define srnamc_1 srnamc_
  17. /* Table of constant values */
  18. static integer c__2 = 2;
  19. static integer c__0 = 0;
  20. static integer c__1 = 1;
  21. static integer c_n1 = -1;
  22. static integer c__18 = 18;
  23. static integer c__3 = 3;
  24. static integer c__32 = 32;
  25. static logical c_true = TRUE_;
  26. static logical c_false = FALSE_;
  27. static integer c__20 = 20;
  28. /* Subroutine */ int derrgg_(char *path, integer *nunit)
  29. {
  30. /* Format strings */
  31. static char fmt_9999[] = "(1x,a3,\002 routines passed the tests of the e"
  32. "rror exits (\002,i3,\002 tests done)\002)";
  33. static char fmt_9998[] = "(\002 *** \002,a3,\002 routines failed the tes"
  34. "ts of the error \002,\002exits ***\002)";
  35. /* Builtin functions */
  36. integer s_wsle(cilist *), e_wsle(void);
  37. /* Subroutine */ int s_copy(char *, char *, ftnlen, ftnlen);
  38. integer s_wsfe(cilist *), do_fio(integer *, char *, ftnlen), e_wsfe(void);
  39. /* Local variables */
  40. static integer info, sdim;
  41. static doublereal anrm, bnrm, tola, tolb;
  42. static integer ifst, ilst;
  43. static doublereal a[9] /* was [3][3] */, b[9] /* was [3][3] */;
  44. static integer i__, j, m;
  45. static doublereal q[9] /* was [3][3] */, u[9] /* was [3][3] */, v[9]
  46. /* was [3][3] */, scale, w[18];
  47. extern /* Subroutine */ int dgges_(char *, char *, char *, L_fp, integer *
  48. , doublereal *, integer *, doublereal *, integer *, integer *,
  49. doublereal *, doublereal *, doublereal *, doublereal *, integer *,
  50. doublereal *, integer *, doublereal *, integer *, logical *,
  51. integer *);
  52. static doublereal z__[9] /* was [3][3] */;
  53. extern /* Subroutine */ int dggev_(char *, char *, integer *, doublereal *
  54. , integer *, doublereal *, integer *, doublereal *, doublereal *,
  55. doublereal *, doublereal *, integer *, doublereal *, integer *,
  56. doublereal *, integer *, integer *);
  57. static char c2[2];
  58. static doublereal r1[3], r2[3], r3[3];
  59. static logical bw[3];
  60. static doublereal ls[3];
  61. static integer iw[3];
  62. extern /* Subroutine */ int dggglm_(integer *, integer *, integer *,
  63. doublereal *, integer *, doublereal *, integer *, doublereal *,
  64. doublereal *, doublereal *, doublereal *, integer *, integer *);
  65. static integer nt;
  66. extern /* Subroutine */ int dgghrd_(char *, char *, integer *, integer *,
  67. integer *, doublereal *, integer *, doublereal *, integer *,
  68. doublereal *, integer *, doublereal *, integer *, integer *);
  69. static doublereal rs[3];
  70. extern /* Subroutine */ int dgglse_(integer *, integer *, integer *,
  71. doublereal *, integer *, doublereal *, integer *, doublereal *,
  72. doublereal *, doublereal *, doublereal *, integer *, integer *),
  73. dggqrf_(integer *, integer *, integer *, doublereal *, integer *,
  74. doublereal *, doublereal *, integer *, doublereal *, doublereal *,
  75. integer *, integer *), dggrqf_(integer *, integer *, integer *,
  76. doublereal *, integer *, doublereal *, doublereal *, integer *,
  77. doublereal *, doublereal *, integer *, integer *);
  78. static integer ncycle;
  79. extern logical dlctes_(), lsamen_(integer *, char *, char *);
  80. extern /* Subroutine */ int dggsvd_(char *, char *, char *, integer *,
  81. integer *, integer *, integer *, integer *, doublereal *, integer
  82. *, doublereal *, integer *, doublereal *, doublereal *,
  83. doublereal *, integer *, doublereal *, integer *, doublereal *,
  84. integer *, doublereal *, integer *, integer *), dggesx_(char *, char *, char *, L_fp, char *, integer *,
  85. doublereal *, integer *, doublereal *, integer *, integer *,
  86. doublereal *, doublereal *, doublereal *, doublereal *, integer *,
  87. doublereal *, integer *, doublereal *, doublereal *, doublereal *
  88. , integer *, integer *, integer *, logical *, integer *), dhgeqz_(char *, char *, char *, integer *
  89. , integer *, integer *, doublereal *, integer *, doublereal *,
  90. integer *, doublereal *, doublereal *, doublereal *, doublereal *,
  91. integer *, doublereal *, integer *, doublereal *, integer *,
  92. integer *), dtgevc_(char *, char *,
  93. logical *, integer *, doublereal *, integer *, doublereal *,
  94. integer *, doublereal *, integer *, doublereal *, integer *,
  95. integer *, integer *, doublereal *, integer *),
  96. chkxer_(char *, integer *, integer *, logical *, logical *), dggevx_(char *, char *, char *, char *, integer *,
  97. doublereal *, integer *, doublereal *, integer *, doublereal *,
  98. doublereal *, doublereal *, doublereal *, integer *, doublereal *,
  99. integer *, integer *, integer *, doublereal *, doublereal *,
  100. doublereal *, doublereal *, doublereal *, doublereal *,
  101. doublereal *, integer *, integer *, logical *, integer *), dtgexc_(logical *, logical *, integer *,
  102. doublereal *, integer *, doublereal *, integer *, doublereal *,
  103. integer *, doublereal *, integer *, integer *, integer *,
  104. doublereal *, integer *, integer *), dtgsen_(integer *, logical *,
  105. logical *, logical *, integer *, doublereal *, integer *,
  106. doublereal *, integer *, doublereal *, doublereal *, doublereal *,
  107. doublereal *, integer *, doublereal *, integer *, integer *,
  108. doublereal *, doublereal *, doublereal *, doublereal *, integer *,
  109. integer *, integer *, integer *), dtgsja_(char *, char *, char *,
  110. integer *, integer *, integer *, integer *, integer *,
  111. doublereal *, integer *, doublereal *, integer *, doublereal *,
  112. doublereal *, doublereal *, doublereal *, doublereal *, integer *,
  113. doublereal *, integer *, doublereal *, integer *, doublereal *,
  114. integer *, integer *), dtgsna_(char *,
  115. char *, logical *, integer *, doublereal *, integer *, doublereal
  116. *, integer *, doublereal *, integer *, doublereal *, integer *,
  117. doublereal *, doublereal *, integer *, integer *, doublereal *,
  118. integer *, integer *, integer *), dggsvp_(char *,
  119. char *, char *, integer *, integer *, integer *, doublereal *,
  120. integer *, doublereal *, integer *, doublereal *, doublereal *,
  121. integer *, integer *, doublereal *, integer *, doublereal *,
  122. integer *, doublereal *, integer *, integer *, doublereal *,
  123. doublereal *, integer *);
  124. extern logical dlctsx_();
  125. static integer dummyk, dummyl;
  126. extern /* Subroutine */ int dtgsyl_(char *, integer *, integer *, integer
  127. *, doublereal *, integer *, doublereal *, integer *, doublereal *,
  128. integer *, doublereal *, integer *, doublereal *, integer *,
  129. doublereal *, integer *, doublereal *, doublereal *, doublereal *,
  130. integer *, integer *, integer *);
  131. static doublereal dif, rce[2];
  132. static logical sel[3];
  133. static doublereal tau[3], rcv[2];
  134. /* Fortran I/O blocks */
  135. static cilist io___1 = { 0, 0, 0, 0, 0 };
  136. static cilist io___38 = { 0, 0, 0, fmt_9999, 0 };
  137. static cilist io___39 = { 0, 0, 0, fmt_9998, 0 };
  138. #define a_ref(a_1,a_2) a[(a_2)*3 + a_1 - 4]
  139. #define b_ref(a_1,a_2) b[(a_2)*3 + a_1 - 4]
  140. /* -- LAPACK test routine (version 3.0) --
  141. Univ. of Tennessee, Univ. of California Berkeley, NAG Ltd.,
  142. Courant Institute, Argonne National Lab, and Rice University
  143. June 30, 1999
  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. infoc_1.nout = *nunit;
  157. io___1.ciunit = infoc_1.nout;
  158. s_wsle(&io___1);
  159. e_wsle();
  160. s_copy(c2, path + 1, (ftnlen)2, (ftnlen)2);
  161. /* Set the variables to innocuous values. */
  162. for (j = 1; j <= 3; ++j) {
  163. sel[j - 1] = TRUE_;
  164. for (i__ = 1; i__ <= 3; ++i__) {
  165. a_ref(i__, j) = 0.;
  166. b_ref(i__, j) = 0.;
  167. /* L10: */
  168. }
  169. /* L20: */
  170. }
  171. for (i__ = 1; i__ <= 3; ++i__) {
  172. a_ref(i__, i__) = 1.;
  173. b_ref(i__, i__) = 1.;
  174. /* L30: */
  175. }
  176. infoc_1.ok = TRUE_;
  177. tola = 1.;
  178. tolb = 1.;
  179. ifst = 1;
  180. ilst = 1;
  181. nt = 0;
  182. /* Test error exits for the GG path. */
  183. if (lsamen_(&c__2, c2, "GG")) {
  184. /* DGGHRD */
  185. s_copy(srnamc_1.srnamt, "DGGHRD", (ftnlen)6, (ftnlen)6);
  186. infoc_1.infot = 1;
  187. dgghrd_("/", "N", &c__0, &c__1, &c__0, a, &c__1, b, &c__1, q, &c__1,
  188. z__, &c__1, &info);
  189. chkxer_("DGGHRD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  190. infoc_1.ok);
  191. infoc_1.infot = 2;
  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 = 3;
  197. dgghrd_("N", "N", &c_n1, &c__0, &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 = 4;
  202. dgghrd_("N", "N", &c__0, &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 = 5;
  207. dgghrd_("N", "N", &c__0, &c__1, &c__1, 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 = 7;
  212. dgghrd_("N", "N", &c__2, &c__1, &c__1, a, &c__1, b, &c__2, 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 = 9;
  217. dgghrd_("N", "N", &c__2, &c__1, &c__1, a, &c__2, 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 = 11;
  222. dgghrd_("V", "N", &c__2, &c__1, &c__1, a, &c__2, b, &c__2, 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 = 13;
  227. dgghrd_("N", "V", &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. nt += 9;
  232. /* DHGEQZ */
  233. s_copy(srnamc_1.srnamt, "DHGEQZ", (ftnlen)6, (ftnlen)6);
  234. infoc_1.infot = 1;
  235. dhgeqz_("/", "N", "N", &c__0, &c__1, &c__0, a, &c__1, b, &c__1, r1,
  236. r2, r3, q, &c__1, z__, &c__1, w, &c__18, &info);
  237. chkxer_("DHGEQZ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  238. infoc_1.ok);
  239. infoc_1.infot = 2;
  240. dhgeqz_("E", "/", "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 = 3;
  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 = 4;
  250. dhgeqz_("E", "N", "N", &c_n1, &c__0, &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 = 5;
  255. dhgeqz_("E", "N", "N", &c__0, &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 = 6;
  260. dhgeqz_("E", "N", "N", &c__0, &c__1, &c__1, 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 = 8;
  265. dhgeqz_("E", "N", "N", &c__2, &c__1, &c__1, a, &c__1, b, &c__2, 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 = 10;
  270. dhgeqz_("E", "N", "N", &c__2, &c__1, &c__1, a, &c__2, 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 = 15;
  275. dhgeqz_("E", "V", "N", &c__2, &c__1, &c__1, a, &c__2, b, &c__2, 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 = 17;
  280. dhgeqz_("E", "N", "V", &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. nt += 10;
  285. /* DTGEVC */
  286. s_copy(srnamc_1.srnamt, "DTGEVC", (ftnlen)6, (ftnlen)6);
  287. infoc_1.infot = 1;
  288. dtgevc_("/", "A", sel, &c__0, a, &c__1, b, &c__1, q, &c__1, z__, &
  289. c__1, &c__0, &m, w, &info);
  290. chkxer_("DTGEVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  291. infoc_1.ok);
  292. infoc_1.infot = 2;
  293. dtgevc_("R", "/", 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 = 4;
  298. dtgevc_("R", "A", sel, &c_n1, 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 = 6;
  303. dtgevc_("R", "A", sel, &c__2, a, &c__1, b, &c__2, q, &c__1, z__, &
  304. c__2, &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 = 8;
  308. dtgevc_("R", "A", sel, &c__2, a, &c__2, b, &c__1, 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 = 10;
  313. dtgevc_("L", "A", sel, &c__2, a, &c__2, b, &c__2, 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 = 12;
  318. dtgevc_("R", "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 = 13;
  323. dtgevc_("R", "A", sel, &c__2, a, &c__2, b, &c__2, q, &c__1, z__, &
  324. c__2, &c__1, &m, w, &info);
  325. chkxer_("DTGEVC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  326. infoc_1.ok);
  327. nt += 8;
  328. /* Test error exits for the GSV path. */
  329. } else if (lsamen_(&c__3, path, "GSV")) {
  330. /* DGGSVD */
  331. s_copy(srnamc_1.srnamt, "DGGSVD", (ftnlen)6, (ftnlen)6);
  332. infoc_1.infot = 1;
  333. dggsvd_("/", "N", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  334. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  335. info);
  336. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  337. infoc_1.ok);
  338. infoc_1.infot = 2;
  339. dggsvd_("N", "/", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  340. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  341. info);
  342. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  343. infoc_1.ok);
  344. infoc_1.infot = 3;
  345. dggsvd_("N", "N", "/", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  346. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  347. info);
  348. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  349. infoc_1.ok);
  350. infoc_1.infot = 4;
  351. dggsvd_("N", "N", "N", &c_n1, &c__0, &c__0, &dummyk, &dummyl, a, &
  352. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  353. info);
  354. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  355. infoc_1.ok);
  356. infoc_1.infot = 5;
  357. dggsvd_("N", "N", "N", &c__0, &c_n1, &c__0, &dummyk, &dummyl, a, &
  358. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  359. info);
  360. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  361. infoc_1.ok);
  362. infoc_1.infot = 6;
  363. dggsvd_("N", "N", "N", &c__0, &c__0, &c_n1, &dummyk, &dummyl, a, &
  364. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  365. info);
  366. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  367. infoc_1.ok);
  368. infoc_1.infot = 10;
  369. dggsvd_("N", "N", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  370. c__0, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  371. info);
  372. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  373. infoc_1.ok);
  374. infoc_1.infot = 12;
  375. dggsvd_("N", "N", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  376. c__1, b, &c__0, r1, r2, u, &c__1, v, &c__1, q, &c__1, w, iw, &
  377. info);
  378. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  379. infoc_1.ok);
  380. infoc_1.infot = 16;
  381. dggsvd_("U", "N", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  382. c__1, b, &c__1, r1, r2, u, &c__0, v, &c__1, q, &c__1, w, iw, &
  383. info);
  384. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  385. infoc_1.ok);
  386. infoc_1.infot = 18;
  387. dggsvd_("N", "V", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  388. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__0, q, &c__1, w, iw, &
  389. info);
  390. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  391. infoc_1.ok);
  392. infoc_1.infot = 20;
  393. dggsvd_("N", "N", "Q", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  394. c__1, b, &c__1, r1, r2, u, &c__1, v, &c__1, q, &c__0, w, iw, &
  395. info);
  396. chkxer_("DGGSVD", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  397. infoc_1.ok);
  398. nt += 11;
  399. /* DGGSVP */
  400. s_copy(srnamc_1.srnamt, "DGGSVP", (ftnlen)6, (ftnlen)6);
  401. infoc_1.infot = 1;
  402. dggsvp_("/", "N", "N", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, &tola,
  403. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  404. tau, w, &info);
  405. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  406. infoc_1.ok);
  407. infoc_1.infot = 2;
  408. dggsvp_("N", "/", "N", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, &tola,
  409. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  410. tau, w, &info);
  411. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  412. infoc_1.ok);
  413. infoc_1.infot = 3;
  414. dggsvp_("N", "N", "/", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, &tola,
  415. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  416. tau, w, &info);
  417. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  418. infoc_1.ok);
  419. infoc_1.infot = 4;
  420. dggsvp_("N", "N", "N", &c_n1, &c__0, &c__0, a, &c__1, b, &c__1, &tola,
  421. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  422. tau, w, &info);
  423. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  424. infoc_1.ok);
  425. infoc_1.infot = 5;
  426. dggsvp_("N", "N", "N", &c__0, &c_n1, &c__0, a, &c__1, b, &c__1, &tola,
  427. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  428. tau, w, &info);
  429. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  430. infoc_1.ok);
  431. infoc_1.infot = 6;
  432. dggsvp_("N", "N", "N", &c__0, &c__0, &c_n1, a, &c__1, b, &c__1, &tola,
  433. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  434. tau, w, &info);
  435. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  436. infoc_1.ok);
  437. infoc_1.infot = 8;
  438. dggsvp_("N", "N", "N", &c__0, &c__0, &c__0, a, &c__0, b, &c__1, &tola,
  439. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  440. tau, w, &info);
  441. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  442. infoc_1.ok);
  443. infoc_1.infot = 10;
  444. dggsvp_("N", "N", "N", &c__0, &c__0, &c__0, a, &c__1, b, &c__0, &tola,
  445. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__1, iw,
  446. tau, w, &info);
  447. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  448. infoc_1.ok);
  449. infoc_1.infot = 16;
  450. dggsvp_("U", "N", "N", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, &tola,
  451. &tolb, &dummyk, &dummyl, u, &c__0, v, &c__1, q, &c__1, iw,
  452. tau, w, &info);
  453. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  454. infoc_1.ok);
  455. infoc_1.infot = 18;
  456. dggsvp_("N", "V", "N", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, &tola,
  457. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__0, q, &c__1, iw,
  458. tau, w, &info);
  459. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  460. infoc_1.ok);
  461. infoc_1.infot = 20;
  462. dggsvp_("N", "N", "Q", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, &tola,
  463. &tolb, &dummyk, &dummyl, u, &c__1, v, &c__1, q, &c__0, iw,
  464. tau, w, &info);
  465. chkxer_("DGGSVP", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  466. infoc_1.ok);
  467. nt += 11;
  468. /* DTGSJA */
  469. s_copy(srnamc_1.srnamt, "DTGSJA", (ftnlen)6, (ftnlen)6);
  470. infoc_1.infot = 1;
  471. dtgsja_("/", "N", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  472. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  473. c__1, w, &ncycle, &info);
  474. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  475. infoc_1.ok);
  476. infoc_1.infot = 2;
  477. dtgsja_("N", "/", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  478. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  479. c__1, w, &ncycle, &info);
  480. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  481. infoc_1.ok);
  482. infoc_1.infot = 3;
  483. dtgsja_("N", "N", "/", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  484. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  485. c__1, w, &ncycle, &info);
  486. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  487. infoc_1.ok);
  488. infoc_1.infot = 4;
  489. dtgsja_("N", "N", "N", &c_n1, &c__0, &c__0, &dummyk, &dummyl, a, &
  490. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  491. c__1, w, &ncycle, &info);
  492. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  493. infoc_1.ok);
  494. infoc_1.infot = 5;
  495. dtgsja_("N", "N", "N", &c__0, &c_n1, &c__0, &dummyk, &dummyl, a, &
  496. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  497. c__1, w, &ncycle, &info);
  498. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  499. infoc_1.ok);
  500. infoc_1.infot = 6;
  501. dtgsja_("N", "N", "N", &c__0, &c__0, &c_n1, &dummyk, &dummyl, a, &
  502. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  503. c__1, w, &ncycle, &info);
  504. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  505. infoc_1.ok);
  506. infoc_1.infot = 10;
  507. dtgsja_("N", "N", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  508. c__0, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  509. c__1, w, &ncycle, &info);
  510. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  511. infoc_1.ok);
  512. infoc_1.infot = 12;
  513. dtgsja_("N", "N", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  514. c__1, b, &c__0, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  515. c__1, w, &ncycle, &info);
  516. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  517. infoc_1.ok);
  518. infoc_1.infot = 18;
  519. dtgsja_("U", "N", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  520. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__0, v, &c__1, q, &
  521. c__1, w, &ncycle, &info);
  522. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  523. infoc_1.ok);
  524. infoc_1.infot = 20;
  525. dtgsja_("N", "V", "N", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  526. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__0, q, &
  527. c__1, w, &ncycle, &info);
  528. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  529. infoc_1.ok);
  530. infoc_1.infot = 22;
  531. dtgsja_("N", "N", "Q", &c__0, &c__0, &c__0, &dummyk, &dummyl, a, &
  532. c__1, b, &c__1, &tola, &tolb, r1, r2, u, &c__1, v, &c__1, q, &
  533. c__0, w, &ncycle, &info);
  534. chkxer_("DTGSJA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  535. infoc_1.ok);
  536. nt += 11;
  537. /* Test error exits for the GLM path. */
  538. } else if (lsamen_(&c__3, path, "GLM")) {
  539. /* DGGGLM */
  540. s_copy(srnamc_1.srnamt, "DGGGLM", (ftnlen)6, (ftnlen)6);
  541. infoc_1.infot = 1;
  542. dggglm_(&c_n1, &c__0, &c__0, a, &c__1, b, &c__1, r1, r2, r3, w, &
  543. c__18, &info);
  544. chkxer_("DGGGLM", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  545. infoc_1.ok);
  546. infoc_1.infot = 2;
  547. dggglm_(&c__0, &c_n1, &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__1, &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 = 3;
  557. dggglm_(&c__0, &c__0, &c_n1, 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__1, &c__0, &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 = 5;
  567. dggglm_(&c__0, &c__0, &c__0, a, &c__0, 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 = 7;
  572. dggglm_(&c__0, &c__0, &c__0, a, &c__1, b, &c__0, 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 = 12;
  577. dggglm_(&c__1, &c__1, &c__1, a, &c__1, b, &c__1, r1, r2, r3, w, &c__1,
  578. &info);
  579. chkxer_("DGGGLM", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  580. infoc_1.ok);
  581. nt += 8;
  582. /* Test error exits for the LSE path. */
  583. } else if (lsamen_(&c__3, path, "LSE")) {
  584. /* DGGLSE */
  585. s_copy(srnamc_1.srnamt, "DGGLSE", (ftnlen)6, (ftnlen)6);
  586. infoc_1.infot = 1;
  587. dgglse_(&c_n1, &c__0, &c__0, a, &c__1, b, &c__1, r1, r2, r3, w, &
  588. c__18, &info);
  589. chkxer_("DGGLSE", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  590. infoc_1.ok);
  591. infoc_1.infot = 2;
  592. dgglse_(&c__0, &c_n1, &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 = 3;
  597. dgglse_(&c__0, &c__0, &c_n1, 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__1, 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__1, &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 = 5;
  612. dgglse_(&c__0, &c__0, &c__0, a, &c__0, 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 = 7;
  617. dgglse_(&c__0, &c__0, &c__0, a, &c__1, b, &c__0, 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 = 12;
  622. dgglse_(&c__1, &c__1, &c__1, a, &c__1, b, &c__1, r1, r2, r3, w, &c__1,
  623. &info);
  624. chkxer_("DGGLSE", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  625. infoc_1.ok);
  626. nt += 8;
  627. /* Test error exits for the GQR path. */
  628. } else if (lsamen_(&c__3, path, "GQR")) {
  629. /* DGGQRF */
  630. s_copy(srnamc_1.srnamt, "DGGQRF", (ftnlen)6, (ftnlen)6);
  631. infoc_1.infot = 1;
  632. dggqrf_(&c_n1, &c__0, &c__0, a, &c__1, r1, b, &c__1, r2, w, &c__18, &
  633. info);
  634. chkxer_("DGGQRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  635. infoc_1.ok);
  636. infoc_1.infot = 2;
  637. dggqrf_(&c__0, &c_n1, &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 = 3;
  642. dggqrf_(&c__0, &c__0, &c_n1, 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 = 5;
  647. dggqrf_(&c__0, &c__0, &c__0, a, &c__0, 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 = 8;
  652. dggqrf_(&c__0, &c__0, &c__0, a, &c__1, r1, b, &c__0, 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 = 11;
  657. dggqrf_(&c__1, &c__1, &c__2, a, &c__1, r1, b, &c__1, r2, w, &c__1, &
  658. info);
  659. chkxer_("DGGQRF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  660. infoc_1.ok);
  661. nt += 6;
  662. /* DGGRQF */
  663. s_copy(srnamc_1.srnamt, "DGGRQF", (ftnlen)6, (ftnlen)6);
  664. infoc_1.infot = 1;
  665. dggrqf_(&c_n1, &c__0, &c__0, a, &c__1, r1, b, &c__1, r2, w, &c__18, &
  666. info);
  667. chkxer_("DGGRQF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  668. infoc_1.ok);
  669. infoc_1.infot = 2;
  670. dggrqf_(&c__0, &c_n1, &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 = 3;
  675. dggrqf_(&c__0, &c__0, &c_n1, 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 = 5;
  680. dggrqf_(&c__0, &c__0, &c__0, a, &c__0, 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 = 8;
  685. dggrqf_(&c__0, &c__0, &c__0, a, &c__1, r1, b, &c__0, 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 = 11;
  690. dggrqf_(&c__1, &c__1, &c__2, a, &c__1, r1, b, &c__1, r2, w, &c__1, &
  691. info);
  692. chkxer_("DGGRQF", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  693. infoc_1.ok);
  694. nt += 6;
  695. /* Test error exits for the DGS, DGV, DGX, and DXV paths. */
  696. } else if (lsamen_(&c__3, path, "DGS") || lsamen_(&
  697. c__3, path, "DGV") || lsamen_(&c__3, path,
  698. "DGX") || lsamen_(&c__3, path, "DXV")) {
  699. /* DGGES */
  700. s_copy(srnamc_1.srnamt, "DGGES ", (ftnlen)6, (ftnlen)6);
  701. infoc_1.infot = 1;
  702. dgges_("/", "N", "S", (L_fp)dlctes_, &c__1, a, &c__1, b, &c__1, &sdim,
  703. r1, r2, r3, q, &c__1, u, &c__1, w, &c__1, bw, &info);
  704. chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  705. infoc_1.ok);
  706. infoc_1.infot = 2;
  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 = 3;
  712. dgges_("N", "V", "/", (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 = 5;
  717. dgges_("N", "V", "S", (L_fp)dlctes_, &c_n1, 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 = 7;
  722. dgges_("N", "V", "S", (L_fp)dlctes_, &c__1, a, &c__0, 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 = 9;
  727. dgges_("N", "V", "S", (L_fp)dlctes_, &c__1, a, &c__1, b, &c__0, &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 = 15;
  732. dgges_("N", "V", "S", (L_fp)dlctes_, &c__1, a, &c__1, b, &c__1, &sdim,
  733. r1, r2, r3, q, &c__0, 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_("V", "V", "S", (L_fp)dlctes_, &c__2, a, &c__2, b, &c__2, &sdim,
  738. r1, r2, r3, q, &c__1, u, &c__2, 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 = 17;
  742. dgges_("N", "V", "S", (L_fp)dlctes_, &c__1, a, &c__1, b, &c__1, &sdim,
  743. r1, r2, r3, q, &c__1, u, &c__0, 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_("V", "V", "S", (L_fp)dlctes_, &c__2, a, &c__2, b, &c__2, &sdim,
  748. r1, r2, r3, q, &c__2, 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 = 19;
  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__2, w, &c__1, bw, &info);
  754. chkxer_("DGGES ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  755. infoc_1.ok);
  756. nt += 11;
  757. /* DGGESX */
  758. s_copy(srnamc_1.srnamt, "DGGESX", (ftnlen)6, (ftnlen)6);
  759. infoc_1.infot = 1;
  760. dggesx_("/", "N", "S", (L_fp)dlctsx_, "N", &c__1, a, &c__1, b, &c__1,
  761. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
  762. &c__1, bw, &info)
  763. ;
  764. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  765. infoc_1.ok);
  766. infoc_1.infot = 2;
  767. dggesx_("N", "/", "S", (L_fp)dlctsx_, "N", &c__1, a, &c__1, b, &c__1,
  768. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
  769. &c__1, bw, &info)
  770. ;
  771. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  772. infoc_1.ok);
  773. infoc_1.infot = 3;
  774. dggesx_("V", "V", "/", (L_fp)dlctsx_, "N", &c__1, a, &c__1, b, &c__1,
  775. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
  776. &c__1, bw, &info)
  777. ;
  778. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  779. infoc_1.ok);
  780. infoc_1.infot = 5;
  781. dggesx_("V", "V", "S", (L_fp)dlctsx_, "/", &c__1, a, &c__1, b, &c__1,
  782. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
  783. &c__1, bw, &info)
  784. ;
  785. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  786. infoc_1.ok);
  787. infoc_1.infot = 6;
  788. dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c_n1, a, &c__1, b, &c__1,
  789. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
  790. &c__1, bw, &info)
  791. ;
  792. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  793. infoc_1.ok);
  794. infoc_1.infot = 8;
  795. dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__1, a, &c__0, b, &c__1,
  796. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
  797. &c__1, bw, &info)
  798. ;
  799. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  800. infoc_1.ok);
  801. infoc_1.infot = 10;
  802. dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__1, a, &c__1, b, &c__0,
  803. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
  804. &c__1, bw, &info)
  805. ;
  806. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  807. infoc_1.ok);
  808. infoc_1.infot = 16;
  809. dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__1, a, &c__1, b, &c__1,
  810. &sdim, r1, r2, r3, q, &c__0, u, &c__1, rce, rcv, w, &c__1, iw,
  811. &c__1, bw, &info)
  812. ;
  813. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  814. infoc_1.ok);
  815. infoc_1.infot = 16;
  816. dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__2, a, &c__2, b, &c__2,
  817. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__1, iw,
  818. &c__1, bw, &info)
  819. ;
  820. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  821. infoc_1.ok);
  822. infoc_1.infot = 18;
  823. dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__1, a, &c__1, b, &c__1,
  824. &sdim, r1, r2, r3, q, &c__1, u, &c__0, rce, rcv, w, &c__1, iw,
  825. &c__1, bw, &info)
  826. ;
  827. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  828. infoc_1.ok);
  829. infoc_1.infot = 18;
  830. dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__2, a, &c__2, b, &c__2,
  831. &sdim, r1, r2, r3, q, &c__2, u, &c__1, rce, rcv, w, &c__1, iw,
  832. &c__1, bw, &info)
  833. ;
  834. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  835. infoc_1.ok);
  836. infoc_1.infot = 22;
  837. dggesx_("V", "V", "S", (L_fp)dlctsx_, "B", &c__2, a, &c__2, b, &c__2,
  838. &sdim, r1, r2, r3, q, &c__2, u, &c__2, rce, rcv, w, &c__1, iw,
  839. &c__1, bw, &info)
  840. ;
  841. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  842. infoc_1.ok);
  843. infoc_1.infot = 24;
  844. dggesx_("V", "V", "S", (L_fp)dlctsx_, "V", &c__1, a, &c__1, b, &c__1,
  845. &sdim, r1, r2, r3, q, &c__1, u, &c__1, rce, rcv, w, &c__32,
  846. iw, &c_n1, bw, &info);
  847. chkxer_("DGGESX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  848. infoc_1.ok);
  849. nt += 13;
  850. /* DGGEV */
  851. s_copy(srnamc_1.srnamt, "DGGEV ", (ftnlen)6, (ftnlen)6);
  852. infoc_1.infot = 1;
  853. dggev_("/", "N", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q, &c__1, u, &
  854. c__1, w, &c__1, &info);
  855. chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  856. infoc_1.ok);
  857. infoc_1.infot = 2;
  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 = 3;
  863. dggev_("V", "V", &c_n1, 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 = 5;
  868. dggev_("V", "V", &c__1, a, &c__0, 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 = 7;
  873. dggev_("V", "V", &c__1, a, &c__1, b, &c__0, 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 = 12;
  878. dggev_("N", "V", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q, &c__0, 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_("V", "V", &c__2, a, &c__2, b, &c__2, r1, r2, r3, q, &c__1, u, &
  884. c__2, w, &c__1, &info);
  885. chkxer_("DGGEV ", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  886. infoc_1.ok);
  887. infoc_1.infot = 14;
  888. dggev_("V", "N", &c__2, a, &c__2, b, &c__2, r1, r2, r3, q, &c__2, u, &
  889. c__0, 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", "V", &c__2, a, &c__2, b, &c__2, r1, r2, r3, q, &c__2, 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 = 16;
  898. dggev_("V", "V", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q, &c__1, 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. nt += 10;
  903. /* DGGEVX */
  904. s_copy(srnamc_1.srnamt, "DGGEVX", (ftnlen)6, (ftnlen)6);
  905. infoc_1.infot = 1;
  906. dggevx_("/", "N", "N", "N", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q,
  907. &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  908. w, &c__1, iw, bw, &info);
  909. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  910. infoc_1.ok);
  911. infoc_1.infot = 2;
  912. dggevx_("N", "/", "N", "N", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q,
  913. &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  914. w, &c__1, iw, bw, &info);
  915. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  916. infoc_1.ok);
  917. infoc_1.infot = 3;
  918. dggevx_("N", "N", "/", "N", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q,
  919. &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  920. w, &c__1, iw, bw, &info);
  921. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  922. infoc_1.ok);
  923. infoc_1.infot = 4;
  924. dggevx_("N", "N", "N", "/", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q,
  925. &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  926. w, &c__1, iw, bw, &info);
  927. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  928. infoc_1.ok);
  929. infoc_1.infot = 5;
  930. dggevx_("N", "N", "N", "N", &c_n1, a, &c__1, b, &c__1, r1, r2, r3, q,
  931. &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  932. w, &c__1, iw, bw, &info);
  933. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  934. infoc_1.ok);
  935. infoc_1.infot = 7;
  936. dggevx_("N", "N", "N", "N", &c__1, a, &c__0, b, &c__1, r1, r2, r3, q,
  937. &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  938. w, &c__1, iw, bw, &info);
  939. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  940. infoc_1.ok);
  941. infoc_1.infot = 9;
  942. dggevx_("N", "N", "N", "N", &c__1, a, &c__1, b, &c__0, r1, r2, r3, q,
  943. &c__1, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  944. w, &c__1, iw, bw, &info);
  945. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  946. infoc_1.ok);
  947. infoc_1.infot = 14;
  948. dggevx_("N", "N", "N", "N", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q,
  949. &c__0, u, &c__1, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  950. w, &c__1, iw, bw, &info);
  951. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  952. infoc_1.ok);
  953. infoc_1.infot = 14;
  954. dggevx_("N", "V", "N", "N", &c__2, a, &c__2, b, &c__2, r1, r2, r3, q,
  955. &c__1, u, &c__2, &c__1, &c__2, ls, rs, &anrm, &bnrm, rce, rcv,
  956. w, &c__1, iw, bw, &info);
  957. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  958. infoc_1.ok);
  959. infoc_1.infot = 16;
  960. dggevx_("N", "N", "N", "N", &c__1, a, &c__1, b, &c__1, r1, r2, r3, q,
  961. &c__1, u, &c__0, &c__1, &c__1, ls, rs, &anrm, &bnrm, rce, rcv,
  962. w, &c__1, iw, bw, &info);
  963. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  964. infoc_1.ok);
  965. infoc_1.infot = 16;
  966. dggevx_("N", "N", "V", "N", &c__2, a, &c__2, b, &c__2, r1, r2, r3, q,
  967. &c__2, u, &c__1, &c__1, &c__2, ls, rs, &anrm, &bnrm, rce, rcv,
  968. w, &c__1, iw, bw, &info);
  969. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  970. infoc_1.ok);
  971. infoc_1.infot = 26;
  972. dggevx_("N", "N", "V", "N", &c__2, a, &c__2, b, &c__2, r1, r2, r3, q,
  973. &c__2, u, &c__2, &c__1, &c__2, ls, rs, &anrm, &bnrm, rce, rcv,
  974. w, &c__1, iw, bw, &info);
  975. chkxer_("DGGEVX", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  976. infoc_1.ok);
  977. nt += 12;
  978. /* DTGEXC */
  979. s_copy(srnamc_1.srnamt, "DTGEXC", (ftnlen)6, (ftnlen)6);
  980. infoc_1.infot = 3;
  981. dtgexc_(&c_true, &c_true, &c_n1, a, &c__1, b, &c__1, q, &c__1, z__, &
  982. c__1, &ifst, &ilst, w, &c__1, &info);
  983. chkxer_("DTGEXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  984. infoc_1.ok);
  985. infoc_1.infot = 5;
  986. dtgexc_(&c_true, &c_true, &c__1, a, &c__0, 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 = 7;
  991. dtgexc_(&c_true, &c_true, &c__1, a, &c__1, b, &c__0, 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 = 9;
  996. dtgexc_(&c_false, &c_true, &c__1, a, &c__1, b, &c__1, q, &c__0, 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_true, &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 = 11;
  1006. dtgexc_(&c_true, &c_false, &c__1, a, &c__1, b, &c__1, q, &c__1, z__, &
  1007. c__0, &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_true, &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 = 15;
  1016. dtgexc_(&c_true, &c_true, &c__1, a, &c__1, b, &c__1, q, &c__1, 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. nt += 8;
  1021. /* DTGSEN */
  1022. s_copy(srnamc_1.srnamt, "DTGSEN", (ftnlen)6, (ftnlen)6);
  1023. infoc_1.infot = 1;
  1024. dtgsen_(&c_n1, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1025. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
  1026. iw, &c__1, &info);
  1027. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1028. infoc_1.ok);
  1029. infoc_1.infot = 5;
  1030. dtgsen_(&c__1, &c_true, &c_true, sel, &c_n1, a, &c__1, b, &c__1, r1,
  1031. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
  1032. iw, &c__1, &info);
  1033. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1034. infoc_1.ok);
  1035. infoc_1.infot = 7;
  1036. dtgsen_(&c__1, &c_true, &c_true, sel, &c__1, a, &c__0, b, &c__1, r1,
  1037. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
  1038. iw, &c__1, &info);
  1039. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1040. infoc_1.ok);
  1041. infoc_1.infot = 9;
  1042. dtgsen_(&c__1, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__0, r1,
  1043. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
  1044. iw, &c__1, &info);
  1045. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1046. infoc_1.ok);
  1047. infoc_1.infot = 14;
  1048. dtgsen_(&c__1, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1049. r2, r3, q, &c__0, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
  1050. iw, &c__1, &info);
  1051. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1052. infoc_1.ok);
  1053. infoc_1.infot = 16;
  1054. dtgsen_(&c__1, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1055. r2, r3, q, &c__1, z__, &c__0, &m, &tola, &tolb, rcv, w, &c__1,
  1056. iw, &c__1, &info);
  1057. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1058. infoc_1.ok);
  1059. infoc_1.infot = 22;
  1060. dtgsen_(&c__0, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1061. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
  1062. iw, &c__1, &info);
  1063. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1064. infoc_1.ok);
  1065. infoc_1.infot = 22;
  1066. dtgsen_(&c__1, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1067. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
  1068. iw, &c__1, &info);
  1069. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1070. infoc_1.ok);
  1071. infoc_1.infot = 22;
  1072. dtgsen_(&c__2, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1073. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &c__1,
  1074. iw, &c__1, &info);
  1075. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1076. infoc_1.ok);
  1077. infoc_1.infot = 24;
  1078. dtgsen_(&c__0, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1079. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &
  1080. c__20, iw, &c__0, &info);
  1081. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1082. infoc_1.ok);
  1083. infoc_1.infot = 24;
  1084. dtgsen_(&c__1, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1085. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &
  1086. c__20, iw, &c__0, &info);
  1087. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1088. infoc_1.ok);
  1089. infoc_1.infot = 24;
  1090. dtgsen_(&c__2, &c_true, &c_true, sel, &c__1, a, &c__1, b, &c__1, r1,
  1091. r2, r3, q, &c__1, z__, &c__1, &m, &tola, &tolb, rcv, w, &
  1092. c__20, iw, &c__1, &info);
  1093. chkxer_("DTGSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1094. infoc_1.ok);
  1095. nt += 12;
  1096. /* DTGSNA */
  1097. s_copy(srnamc_1.srnamt, "DTGSNA", (ftnlen)6, (ftnlen)6);
  1098. infoc_1.infot = 1;
  1099. dtgsna_("/", "A", sel, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &c__1,
  1100. r1, r2, &c__1, &m, w, &c__1, iw, &info);
  1101. chkxer_("DTGSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1102. infoc_1.ok);
  1103. infoc_1.infot = 2;
  1104. dtgsna_("B", "/", 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 = 4;
  1109. dtgsna_("B", "A", sel, &c_n1, 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 = 6;
  1114. dtgsna_("B", "A", sel, &c__1, a, &c__0, 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 = 8;
  1119. dtgsna_("B", "A", sel, &c__1, a, &c__1, b, &c__0, 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 = 10;
  1124. dtgsna_("E", "A", sel, &c__1, a, &c__1, b, &c__1, q, &c__0, 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 = 12;
  1129. dtgsna_("E", "A", sel, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &c__0,
  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 = 15;
  1134. dtgsna_("E", "A", sel, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &c__1,
  1135. r1, r2, &c__0, &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 = 18;
  1139. dtgsna_("E", "A", sel, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &c__1,
  1140. r1, r2, &c__1, &m, w, &c__0, iw, &info);
  1141. chkxer_("DTGSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1142. infoc_1.ok);
  1143. nt += 9;
  1144. /* DTGSYL */
  1145. s_copy(srnamc_1.srnamt, "DTGSYL", (ftnlen)6, (ftnlen)6);
  1146. infoc_1.infot = 1;
  1147. dtgsyl_("/", &c__0, &c__1, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &
  1148. c__1, v, &c__1, z__, &c__1, &scale, &dif, w, &c__1, iw, &info);
  1149. chkxer_("DTGSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
  1150. infoc_1.ok);
  1151. infoc_1.infot = 2;
  1152. dtgsyl_("N", &c_n1, &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 = 3;
  1157. dtgsyl_("N", &c__0, &c__0, &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 = 4;
  1162. dtgsyl_("N", &c__0, &c__1, &c__0, 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 = 6;
  1167. dtgsyl_("N", &c__0, &c__1, &c__1, a, &c__0, 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 = 8;
  1172. dtgsyl_("N", &c__0, &c__1, &c__1, a, &c__1, b, &c__0, 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 = 10;
  1177. dtgsyl_("N", &c__0, &c__1, &c__1, a, &c__1, b, &c__1, q, &c__0, 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 = 12;
  1182. dtgsyl_("N", &c__0, &c__1, &c__1, a, &c__1, b, &c__1, q, &c__1, u, &
  1183. c__0, 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 = 14;
  1187. dtgsyl_

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