Lines Matching +full:- +full:a
6 * Developed at SunPro, a Sun Microsystems, Inc. business.
19 * y0(0)=y1(0)=yn(n,0) = -inf with division by zero signal;
20 * y0(-ve)=y1(-ve)=yn(n,-ve) are NaN with invalid signal.
26 * for n>x, a continued fraction approximation to
27 * j(n,x)/j(n-1,x) is evaluated and then backward
28 * recursion is used starting from a supposed value
44 invsqrtpi = _F_64(5.64189583547756279280e-01), /* 0x3FE20DD7, 0x50429B6D */
54 __float64 a, b, temp, di; in jn64() local
63 /* J(-n,x) = (-1)^n * J(n, x), J(n, -x) = (-1)^n * J(n, x) in jn64()
64 * Thus, J(-n,x) = J(n,-x) in jn64()
70 n = -n; in jn64()
71 x = -x; in jn64()
78 sgn = (n & 1) & (hx >> 31); /* even n -- 0, odd n -- sign(x) */ in jn64()
83 /* Safe to use J(n+1,x)=2n/x *J(n,x)-J(n-1,x) */ in jn64()
86 * Jn(x) = cos(x-(2n+1)*pi/4)*sqrt(2/x*pi) in jn64()
87 * Yn(x) = sin(x-(2n+1)*pi/4)*sqrt(2/x*pi) in jn64()
89 * xn=x-(2n+1)*pi/4, sqt2 = sqrt(2),then in jn64()
92 * ---------------------------------- in jn64()
93 * 0 s-c c+s in jn64()
94 * 1 -s-c -c+s in jn64()
95 * 2 -s+c -c-s in jn64()
96 * 3 s+c c-s in jn64()
104 temp = -cos64(x) + sin64(x); in jn64()
107 temp = -cos64(x) - sin64(x); in jn64()
110 temp = cos64(x) - sin64(x); in jn64()
115 a = j064(x); in jn64()
119 b = b * ((__float64)(i + i) / x) - a; /* avoid underflow */ in jn64()
120 a = temp; in jn64()
124 if (ix < 0x3e100000) { /* x < 2**-29 */ in jn64()
126 * J(n,x) = 1/n!*(x/2)^n - ... in jn64()
133 for (a = one, i = 2; i <= n; i++) { in jn64()
134 a *= (__float64)i; /* a = n! */ in jn64()
137 b = b / a; in jn64()
142 * J(n,x)/J(n-1,x) = ---- ------ ------ ..... in jn64()
143 * 2n - 2(n+1) - 2(n+2) in jn64()
146 * (for large x) = ---- ------ ------ ..... in jn64()
148 * -- - ------ - ------ - in jn64()
154 * = ----------------------- in jn64()
156 * w - ----------------- in jn64()
158 * w+h - --------- in jn64()
159 * w+2h - ... in jn64()
162 * Q(0) = w, Q(1) = w(w+h) - 1, in jn64()
163 * Q(k) = (w+k*h)*Q(k-1) - Q(k-2), in jn64()
176 q1 = w * z - _F_64(1.0); in jn64()
181 tmp = z * q1 - q0; in jn64()
186 for (t = zero, i = 2 * (n + k); i >= m; i -= 2) in jn64()
187 t = one / (i / x - t); in jn64()
188 a = t; in jn64()
202 for (i = n - 1, di = (__float64)(i + i); i > 0; i--) { in jn64()
205 b = b / x - a; in jn64()
206 a = temp; in jn64()
207 di -= two; in jn64()
210 for (i = n - 1, di = (__float64)(i + i); i > 0; i--) { in jn64()
213 b = b / x - a; in jn64()
214 a = temp; in jn64()
215 di -= two; in jn64()
218 a /= b; in jn64()
228 return -b; in jn64()
240 __float64 a, b, temp; in _MATH_ALIAS_d_id() local
260 n = -n; in _MATH_ALIAS_d_id()
261 sign = 1 - ((n & 1) << 1); in _MATH_ALIAS_d_id()
270 * Jn(x) = cos(x-(2n+1)*pi/4)*sqrt(2/x*pi) in _MATH_ALIAS_d_id()
271 * Yn(x) = sin(x-(2n+1)*pi/4)*sqrt(2/x*pi) in _MATH_ALIAS_d_id()
273 * xn=x-(2n+1)*pi/4, sqt2 = sqrt(2),then in _MATH_ALIAS_d_id()
276 * ---------------------------------- in _MATH_ALIAS_d_id()
277 * 0 s-c c+s in _MATH_ALIAS_d_id()
278 * 1 -s-c -c+s in _MATH_ALIAS_d_id()
279 * 2 -s+c -c-s in _MATH_ALIAS_d_id()
280 * 3 s+c c-s in _MATH_ALIAS_d_id()
285 temp = sin64(x) - cos64(x); in _MATH_ALIAS_d_id()
288 temp = -sin64(x) - cos64(x); in _MATH_ALIAS_d_id()
291 temp = -sin64(x) + cos64(x); in _MATH_ALIAS_d_id()
300 a = y064(x); in _MATH_ALIAS_d_id()
302 /* quit if b is -inf */ in _MATH_ALIAS_d_id()
306 b = ((__float64)(i + i) / x) * b - a; in _MATH_ALIAS_d_id()
308 a = temp; in _MATH_ALIAS_d_id()
314 return -b; in _MATH_ALIAS_d_id()