blob: cc8e39fbcb15a5d3911857ee641b6a0d8caa4256 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
|
// Copyright 2009 The Go Authors. All rights reserved.
// Use of this source code is governed by a BSD-style
// license that can be found in the LICENSE file.
package math
const
(
p0 = .1357884097877375669092680e8;
p1 = -.4942908100902844161158627e7;
p2 = .4401030535375266501944918e6;
p3 = -.1384727249982452873054457e5;
p4 = .1459688406665768722226959e3;
q0 = .8644558652922534429915149e7;
q1 = .4081792252343299749395779e6;
q2 = .9463096101538208180571257e4;
q3 = .1326534908786136358911494e3;
piu2 = .6366197723675813430755350e0; // 2/pi
)
func
sinus(arg float64, quad int) float64
{
var e, f, ysq, x, y, temp1, temp2 float64;
var k int32;
x = arg;
if(x < 0) {
x = -x;
quad = quad+2;
}
x = x * piu2; /* underflow? */
if x > 32764 {
e,y = sys.modf(x);
e = e + float64(quad);
temp1,f = sys.modf(0.25*e);
quad = int(e - 4*f);
} else {
k = int32(x);
y = x - float64(k);
quad = (quad + int(k)) & 3;
}
if quad&1 != 0 {
y = 1-y;
}
if quad > 1 {
y = -y;
}
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;
}
export func
cos(arg float64) float64
{
if arg < 0 {
arg = -arg;
}
return sinus(arg, 1);
}
export func
sin(arg float64) float64
{
return sinus(arg, 0);
}
|
