![[BACK]](/cvsweb/icons/back.gif){kind=icon}
Return to ksinfcosf.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 ksinfcosf.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: * @(#)Ksinfcosf.s 7.1 (Berkeley) 12/6/90
37: */
38:
39: #include "../tahoe/SYS.h"
40:
41: .text
42: ENTRY(Kcosf, 0)
43: tstl 4(fp) # if (arg < 0)
44: jgeq L23
45: lnd 4(fp)
46: std 4(fp) # arg = -arg;
47: L23: pushl 20(fp) # hfs
48: pushl $1
49: pushl 8(fp)
50: pushl 4(fp)
51: callf $20,_sinus
52: ret # return(sinus(arg,1));
53:
54: ENTRY(Ksinf, 0)
55: pushl 20(fp) # hfs
56: pushl $0
57: pushl 8(fp)
58: pushl 4(fp)
59: callf $20,_sinus
60: ret # return(sinus(arg, 0));
61:
62: ENTRY(sinus, 0)
63: subl3 $112,fp,sp
64: movl 8(fp),-80(fp)
65: movl 4(fp),-84(fp) # x = arg;
66: jgeq L34 # if (x < 0) {
67: lnd -84(fp)
68: std -84(fp) # x = -x;
69: addl2 $2,12(fp) # quad = quad + 2;}
70: L34: pushl 16(fp)
71: ldd _twoopi
72: pushd
73: ldd -84(fp)
74: pushd
75: callf $24,_Kmuld
76: movl r1,-88(fp)
77: movl r0,-84(fp) # std -84(fp) x = x * twoopi;
78: cmpd2 -84(fp),L36 # if (x > 32764) {
79: jleq L35
80: pushl 16(fp) # hfs
81: subl3 $60,fp,-(sp) # &e
82: pushl -80(fp)
83: pushl -84(fp) # x
84: callf $20,_Kmodf
85: clrl -92+4(fp)
86: movl r0,-92(fp) # y = modf(x, &e);
87: pushl 16(fp) # hfs
88: pushl $0
89: pushl 12(fp) # quad
90: pushl $0
91: pushl $0 # dummy acc
92: callf $24,_Kcvtld
93: pushl 16(fp) # hfs
94: pushl r1
95: pushl r0 # pushd (double)quad
96: ldd -60(fp)
97: pushd # &e
98: callf $24,_Kaddd
99: movl r1,-56(fp)
100: movl r0,-60(fp) # std -60(fp) e = e + quad;
101: pushl 16(fp) # hfs
102: subl3 $68,fp,-(sp) # &f
103:
104: pushl 16(fp) # hfs
105: ldd -60(fp)
106: pushd # &e
107: ldd L37
108: pushd # 0.25
109: callf $24,_Kmuld
110:
111: pushl r1
112: pushl r0 # pushd 0.25 * e
113: callf $20,_Kmodf # modf(0.25 * e, &f);
114:
115: pushl 16(fp) # hfs
116: movl r1,-64(fp)
117: movl r0,-68(fp) # ldd -68(fp) &f
118: pushl r1
119: pushl r0 # pushd
120: ldd L38
121: pushd
122: callf $24,_Kmuld # 4*f
123: # std r0
124: pushl 16(fp) # hfs
125: pushl r1
126: pushl r0 # pushd
127: ldd -60(fp)
128: pushd
129: callf $24,_Ksubd # e - 4 * f
130: ldd r0 # load acc with result of Ksubd
131: cvdl 12(fp) # quad = (e - 4 * f); }
132: jbr L39
133: L35: ldd -84(fp) # else {
134: cvdl -96(fp) # k = x;
135: pushl 16(fp) # hfs
136: pushl $0
137: pushl -96(fp)
138: pushl $0
139: pushl $0 # acc
140: callf $24,_Kcvtld
141: # std r0
142: pushl 16(fp) # hfs
143: pushl r1
144: pushl r0 # pushd (double)k
145: ldd -84(fp)
146: pushd # x
147: callf $24,_Ksubd
148: movl r1,-88(fp)
149: movl r0,-92(fp) # std -92(fp) y = x - k;
150: addl3 -96(fp),12(fp),r0
151: andl3 $3,r0,12(fp) # quad = (quad + k) & 03; }
152: L39: bitl $1,12(fp) # if (quad & 01)
153: jeql L40
154: pushl 16(fp) # hfs
155: ldd -92(fp)
156: pushd # y
157: ldd L41
158: pushd # 1.0
159: callf $24,_Ksubd
160: movl r1,-88(fp)
161: movl r0,-92(fp) # std -92(fp) y = 1 - y;
162: L40: cmpl 12(fp),$1 # if (quad > 1)
163: jleq L42
164: lnd -92(fp)
165: std -92(fp) # y = -y;
166: L42: pushl 16(fp) # hfs
167: ldd -92(fp)
168: pushd
169: pushd
170: callf $24,_Kmuld # y*y
171: movl r1,-72(fp)
172: movl r0,-76(fp) # std -76(fp) ysq=y*y;
173:
174: pushl 16(fp) # hfs
175: ldd -76(fp)
176: pushd
177: ldd _p4
178: pushd
179: callf $24,_Kmuld # p4*ysq
180:
181: pushl 16(fp) # hfs
182: pushl _p3+4
183: pushl _p3
184: pushl r1
185: pushl r0 # pushd
186: callf $24,_Kaddd # p4*ysq+p3
187:
188: pushl 16(fp) # hfs
189: pushl -72(fp)
190: pushl -76(fp)
191: pushl r1
192: pushl r0 # pushd
193: callf $24,_Kmuld # (p4*ysq+p3)*ysq
194:
195: pushl 16(fp) # hfs
196: pushl _p2+4
197: pushl _p2
198: pushl r1
199: pushl r0 # pushd
200: callf $24,_Kaddd # (p4*ysq+p3)*ysq+p2
201:
202: pushl 16(fp) # hfs
203: pushl -72(fp)
204: pushl -76(fp)
205: pushl r1
206: pushl r0 # pushd
207: callf $24,_Kmuld # ((p4*ysq+p3)*ysq+p2)*ysq
208:
209: pushl 16(fp) # hfs
210: pushl _p1+4
211: pushl _p1
212: pushl r1
213: pushl r0 # pushd
214: callf $24,_Kaddd # ((p4*ysq+p3)*ysq+p2)*ysq+p1
215:
216: pushl 16(fp) # hfs
217: pushl -72(fp)
218: pushl -76(fp)
219: pushl r1
220: pushl r0 # pushd
221: callf $24,_Kmuld # (((p4*ysq+p3)*ysq+p2)*ysq+p1)*ysq
222:
223: pushl 16(fp) # hfs
224: pushl _p0+4
225: pushl _p0
226: pushl r1
227: pushl r0 # pushd
228: callf $24,_Kaddd # (((p4*ysq+p3)*ysq+p2)*ysq+p1)*ysq+p0
229:
230: pushl 16(fp) # hfs
231: pushl -88(fp)
232: pushl -92(fp)
233: pushl r1
234: pushl r0 # pushd
235: callf $24,_Kmuld # ((((p4*ysq+p3)*ysq+p2)*ysq+p1)*ysq+p0)*y;
236:
237: movl r1,-100(fp)
238: movl r0,-104(fp) # std -104(fp) temp1
239: ldd r0
240: cvdf
241: stf r0
242: ret#1
243:
244: .data
245: .align 2
246: _twoopi:.long 0x4022F983, 0x6E4E4415 # .double .63661977236758134
247: _p0: .long 0x40C90FDA, 0x90304197 # .double 1.57079631844
248: _p1: .long 0xC0255DE0, 0xB36CEE75 # .double -.645963710599
249: _p2: .long 0x3EA33457, 0xA736E807 # .double .079689678945999999
250: _p3: .long 0xBC992665, 0x5E9A6554 # .double -.0046737666099999999
251: _p4: .long 0x3A1ED7FA, 0xCC54924E # .double .000151485129
252: L36: .long 0x47FFF800, 0x00000000 # .double 32764
253: L37: .long 0x3F800000, 0x00000000 # .double .25
254: L38: .long 0x41800000, 0x00000000 # .double 4
255: L41: .long 0x40800000, 0x00000000 # .double 1
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.