/clapack/CLAPACK-3.2.1/TESTING/EIG/derrec.c

https://github.com/zye/cmvs_dep · C · 359 lines · 263 code · 39 blank · 57 comment · 5 complexity · 9dc9328f9ae27b670d2a9fc45134eed0 MD5 · raw file 

    

  1. /* derrec.f -- translated by f2c (version 20061008).
    2. 

  2.    You must link the resulting object file with libf2c:
    3. 

  3. 	on Microsoft Windows system, link with libf2c.lib;
    4. 

  4. 	on Linux or Unix systems, link with .../path/to/libf2c.a -lm
    5. 

  5. 	or, if you install libf2c.a in a standard place, with -lf2c -lm
    6. 

  6. 	-- in that order, at the end of the command line, as in
    7. 

  7. 		cc *.o -lf2c -lm
    8. 

  8. 	Source for libf2c is in /netlib/f2c/libf2c.zip, e.g.,
    9. 

    10. 

  10. 		http://www.netlib.org/f2c/libf2c.zip
    11. 

  11. */
    12. 

  12. 

  13. #include "f2c.h"
    14. 

  14. #include "blaswrap.h"
    15. 

  15. 

  16. /* Common Block Declarations */
    17. 

  17. 

  18. struct {
    19. 

  19.     integer infot, nout;
    20. 

  20.     logical ok, lerr;
    21. 

  21. } infoc_;
    22. 

  22. 

  23. #define infoc_1 infoc_
    24. 

  24. 

  25. struct {
    26. 

  26.     char srnamt[32];
    27. 

  27. } srnamc_;
    28. 

  28. 

  29. #define srnamc_1 srnamc_
    30. 

  30. 

  31. /* Table of constant values */
    32. 

  32. 

  33. static integer c__1 = 1;
    34. 

  34. static integer c__0 = 0;
    35. 

  35. static integer c_n1 = -1;
    36. 

  36. static integer c__2 = 2;
    37. 

  37. static integer c__3 = 3;
    38. 

  38. static integer c__4 = 4;
    39. 

  39. 

  40. /* Subroutine */ int derrec_(char *path, integer *nunit)
    41. 

  41. {
    42. 

  42.     /* Format strings */
    43. 

  43.     static char fmt_9999[] = "(1x,a3,\002 routines passed the tests of the e"
    44. 

  44. 	    "rror exits (\002,i3,\002 tests done)\002)";
    45. 

  45.     static char fmt_9998[] = "(\002 *** \002,a3,\002 routines failed the tes"
    46. 

  46. 	    "ts of the error ex\002,\002its ***\002)";
    47. 

  47. 

  48.     /* Builtin functions */
    49. 

  49.     /* Subroutine */ int s_copy(char *, char *, ftnlen, ftnlen);
    50. 

  50.     integer s_wsfe(cilist *), do_fio(integer *, char *, ftnlen), e_wsfe(void);
    51. 

  51. 

  52.     /* Local variables */
    53. 

  53.     doublereal a[16]	/* was [4][4] */, b[16]	/* was [4][4] */, c__[16]	
    54. 

  54. 	    /* was [4][4] */;
    55. 

  55.     integer i__, j, m;
    56. 

  56.     doublereal s[4], wi[4];
    57. 

  57.     integer nt;
    58. 

  58.     doublereal wr[4];
    59. 

  59.     logical sel[4];
    60. 

  60.     doublereal sep[4];
    61. 

  61.     integer info, ifst, ilst;
    62. 

  62.     doublereal work[4], scale;
    63. 

  63.     integer iwork[4];
    64. 

  64.     extern /* Subroutine */ int chkxer_(char *, integer *, integer *, logical 
    65. 

  65. 	    *, logical *), dtrexc_(char *, integer *, doublereal *, 
    66. 

  66. 	    integer *, doublereal *, integer *, integer *, integer *, 
    67. 

  67. 	    doublereal *, integer *), dtrsna_(char *, char *, logical 
    68. 

  68. 	    *, integer *, doublereal *, integer *, doublereal *, integer *, 
    69. 

  69. 	    doublereal *, integer *, doublereal *, doublereal *, integer *, 
    70. 

  70. 	    integer *, doublereal *, integer *, integer *, integer *), dtrsen_(char *, char *, logical *, integer *, doublereal 
    71. 

  71. 	    *, integer *, doublereal *, integer *, doublereal *, doublereal *, 
    72. 

  72. 	     integer *, doublereal *, doublereal *, doublereal *, integer *, 
    73. 

  73. 	    integer *, integer *, integer *), dtrsyl_(char *, 
    74. 

  74. 	    char *, integer *, integer *, integer *, doublereal *, integer *, 
    75. 

  75. 	    doublereal *, integer *, doublereal *, integer *, doublereal *, 
    76. 

  76. 	    integer *);
    77. 

  77. 

  78.     /* Fortran I/O blocks */
    79. 

  79.     static cilist io___19 = { 0, 0, 0, fmt_9999, 0 };
    80. 

  80.     static cilist io___20 = { 0, 0, 0, fmt_9998, 0 };
    81. 

  81. 

  82. 

  83. 

  84. /*  -- LAPACK test routine (version 3.1) -- */
    85. 

  85. /*     Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. */
    86. 

  86. /*     November 2006 */
    87. 

  87. 

  88. /*     .. Scalar Arguments .. */
    89. 

  89. /*     .. */
    90. 

  90. 

  91. /*  Purpose */
    92. 

  92. /*  ======= */
    93. 

  93. 

  94. /*  DERREC tests the error exits for the routines for eigen- condition */
    95. 

  95. /*  estimation for DOUBLE PRECISION matrices: */
    96. 

  96. /*     DTRSYL, STREXC, STRSNA and STRSEN. */
    97. 

  97. 

  98. /*  Arguments */
    99. 

  99. /*  ========= */
    100. 

  100. 

  101. /*  PATH    (input) CHARACTER*3 */
    102. 

  102. /*          The LAPACK path name for the routines to be tested. */
    103. 

  103. 

  104. /*  NUNIT   (input) INTEGER */
    105. 

  105. /*          The unit number for output. */
    106. 

  106. 

  107. /*  ===================================================================== */
    108. 

  108. 

  109. /*     .. Parameters .. */
    110. 

  110. /*     .. */
    111. 

  111. /*     .. Local Scalars .. */
    112. 

  112. /*     .. */
    113. 

  113. /*     .. Local Arrays .. */
    114. 

  114. /*     .. */
    115. 

  115. /*     .. External Subroutines .. */
    116. 

  116. /*     .. */
    117. 

  117. /*     .. Scalars in Common .. */
    118. 

  118. /*     .. */
    119. 

  119. /*     .. Common blocks .. */
    120. 

  120. /*     .. */
    121. 

  121. /*     .. Executable Statements .. */
    122. 

  122. 

  123.     infoc_1.nout = *nunit;
    124. 

  124.     infoc_1.ok = TRUE_;
    125. 

  125.     nt = 0;
    126. 

  126. 

  127. /*     Initialize A, B and SEL */
    128. 

  128. 

  129.     for (j = 1; j <= 4; ++j) {
    130. 

  130. 	for (i__ = 1; i__ <= 4; ++i__) {
    131. 

  131. 	    a[i__ + (j << 2) - 5] = 0.;
    132. 

  132. 	    b[i__ + (j << 2) - 5] = 0.;
    133. 

  133. /* L10: */
    134. 

  134. 	}
    135. 

  135. /* L20: */
    136. 

  136.     }
    137. 

  137.     for (i__ = 1; i__ <= 4; ++i__) {
    138. 

  138. 	a[i__ + (i__ << 2) - 5] = 1.;
    139. 

  139. 	sel[i__ - 1] = TRUE_;
    140. 

  140. /* L30: */
    141. 

  141.     }
    142. 

  142. 

  143. /*     Test DTRSYL */
    144. 

  144. 

  145.     s_copy(srnamc_1.srnamt, "DTRSYL", (ftnlen)32, (ftnlen)6);
    146. 

  146.     infoc_1.infot = 1;
    147. 

  147.     dtrsyl_("X", "N", &c__1, &c__0, &c__0, a, &c__1, b, &c__1, c__, &c__1, &
    148. 

  148. 	    scale, &info);
    149. 

  149.     chkxer_("DTRSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    150. 

  150. 	    infoc_1.ok);
    151. 

  151.     infoc_1.infot = 2;
    152. 

  152.     dtrsyl_("N", "X", &c__1, &c__0, &c__0, a, &c__1, b, &c__1, c__, &c__1, &
    153. 

  153. 	    scale, &info);
    154. 

  154.     chkxer_("DTRSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    155. 

  155. 	    infoc_1.ok);
    156. 

  156.     infoc_1.infot = 3;
    157. 

  157.     dtrsyl_("N", "N", &c__0, &c__0, &c__0, a, &c__1, b, &c__1, c__, &c__1, &
    158. 

  158. 	    scale, &info);
    159. 

  159.     chkxer_("DTRSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    160. 

  160. 	    infoc_1.ok);
    161. 

  161.     infoc_1.infot = 4;
    162. 

  162.     dtrsyl_("N", "N", &c__1, &c_n1, &c__0, a, &c__1, b, &c__1, c__, &c__1, &
    163. 

  163. 	    scale, &info);
    164. 

  164.     chkxer_("DTRSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    165. 

  165. 	    infoc_1.ok);
    166. 

  166.     infoc_1.infot = 5;
    167. 

  167.     dtrsyl_("N", "N", &c__1, &c__0, &c_n1, a, &c__1, b, &c__1, c__, &c__1, &
    168. 

  168. 	    scale, &info);
    169. 

  169.     chkxer_("DTRSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    170. 

  170. 	    infoc_1.ok);
    171. 

  171.     infoc_1.infot = 7;
    172. 

  172.     dtrsyl_("N", "N", &c__1, &c__2, &c__0, a, &c__1, b, &c__1, c__, &c__2, &
    173. 

  173. 	    scale, &info);
    174. 

  174.     chkxer_("DTRSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    175. 

  175. 	    infoc_1.ok);
    176. 

  176.     infoc_1.infot = 9;
    177. 

  177.     dtrsyl_("N", "N", &c__1, &c__0, &c__2, a, &c__1, b, &c__1, c__, &c__1, &
    178. 

  178. 	    scale, &info);
    179. 

  179.     chkxer_("DTRSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    180. 

  180. 	    infoc_1.ok);
    181. 

  181.     infoc_1.infot = 11;
    182. 

  182.     dtrsyl_("N", "N", &c__1, &c__2, &c__0, a, &c__2, b, &c__1, c__, &c__1, &
    183. 

  183. 	    scale, &info);
    184. 

  184.     chkxer_("DTRSYL", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    185. 

  185. 	    infoc_1.ok);
    186. 

  186.     nt += 8;
    187. 

  187. 

  188. /*     Test DTREXC */
    189. 

  189. 

  190.     s_copy(srnamc_1.srnamt, "DTREXC", (ftnlen)32, (ftnlen)6);
    191. 

  191.     ifst = 1;
    192. 

  192.     ilst = 1;
    193. 

  193.     infoc_1.infot = 1;
    194. 

  194.     dtrexc_("X", &c__1, a, &c__1, b, &c__1, &ifst, &ilst, work, &info);
    195. 

  195.     chkxer_("DTREXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    196. 

  196. 	    infoc_1.ok);
    197. 

  197.     infoc_1.infot = 7;
    198. 

  198.     dtrexc_("N", &c__0, a, &c__1, b, &c__1, &ifst, &ilst, work, &info);
    199. 

  199.     chkxer_("DTREXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    200. 

  200. 	    infoc_1.ok);
    201. 

  201.     infoc_1.infot = 4;
    202. 

  202.     ilst = 2;
    203. 

  203.     dtrexc_("N", &c__2, a, &c__1, b, &c__1, &ifst, &ilst, work, &info);
    204. 

  204.     chkxer_("DTREXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    205. 

  205. 	    infoc_1.ok);
    206. 

  206.     infoc_1.infot = 6;
    207. 

  207.     dtrexc_("V", &c__2, a, &c__2, b, &c__1, &ifst, &ilst, work, &info);
    208. 

  208.     chkxer_("DTREXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    209. 

  209. 	    infoc_1.ok);
    210. 

  210.     infoc_1.infot = 7;
    211. 

  211.     ifst = 0;
    212. 

  212.     ilst = 1;
    213. 

  213.     dtrexc_("V", &c__1, a, &c__1, b, &c__1, &ifst, &ilst, work, &info);
    214. 

  214.     chkxer_("DTREXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    215. 

  215. 	    infoc_1.ok);
    216. 

  216.     infoc_1.infot = 7;
    217. 

  217.     ifst = 2;
    218. 

  218.     dtrexc_("V", &c__1, a, &c__1, b, &c__1, &ifst, &ilst, work, &info);
    219. 

  219.     chkxer_("DTREXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    220. 

  220. 	    infoc_1.ok);
    221. 

  221.     infoc_1.infot = 8;
    222. 

  222.     ifst = 1;
    223. 

  223.     ilst = 0;
    224. 

  224.     dtrexc_("V", &c__1, a, &c__1, b, &c__1, &ifst, &ilst, work, &info);
    225. 

  225.     chkxer_("DTREXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    226. 

  226. 	    infoc_1.ok);
    227. 

  227.     infoc_1.infot = 8;
    228. 

  228.     ilst = 2;
    229. 

  229.     dtrexc_("V", &c__1, a, &c__1, b, &c__1, &ifst, &ilst, work, &info);
    230. 

  230.     chkxer_("DTREXC", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    231. 

  231. 	    infoc_1.ok);
    232. 

  232.     nt += 8;
    233. 

  233. 

  234. /*     Test DTRSNA */
    235. 

  235. 

  236.     s_copy(srnamc_1.srnamt, "DTRSNA", (ftnlen)32, (ftnlen)6);
    237. 

  237.     infoc_1.infot = 1;
    238. 

  238.     dtrsna_("X", "A", sel, &c__0, a, &c__1, b, &c__1, c__, &c__1, s, sep, &
    239. 

  239. 	    c__1, &m, work, &c__1, iwork, &info);
    240. 

  240.     chkxer_("DTRSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    241. 

  241. 	    infoc_1.ok);
    242. 

  242.     infoc_1.infot = 2;
    243. 

  243.     dtrsna_("B", "X", sel, &c__0, a, &c__1, b, &c__1, c__, &c__1, s, sep, &
    244. 

  244. 	    c__1, &m, work, &c__1, iwork, &info);
    245. 

  245.     chkxer_("DTRSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    246. 

  246. 	    infoc_1.ok);
    247. 

  247.     infoc_1.infot = 4;
    248. 

  248.     dtrsna_("B", "A", sel, &c_n1, a, &c__1, b, &c__1, c__, &c__1, s, sep, &
    249. 

  249. 	    c__1, &m, work, &c__1, iwork, &info);
    250. 

  250.     chkxer_("DTRSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    251. 

  251. 	    infoc_1.ok);
    252. 

  252.     infoc_1.infot = 6;
    253. 

  253.     dtrsna_("V", "A", sel, &c__2, a, &c__1, b, &c__1, c__, &c__1, s, sep, &
    254. 

  254. 	    c__2, &m, work, &c__2, iwork, &info);
    255. 

  255.     chkxer_("DTRSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    256. 

  256. 	    infoc_1.ok);
    257. 

  257.     infoc_1.infot = 8;
    258. 

  258.     dtrsna_("B", "A", sel, &c__2, a, &c__2, b, &c__1, c__, &c__2, s, sep, &
    259. 

  259. 	    c__2, &m, work, &c__2, iwork, &info);
    260. 

  260.     chkxer_("DTRSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    261. 

  261. 	    infoc_1.ok);
    262. 

  262.     infoc_1.infot = 10;
    263. 

  263.     dtrsna_("B", "A", sel, &c__2, a, &c__2, b, &c__2, c__, &c__1, s, sep, &
    264. 

  264. 	    c__2, &m, work, &c__2, iwork, &info);
    265. 

  265.     chkxer_("DTRSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    266. 

  266. 	    infoc_1.ok);
    267. 

  267.     infoc_1.infot = 13;
    268. 

  268.     dtrsna_("B", "A", sel, &c__1, a, &c__1, b, &c__1, c__, &c__1, s, sep, &
    269. 

  269. 	    c__0, &m, work, &c__1, iwork, &info);
    270. 

  270.     chkxer_("DTRSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    271. 

  271. 	    infoc_1.ok);
    272. 

  272.     infoc_1.infot = 13;
    273. 

  273.     dtrsna_("B", "S", sel, &c__2, a, &c__2, b, &c__2, c__, &c__2, s, sep, &
    274. 

  274. 	    c__1, &m, work, &c__2, iwork, &info);
    275. 

  275.     chkxer_("DTRSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    276. 

  276. 	    infoc_1.ok);
    277. 

  277.     infoc_1.infot = 16;
    278. 

  278.     dtrsna_("B", "A", sel, &c__2, a, &c__2, b, &c__2, c__, &c__2, s, sep, &
    279. 

  279. 	    c__2, &m, work, &c__1, iwork, &info);
    280. 

  280.     chkxer_("DTRSNA", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    281. 

  281. 	    infoc_1.ok);
    282. 

  282.     nt += 9;
    283. 

  283. 

  284. /*     Test DTRSEN */
    285. 

  285. 

  286.     sel[0] = FALSE_;
    287. 

  287.     s_copy(srnamc_1.srnamt, "DTRSEN", (ftnlen)32, (ftnlen)6);
    288. 

  288.     infoc_1.infot = 1;
    289. 

  289.     dtrsen_("X", "N", sel, &c__0, a, &c__1, b, &c__1, wr, wi, &m, s, sep, 
    290. 

  290. 	    work, &c__1, iwork, &c__1, &info);
    291. 

  291.     chkxer_("DTRSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    292. 

  292. 	    infoc_1.ok);
    293. 

  293.     infoc_1.infot = 2;
    294. 

  294.     dtrsen_("N", "X", sel, &c__0, a, &c__1, b, &c__1, wr, wi, &m, s, sep, 
    295. 

  295. 	    work, &c__1, iwork, &c__1, &info);
    296. 

  296.     chkxer_("DTRSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    297. 

  297. 	    infoc_1.ok);
    298. 

  298.     infoc_1.infot = 4;
    299. 

  299.     dtrsen_("N", "N", sel, &c_n1, a, &c__1, b, &c__1, wr, wi, &m, s, sep, 
    300. 

  300. 	    work, &c__1, iwork, &c__1, &info);
    301. 

  301.     chkxer_("DTRSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    302. 

  302. 	    infoc_1.ok);
    303. 

  303.     infoc_1.infot = 6;
    304. 

  304.     dtrsen_("N", "N", sel, &c__2, a, &c__1, b, &c__1, wr, wi, &m, s, sep, 
    305. 

  305. 	    work, &c__2, iwork, &c__1, &info);
    306. 

  306.     chkxer_("DTRSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    307. 

  307. 	    infoc_1.ok);
    308. 

  308.     infoc_1.infot = 8;
    309. 

  309.     dtrsen_("N", "V", sel, &c__2, a, &c__2, b, &c__1, wr, wi, &m, s, sep, 
    310. 

  310. 	    work, &c__1, iwork, &c__1, &info);
    311. 

  311.     chkxer_("DTRSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    312. 

  312. 	    infoc_1.ok);
    313. 

  313.     infoc_1.infot = 15;
    314. 

  314.     dtrsen_("N", "V", sel, &c__2, a, &c__2, b, &c__2, wr, wi, &m, s, sep, 
    315. 

  315. 	    work, &c__0, iwork, &c__1, &info);
    316. 

  316.     chkxer_("DTRSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    317. 

  317. 	    infoc_1.ok);
    318. 

  318.     infoc_1.infot = 15;
    319. 

  319.     dtrsen_("E", "V", sel, &c__3, a, &c__3, b, &c__3, wr, wi, &m, s, sep, 
    320. 

  320. 	    work, &c__1, iwork, &c__1, &info);
    321. 

  321.     chkxer_("DTRSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    322. 

  322. 	    infoc_1.ok);
    323. 

  323.     infoc_1.infot = 15;
    324. 

  324.     dtrsen_("V", "V", sel, &c__3, a, &c__3, b, &c__3, wr, wi, &m, s, sep, 
    325. 

  325. 	    work, &c__3, iwork, &c__2, &info);
    326. 

  326.     chkxer_("DTRSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    327. 

  327. 	    infoc_1.ok);
    328. 

  328.     infoc_1.infot = 17;
    329. 

  329.     dtrsen_("E", "V", sel, &c__2, a, &c__2, b, &c__2, wr, wi, &m, s, sep, 
    330. 

  330. 	    work, &c__1, iwork, &c__0, &info);
    331. 

  331.     chkxer_("DTRSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    332. 

  332. 	    infoc_1.ok);
    333. 

  333.     infoc_1.infot = 17;
    334. 

  334.     dtrsen_("V", "V", sel, &c__3, a, &c__3, b, &c__3, wr, wi, &m, s, sep, 
    335. 

  335. 	    work, &c__4, iwork, &c__1, &info);
    336. 

  336.     chkxer_("DTRSEN", &infoc_1.infot, &infoc_1.nout, &infoc_1.lerr, &
    337. 

  337. 	    infoc_1.ok);
    338. 

  338.     nt += 10;
    339. 

  339. 

  340. /*     Print a summary line. */
    341. 

  341. 

  342.     if (infoc_1.ok) {
    343. 

  343. 	io___19.ciunit = infoc_1.nout;
    344. 

  344. 	s_wsfe(&io___19);
    345. 

  345. 	do_fio(&c__1, path, (ftnlen)3);
    346. 

  346. 	do_fio(&c__1, (char *)&nt, (ftnlen)sizeof(integer));
    347. 

  347. 	e_wsfe();
    348. 

  348.     } else {
    349. 

  349. 	io___20.ciunit = infoc_1.nout;
    350. 

  350. 	s_wsfe(&io___20);
    351. 

  351. 	do_fio(&c__1, path, (ftnlen)3);
    352. 

  352. 	e_wsfe();
    353. 

  353.     }
    354. 

  354. 

  355.     return 0;
    356. 

  356. 

  357. /*     End of DERREC */
    358. 

  358. 

  359. } /* derrec_ */
    360.