summaryrefslogtreecommitdiff
path: root/src/lib/math/sin.go
diff options
context:
space:
mode:
Diffstat (limited to 'src/lib/math/sin.go')
-rw-r--r--src/lib/math/sin.go36
1 files changed, 18 insertions, 18 deletions
diff --git a/src/lib/math/sin.go b/src/lib/math/sin.go
index 57de55913..077506b30 100644
--- a/src/lib/math/sin.go
+++ b/src/lib/math/sin.go
@@ -9,32 +9,32 @@ package math
*/
const
(
- p0 = .1357884097877375669092680e8;
- p1 = -.4942908100902844161158627e7;
- p2 = .4401030535375266501944918e6;
- p3 = -.1384727249982452873054457e5;
- p4 = .1459688406665768722226959e3;
- q0 = .8644558652922534429915149e7;
- q1 = .4081792252343299749395779e6;
- q2 = .9463096101538208180571257e4;
- q3 = .1326534908786136358911494e3;
+ sp0 = .1357884097877375669092680e8;
+ sp1 = -.4942908100902844161158627e7;
+ sp2 = .4401030535375266501944918e6;
+ sp3 = -.1384727249982452873054457e5;
+ sp4 = .1459688406665768722226959e3;
+ sq0 = .8644558652922534429915149e7;
+ sq1 = .4081792252343299749395779e6;
+ sq2 = .9463096101538208180571257e4;
+ sq3 = .1326534908786136358911494e3;
- piu2 = .6366197723675813430755350e0; // 2/pi
+ spiu2 = .6366197723675813430755350e0; // 2/pi
)
-func Sinus(arg float64, quad int) float64 {
+func sinus(arg float64, quad int) float64 {
x := arg;
if(x < 0) {
x = -x;
quad = quad+2;
}
- x = x * piu2; /* underflow? */
+ x = x * spiu2; /* underflow? */
var y float64;
if x > 32764 {
var e float64;
e, y = sys.modf(x);
e = e + float64(quad);
- temp1, f := sys.modf(0.25*e);
+ temsp1, f := sys.modf(0.25*e);
quad = int(e - 4*f);
} else {
k := int32(x);
@@ -50,18 +50,18 @@ func Sinus(arg float64, quad int) float64 {
}
ysq := y*y;
- temp1 := ((((p4*ysq+p3)*ysq+p2)*ysq+p1)*ysq+p0)*y;
- temp2 := ((((ysq+q3)*ysq+q2)*ysq+q1)*ysq+q0);
- return temp1/temp2;
+ temsp1 := ((((sp4*ysq+sp3)*ysq+sp2)*ysq+sp1)*ysq+sp0)*y;
+ temsp2 := ((((ysq+sq3)*ysq+sq2)*ysq+sq1)*ysq+sq0);
+ return temsp1/temsp2;
}
export func Cos(arg float64) float64 {
if arg < 0 {
arg = -arg;
}
- return Sinus(arg, 1);
+ return sinus(arg, 1);
}
export func Sin(arg float64) float64 {
- return Sinus(arg, 0);
+ return sinus(arg, 0);
}