| dnl PowerPC-64 mpn_gcd_1. |
| |
| dnl Copyright 2000-2002, 2005, 2009, 2011-2013 Free Software Foundation, Inc. |
| |
| dnl This file is part of the GNU MP Library. |
| dnl |
| dnl The GNU MP Library is free software; you can redistribute it and/or modify |
| dnl it under the terms of either: |
| dnl |
| dnl * the GNU Lesser General Public License as published by the Free |
| dnl Software Foundation; either version 3 of the License, or (at your |
| dnl option) any later version. |
| dnl |
| dnl or |
| dnl |
| dnl * the GNU General Public License as published by the Free Software |
| dnl Foundation; either version 2 of the License, or (at your option) any |
| dnl later version. |
| dnl |
| dnl or both in parallel, as here. |
| dnl |
| dnl The GNU MP Library is distributed in the hope that it will be useful, but |
| dnl WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY |
| dnl or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License |
| dnl for more details. |
| dnl |
| dnl You should have received copies of the GNU General Public License and the |
| dnl GNU Lesser General Public License along with the GNU MP Library. If not, |
| dnl see https://www.gnu.org/licenses/. |
| |
| include(`../config.m4') |
| |
| C cycles/bit (approx) |
| C POWER3/PPC630 ? |
| C POWER4/PPC970 8.5 |
| C POWER5 ? |
| C POWER6 10.1 |
| C POWER7 9.4 |
| C Numbers measured with: speed -CD -s16-64 -t48 mpn_gcd_1 |
| |
| C INPUT PARAMETERS |
| define(`up', `r3') |
| define(`n', `r4') |
| define(`v0', `r5') |
| |
| ifdef(`BMOD_1_TO_MOD_1_THRESHOLD',, |
| `define(`BMOD_1_TO_MOD_1_THRESHOLD',30)') |
| |
| EXTERN_FUNC(mpn_mod_1) |
| EXTERN_FUNC(mpn_modexact_1c_odd) |
| |
| ASM_START() |
| PROLOGUE(mpn_gcd_1,toc) |
| mflr r0 |
| std r30, -16(r1) |
| std r31, -8(r1) |
| std r0, 16(r1) |
| stdu r1, -128(r1) |
| |
| ld r7, 0(up) C U low limb |
| or r0, r5, r7 C x | y |
| |
| neg r6, r0 |
| and r6, r6, r0 |
| cntlzd r31, r6 C common twos |
| subfic r31, r31, 63 |
| |
| neg r6, r5 |
| and r6, r6, r5 |
| cntlzd r8, r6 |
| subfic r8, r8, 63 |
| srd r5, r5, r8 |
| mr r30, r5 C v0 saved |
| |
| cmpdi r4, BMOD_1_TO_MOD_1_THRESHOLD |
| blt L(bmod) |
| CALL( mpn_mod_1) |
| b L(reduced) |
| L(bmod): |
| li r6, 0 |
| CALL( mpn_modexact_1c_odd) |
| L(reduced): |
| |
| define(`mask', `r0')dnl |
| define(`a1', `r4')dnl |
| define(`a2', `r5')dnl |
| define(`d1', `r6')dnl |
| define(`d2', `r7')dnl |
| define(`cnt', `r9')dnl |
| |
| neg. r6, r3 |
| and r6, r6, r3 |
| cntlzd cnt, r6 |
| subfic cnt, cnt, 63 |
| li r12, 63 |
| bne L(mid) |
| b L(end) |
| |
| ALIGN(16) |
| L(top): |
| and a1, r10, mask C d - a |
| andc a2, r11, mask C a - d |
| and d1, r3, mask C a |
| andc d2, r30, mask C d |
| or r3, a1, a2 C new a |
| subf cnt, cnt, r12 |
| or r30, d1, d2 C new d |
| L(mid): srd r3, r3, cnt |
| sub. r10, r30, r3 C r10 = d - a |
| subc r11, r3, r30 C r11 = a - d |
| neg r8, r10 |
| and r8, r8, r10 |
| subfe mask, mask, mask |
| cntlzd cnt, r8 |
| bne L(top) |
| |
| L(end): sld r3, r30, r31 |
| |
| addi r1, r1, 128 |
| ld r0, 16(r1) |
| ld r30, -16(r1) |
| ld r31, -8(r1) |
| mtlr r0 |
| blr |
| EPILOGUE() |
