Return to add_n.s CVS log

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

Annotation of OpenXM_contrib/gmp/mpn/m88k/add_n.s, Revision 1.1

1.1     ! maekawa     1: ; mc88100 __mpn_add -- Add two limb vectors of the same length > 0 and store
        !             2: ; sum in a third limb vector.
        !             3:
        !             4: ; Copyright (C) 1992, 1994, 1995 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      r2
        !            26: ; s1_ptr       r3
        !            27: ; s2_ptr       r4
        !            28: ; size         r5
        !            29:
        !            30: ; This code has been optimized to run one instruction per clock, avoiding
        !            31: ; load stalls and writeback contention.  As a result, the instruction
        !            32: ; order is not always natural.
        !            33:
        !            34: ; The speed is about 4.6 clocks/limb + 18 clocks/limb-vector on an 88100,
        !            35: ; but on the 88110, it seems to run much slower, 6.6 clocks/limb.
        !            36:
        !            37:        text
        !            38:        align    16
        !            39:        global   ___mpn_add_n
        !            40: ___mpn_add_n:
        !            41:        ld      r6,r3,0                 ; read first limb from s1_ptr
        !            42:        extu    r10,r5,3
        !            43:        ld      r7,r4,0                 ; read first limb from s2_ptr
        !            44:
        !            45:        subu.co r5,r0,r5                ; (clear carry as side effect)
        !            46:        mak     r5,r5,3<4>
        !            47:        bcnd    eq0,r5,Lzero
        !            48:
        !            49:        or      r12,r0,lo16(Lbase)
        !            50:        or.u    r12,r12,hi16(Lbase)
        !            51:        addu    r12,r12,r5              ; r12 is address for entering in loop
        !            52:
        !            53:        extu    r5,r5,2                 ; divide by 4
        !            54:        subu    r2,r2,r5                ; adjust res_ptr
        !            55:        subu    r3,r3,r5                ; adjust s1_ptr
        !            56:        subu    r4,r4,r5                ; adjust s2_ptr
        !            57:
        !            58:        or      r8,r6,r0
        !            59:
        !            60:        jmp.n   r12
        !            61:         or     r9,r7,r0
        !            62:
        !            63: Loop:  addu    r3,r3,32
        !            64:        st      r8,r2,28
        !            65:        addu    r4,r4,32
        !            66:        ld      r6,r3,0
        !            67:        addu    r2,r2,32
        !            68:        ld      r7,r4,0
        !            69: Lzero: subu    r10,r10,1               ; add 0 + 8r limbs (adj loop cnt)
        !            70: Lbase: ld      r8,r3,4
        !            71:        addu.cio r6,r6,r7
        !            72:        ld      r9,r4,4
        !            73:        st      r6,r2,0
        !            74:        ld      r6,r3,8                 ; add 7 + 8r limbs
        !            75:        addu.cio r8,r8,r9
        !            76:        ld      r7,r4,8
        !            77:        st      r8,r2,4
        !            78:        ld      r8,r3,12                ; add 6 + 8r limbs
        !            79:        addu.cio r6,r6,r7
        !            80:        ld      r9,r4,12
        !            81:        st      r6,r2,8
        !            82:        ld      r6,r3,16                ; add 5 + 8r limbs
        !            83:        addu.cio r8,r8,r9
        !            84:        ld      r7,r4,16
        !            85:        st      r8,r2,12
        !            86:        ld      r8,r3,20                ; add 4 + 8r limbs
        !            87:        addu.cio r6,r6,r7
        !            88:        ld      r9,r4,20
        !            89:        st      r6,r2,16
        !            90:        ld      r6,r3,24                ; add 3 + 8r limbs
        !            91:        addu.cio r8,r8,r9
        !            92:        ld      r7,r4,24
        !            93:        st      r8,r2,20
        !            94:        ld      r8,r3,28                ; add 2 + 8r limbs
        !            95:        addu.cio r6,r6,r7
        !            96:        ld      r9,r4,28
        !            97:        st      r6,r2,24
        !            98:        bcnd.n  ne0,r10,Loop            ; add 1 + 8r limbs
        !            99:         addu.cio r8,r8,r9
        !           100:
        !           101:        st      r8,r2,28                ; store most significant limb
        !           102:
        !           103:        jmp.n    r1
        !           104:         addu.ci r2,r0,r0               ; return carry-out from most sign. limb