diff options
Diffstat (limited to 'gmp/mpn/powerpc32/sqr_diagonal.asm')
| -rw-r--r-- | gmp/mpn/powerpc32/sqr_diagonal.asm | 103 |
1 files changed, 103 insertions, 0 deletions
diff --git a/gmp/mpn/powerpc32/sqr_diagonal.asm b/gmp/mpn/powerpc32/sqr_diagonal.asm new file mode 100644 index 0000000000..d315349f63 --- /dev/null +++ b/gmp/mpn/powerpc32/sqr_diagonal.asm @@ -0,0 +1,103 @@ +dnl PowerPC-32 mpn_sqr_diagonal. + +dnl Copyright 2001, 2002, 2003, 2005 Free Software Foundation, Inc. + +dnl This file is part of the GNU MP Library. + +dnl The GNU MP Library is free software; you can redistribute it and/or modify +dnl it under the terms of the GNU Lesser General Public License as published +dnl by the Free Software Foundation; either version 3 of the License, or (at +dnl your option) any later version. + +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 Lesser General Public +dnl License for more details. + +dnl You should have received a copy of the GNU Lesser General Public License +dnl along with the GNU MP Library. If not, see http://www.gnu.org/licenses/. + +include(`../config.m4') + +C cycles/limb +C 603e: ? +C 604e: 4.0 +C 75x (G3): 10.5 +C 7400,7410 (G4): 10.5 +C 744x,745x (G4+): 4.0 +C power4/ppc970: 8.6 +C power5: 7.0 + +C INPUT PARAMETERS +C rp r3 +C up r4 +C n r5 + +ASM_START() +PROLOGUE(mpn_sqr_diagonal) + lwz r6,0(r4) + mtctr r5 + + addi r3,r3,-4 + bdz L(end1) + + lwzu r7,4(r4) + mullw r9,r6,r6 + mulhwu r11,r6,r6 + bdz L(end2) + + lwzu r6,4(r4) + mullw r8,r7,r7 + mulhwu r10,r7,r7 + bdz L(ende) + +L(loop): + lwzu r7,4(r4) + stw r9,4(r3) + mullw r9,r6,r6 + stwu r11,8(r3) + mulhwu r11,r6,r6 + bdz L(endo) + lwzu r6,4(r4) + stw r8,4(r3) + mullw r8,r7,r7 + stwu r10,8(r3) + mulhwu r10,r7,r7 + bdnz L(loop) + +L(ende): + stw r9,4(r3) + mullw r9,r6,r6 + stw r11,8(r3) + mulhwu r11,r6,r6 + stw r8,12(r3) + stw r10,16(r3) + stw r9,20(r3) + stw r11,24(r3) + blr +L(endo): + stw r8,4(r3) + mullw r8,r7,r7 + stw r10,8(r3) + mulhwu r10,r7,r7 + stw r9,12(r3) + stw r11,16(r3) + stw r8,20(r3) + stw r10,24(r3) + blr + +L(end2): + mullw r8,r7,r7 + stw r9,4(r3) + mulhwu r10,r7,r7 + stw r11,8(r3) + stw r8,12(r3) + stw r10,16(r3) + blr +L(end1): + mullw r9,r6,r6 + mulhwu r11,r6,r6 + stw r9,4(r3) + stw r11,8(r3) + blr +EPILOGUE(mpn_sqr_diagonal) |