summaryrefslogtreecommitdiff
path: root/src/pkg/math/sincos.go
diff options
context:
space:
mode:
Diffstat (limited to 'src/pkg/math/sincos.go')
-rw-r--r--src/pkg/math/sincos.go69
1 files changed, 0 insertions, 69 deletions
diff --git a/src/pkg/math/sincos.go b/src/pkg/math/sincos.go
deleted file mode 100644
index 718030319..000000000
--- a/src/pkg/math/sincos.go
+++ /dev/null
@@ -1,69 +0,0 @@
-// Copyright 2010 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
-
-// Coefficients _sin[] and _cos[] are found in pkg/math/sin.go.
-
-// Sincos returns Sin(x), Cos(x).
-//
-// Special cases are:
-// Sincos(±0) = ±0, 1
-// Sincos(±Inf) = NaN, NaN
-// Sincos(NaN) = NaN, NaN
-func Sincos(x float64) (sin, cos float64)
-
-func sincos(x float64) (sin, cos float64) {
- const (
- PI4A = 7.85398125648498535156E-1 // 0x3fe921fb40000000, Pi/4 split into three parts
- PI4B = 3.77489470793079817668E-8 // 0x3e64442d00000000,
- PI4C = 2.69515142907905952645E-15 // 0x3ce8469898cc5170,
- M4PI = 1.273239544735162542821171882678754627704620361328125 // 4/pi
- )
- // special cases
- switch {
- case x == 0:
- return x, 1 // return ±0.0, 1.0
- case IsNaN(x) || IsInf(x, 0):
- return NaN(), NaN()
- }
-
- // make argument positive
- sinSign, cosSign := false, false
- if x < 0 {
- x = -x
- sinSign = true
- }
-
- j := int64(x * M4PI) // integer part of x/(Pi/4), as integer for tests on the phase angle
- y := float64(j) // integer part of x/(Pi/4), as float
-
- if j&1 == 1 { // map zeros to origin
- j += 1
- y += 1
- }
- j &= 7 // octant modulo 2Pi radians (360 degrees)
- if j > 3 { // reflect in x axis
- j -= 4
- sinSign, cosSign = !sinSign, !cosSign
- }
- if j > 1 {
- cosSign = !cosSign
- }
-
- z := ((x - y*PI4A) - y*PI4B) - y*PI4C // Extended precision modular arithmetic
- zz := z * z
- cos = 1.0 - 0.5*zz + zz*zz*((((((_cos[0]*zz)+_cos[1])*zz+_cos[2])*zz+_cos[3])*zz+_cos[4])*zz+_cos[5])
- sin = z + z*zz*((((((_sin[0]*zz)+_sin[1])*zz+_sin[2])*zz+_sin[3])*zz+_sin[4])*zz+_sin[5])
- if j == 1 || j == 2 {
- sin, cos = cos, sin
- }
- if cosSign {
- cos = -cos
- }
- if sinSign {
- sin = -sin
- }
- return
-}