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