![[BACK]](/cvsweb/icons/back.gif){kind=icon}
Return to atan.s CVS log ![[TXT]](/cvsweb/icons/text.gif){kind=icon}
![[DIR]](/cvsweb/icons/dir.gif){kind=icon}
Up to [Power 6/32 Unix Tahoe 4.2BSD] / cci / usr / src / usr.lib / libm|
Return to atan.s CVS log | Up to [Power 6/32 Unix Tahoe 4.2BSD] / cci / usr / src / usr.lib / libm |
1.1 root 1: /* @(#)atan.c 4.1/4.2 10/31/84 CCI-CPG */
2:
3:
4: /*
5: Double-precision arctangent function recoded with
6: new algorithm polynomial at CCI-CPG.
7: Coefficients are #4969 from Hart & etc. (17.92D)
8:
9: atan returns the value of the arctangent of its
10: argument in the range [-pi/2,pi/2].
11:
12: atan2 returns the arctangent of arg1/arg2
13: in the range [-pi,pi].
14:
15: there are no error returns.
16:
17: */
18: .data
19: .align 2
20: _sq2p1: .long 0x411a8279, 0x99fcef32 # .double 2.414213562373095
21: .align 2
22: _sq2m1: .long 0x3fd413cc, 0xcfe77992 # .double 0.41421356237309505
23: .align 2
24: _pio2: .long 0x40c90fda, 0xa22168c2 # .double 1.5707963267948966
25: .align 2
26: _pio4: .long 0x40490fda, 0xa22168c2 # .double 0.78539816339744831
27: .align 2
28: _p0: .long 0x40800000, 0x00000000 # .double 1
29: .align 2
30: _p1: .long 0xbfaaaaaa, 0xaaaaa984 # .double -0.33333333333333129
31: .align 2
32: _p2: .long 0x3f4ccccc, 0xccca47cb # .double 0.19999999999942712
33: .align 2
34: _p3: .long 0xbf124924, 0x91341847 # .double -0.14285714279414943
35: .align 2
36: _p4: .long 0x3ee38e38, 0x685c2f79 # .double 0.11111110752565144
37: .align 2
38: _p5: .long 0xbeba2e7b, 0x6e11c9ad # .double -0.090908970161966385
39: .align 2
40: _p6: .long 0x3e9d887d, 0x2d47ecb2 # .double 0.076920488304486365
41: .align 2
42: _p7: .long 0xbe887541, 0x84d86747 # .double -0.066629897944085704
43: .align 2
44: _p8: .long 0x3e6f7fbd, 0x60c7567d # .double 0.05847143150136491
45: .align 2
46: _p9: .long 0xbe4e48e7, 0x446d54d6 # .double -0.0503624948385852
47: .align 2
48: _p10: .long 0x3e1b9ee6, 0x004c28d5 # .double 0.03799333423809
49: .align 2
50: _p11: .long 0xbd921adf, 0x2d9e4583 # .double -0.017835079092898
51: .text
52: LL0: .align 1
53: .globl _atan
54: .set L28,0x0
55: .data
56: .text
57: _atan: .word L28
58: tstl 4(fp)
59: jleq L33
60: pushl 8(fp)
61: pushl 4(fp)
62: callf $12,_satan
63: ret
64: L33: lnd 4(fp)
65: pushd
66: callf $12,_satan
67: lnd r0
68: std r0
69: ret
70: .align 1
71: .globl _atan2
72: .set L35,0xc
73: .data
74: .text
75: .data
76: .align 2
77: L51: .long 0x40800000, 0x00000000 # .double 1
78: .text
79: .set L44,0xc
80: .data
81: .text
82: .set L54,0x0
83: .data
84: .text
85: _atan2: .word L35
86: subl3 $60,fp,sp
87: ldd 4(fp)
88: addd 12(fp)
89: cmpd 4(fp)
90: jneq L39
91: tstl 4(fp)
92: jlss L40
93: movl _pio2+4,r1
94: movl _pio2,r0
95: ret
96: L40: lnd _pio2
97: jbr L2000000
98: L2000003:
99: lnd 4(fp)
100: divd 12(fp)
101: pushd
102: callf $12,_satan
103: movl r1,-56(fp)
104: movl r0,-60(fp)
105: ldd _pio2
106: addd _pio2
107: L2000001:
108: subd -60(fp)
109: jbr L2000000
110: L2000005:
111: tstl 4(fp)
112: jgeq L2000003
113: ldd 4(fp)
114: divd 12(fp)
115: pushd
116: callf $12,_satan
117: lnd _pio2
118: subd _pio2
119: addd r0
120: L2000000:
121: std r0
122: ret
123: L39: tstl 12(fp)
124: jlss L2000005
125: tstl 4(fp)
126: jleq L43
127: ldd 4(fp)
128: divd 12(fp)
129: pushd
130: callf $12,_satan
131: ret
132: L43: lnd 4(fp)
133: divd 12(fp)
134: pushd
135: callf $12,_satan
136: lnd r0
137: jbr L2000000
138: .align 1
139: _satan: .word L44
140: subl3 $60,fp,sp
141: cmpd2 4(fp),_sq2m1
142: jgeq L49
143: pushl 8(fp)
144: pushl 4(fp)
145: callf $12,_xatan
146: ret
147: L49: cmpd2 4(fp),_sq2p1
148: jleq L50
149: ldd L51
150: divd 4(fp)
151: pushd
152: callf $12,_xatan
153: movl r1,-56(fp)
154: movl r0,-60(fp)
155: ldd _pio2
156: jbr L2000001
157: L50: ldd 4(fp)
158: subd L51
159: std r0
160: ldd 4(fp)
161: addd L51
162: std r2
163: ldd r0
164: divd r2
165: pushd
166: callf $12,_xatan
167: ldd r0
168: addd _pio4
169: jbr L2000000
170: .align 1
171: _xatan: .word L54
172: subl3 $68,fp,sp
173: ldd 4(fp)
174: muld 4(fp)
175: std -60(fp)
176: ldd _p11
177: muld -60(fp)
178: addd _p10
179: muld -60(fp)
180: addd _p9
181: muld -60(fp)
182: addd _p8
183: muld -60(fp)
184: addd _p7
185: muld -60(fp)
186: addd _p6
187: muld -60(fp)
188: addd _p5
189: muld -60(fp)
190: addd _p4
191: muld -60(fp)
192: addd _p3
193: muld -60(fp)
194: addd _p2
195: muld -60(fp)
196: addd _p1
197: muld -60(fp)
198: addd _p0
199: muld 4(fp)
200: std -68(fp)
201: movl -64(fp),r1
202: movl -68(fp),r0
203: ret
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.