![[BACK]](/cvsweb/icons/back.gif){kind=icon}
Return to katanf.s CVS log ![[TXT]](/cvsweb/icons/text.gif){kind=icon}
![[DIR]](/cvsweb/icons/dir.gif){kind=icon}
Up to [Net/2] / Net2 / arch / tahoe / math|
Return to katanf.s CVS log | Up to [Net/2] / Net2 / arch / tahoe / math |
1.1 root 1: /*-
2: * Copyright (c) 1985 The Regents of the University of California.
3: * All rights reserved.
4: *
5: * This code is derived from software contributed to Berkeley by
6: * Computer Consoles Inc.
7: *
8: * Redistribution and use in source and binary forms, with or without
9: * modification, are permitted provided that the following conditions
10: * are met:
11: * 1. Redistributions of source code must retain the above copyright
12: * notice, this list of conditions and the following disclaimer.
13: * 2. Redistributions in binary form must reproduce the above copyright
14: * notice, this list of conditions and the following disclaimer in the
15: * documentation and/or other materials provided with the distribution.
16: * 3. All advertising materials mentioning features or use of this software
17: * must display the following acknowledgement:
18: * This product includes software developed by the University of
19: * California, Berkeley and its contributors.
20: * 4. Neither the name of the University nor the names of its contributors
21: * may be used to endorse or promote products derived from this software
22: * without specific prior written permission.
23: *
24: * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND
25: * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
26: * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
27: * ARE DISCLAIMED. IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
28: * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
29: * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
30: * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
31: * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
32: * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
33: * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
34: * SUCH DAMAGE.
35: *
36: * @(#)Katanf.s 7.1 (Berkeley) 12/6/90
37: */
38:
39: #include "../tahoe/SYS.h"
40:
41: .text
42: ENTRY(Katanf, 0)
43: tstl 4(fp) # if (arg > 0)
44: jleq 1f
45: pushl 20(fp) # hfs
46: pushl 8(fp)
47: pushl 4(fp)
48: callf $16,satan
49: ret # return(satan(arg));
50: 1: # else
51: pushl 20(fp) # hfs
52: lnd 4(fp)
53: pushd
54: callf $16,satan
55: lnf r0
56: stf r0
57: ret # return(-satan(-arg));
58:
59: ASENTRY(satan, R2|R3)
60: subl3 $60,fp,sp
61: cmpd2 4(fp),_sq2m1 # if (arg < sq2m1)
62: jgeq 1f
63: pushl 12(fp) # hfs
64: pushl 8(fp)
65: pushl 4(fp)
66: callf $16,xatan
67: ret # return(xatan(arg));
68: 1:
69: cmpd2 4(fp),_sq2p1 # else if (arg > sq2p1)
70: jgtr 9f
71: pushl 12(fp) # hfs
72: ldd one; pushd
73: ldd 4(fp); pushd
74: callf $24,_Kaddd # (arg+1.0)
75: pushl 12(fp) # hfs of _Kdivd
76: pushl r1
77: pushl r0
78: pushl 12(fp) # hfs
79: ldd one; pushd
80: ldd 4(fp); pushd
81: callf $24,_Ksubd # (arg-1.0)
82: pushl r1
83: pushl r0
84: callf $24,_Kdivd # (arg-1.0)/(arg+1.0)
85: pushl 12(fp) # hfs
86: pushl r1
87: pushl r0
88: callf $16,xatan # xatan((ag-1.0)/(arg+1.0))
89: pushl 12(fp) # hfs
90: pushl r1
91: pushl r0
92: ldd _pio4; pushd
93: callf $24,_Kaddd
94: ldd r0; cvdf; stf r0
95: ret # return(pio4+xatan((xatan(...)));
96:
97: 9:
98: pushl 12(fp)
99: ldd 4(fp); pushd
100: ldd one; pushd
101: callf $24,_Kdivd # (1.0/arg)
102: pushl 12(fp) # hfs
103: ldd r0; pushd
104: callf $16,xatan
105: # clrl -60+4(fp)
106: # movl r0,-60(fp)
107: pushl 12(fp) # hfs
108: ldd r0; pushd
109: ldd _pio2; pushd
110: callf $24,_Ksubd
111: ldd r0; cvdf; stf r0
112: ret
113:
114: ASENTRY(xatan, 0)
115: subl3 $68,fp,sp
116: pushl 12(fp) # hfs
117: ldd 4(fp); pushd; pushd
118: callf $24,_Kmuld # argsq = arg*arg;
119: ldd r0
120: std -60(fp) # argsq
121: pushl 12(fp) # hfs
122: pushd
123: ldd _p5; pushd
124: callf $24,_Kmuld # p5*argsq
125: pushl 12(fp) # hfs
126: ldd _p4; pushd
127: ldd r0; pushd
128: callf $24,_Kaddd # (p5*argsq+p4)
129: pushl 12(fp) # hfs
130: ldd -60(fp); pushd
131: ldd r0; pushd
132: callf $24,_Kmuld # (p5*argsq+p4)*argsq
133: pushl 12(fp) # hfs
134: ldd _p3; pushd
135: ldd r0; pushd
136: callf $24,_Kaddd # ((p5*argsq+p4)*argsq+p3)
137: pushl 12(fp) # hfs
138: ldd -60(fp); pushd
139: ldd r0; pushd
140: callf $24,_Kmuld # (..)*argsq
141: pushl 12(fp) # hfs
142: ldd _p2; pushd
143: ldd r0; pushd
144: callf $24,_Kaddd # (..)*argsq+p2)
145: pushl 12(fp) # hfs
146: ldd -60(fp); pushd
147: ldd r0; pushd
148: callf $24,_Kmuld # ((..)*argsq+p2)*argsq
149: pushl 12(fp) # hfs
150: ldd _p1; pushd
151: ldd r0; pushd
152: callf $24,_Kaddd # ((..)*argsq+p2)*argsq+p1)
153: pushl 12(fp) # hfs
154: ldd -60(fp); pushd
155: ldd r0; pushd
156: callf $24,_Kmuld # (..)*argsq
157: pushl 12(fp) # hfs
158: ldd _p0; pushd
159: ldd r0; pushd
160: callf $24,_Kaddd # ((..)*argsq+p1)*argsq+p0)
161: pushl 12(fp) # hfs
162: ldd 4(fp); pushd
163: ldd r0; pushd
164: callf $24,_Kmuld # (..)*arg
165: ldd r0
166: std -68(fp); cvdf; stf r0 # value
167: ret
168:
169: .data
170: .align 2
171: _sq2p1: .long 0x411A8279, 0x99FCEF31 # .double 2.414213562373095
172: _sq2m1: .long 0x3FD413CC, 0xCFE77990 # .double .41421356237309503
173: _pio2: .long 0x40C90FDA, 0xA22168C1 # .double 1.5707963267948966
174: _pio4: .long 0x40490FDA, 0xA22168C1 # .double .78539816339744829
175: _p0: .long 0x407FFFFF, 0xFD687A4B # .double .99999999939652999
176: _p1: .long 0xBFAAAAA2, 0x09F9DBF2 # .double -.3333330762079
177: _p2: .long 0x3F4CC820, 0x0670059B # .double .199982166665
178: _p3: .long 0xBF11D182, 0x6601878B # .double -.142400777317
179: _p4: .long 0x3ED88B47, 0x4EFC9AF9 # .double .10573440275
180: _p5: .long 0xBE772E4B, 0x0E689AEB # .double -.060346883
181: one: .long 0x40800000, 0x00000000 # .double 1
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.