summaryrefslogtreecommitdiff
path: root/src/pkg/math/sin.go
diff options
context:
space:
mode:
Diffstat (limited to 'src/pkg/math/sin.go')
-rw-r--r--src/pkg/math/sin.go54
1 files changed, 27 insertions, 27 deletions
diff --git a/src/pkg/math/sin.go b/src/pkg/math/sin.go
index 25f676355..e17daf688 100644
--- a/src/pkg/math/sin.go
+++ b/src/pkg/math/sin.go
@@ -8,32 +8,32 @@ package math
func sinus(x float64, quad int) float64 {
// Coefficients are #3370 from Hart & Cheney (18.80D).
const (
- P0 = .1357884097877375669092680e8;
- P1 = -.4942908100902844161158627e7;
- P2 = .4401030535375266501944918e6;
- P3 = -.1384727249982452873054457e5;
- P4 = .1459688406665768722226959e3;
- Q0 = .8644558652922534429915149e7;
- Q1 = .4081792252343299749395779e6;
- Q2 = .9463096101538208180571257e4;
- Q3 = .1326534908786136358911494e3;
+ P0 = .1357884097877375669092680e8
+ P1 = -.4942908100902844161158627e7
+ P2 = .4401030535375266501944918e6
+ P3 = -.1384727249982452873054457e5
+ P4 = .1459688406665768722226959e3
+ Q0 = .8644558652922534429915149e7
+ Q1 = .4081792252343299749395779e6
+ Q2 = .9463096101538208180571257e4
+ Q3 = .1326534908786136358911494e3
)
if x < 0 {
- x = -x;
- quad = quad + 2;
+ x = -x
+ quad = quad + 2
}
- x = x * (2 / Pi); /* underflow? */
- var y float64;
+ x = x * (2 / Pi) /* underflow? */
+ var y float64
if x > 32764 {
- var e float64;
- e, y = Modf(x);
- e = e + float64(quad);
- f, _ := Modf(0.25 * e);
- quad = int(e - 4*f);
+ var e float64
+ e, y = Modf(x)
+ e = e + float64(quad)
+ f, _ := Modf(0.25 * e)
+ quad = int(e - 4*f)
} else {
- k := int32(x);
- y = x - float64(k);
- quad = (quad + int(k)) & 3;
+ k := int32(x)
+ y = x - float64(k)
+ quad = (quad + int(k)) & 3
}
if quad&1 != 0 {
@@ -43,10 +43,10 @@ func sinus(x float64, quad int) float64 {
y = -y
}
- yy := y * y;
- temp1 := ((((P4*yy+P3)*yy+P2)*yy+P1)*yy + P0) * y;
- temp2 := ((((yy+Q3)*yy+Q2)*yy+Q1)*yy + Q0);
- return temp1 / temp2;
+ yy := y * y
+ temp1 := ((((P4*yy+P3)*yy+P2)*yy+P1)*yy + P0) * y
+ temp2 := ((((yy+Q3)*yy+Q2)*yy+Q1)*yy + Q0)
+ return temp1 / temp2
}
// Cos returns the cosine of x.
@@ -54,8 +54,8 @@ func Cos(x float64) float64 {
if x < 0 {
x = -x
}
- return sinus(x, 1);
+ return sinus(x, 1)
}
// Sin returns the sine of x.
-func Sin(x float64) float64 { return sinus(x, 0) }
+func Sin(x float64) float64 { return sinus(x, 0) }