Return to rshift.asm CVS log

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

Annotation of OpenXM_contrib/gmp/mpn/x86/rshift.asm, Revision 1.1

1.1     ! maekawa     1: dnl  x86 mpn_rshift -- mpn right shift.
        !             2:
        !             3: dnl  Copyright (C) 1992, 1994, 1996, 1999, 2000 Free Software Foundation,
        !             4: dnl  Inc.
        !             5: dnl
        !             6: dnl  This file is part of the GNU MP Library.
        !             7: dnl
        !             8: dnl  The GNU MP Library is free software; you can redistribute it and/or
        !             9: dnl  modify it under the terms of the GNU Lesser General Public License as
        !            10: dnl  published by the Free Software Foundation; either version 2.1 of the
        !            11: dnl  License, or (at your option) any later version.
        !            12: dnl
        !            13: dnl  The GNU MP Library is distributed in the hope that it will be useful,
        !            14: dnl  but WITHOUT ANY WARRANTY; without even the implied warranty of
        !            15: dnl  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
        !            16: dnl  Lesser General Public License for more details.
        !            17: dnl
        !            18: dnl  You should have received a copy of the GNU Lesser General Public
        !            19: dnl  License along with the GNU MP Library; see the file COPYING.LIB.  If
        !            20: dnl  not, write to the Free Software Foundation, Inc., 59 Temple Place -
        !            21: dnl  Suite 330, Boston, MA 02111-1307, USA.
        !            22:
        !            23:
        !            24: include(`../config.m4')
        !            25:
        !            26:
        !            27: C mp_limb_t mpn_rshift (mp_ptr dst, mp_srcptr src, mp_size_t size,
        !            28: C                       unsigned shift);
        !            29:
        !            30: defframe(PARAM_SHIFT,16)
        !            31: defframe(PARAM_SIZE, 12)
        !            32: defframe(PARAM_SRC,  8)
        !            33: defframe(PARAM_DST,  4)
        !            34:
        !            35:        .text
        !            36:        ALIGN(8)
        !            37: PROLOGUE(mpn_rshift)
        !            38:
        !            39:        pushl   %edi
        !            40:        pushl   %esi
        !            41:        pushl   %ebx
        !            42: deflit(`FRAME',12)
        !            43:
        !            44:        movl    PARAM_DST,%edi
        !            45:        movl    PARAM_SRC,%esi
        !            46:        movl    PARAM_SIZE,%edx
        !            47:        movl    PARAM_SHIFT,%ecx
        !            48:
        !            49:        leal    -4(%edi,%edx,4),%edi
        !            50:        leal    (%esi,%edx,4),%esi
        !            51:        negl    %edx
        !            52:
        !            53:        movl    (%esi,%edx,4),%ebx      C read least significant limb
        !            54:        xorl    %eax,%eax
        !            55:        shrdl(  %cl, %ebx, %eax)        C compute carry limb
        !            56:        incl    %edx
        !            57:        jz      L(end)
        !            58:        pushl   %eax                    C push carry limb onto stack
        !            59:        testb   $1,%dl
        !            60:        jnz     L(1)                    C enter loop in the middle
        !            61:        movl    %ebx,%eax
        !            62:
        !            63:        ALIGN(8)
        !            64: L(oop):        movl    (%esi,%edx,4),%ebx      C load next higher limb
        !            65:        shrdl(  %cl, %ebx, %eax)        C compute result limb
        !            66:        movl    %eax,(%edi,%edx,4)      C store it
        !            67:        incl    %edx
        !            68: L(1):  movl    (%esi,%edx,4),%eax
        !            69:        shrdl(  %cl, %eax, %ebx)
        !            70:        movl    %ebx,(%edi,%edx,4)
        !            71:        incl    %edx
        !            72:        jnz     L(oop)
        !            73:
        !            74:        shrl    %cl,%eax                C compute most significant limb
        !            75:        movl    %eax,(%edi)             C store it
        !            76:
        !            77:        popl    %eax                    C pop carry limb
        !            78:
        !            79:        popl    %ebx
        !            80:        popl    %esi
        !            81:        popl    %edi
        !            82:        ret
        !            83:
        !            84: L(end):        shrl    %cl,%ebx                C compute most significant limb
        !            85:        movl    %ebx,(%edi)             C store it
        !            86:
        !            87:        popl    %ebx
        !            88:        popl    %esi
        !            89:        popl    %edi
        !            90:        ret
        !            91:
        !            92: EPILOGUE()