Return to addmul_1.S CVS log

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

Annotation of OpenXM_contrib/gmp/mpn/x86/addmul_1.S, Revision 1.1

1.1     ! maekawa     1: /* i80386 __mpn_addmul_1 -- Multiply a limb vector with a limb and add
        !             2:    the result to 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_addmul_1)
        !            42: C_SYMBOL_NAME(__mpn_addmul_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(add,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