Return to mul_1.S CVS log

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

Annotation of OpenXM_contrib/gmp/mpn/arm/mul_1.S, Revision 1.1.1.1

1.1       maekawa     1: @ ARM mpn_addmul_1 -- Multiply a limb vector with a limb and store the result
                      2: @ in a second limb vector.
                      3: @ Contributed by Robert Harley.
                      4:
                      5: @ Copyright (C) 1998, 2000 Free Software Foundation, Inc.
                      6:
                      7: @ This file is part of the GNU MP Library.
                      8:
                      9: @ The GNU MP Library is free software; you can redistribute it and/or modify
                     10: @ it under the terms of the GNU Lesser General Public License as published by
                     11: @ the Free Software Foundation; either version 2.1 of the License, or (at your
                     12: @ option) any later version.
                     13:
                     14: @ The GNU MP Library is distributed in the hope that it will be useful, but
                     15: @ WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
                     16: @ or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
                     17: @ License for more details.
                     18:
                     19: @ You should have received a copy of the GNU Lesser General Public License
                     20: @ along with the GNU MP Library; see the file COPYING.LIB.  If not, write to
                     21: @ the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
                     22: @ MA 02111-1307, USA.
                     23:
                     24: #define p r0
                     25: #define a r1
                     26: #define n r2
                     27: #define w r3
                     28:
                     29: #define sl r10
                     30: #define fp r11
                     31: #define ip r12
                     32: #define sp r13
                     33: #define lr r14
                     34: #define pc r15
                     35:
                     36: .text
                     37:        .align  0
                     38:        .global __gmpn_mul_1
                     39:        .type   __gmpn_mul_1,%function
                     40: __gmpn_mul_1:
                     41:        stmfd   sp!, { r8, r9, lr }
                     42:        ands    ip, n, #1
                     43:        beq     skip1
                     44:        ldr     lr, [a], #4
                     45:        umull   r9, ip, w, lr
                     46:        str     r9, [p], #4
                     47: skip1:
                     48:        tst     n, #2
                     49:        beq     skip2
                     50:        mov     r8, ip
                     51:        ldmia   a!, { ip, lr }
                     52:        mov     r9, #0
                     53:        umlal   r8, r9, w, ip
                     54:        mov     ip, #0
                     55:        umlal   r9, ip, w, lr
                     56:        stmia   p!, { r8, r9 }
                     57: skip2:
                     58:        bics    n, n, #3
                     59:        beq     return
                     60:        stmfd   sp!, { r6, r7 }
                     61: mul_1_loop:
                     62:        mov     r6, ip
                     63:        ldmia   a!, { r8, r9, ip, lr }
                     64:        ldr     r7, [p] /* Bring stuff into cache. */
                     65:        mov     r7, #0
                     66:        umlal   r6, r7, w, r8
                     67:        mov     r8, #0
                     68:        umlal   r7, r8, w, r9
                     69:        mov     r9, #0
                     70:        umlal   r8, r9, w, ip
                     71:        mov     ip, #0
                     72:        umlal   r9, ip, w, lr
                     73:        subs    n, n, #4
                     74:        stmia   p!, { r6, r7, r8, r9 }
                     75:        bne     mul_1_loop
                     76:        ldmfd   sp!, { r6, r7 }
                     77: return:
                     78:        mov     r0, ip
                     79:        ldmfd   sp!, { r8, r9, pc }
                     80: end:
                     81:        .size   __gmpn_mul_1, end - __gmpn_mul_1