diff options
Diffstat (limited to 'src/lib/math/sin.go')
-rw-r--r-- | src/lib/math/sin.go | 36 |
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); } |