d0p1/ack
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

fixed some errors

Browse source
This commit is contained in:
ceriel 1989-05-18 16:48:28 +00:00
parent 9f7ee118f7
commit e431cec90a
2 changed files with 10 additions and 10 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

18
lang/cem/libcc/math/asin.c
View file

@ -45,25 +45,25 @@ asin_acos(x, cosfl)
errno = EDOM; errno = EDOM;
return 0; return 0;
} }
g = 0.5 - 0.5 * y; g = 0.5 - 0.5 * x;
y = - sqrt(g); x = - sqrt(g);
y += y; x += x;
} }
else { else {
/* ??? avoid underflow ??? */ /* ??? avoid underflow ??? */
g = y * y; g = x * x;
} }
y += y * g * POLYNOM4(g, x) / POLYNOM5(g, y); x += x * g * POLYNOM4(g, p) / POLYNOM5(g, q);
if (i == 1) { if (i == 1) {
if (cosfl == 0 || ! negative) { if (cosfl == 0 || ! negative) {
y = (y + M_PI_4) + M_PI_4; x = (x + M_PI_4) + M_PI_4;
} }
else if (cosfl && negative) { else if (cosfl && negative) {
y = (y + M_PI_2) + M_PI_2; x = (x + M_PI_2) + M_PI_2;
} }
} }
if (! cosfl && negative) y = -y; if (! cosfl && negative) x = -x;
return y; return x;
} }
double double

2
lang/cem/libcc/math/tan.c
View file

@ -68,6 +68,6 @@ tan(x)
y = x * x; y = x * x;
x += x * y * POLYNOM2(y, p+1); x += x * y * POLYNOM2(y, p+1);
y = POLYNOM4(y, q); y = POLYNOM4(y, q);
if (neg) x = -x; if (negative) x = -x;
return invert ? -y/x : x/y; return invert ? -y/x : x/y;
} }