Mail Archives: cygwin-developers/1998/11/07/16:27:42
Were some of the bessel functions left out on purpose? I'm curious since
the non-standard float versions are there, but not the "more standard"
double precision versions.
Here's a trivial patch if not. I'm appending a simple test case as well.
/* bessel function test case. */
#include <stdio.h>
#include <math.h>
static int
dbl_equal (double d1, double d2)
{
const double eps = 1.0e-15L;
return (fabs (d1 - d2) <= eps);
}
int
main ()
{
double t1 = j1 (5.0);
double t2 = jn (1, 5.0);
double t3 = y1 (5.0);
double t4 = yn (1, 5.0);
printf ("j1 (5.0) = %f\n", t1);
printf ("jn (1, 5.0) = %f\n", t2);
printf ("y1 (5.0) = %f\n", t3);
printf ("yn (1, 5.0) = %f\n", t4);
printf ("Testing `j1 (x) == jn (1, x)' ... %s\n",
(dbl_equal (t1, t2)) ? "passed" : "failed");
printf ("Testing `y1 (x) == yn (1, x)' ... %s\n",
(dbl_equal (t3, t4)) ? "passed" : "failed");
return 0;
}
Sat Nov 7 16:54:34 1998 Mumit Khan <khan AT xraylith DOT wisc DOT edu>
* cygwin.din ({j1,jn,y1,yn}): Add missing exports.
Index: cygwin.din
===================================================================
RCS file: /mounts/sdb5/src/cdk-b20/CVSROOT/b20/winsup/cygwin.din,v
retrieving revision 1.1.1.1
diff -u -3 -p -r1.1.1.1 cygwin.din
--- cygwin.din 1998/11/05 21:47:13 1.1.1.1
+++ cygwin.din 1998/11/07 22:54:17
@@ -370,7 +370,9 @@ isxdigit
_isxdigit = isxdigit
j0
j0f
+j1
j1f
+jn
jnf
kill
_kill = kill
@@ -781,7 +783,9 @@ writev
_writev = writev
y0
y0f
+y1
y1f
+yn
ynf
setmode
_setmode = setmode
Regards,
Mumit -- khan AT xraylith DOT wisc DOT edu
- Raw text -