![[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 [CSRG BSD Unix] / 43BSDReno / sys / tahoemath|
Return to Ksinfcosf.s CVS log | Up to [CSRG BSD Unix] / 43BSDReno / sys / tahoemath |
1.1 root 1: /* Ksinfcosf.s 1.3 86/01/05 */
2:
3: #include "../tahoe/SYS.h"
4:
5: .text
6: ENTRY(Kcosf, 0)
7: tstl 4(fp) # if (arg < 0)
8: jgeq L23
9: lnd 4(fp)
10: std 4(fp) # arg = -arg;
11: L23: pushl 20(fp) # hfs
12: pushl $1
13: pushl 8(fp)
14: pushl 4(fp)
15: callf $20,_sinus
16: ret # return(sinus(arg,1));
17:
18: ENTRY(Ksinf, 0)
19: pushl 20(fp) # hfs
20: pushl $0
21: pushl 8(fp)
22: pushl 4(fp)
23: callf $20,_sinus
24: ret # return(sinus(arg, 0));
25:
26: ENTRY(sinus, 0)
27: subl3 $112,fp,sp
28: movl 8(fp),-80(fp)
29: movl 4(fp),-84(fp) # x = arg;
30: jgeq L34 # if (x < 0) {
31: lnd -84(fp)
32: std -84(fp) # x = -x;
33: addl2 $2,12(fp) # quad = quad + 2;}
34: L34: pushl 16(fp)
35: ldd _twoopi
36: pushd
37: ldd -84(fp)
38: pushd
39: callf $24,_Kmuld
40: movl r1,-88(fp)
41: movl r0,-84(fp) # std -84(fp) x = x * twoopi;
42: cmpd2 -84(fp),L36 # if (x > 32764) {
43: jleq L35
44: pushl 16(fp) # hfs
45: subl3 $60,fp,-(sp) # &e
46: pushl -80(fp)
47: pushl -84(fp) # x
48: callf $20,_Kmodf
49: clrl -92+4(fp)
50: movl r0,-92(fp) # y = modf(x, &e);
51: pushl 16(fp) # hfs
52: pushl $0
53: pushl 12(fp) # quad
54: pushl $0
55: pushl $0 # dummy acc
56: callf $24,_Kcvtld
57: pushl 16(fp) # hfs
58: pushl r1
59: pushl r0 # pushd (double)quad
60: ldd -60(fp)
61: pushd # &e
62: callf $24,_Kaddd
63: movl r1,-56(fp)
64: movl r0,-60(fp) # std -60(fp) e = e + quad;
65: pushl 16(fp) # hfs
66: subl3 $68,fp,-(sp) # &f
67:
68: pushl 16(fp) # hfs
69: ldd -60(fp)
70: pushd # &e
71: ldd L37
72: pushd # 0.25
73: callf $24,_Kmuld
74:
75: pushl r1
76: pushl r0 # pushd 0.25 * e
77: callf $20,_Kmodf # modf(0.25 * e, &f);
78:
79: pushl 16(fp) # hfs
80: movl r1,-64(fp)
81: movl r0,-68(fp) # ldd -68(fp) &f
82: pushl r1
83: pushl r0 # pushd
84: ldd L38
85: pushd
86: callf $24,_Kmuld # 4*f
87: # std r0
88: pushl 16(fp) # hfs
89: pushl r1
90: pushl r0 # pushd
91: ldd -60(fp)
92: pushd
93: callf $24,_Ksubd # e - 4 * f
94: ldd r0 # load acc with result of Ksubd
95: cvdl 12(fp) # quad = (e - 4 * f); }
96: jbr L39
97: L35: ldd -84(fp) # else {
98: cvdl -96(fp) # k = x;
99: pushl 16(fp) # hfs
100: pushl $0
101: pushl -96(fp)
102: pushl $0
103: pushl $0 # acc
104: callf $24,_Kcvtld
105: # std r0
106: pushl 16(fp) # hfs
107: pushl r1
108: pushl r0 # pushd (double)k
109: ldd -84(fp)
110: pushd # x
111: callf $24,_Ksubd
112: movl r1,-88(fp)
113: movl r0,-92(fp) # std -92(fp) y = x - k;
114: addl3 -96(fp),12(fp),r0
115: andl3 $3,r0,12(fp) # quad = (quad + k) & 03; }
116: L39: bitl $1,12(fp) # if (quad & 01)
117: jeql L40
118: pushl 16(fp) # hfs
119: ldd -92(fp)
120: pushd # y
121: ldd L41
122: pushd # 1.0
123: callf $24,_Ksubd
124: movl r1,-88(fp)
125: movl r0,-92(fp) # std -92(fp) y = 1 - y;
126: L40: cmpl 12(fp),$1 # if (quad > 1)
127: jleq L42
128: lnd -92(fp)
129: std -92(fp) # y = -y;
130: L42: pushl 16(fp) # hfs
131: ldd -92(fp)
132: pushd
133: pushd
134: callf $24,_Kmuld # y*y
135: movl r1,-72(fp)
136: movl r0,-76(fp) # std -76(fp) ysq=y*y;
137:
138: pushl 16(fp) # hfs
139: ldd -76(fp)
140: pushd
141: ldd _p4
142: pushd
143: callf $24,_Kmuld # p4*ysq
144:
145: pushl 16(fp) # hfs
146: pushl _p3+4
147: pushl _p3
148: pushl r1
149: pushl r0 # pushd
150: callf $24,_Kaddd # p4*ysq+p3
151:
152: pushl 16(fp) # hfs
153: pushl -72(fp)
154: pushl -76(fp)
155: pushl r1
156: pushl r0 # pushd
157: callf $24,_Kmuld # (p4*ysq+p3)*ysq
158:
159: pushl 16(fp) # hfs
160: pushl _p2+4
161: pushl _p2
162: pushl r1
163: pushl r0 # pushd
164: callf $24,_Kaddd # (p4*ysq+p3)*ysq+p2
165:
166: pushl 16(fp) # hfs
167: pushl -72(fp)
168: pushl -76(fp)
169: pushl r1
170: pushl r0 # pushd
171: callf $24,_Kmuld # ((p4*ysq+p3)*ysq+p2)*ysq
172:
173: pushl 16(fp) # hfs
174: pushl _p1+4
175: pushl _p1
176: pushl r1
177: pushl r0 # pushd
178: callf $24,_Kaddd # ((p4*ysq+p3)*ysq+p2)*ysq+p1
179:
180: pushl 16(fp) # hfs
181: pushl -72(fp)
182: pushl -76(fp)
183: pushl r1
184: pushl r0 # pushd
185: callf $24,_Kmuld # (((p4*ysq+p3)*ysq+p2)*ysq+p1)*ysq
186:
187: pushl 16(fp) # hfs
188: pushl _p0+4
189: pushl _p0
190: pushl r1
191: pushl r0 # pushd
192: callf $24,_Kaddd # (((p4*ysq+p3)*ysq+p2)*ysq+p1)*ysq+p0
193:
194: pushl 16(fp) # hfs
195: pushl -88(fp)
196: pushl -92(fp)
197: pushl r1
198: pushl r0 # pushd
199: callf $24,_Kmuld # ((((p4*ysq+p3)*ysq+p2)*ysq+p1)*ysq+p0)*y;
200:
201: movl r1,-100(fp)
202: movl r0,-104(fp) # std -104(fp) temp1
203: ldd r0
204: cvdf
205: stf r0
206: ret#1
207:
208: .data
209: .align 2
210: _twoopi:.long 0x4022F983, 0x6E4E4415 # .double .63661977236758134
211: _p0: .long 0x40C90FDA, 0x90304197 # .double 1.57079631844
212: _p1: .long 0xC0255DE0, 0xB36CEE75 # .double -.645963710599
213: _p2: .long 0x3EA33457, 0xA736E807 # .double .079689678945999999
214: _p3: .long 0xBC992665, 0x5E9A6554 # .double -.0046737666099999999
215: _p4: .long 0x3A1ED7FA, 0xCC54924E # .double .000151485129
216: L36: .long 0x47FFF800, 0x00000000 # .double 32764
217: L37: .long 0x3F800000, 0x00000000 # .double .25
218: L38: .long 0x41800000, 0x00000000 # .double 4
219: 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.