Print this page
@@ -56,10 +56,11 @@
* values of n>1.
*
*/
#include "libm.h"
+#include "longdouble.h"
#include <float.h> /* LDBL_MAX */
#define GENERIC long double
static const GENERIC
@@ -69,11 +70,11 @@
one = 1.0L;
GENERIC
jnl(n,x) int n; GENERIC x;{
int i, sgn;
- GENERIC a, b, temp, z, w;
+ GENERIC a, b, temp = 0, z, w;
/* J(-n,x) = (-1)^n * J(n, x), J(n, -x) = (-1)^n * J(n, x)
* Thus, J(-n,x) = J(n,-x)
*/
if(n<0){
@@ -209,18 +210,20 @@
GENERIC ynl(n,x)
int n; GENERIC x;{
int i;
int sign;
- GENERIC a, b, temp;
+ GENERIC a, b, temp = 0;
- if(x!=x) return x+x;
- if (x <= zero)
+ if(x!=x)
+ return x+x;
+ if (x <= zero) {
if(x==zero)
return -one/zero;
else
return zero/zero;
+ }
sign = 1;
if(n<0){
n = -n;
if((n&1) == 1) sign = -1;
}