1 #Region "Microsoft.VisualBasic::fbb2ea48f2dd4d815304d189f668b4b9, Microsoft.VisualBasic.Core\Extensions\Math\Trigonometric\Arctan.vb"
2
3     ' Author:
4     
5     '       asuka (amethyst.asuka@gcmodeller.org)
6     '       xie (genetics@smrucc.org)
7     '       xieguigang (xie.guigang@live.com)
8     
9     ' Copyright (c) 2018 GPL3 Licensed
10     
11     
12     ' GNU GENERAL PUBLIC LICENSE (GPL3)
13     
14     
15     ' This program is free software: you can redistribute it and/or modify
16     ' it under the terms of the GNU General Public License as published by
17     ' the Free Software Foundation, either version 3 of the License, or
18     ' (at your option) any later version.
19     
20     ' This program is distributed in the hope that it will be useful,
21     ' but WITHOUT ANY WARRANTY; without even the implied warranty of
22     ' MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
23     ' GNU General Public License for more details.
24     
25     ' You should have received a copy of the GNU General Public License
26     ' along with this program. If not, see <http://www.gnu.org/licenses/>.
27
28
29
30     ' /********************************************************************************/
31
32     ' Summaries:
33
34     '     Module Trigonometric
35     
36     '         Function: Atn
37     
38     
39     ' /********************************************************************************/
40
41 #End Region
42
43 Imports sys = System.Math
44
45 Namespace Math
46
47     Partial Module Trigonometric
48
49         Const HalfPI# = PI / 2
50         Const n1# = 1.0R
51         Const n05# = 0.5R
52
53         ''' <summary>
54         ''' 通过这个参数来控制计算精度,这个参数值越大,计算精度越高
55         ''' </summary>
56         Public AtanPrecise% = 500
57
58         ''' <summary>
59         ''' Taylor Atan
60         ''' </summary>
61         ''' <param name="x"></param>
62         ''' <returns></returns>
63         ''' <remarks>Atan测试没有问题</remarks>>
64         Public Function Atn(x#) As Double
65             If x# = 1.0 Then
66                 Return PI / 4
67             End If
68             If sys.Sign(x) = sys.Sign(-1) Then
69                 Return -Atn(-x)
70             End If
71             If x > n1 Then
72                 Return HalfPI - Atn(1 / x)
73             End If
74             If x > n05 Then
75                 Return Atn(n1) + Atn((x - n1) / (1 + x))
76             Else
77                 Dim xPow2 As Double = x * x
78                 Dim n__1 As Integer = AtanPrecise
79                 Dim y As Double = 1 / (2 * n__1 + 1)
80                 Dim i As Integer = AtanPrecise
81
82                 While i > 0
83                     y = (1 / (2 * n__1 - 1)) - (xPow2 * y)
84                     i -= 1
85                     n__1 -= 1
86                 End While
87
88                 Return x * y
89             End If
90         End Function
91     End Module
92 End Namespace