Return to invert_limb.asm CVS log

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

Annotation of OpenXM_contrib/gmp/mpn/arm/invert_limb.asm, Revision 1.1.1.1

1.1       ohara       1: dnl  ARM mpn_invert_limb -- Invert a normalized limb.
                      2:
                      3: dnl  Copyright 2001 Free Software Foundation, Inc.
                      4:
                      5: dnl  This file is part of the GNU MP Library.
                      6:
                      7: dnl  The GNU MP Library is free software; you can redistribute it and/or modify
                      8: dnl  it under the terms of the GNU Lesser General Public License as published
                      9: dnl  by the Free Software Foundation; either version 2.1 of the License, or (at
                     10: dnl  your option) any later version.
                     11:
                     12: dnl  The GNU MP Library is distributed in the hope that it will be useful, but
                     13: dnl  WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
                     14: dnl  or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
                     15: dnl  License for more details.
                     16:
                     17: dnl  You should have received a copy of the GNU Lesser General Public License
                     18: dnl  along with the GNU MP Library; see the file COPYING.LIB.  If not, write to
                     19: dnl  the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
                     20: dnl  MA 02111-1307, USA.
                     21:
                     22: include(`../config.m4')
                     23:
                     24: C INPUT PARAMETERS
                     25: define(`d',`r0')       C number to be inverted
                     26:
                     27:
                     28: PROLOGUE(mpn_invert_limb)
                     29:        stmfd   sp!, {r4, lr}
                     30:        mov     r3, d, lsr #23
                     31:        sub     r3, r3, #256
                     32:        add     r2, pc, #invtab-.-8
                     33:        mov     r3, r3, lsl #1
                     34:        ldrh    r1, [r2, r3]            C get initial approximation from table
                     35:        mov     r2, r1, lsl #6          C start iteration 1
                     36:        mul     ip, r2, r2
                     37:        umull   lr, r4, ip, d
                     38:        mov     r2, r4, lsl #1
                     39:        rsb     r2, r2, r1, lsl #23     C iteration 1 complete
                     40:        umull   ip, r3, r2, r2          C start iteration 2
                     41:        umull   lr, r4, r3, d
                     42:        umull   r3, r1, ip, d
                     43:        adds    lr, lr, r1
                     44:        addcs   r4, r4, #1
                     45:        mov     r3, lr, lsr #30
                     46:        orr     r4, r3, r4, lsl #2
                     47:        mov     lr, lr, lsl #2
                     48:        cmn     lr, #1
                     49:        rsc     r2, r4, r2, lsl #2      C iteration 2 complete
                     50:        umull   ip, r1, d, r2           C start adjustment step
                     51:        add     r1, r1, d
                     52:        cmn     r1, #1
                     53:        beq     L(1)
                     54:        adds    ip, ip, d
                     55:        adc     r1, r1, #0
                     56:        add     r2, r2, #1
                     57: L(1):
                     58:        adds    r3, ip, d
                     59:        adcs    r1, r1, #0
                     60:        moveq   r0, r2
                     61:        addne   r0, r2, #1
                     62:        ldmfd   sp!, {r4, pc}
                     63:
                     64: invtab:
                     65:        .short  1023,1020,1016,1012,1008,1004,1000,996
                     66:        .short  992,989,985,981,978,974,970,967
                     67:        .short  963,960,956,953,949,946,942,939
                     68:        .short  936,932,929,926,923,919,916,913
                     69:        .short  910,907,903,900,897,894,891,888
                     70:        .short  885,882,879,876,873,870,868,865
                     71:        .short  862,859,856,853,851,848,845,842
                     72:        .short  840,837,834,832,829,826,824,821
                     73:        .short  819,816,814,811,809,806,804,801
                     74:        .short  799,796,794,791,789,787,784,782
                     75:        .short  780,777,775,773,771,768,766,764
                     76:        .short  762,759,757,755,753,751,748,746
                     77:        .short  744,742,740,738,736,734,732,730
                     78:        .short  728,726,724,722,720,718,716,714
                     79:        .short  712,710,708,706,704,702,700,699
                     80:        .short  697,695,693,691,689,688,686,684
                     81:        .short  682,680,679,677,675,673,672,670
                     82:        .short  668,667,665,663,661,660,658,657
                     83:        .short  655,653,652,650,648,647,645,644
                     84:        .short  642,640,639,637,636,634,633,631
                     85:        .short  630,628,627,625,624,622,621,619
                     86:        .short  618,616,615,613,612,611,609,608
                     87:        .short  606,605,604,602,601,599,598,597
                     88:        .short  595,594,593,591,590,589,587,586
                     89:        .short  585,583,582,581,579,578,577,576
                     90:        .short  574,573,572,571,569,568,567,566
                     91:        .short  564,563,562,561,560,558,557,556
                     92:        .short  555,554,553,551,550,549,548,547
                     93:        .short  546,544,543,542,541,540,539,538
                     94:        .short  537,536,534,533,532,531,530,529
                     95:        .short  528,527,526,525,524,523,522,521
                     96:        .short  520,519,518,517,516,515,514,513
                     97: EPILOGUE(mpn_invert_limb)