[BACK]Return to submul_1.s CVS log [TXT][DIR] Up to [local] / OpenXM_contrib / gmp / mpn / a29k

Annotation of OpenXM_contrib/gmp/mpn/a29k/submul_1.s, Revision 1.1.1.1

1.1       maekawa     1: ; 29000 __mpn_submul_1 -- Multiply a limb vector with a single limb and
                      2: ; subtract the product 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      lr2
                     26: ; s1_ptr       lr3
                     27: ; size         lr4
                     28: ; s2_limb      lr5
                     29:
                     30:        .cputype 29050
                     31:        .sect .lit,lit
                     32:        .text
                     33:        .align  4
                     34:        .global ___mpn_submul_1
                     35:        .word   0x60000
                     36: ___mpn_submul_1:
                     37:        sub     lr4,lr4,8
                     38:        jmpt    lr4,Ltail
                     39:         const  gr120,0                 ; init cylimb reg
                     40:
                     41:        srl     gr117,lr4,3             ; divide by 8
                     42:        sub     gr117,gr117,1           ; count for jmpfdec
                     43:
                     44: Loop:  mtsrim  cr,(8-1)
                     45:        loadm   0,0,gr96,lr3
                     46:        add     lr3,lr3,32
                     47:
                     48:        multiplu gr104,gr96,lr5
                     49:        multmu   gr96,gr96,lr5
                     50:        multiplu gr105,gr97,lr5
                     51:        multmu   gr97,gr97,lr5
                     52:        multiplu gr106,gr98,lr5
                     53:        multmu   gr98,gr98,lr5
                     54:        multiplu gr107,gr99,lr5
                     55:        multmu   gr99,gr99,lr5
                     56:        multiplu gr108,gr100,lr5
                     57:        multmu   gr100,gr100,lr5
                     58:        multiplu gr109,gr101,lr5
                     59:        multmu   gr101,gr101,lr5
                     60:        multiplu gr110,gr102,lr5
                     61:        multmu   gr102,gr102,lr5
                     62:        multiplu gr111,gr103,lr5
                     63:        multmu   gr103,gr103,lr5
                     64:
                     65:        add     gr104,gr104,gr120
                     66:        addc    gr105,gr105,gr96
                     67:        addc    gr106,gr106,gr97
                     68:        addc    gr107,gr107,gr98
                     69:        addc    gr108,gr108,gr99
                     70:        addc    gr109,gr109,gr100
                     71:        addc    gr110,gr110,gr101
                     72:        addc    gr111,gr111,gr102
                     73:        addc    gr120,gr103,0
                     74:
                     75:        mtsrim  cr,(8-1)
                     76:        loadm   0,0,gr96,lr2
                     77:
                     78:        sub     gr96,gr96,gr104
                     79:        subc    gr97,gr97,gr105
                     80:        subc    gr98,gr98,gr106
                     81:        subc    gr99,gr99,gr107
                     82:        subc    gr100,gr100,gr108
                     83:        subc    gr101,gr101,gr109
                     84:        subc    gr102,gr102,gr110
                     85:        subc    gr103,gr103,gr111
                     86:
                     87:        add     gr104,gr103,gr111       ; invert carry from previus sub
                     88:        addc    gr120,gr120,0
                     89:
                     90:        mtsrim  cr,(8-1)
                     91:        storem  0,0,gr96,lr2
                     92:        jmpfdec gr117,Loop
                     93:         add    lr2,lr2,32
                     94:
                     95: Ltail: and     lr4,lr4,(8-1)
                     96:        sub     gr118,lr4,1             ; count for CR
                     97:        jmpt    gr118,Lend
                     98:         sub    lr4,lr4,2
                     99:        sub     lr2,lr2,4               ; offset res_ptr by one limb
                    100:
                    101: Loop2: load    0,0,gr116,lr3
                    102:        add     lr3,lr3,4
                    103:        multiplu gr117,gr116,lr5
                    104:        multmu  gr118,gr116,lr5
                    105:        add     lr2,lr2,4
                    106:        load    0,0,gr119,lr2
                    107:        add     gr117,gr117,gr120
                    108:        addc    gr118,gr118,0
                    109:        sub     gr119,gr119,gr117
                    110:        add     gr104,gr119,gr117       ; invert carry from previus sub
                    111:        store   0,0,gr119,lr2
                    112:        jmpfdec lr4,Loop2
                    113:         addc   gr120,gr118,0
                    114:
                    115: Lend:  jmpi    lr0
                    116:         or     gr96,gr120,0            ; copy

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>