Annotation of OpenXM_contrib/pari/src/kernel/ix86/l0asm.c, Revision 1.1
1.1 ! maekawa 1: /***********************************************************************/
! 2: /** **/
! 3: /** Low level arithmetic for PARI **/
! 4: /** written by Bruno Haible 14.11.1992 **/
! 5: /** macroified 21.1.1998 **/
! 6: /** **/
! 7: /***********************************************************************/
! 8: /* $Id: l0asm.c,v 1.1.1.1 1999/09/16 13:47:45 karim Exp $
! 9:
! 10: /* processor: Intel ix86 in native mode
! 11: * assembler syntax: GNU or SUN, moves go from left to right
! 12: * compiler: GNU gcc or SUN cc
! 13: * parameter passing convention: on the stack 4(%esp),8(%esp),...
! 14: * registers: %eax,%edx,%ecx may be modified,
! 15: * everything else must be saved and restored
! 16: * result: passed in %eax
! 17: * word length: 32 bits
! 18: */
! 19:
! 20: /* If the C compiler is GNU C, better stuff is contained in asmix86inline.h,
! 21: but we compile and link this file nevertheless because it defines the
! 22: global variable `hiremainder'. */
! 23:
! 24: /* This should ideally be determined at configure time. */
! 25: #if defined(__EMX__) || defined(__DJGCC__) || defined(__GO32__) || (defined(linux) && !defined(__ELF__)) || defined(__386BSD__) || defined(__NetBSD__) || (defined(__FreeBSD__) && !defined(__ELF__)) || defined(NeXT) || defined(__CYGWIN32__) || defined(__MINGW32__)
! 26: # define ASM_UNDERSCORE
! 27: #endif
! 28:
! 29: #include "l0asm.h"
! 30:
! 31: #if defined(__EMX__) && defined(__NO_AOUT) /* No underscores, rest as is */
! 32: # undef C
! 33: # define C(entrypoint) /**/entrypoint
! 34: #endif
! 35:
! 36:
! 37: #undef ALIGN
! 38: #ifdef _MSC_VER
! 39: #define ALIGN
! 40: #else
! 41: #if defined(ASM_UNDERSCORE) && !(defined(__CYGWIN32__) || defined(__MINGW32__))
! 42: /* BSD syntax assembler */
! 43: #define ALIGN .align 3
! 44: #else
! 45: /* ELF syntax assembler */
! 46: #define ALIGN .align 8
! 47: #endif
! 48: #endif
! 49:
! 50: GLOBL(C(addll))
! 51: GLOBL(C(subll))
! 52: GLOBL(C(addllx))
! 53: GLOBL(C(subllx))
! 54: GLOBL(C(shiftl))
! 55: GLOBL(C(shiftlr))
! 56: GLOBL(C(bfffo))
! 57: GLOBL(C(mulll))
! 58: GLOBL(C(addmul))
! 59: GLOBL(C(divll))
! 60: GLOBL(C(overflow))
! 61: GLOBL(C(hiremainder))
! 62:
! 63: #ifdef _MSC_VER
! 64: unsigned long overflow;
! 65: unsigned long hiremainder;
! 66: #else
! 67: #ifdef BSD_SYNTAX
! 68: # if defined(__EMX__) && defined(__NO_AOUT) /* Otherwise IBM linker will
! 69: ignore this module */
! 70: .data
! 71: .align 2
! 72: C(overflow):
! 73: .long 0
! 74: C(hiremainder):
! 75: .long 0
! 76: .text
! 77: # else
! 78: .comm C(overflow),4
! 79: .comm C(hiremainder),4
! 80: # endif
! 81: #endif
! 82: #ifdef ELF_SYNTAX
! 83: .comm C(overflow),4,4
! 84: .comm C(hiremainder),4,4
! 85: #endif
! 86: #endif
! 87:
! 88: TEXT()
! 89:
! 90: ALIGN
! 91: FUNBEGIN(addll)
! 92: INSN2(xor,l ,R(edx),R(edx)) /* clear %edx */
! 93: INSN2(mov,l ,X4 MEM_DISP(esp,4),R(eax)) /* get x */
! 94: INSN2(add,l ,X4 MEM_DISP(esp,8),R(eax)) /* add y */
! 95: INSN2(adc,l ,R(edx),R(edx)) /* %edx := carry */
! 96: INSN2(mov,l ,R(edx),C(overflow)) /* set overflow */
! 97: ret /* return %eax */
! 98: FUNEND()
! 99:
! 100: ALIGN
! 101: FUNBEGIN(addllx)
! 102: INSN2(xor,l ,R(edx),R(edx)) /* clear %edx */
! 103: INSN2(xor,l ,R(eax),R(eax)) /* clear %eax */
! 104: INSN2(sub,l ,C(overflow),R(eax)) /* set carry */
! 105: INSN2(mov,l ,X4 MEM_DISP(esp,4),R(eax)) /* get x */
! 106: INSN2(adc,l ,X4 MEM_DISP(esp,8),R(eax)) /* add y and carry */
! 107: INSN2(adc,l ,R(edx),R(edx)) /* %edx := carry */
! 108: INSN2(mov,l ,R(edx),C(overflow)) /* set overflow */
! 109: ret /* return %eax */
! 110: FUNEND()
! 111:
! 112: ALIGN
! 113: FUNBEGIN(subll)
! 114: INSN2(xor,l ,R(edx),R(edx)) /* clear %edx */
! 115: INSN2(mov,l ,X4 MEM_DISP(esp,4),R(eax)) /* get x */
! 116: INSN2(sub,l ,X4 MEM_DISP(esp,8),R(eax)) /* subtract y */
! 117: INSN2(adc,l ,R(edx),R(edx)) /* %edx := carry */
! 118: INSN2(mov,l ,R(edx),C(overflow)) /* set overflow */
! 119: ret /* return %eax */
! 120: FUNEND()
! 121:
! 122: ALIGN
! 123: FUNBEGIN(subllx)
! 124: INSN2(xor,l ,R(edx),R(edx)) /* clear %edx */
! 125: INSN2(xor,l ,R(eax),R(eax)) /* clear %eax */
! 126: INSN2(sub,l ,C(overflow),R(eax)) /* set carry */
! 127: INSN2(mov,l ,X4 MEM_DISP(esp,4),R(eax)) /* get x */
! 128: INSN2(sbb,l ,X4 MEM_DISP(esp,8),R(eax)) /* subtract y and carry */
! 129: INSN2(adc,l ,R(edx),R(edx)) /* %edx := carry */
! 130: INSN2(mov,l ,R(edx),C(overflow)) /* set overflow */
! 131: ret /* return %eax */
! 132: FUNEND()
! 133:
! 134: ALIGN
! 135: FUNBEGIN(shiftl)
! 136: INSN2(mov,l ,X4 MEM_DISP(esp,4),R(eax)) /* get x */
! 137: INSN2(mov,b ,X1 MEM_DISP(esp,8),R(cl)) /* get shift count i */
! 138: INSN2(xor,l ,R(edx),R(edx)) /* clear %edx */
! 139: INSN2SHCL(shld,l,R(eax),R(edx)) /* shift %edx left by i bits,
! 140: feeding in %eax from the right */
! 141: INSN2(shl,l ,R(cl),R(eax)) /* shift %eax left by i bits */
! 142: INSN2(mov,l ,R(edx),C(hiremainder)) /* set hiremainder */
! 143: ret /* return %eax */
! 144: FUNEND()
! 145:
! 146: ALIGN
! 147: FUNBEGIN(shiftlr)
! 148: INSN2(mov,l ,X4 MEM_DISP(esp,4),R(eax)) /* get x */
! 149: INSN2(mov,b ,X1 MEM_DISP(esp,8),R(cl)) /* get shift count i */
! 150: INSN2(xor,l ,R(edx),R(edx)) /* clear %edx */
! 151: INSN2SHCL(shrd,l,R(eax),R(edx)) /* shift %edx right by i bits,
! 152: feeding in %eax from the left */
! 153: INSN2(shr,l ,R(cl),R(eax)) /* shift %eax right by i bits */
! 154: INSN2(mov,l ,R(edx),C(hiremainder)) /* set hiremainder */
! 155: ret /* return %eax */
! 156: FUNEND()
! 157:
! 158: ALIGN
! 159: FUNBEGIN(bfffo)
! 160: INSN2(mov,l ,X4 MEM_DISP(esp,4),R(eax)) /* get x */
! 161: INSN2(bsr,l ,R(eax),R(edx)) /* %edx := number of leading bit */
! 162: INSN2(mov,l ,NUM(31),R(eax))
! 163: INSN2(sub,l ,R(edx),R(eax)) /* result is 31 - %edx */
! 164: ret /* return %eax */
! 165: FUNEND()
! 166:
! 167: ALIGN
! 168: FUNBEGIN(mulll)
! 169: INSN2(mov,l ,X4 MEM_DISP(esp,4),R(eax)) /* get x */
! 170: INSN1(mul,l ,X4 MEM_DISP(esp,8)) /* %edx|%eax := x * y */
! 171: INSN2(mov,l ,R(edx),C(hiremainder)) /* store high word */
! 172: ret /* return low word %eax */
! 173: FUNEND()
! 174:
! 175: ALIGN
! 176: FUNBEGIN(addmul)
! 177: INSN2(xor,l ,R(ecx),R(ecx)) /* clear %ecx */
! 178: INSN2(mov,l ,X4 MEM_DISP(esp,4),R(eax)) /* get x */
! 179: INSN1(mul,l ,X4 MEM_DISP(esp,8)) /* %edx|%eax := x * y */
! 180: INSN2(add,l ,C(hiremainder),R(eax)) /* add 0|hiremainder */
! 181: INSN2(adc,l ,R(ecx),R(edx))
! 182: INSN2(mov,l ,R(edx),C(hiremainder)) /* store high word */
! 183: ret /* return low word %eax */
! 184: FUNEND()
! 185:
! 186: ALIGN
! 187: FUNBEGIN(divll)
! 188: INSN2(mov,l ,X4 MEM_DISP(esp,4),R(eax)) /* get low word x */
! 189: INSN2(mov,l ,C(hiremainder),R(edx)) /* get high word */
! 190: INSN1(div,l ,X4 MEM_DISP(esp,8)) /* divide %edx|%eax by y */
! 191: INSN2(mov,l ,R(edx), C(hiremainder)) /* store remainder */
! 192: ret /* return quotient %eax */
! 193: FUNEND()
! 194:
! 195: ALIGN
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to l0asm.c CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib / pari / src / kernel / ix86