123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237 |
- package math
- var _sin = [...]float64{
- 1.58962301576546568060E-10,
- -2.50507477628578072866E-8,
- 2.75573136213857245213E-6,
- -1.98412698295895385996E-4,
- 8.33333333332211858878E-3,
- -1.66666666666666307295E-1,
- }
- var _cos = [...]float64{
- -1.13585365213876817300E-11,
- 2.08757008419747316778E-9,
- -2.75573141792967388112E-7,
- 2.48015872888517045348E-5,
- -1.38888888888730564116E-3,
- 4.16666666666665929218E-2,
- }
- func libc_cos(float64) float64
- func Cos(x float64) float64 {
- return libc_cos(x)
- }
- func cos(x float64) float64 {
- const (
- PI4A = 7.85398125648498535156E-1
- PI4B = 3.77489470793079817668E-8
- PI4C = 2.69515142907905952645E-15
- M4PI = 1.273239544735162542821171882678754627704620361328125
- )
-
- switch {
- case IsNaN(x) || IsInf(x, 0):
- return NaN()
- }
-
- sign := false
- if x < 0 {
- x = -x
- }
- j := int64(x * M4PI)
- y := float64(j)
-
- if j&1 == 1 {
- j += 1
- y += 1
- }
- j &= 7
- if j > 3 {
- j -= 4
- sign = !sign
- }
- if j > 1 {
- sign = !sign
- }
- z := ((x - y*PI4A) - y*PI4B) - y*PI4C
- zz := z * z
- if j == 1 || j == 2 {
- y = z + z*zz*((((((_sin[0]*zz)+_sin[1])*zz+_sin[2])*zz+_sin[3])*zz+_sin[4])*zz+_sin[5])
- } else {
- y = 1.0 - 0.5*zz + zz*zz*((((((_cos[0]*zz)+_cos[1])*zz+_cos[2])*zz+_cos[3])*zz+_cos[4])*zz+_cos[5])
- }
- if sign {
- y = -y
- }
- return y
- }
- func libc_sin(float64) float64
- func Sin(x float64) float64 {
- return libc_sin(x)
- }
- func sin(x float64) float64 {
- const (
- PI4A = 7.85398125648498535156E-1
- PI4B = 3.77489470793079817668E-8
- PI4C = 2.69515142907905952645E-15
- M4PI = 1.273239544735162542821171882678754627704620361328125
- )
-
- switch {
- case x == 0 || IsNaN(x):
- return x
- case IsInf(x, 0):
- return NaN()
- }
-
- sign := false
- if x < 0 {
- x = -x
- sign = true
- }
- j := int64(x * M4PI)
- y := float64(j)
-
- if j&1 == 1 {
- j += 1
- y += 1
- }
- j &= 7
-
- if j > 3 {
- sign = !sign
- j -= 4
- }
- z := ((x - y*PI4A) - y*PI4B) - y*PI4C
- zz := z * z
- if j == 1 || j == 2 {
- y = 1.0 - 0.5*zz + zz*zz*((((((_cos[0]*zz)+_cos[1])*zz+_cos[2])*zz+_cos[3])*zz+_cos[4])*zz+_cos[5])
- } else {
- y = z + z*zz*((((((_sin[0]*zz)+_sin[1])*zz+_sin[2])*zz+_sin[3])*zz+_sin[4])*zz+_sin[5])
- }
- if sign {
- y = -y
- }
- return y
- }
|