![[BACK]](/cvsweb/icons/back.gif){kind=icon}
Return to math.s CVS log ![[TXT]](/cvsweb/icons/text.gif){kind=icon}
![[DIR]](/cvsweb/icons/dir.gif){kind=icon}
Up to [Quake II] / quake2 / linux|
Return to math.s CVS log | Up to [Quake II] / quake2 / linux |
1.1 root 1: //
2: // math.s
3: // x86 assembly-language math routines.
4:
5: #define GLQUAKE 1 // don't include unneeded defs
6: #include "qasm.h"
7:
8:
9: #if id386
10:
11: .data
12:
13: #if 0
14: .align 4
15: Ljmptab: .long Lcase0, Lcase1, Lcase2, Lcase3
16: .long Lcase4, Lcase5, Lcase6, Lcase7
17: #endif
18:
19: .text
20:
21: // TODO: rounding needed?
22: // stack parameter offset
23: #define val 4
24:
25: .globl C(Invert24To16)
26: C(Invert24To16):
27:
28: movl val(%esp),%ecx
29: movl $0x100,%edx // 0x10000000000 as dividend
30: cmpl %edx,%ecx
31: jle LOutOfRange
32:
33: subl %eax,%eax
34: divl %ecx
35:
36: ret
37:
38: LOutOfRange:
39: movl $0xFFFFFFFF,%eax
40: ret
41:
42: #define in 4
43: #define out 8
44:
45: .align 2
46: .globl C(TransformVector)
47: C(TransformVector):
48: movl in(%esp),%eax
49: movl out(%esp),%edx
50:
51: flds (%eax) // in[0]
52: fmuls C(vright) // in[0]*vright[0]
53: flds (%eax) // in[0] | in[0]*vright[0]
54: fmuls C(vup) // in[0]*vup[0] | in[0]*vright[0]
55: flds (%eax) // in[0] | in[0]*vup[0] | in[0]*vright[0]
56: fmuls C(vpn) // in[0]*vpn[0] | in[0]*vup[0] | in[0]*vright[0]
57:
58: flds 4(%eax) // in[1] | ...
59: fmuls C(vright)+4 // in[1]*vright[1] | ...
60: flds 4(%eax) // in[1] | in[1]*vright[1] | ...
61: fmuls C(vup)+4 // in[1]*vup[1] | in[1]*vright[1] | ...
62: flds 4(%eax) // in[1] | in[1]*vup[1] | in[1]*vright[1] | ...
63: fmuls C(vpn)+4 // in[1]*vpn[1] | in[1]*vup[1] | in[1]*vright[1] | ...
64: fxch %st(2) // in[1]*vright[1] | in[1]*vup[1] | in[1]*vpn[1] | ...
65:
66: faddp %st(0),%st(5) // in[1]*vup[1] | in[1]*vpn[1] | ...
67: faddp %st(0),%st(3) // in[1]*vpn[1] | ...
68: faddp %st(0),%st(1) // vpn_accum | vup_accum | vright_accum
69:
70: flds 8(%eax) // in[2] | ...
71: fmuls C(vright)+8 // in[2]*vright[2] | ...
72: flds 8(%eax) // in[2] | in[2]*vright[2] | ...
73: fmuls C(vup)+8 // in[2]*vup[2] | in[2]*vright[2] | ...
74: flds 8(%eax) // in[2] | in[2]*vup[2] | in[2]*vright[2] | ...
75: fmuls C(vpn)+8 // in[2]*vpn[2] | in[2]*vup[2] | in[2]*vright[2] | ...
76: fxch %st(2) // in[2]*vright[2] | in[2]*vup[2] | in[2]*vpn[2] | ...
77:
78: faddp %st(0),%st(5) // in[2]*vup[2] | in[2]*vpn[2] | ...
79: faddp %st(0),%st(3) // in[2]*vpn[2] | ...
80: faddp %st(0),%st(1) // vpn_accum | vup_accum | vright_accum
81:
82: fstps 8(%edx) // out[2]
83: fstps 4(%edx) // out[1]
84: fstps (%edx) // out[0]
85:
86: ret
87:
88: #if 0 // in C
89:
90: #define EMINS 4+4
91: #define EMAXS 4+8
92: #define P 4+12
93:
94: .align 2
95: .globl C(BoxOnPlaneSide)
96: C(BoxOnPlaneSide):
97: pushl %ebx
98:
99: movl P(%esp),%edx
100: movl EMINS(%esp),%ecx
101: xorl %eax,%eax
102: movl EMAXS(%esp),%ebx
103: movb pl_signbits(%edx),%al
104: cmpb $8,%al
105: jge Lerror
106: flds pl_normal(%edx) // p->normal[0]
107: fld %st(0) // p->normal[0] | p->normal[0]
108: jmp Ljmptab(,%eax,4)
109:
110:
111: //dist1= p->normal[0]*emaxs[0] + p->normal[1]*emaxs[1] + p->normal[2]*emaxs[2];
112: //dist2= p->normal[0]*emins[0] + p->normal[1]*emins[1] + p->normal[2]*emins[2];
113: Lcase0:
114: fmuls (%ebx) // p->normal[0]*emaxs[0] | p->normal[0]
115: flds pl_normal+4(%edx) // p->normal[1] | p->normal[0]*emaxs[0] |
116: // p->normal[0]
117: fxch %st(2) // p->normal[0] | p->normal[0]*emaxs[0] |
118: // p->normal[1]
119: fmuls (%ecx) // p->normal[0]*emins[0] |
120: // p->normal[0]*emaxs[0] | p->normal[1]
121: fxch %st(2) // p->normal[1] | p->normal[0]*emaxs[0] |
122: // p->normal[0]*emins[0]
123: fld %st(0) // p->normal[1] | p->normal[1] |
124: // p->normal[0]*emaxs[0] |
125: // p->normal[0]*emins[0]
126: fmuls 4(%ebx) // p->normal[1]*emaxs[1] | p->normal[1] |
127: // p->normal[0]*emaxs[0] |
128: // p->normal[0]*emins[0]
129: flds pl_normal+8(%edx) // p->normal[2] | p->normal[1]*emaxs[1] |
130: // p->normal[1] | p->normal[0]*emaxs[0] |
131: // p->normal[0]*emins[0]
132: fxch %st(2) // p->normal[1] | p->normal[1]*emaxs[1] |
133: // p->normal[2] | p->normal[0]*emaxs[0] |
134: // p->normal[0]*emins[0]
135: fmuls 4(%ecx) // p->normal[1]*emins[1] |
136: // p->normal[1]*emaxs[1] |
137: // p->normal[2] | p->normal[0]*emaxs[0] |
138: // p->normal[0]*emins[0]
139: fxch %st(2) // p->normal[2] | p->normal[1]*emaxs[1] |
140: // p->normal[1]*emins[1] |
141: // p->normal[0]*emaxs[0] |
142: // p->normal[0]*emins[0]
143: fld %st(0) // p->normal[2] | p->normal[2] |
144: // p->normal[1]*emaxs[1] |
145: // p->normal[1]*emins[1] |
146: // p->normal[0]*emaxs[0] |
147: // p->normal[0]*emins[0]
148: fmuls 8(%ebx) // p->normal[2]*emaxs[2] |
149: // p->normal[2] |
150: // p->normal[1]*emaxs[1] |
151: // p->normal[1]*emins[1] |
152: // p->normal[0]*emaxs[0] |
153: // p->normal[0]*emins[0]
154: fxch %st(5) // p->normal[0]*emins[0] |
155: // p->normal[2] |
156: // p->normal[1]*emaxs[1] |
157: // p->normal[1]*emins[1] |
158: // p->normal[0]*emaxs[0] |
159: // p->normal[2]*emaxs[2]
160: faddp %st(0),%st(3) //p->normal[2] |
161: // p->normal[1]*emaxs[1] |
162: // p->normal[1]*emins[1]+p->normal[0]*emins[0]|
163: // p->normal[0]*emaxs[0] |
164: // p->normal[2]*emaxs[2]
165: fmuls 8(%ecx) //p->normal[2]*emins[2] |
166: // p->normal[1]*emaxs[1] |
167: // p->normal[1]*emins[1]+p->normal[0]*emins[0]|
168: // p->normal[0]*emaxs[0] |
169: // p->normal[2]*emaxs[2]
170: fxch %st(1) //p->normal[1]*emaxs[1] |
171: // p->normal[2]*emins[2] |
172: // p->normal[1]*emins[1]+p->normal[0]*emins[0]|
173: // p->normal[0]*emaxs[0] |
174: // p->normal[2]*emaxs[2]
175: faddp %st(0),%st(3) //p->normal[2]*emins[2] |
176: // p->normal[1]*emins[1]+p->normal[0]*emins[0]|
177: // p->normal[0]*emaxs[0]+p->normal[1]*emaxs[1]|
178: // p->normal[2]*emaxs[2]
179: fxch %st(3) //p->normal[2]*emaxs[2] +
180: // p->normal[1]*emins[1]+p->normal[0]*emins[0]|
181: // p->normal[0]*emaxs[0]+p->normal[1]*emaxs[1]|
182: // p->normal[2]*emins[2]
183: faddp %st(0),%st(2) //p->normal[1]*emins[1]+p->normal[0]*emins[0]|
184: // dist1 | p->normal[2]*emins[2]
185:
186: jmp LSetSides
187:
188: //dist1= p->normal[0]*emins[0] + p->normal[1]*emaxs[1] + p->normal[2]*emaxs[2];
189: //dist2= p->normal[0]*emaxs[0] + p->normal[1]*emins[1] + p->normal[2]*emins[2];
190: Lcase1:
191: fmuls (%ecx) // emins[0]
192: flds pl_normal+4(%edx)
193: fxch %st(2)
194: fmuls (%ebx) // emaxs[0]
195: fxch %st(2)
196: fld %st(0)
197: fmuls 4(%ebx) // emaxs[1]
198: flds pl_normal+8(%edx)
199: fxch %st(2)
200: fmuls 4(%ecx) // emins[1]
201: fxch %st(2)
202: fld %st(0)
203: fmuls 8(%ebx) // emaxs[2]
204: fxch %st(5)
205: faddp %st(0),%st(3)
206: fmuls 8(%ecx) // emins[2]
207: fxch %st(1)
208: faddp %st(0),%st(3)
209: fxch %st(3)
210: faddp %st(0),%st(2)
211:
212: jmp LSetSides
213:
214: //dist1= p->normal[0]*emaxs[0] + p->normal[1]*emins[1] + p->normal[2]*emaxs[2];
215: //dist2= p->normal[0]*emins[0] + p->normal[1]*emaxs[1] + p->normal[2]*emins[2];
216: Lcase2:
217: fmuls (%ebx) // emaxs[0]
218: flds pl_normal+4(%edx)
219: fxch %st(2)
220: fmuls (%ecx) // emins[0]
221: fxch %st(2)
222: fld %st(0)
223: fmuls 4(%ecx) // emins[1]
224: flds pl_normal+8(%edx)
225: fxch %st(2)
226: fmuls 4(%ebx) // emaxs[1]
227: fxch %st(2)
228: fld %st(0)
229: fmuls 8(%ebx) // emaxs[2]
230: fxch %st(5)
231: faddp %st(0),%st(3)
232: fmuls 8(%ecx) // emins[2]
233: fxch %st(1)
234: faddp %st(0),%st(3)
235: fxch %st(3)
236: faddp %st(0),%st(2)
237:
238: jmp LSetSides
239:
240: //dist1= p->normal[0]*emins[0] + p->normal[1]*emins[1] + p->normal[2]*emaxs[2];
241: //dist2= p->normal[0]*emaxs[0] + p->normal[1]*emaxs[1] + p->normal[2]*emins[2];
242: Lcase3:
243: fmuls (%ecx) // emins[0]
244: flds pl_normal+4(%edx)
245: fxch %st(2)
246: fmuls (%ebx) // emaxs[0]
247: fxch %st(2)
248: fld %st(0)
249: fmuls 4(%ecx) // emins[1]
250: flds pl_normal+8(%edx)
251: fxch %st(2)
252: fmuls 4(%ebx) // emaxs[1]
253: fxch %st(2)
254: fld %st(0)
255: fmuls 8(%ebx) // emaxs[2]
256: fxch %st(5)
257: faddp %st(0),%st(3)
258: fmuls 8(%ecx) // emins[2]
259: fxch %st(1)
260: faddp %st(0),%st(3)
261: fxch %st(3)
262: faddp %st(0),%st(2)
263:
264: jmp LSetSides
265:
266: //dist1= p->normal[0]*emaxs[0] + p->normal[1]*emaxs[1] + p->normal[2]*emins[2];
267: //dist2= p->normal[0]*emins[0] + p->normal[1]*emins[1] + p->normal[2]*emaxs[2];
268: Lcase4:
269: fmuls (%ebx) // emaxs[0]
270: flds pl_normal+4(%edx)
271: fxch %st(2)
272: fmuls (%ecx) // emins[0]
273: fxch %st(2)
274: fld %st(0)
275: fmuls 4(%ebx) // emaxs[1]
276: flds pl_normal+8(%edx)
277: fxch %st(2)
278: fmuls 4(%ecx) // emins[1]
279: fxch %st(2)
280: fld %st(0)
281: fmuls 8(%ecx) // emins[2]
282: fxch %st(5)
283: faddp %st(0),%st(3)
284: fmuls 8(%ebx) // emaxs[2]
285: fxch %st(1)
286: faddp %st(0),%st(3)
287: fxch %st(3)
288: faddp %st(0),%st(2)
289:
290: jmp LSetSides
291:
292: //dist1= p->normal[0]*emins[0] + p->normal[1]*emaxs[1] + p->normal[2]*emins[2];
293: //dist2= p->normal[0]*emaxs[0] + p->normal[1]*emins[1] + p->normal[2]*emaxs[2];
294: Lcase5:
295: fmuls (%ecx) // emins[0]
296: flds pl_normal+4(%edx)
297: fxch %st(2)
298: fmuls (%ebx) // emaxs[0]
299: fxch %st(2)
300: fld %st(0)
301: fmuls 4(%ebx) // emaxs[1]
302: flds pl_normal+8(%edx)
303: fxch %st(2)
304: fmuls 4(%ecx) // emins[1]
305: fxch %st(2)
306: fld %st(0)
307: fmuls 8(%ecx) // emins[2]
308: fxch %st(5)
309: faddp %st(0),%st(3)
310: fmuls 8(%ebx) // emaxs[2]
311: fxch %st(1)
312: faddp %st(0),%st(3)
313: fxch %st(3)
314: faddp %st(0),%st(2)
315:
316: jmp LSetSides
317:
318: //dist1= p->normal[0]*emaxs[0] + p->normal[1]*emins[1] + p->normal[2]*emins[2];
319: //dist2= p->normal[0]*emins[0] + p->normal[1]*emaxs[1] + p->normal[2]*emaxs[2];
320: Lcase6:
321: fmuls (%ebx) // emaxs[0]
322: flds pl_normal+4(%edx)
323: fxch %st(2)
324: fmuls (%ecx) // emins[0]
325: fxch %st(2)
326: fld %st(0)
327: fmuls 4(%ecx) // emins[1]
328: flds pl_normal+8(%edx)
329: fxch %st(2)
330: fmuls 4(%ebx) // emaxs[1]
331: fxch %st(2)
332: fld %st(0)
333: fmuls 8(%ecx) // emins[2]
334: fxch %st(5)
335: faddp %st(0),%st(3)
336: fmuls 8(%ebx) // emaxs[2]
337: fxch %st(1)
338: faddp %st(0),%st(3)
339: fxch %st(3)
340: faddp %st(0),%st(2)
341:
342: jmp LSetSides
343:
344: //dist1= p->normal[0]*emins[0] + p->normal[1]*emins[1] + p->normal[2]*emins[2];
345: //dist2= p->normal[0]*emaxs[0] + p->normal[1]*emaxs[1] + p->normal[2]*emaxs[2];
346: Lcase7:
347: fmuls (%ecx) // emins[0]
348: flds pl_normal+4(%edx)
349: fxch %st(2)
350: fmuls (%ebx) // emaxs[0]
351: fxch %st(2)
352: fld %st(0)
353: fmuls 4(%ecx) // emins[1]
354: flds pl_normal+8(%edx)
355: fxch %st(2)
356: fmuls 4(%ebx) // emaxs[1]
357: fxch %st(2)
358: fld %st(0)
359: fmuls 8(%ecx) // emins[2]
360: fxch %st(5)
361: faddp %st(0),%st(3)
362: fmuls 8(%ebx) // emaxs[2]
363: fxch %st(1)
364: faddp %st(0),%st(3)
365: fxch %st(3)
366: faddp %st(0),%st(2)
367:
368: LSetSides:
369:
370: // sides = 0;
371: // if (dist1 >= p->dist)
372: // sides = 1;
373: // if (dist2 < p->dist)
374: // sides |= 2;
375:
376: faddp %st(0),%st(2) // dist1 | dist2
377: fcomps pl_dist(%edx)
378: xorl %ecx,%ecx
379: fnstsw %ax
380: fcomps pl_dist(%edx)
381: andb $1,%ah
382: xorb $1,%ah
383: addb %ah,%cl
384:
385: fnstsw %ax
386: andb $1,%ah
387: addb %ah,%ah
388: addb %ah,%cl
389:
390: // return sides;
391:
392: popl %ebx
393: movl %ecx,%eax // return status
394:
395: ret
396:
397:
398: Lerror:
399: call C(BOPS_Error)
400:
401: #endif
402:
403: #endif // id386
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.