1dnl PowerPC-32 mpn_sqr_diag_addlsh1. 2 3dnl Copyright 2013 Free Software Foundation, Inc. 4 5dnl This file is part of the GNU MP Library. 6dnl 7dnl The GNU MP Library is free software; you can redistribute it and/or modify 8dnl it under the terms of either: 9dnl 10dnl * the GNU Lesser General Public License as published by the Free 11dnl Software Foundation; either version 3 of the License, or (at your 12dnl option) any later version. 13dnl 14dnl or 15dnl 16dnl * the GNU General Public License as published by the Free Software 17dnl Foundation; either version 2 of the License, or (at your option) any 18dnl later version. 19dnl 20dnl or both in parallel, as here. 21dnl 22dnl The GNU MP Library is distributed in the hope that it will be useful, but 23dnl WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 24dnl or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License 25dnl for more details. 26dnl 27dnl You should have received copies of the GNU General Public License and the 28dnl GNU Lesser General Public License along with the GNU MP Library. If not, 29dnl see https://www.gnu.org/licenses/. 30 31include(`../config.m4') 32 33C cycles/limb 34C 603e ? 35C 604e ? 36C 75x (G3) ? 37C 7400,7410 (G4) ? 38C 744x,745x (G4+) 6 39C power4/ppc970 ? 40C power5 ? 41 42C This has been feebly optimised for 7447 but not for any other CPU. 43 44define(`rp', r3) 45define(`tp', r4) 46define(`up', r5) 47define(`n', r6) 48 49ASM_START() 50PROLOGUE(mpn_sqr_diag_addlsh1) 51 addi n, n, -1 52 addi tp, tp, -4 53 mtctr n 54 lwz r0, 0(up) 55 li r10, 0 56 mullw r7, r0, r0 57 stw r7, 0(rp) 58 mulhwu r6, r0, r0 59 addic r31, r31, 0 C clear CF 60 61 ALIGN(16) 62L(top): lwzu r0, 4(up) 63 mullw r7, r0, r0 64 lwz r8, 4(tp) 65 lwzu r9, 8(tp) 66 rlwimi r10, r8, 1,0,30 67 srwi r11, r8, 31 68 rlwimi r11, r9, 1,0,30 69 adde r10, r10, r6 70 adde r11, r11, r7 71 stw r10, 4(rp) 72 srwi r10, r9, 31 73 mulhwu r6, r0, r0 74 stwu r11, 8(rp) 75 bdnz L(top) 76 77 adde r10, r10, r6 78 stw r10, 4(rp) 79 blr 80EPILOGUE() 81