2008-03-28 14:56:47 -06:00
|
|
|
// 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.
|
|
|
|
|
2008-06-27 18:06:23 -06:00
|
|
|
package math
|
2008-03-28 14:56:47 -06:00
|
|
|
|
2009-01-15 20:11:32 -07:00
|
|
|
|
2008-03-28 14:56:47 -06:00
|
|
|
/*
|
2008-07-08 21:48:41 -06:00
|
|
|
* floating point tangent
|
2008-03-28 14:56:47 -06:00
|
|
|
*/
|
|
|
|
|
2009-03-05 14:31:01 -07:00
|
|
|
// Tan returns the tangent of x.
|
|
|
|
func Tan(x float64) float64 {
|
2009-01-15 20:11:32 -07:00
|
|
|
// Coefficients are #4285 from Hart & Cheney. (19.74D)
|
2009-10-06 20:40:35 -06:00
|
|
|
const (
|
2009-12-15 16:35:38 -07:00
|
|
|
P0 = -.1306820264754825668269611177e+5
|
|
|
|
P1 = .1055970901714953193602353981e+4
|
|
|
|
P2 = -.1550685653483266376941705728e+2
|
|
|
|
P3 = .3422554387241003435328470489e-1
|
|
|
|
P4 = .3386638642677172096076369e-4
|
|
|
|
Q0 = -.1663895238947119001851464661e+5
|
|
|
|
Q1 = .4765751362916483698926655581e+4
|
|
|
|
Q2 = -.1555033164031709966900124574e+3
|
2009-01-15 20:11:32 -07:00
|
|
|
)
|
|
|
|
|
2009-12-15 16:35:38 -07:00
|
|
|
flag := false
|
|
|
|
sign := false
|
2009-06-24 21:12:50 -06:00
|
|
|
if x < 0 {
|
2009-12-15 16:35:38 -07:00
|
|
|
x = -x
|
|
|
|
sign = true
|
2008-03-28 14:56:47 -06:00
|
|
|
}
|
2009-12-15 16:35:38 -07:00
|
|
|
x = x * (4 / Pi) /* overflow? */
|
|
|
|
var e float64
|
|
|
|
e, x = Modf(x)
|
|
|
|
i := int32(e)
|
2008-03-28 14:56:47 -06:00
|
|
|
|
2009-11-09 22:23:52 -07:00
|
|
|
switch i & 3 {
|
2008-03-28 14:56:47 -06:00
|
|
|
case 1:
|
2009-12-15 16:35:38 -07:00
|
|
|
x = 1 - x
|
|
|
|
flag = true
|
2008-03-28 14:56:47 -06:00
|
|
|
|
|
|
|
case 2:
|
2009-12-15 16:35:38 -07:00
|
|
|
sign = !sign
|
|
|
|
flag = true
|
2008-03-28 14:56:47 -06:00
|
|
|
|
|
|
|
case 3:
|
2009-12-15 16:35:38 -07:00
|
|
|
x = 1 - x
|
|
|
|
sign = !sign
|
2008-03-28 14:56:47 -06:00
|
|
|
}
|
|
|
|
|
2009-12-15 16:35:38 -07:00
|
|
|
xsq := x * x
|
|
|
|
temp := ((((P4*xsq+P3)*xsq+P2)*xsq+P1)*xsq + P0) * x
|
|
|
|
temp = temp / (((xsq+Q2)*xsq+Q1)*xsq + Q0)
|
2008-03-28 14:56:47 -06:00
|
|
|
|
|
|
|
if flag {
|
2009-06-24 21:12:50 -06:00
|
|
|
if temp == 0 {
|
2009-11-09 13:07:39 -07:00
|
|
|
panic(NaN())
|
2008-03-28 14:56:47 -06:00
|
|
|
}
|
2009-12-15 16:35:38 -07:00
|
|
|
temp = 1 / temp
|
2008-03-28 14:56:47 -06:00
|
|
|
}
|
|
|
|
if sign {
|
2009-11-09 13:07:39 -07:00
|
|
|
temp = -temp
|
2008-03-28 14:56:47 -06:00
|
|
|
}
|
2009-12-15 16:35:38 -07:00
|
|
|
return temp
|
2008-03-28 14:56:47 -06:00
|
|
|
}
|