|
Line 0
Link Here
|
|
|
1 |
--- src/lib/dev/ltra/ltramisc.c.orig Wed Jan 5 03:25:41 2005 |
| 2 |
+++ src/lib/dev/ltra/ltramisc.c Thu Jul 20 23:59:44 2006 |
| 3 |
@@ -965,10 +965,12 @@ |
| 4 |
double ax = fabs(x); |
| 5 |
if (ax < 3.75) { |
| 6 |
double y = x/3.75; |
| 7 |
+ double ans1; |
| 8 |
+ double ans2; |
| 9 |
y *= y; |
| 10 |
- double ans1 = 1.0+y*(3.5156229+y*(3.0899424+y*(1.2067492 |
| 11 |
+ ans1 = 1.0+y*(3.5156229+y*(3.0899424+y*(1.2067492 |
| 12 |
+y*(0.2659732+y*(0.360768e-1+y*0.45813e-2))))); |
| 13 |
- double ans2 = ax*(0.5+y*(0.87890594+y*(0.51498869+y*(0.15084934 |
| 14 |
+ ans2 = ax*(0.5+y*(0.87890594+y*(0.51498869+y*(0.15084934 |
| 15 |
+y*(0.2658733e-1+y*(0.301532e-2+y*0.32411e-3)))))); |
| 16 |
if (x > 0) |
| 17 |
ans1 += ans2; |
| 18 |
@@ -978,11 +980,12 @@ |
| 19 |
} |
| 20 |
else { |
| 21 |
double y = 3.75/ax; |
| 22 |
+ double ans2; |
| 23 |
double ans1 = (0.39894228+y*(0.1328592e-1 |
| 24 |
+y*(0.225319e-2+y*(-0.157565e-2+y*(0.916281e-2 |
| 25 |
+y*(-0.2057706e-1+y*(0.2635537e-1+y*(-0.1647633e-1 |
| 26 |
+y*0.392377e-2)))))))); |
| 27 |
- double ans2 = 0.2282967e-1+y*(-0.2895312e-1+y*(0.1787654e-1 |
| 28 |
+ ans2 = 0.2282967e-1+y*(-0.2895312e-1+y*(0.1787654e-1 |
| 29 |
-y*0.420059e-2)); |
| 30 |
ans2 = 0.39894228+y*(-0.3988024e-1+y*(-0.362018e-2 |
| 31 |
+y*(0.163801e-2+y*(-0.1031555e-1+y*ans2)))); |