Return to submul_1.S CVS log

Up to [local] / OpenXM_contrib / gmp / mpn / x86

Annotation of OpenXM_contrib/gmp/mpn/x86/submul_1.S, Revision 1.1.1.1

1.1       maekawa     1: /* i80386 __mpn_submul_1 -- Multiply a limb vector with a limb and subtract
                      2:    the result from a second limb vector.
                      3:
                      4: Copyright (C) 1992, 1994 Free Software Foundation, Inc.
                      5:
                      6: This file is part of the GNU MP Library.
                      7:
                      8: The GNU MP Library is free software; you can redistribute it and/or modify
                      9: it under the terms of the GNU Library General Public License as published by
                     10: the Free Software Foundation; either version 2 of the License, or (at your
                     11: option) any later version.
                     12:
                     13: The GNU MP Library is distributed in the hope that it will be useful, but
                     14: WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
                     15: or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Library General Public
                     16: License for more details.
                     17:
                     18: You should have received a copy of the GNU Library General Public License
                     19: along with the GNU MP Library; see the file COPYING.LIB.  If not, write to
                     20: the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
                     21: MA 02111-1307, USA. */
                     22:
                     23: /*
                     24:    INPUT PARAMETERS
                     25:    res_ptr     (sp + 4)
                     26:    s1_ptr      (sp + 8)
                     27:    size                (sp + 12)
                     28:    s2_limb     (sp + 16)
                     29: */
                     30:
                     31: #include "sysdep.h"
                     32: #include "asm-syntax.h"
                     33:
                     34: #define res_ptr edi
                     35: #define s1_ptr esi
                     36: #define size ecx
                     37: #define s2_limb ebp
                     38:
                     39:        TEXT
                     40:        ALIGN (3)
                     41:        GLOBL   C_SYMBOL_NAME(__mpn_submul_1)
                     42: C_SYMBOL_NAME(__mpn_submul_1:)
                     43:
                     44:        INSN1(push,l    ,R(edi))
                     45:        INSN1(push,l    ,R(esi))
                     46:        INSN1(push,l    ,R(ebx))
                     47:        INSN1(push,l    ,R(ebp))
                     48:
                     49:        INSN2(mov,l     ,R(res_ptr),MEM_DISP(esp,20))
                     50:        INSN2(mov,l     ,R(s1_ptr),MEM_DISP(esp,24))
                     51:        INSN2(mov,l     ,R(size),MEM_DISP(esp,28))
                     52:        INSN2(mov,l     ,R(s2_limb),MEM_DISP(esp,32))
                     53:
                     54:        INSN2(lea,l     ,R(res_ptr),MEM_INDEX(res_ptr,size,4))
                     55:        INSN2(lea,l     ,R(s1_ptr),MEM_INDEX(s1_ptr,size,4))
                     56:        INSN1(neg,l     ,R(size))
                     57:        INSN2(xor,l     ,R(ebx),R(ebx))
                     58:        ALIGN (3)
                     59: Loop:
                     60:        INSN2(mov,l     ,R(eax),MEM_INDEX(s1_ptr,size,4))
                     61:        INSN1(mul,l     ,R(s2_limb))
                     62:        INSN2(add,l     ,R(eax),R(ebx))
                     63:        INSN2(adc,l     ,R(edx),$0)
                     64:        INSN2(sub,l     ,MEM_INDEX(res_ptr,size,4),R(eax))
                     65:        INSN2(adc,l     ,R(edx),$0)
                     66:        INSN2(mov,l     ,R(ebx),R(edx))
                     67:
                     68:        INSN1(inc,l     ,R(size))
                     69:        INSN1(jnz,      ,Loop)
                     70:        INSN2(mov,l     ,R(eax),R(ebx))
                     71:
                     72:        INSN1(pop,l     ,R(ebp))
                     73:        INSN1(pop,l     ,R(ebx))
                     74:        INSN1(pop,l     ,R(esi))
                     75:        INSN1(pop,l     ,R(edi))
                     76:        ret