Imported gcc-4.4.3

[msp430-gcc.git] / gmp / mpn / mips64 / sqr_diagonal.asm

diff --git a/gmp/mpn/mips64/sqr_diagonal.asm b/gmp/mpn/mips64/sqr_diagonal.asm
new file mode 100644 (file)
index 0000000..511a755
--- /dev/null
+++ b/gmp/mpn/mips64/sqr_diagonal.asm
@@ -0,0 +1,66 @@
+dnl  MIPS64 mpn_sqr_diagonal.
+
+dnl  Copyright 2001, 2002 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/.
+
+
+dnl  INPUT PARAMETERS
+dnl  rp                $4
+dnl  up                $5
+dnl  n         $6
+
+include(`../config.m4')
+
+ASM_START()
+PROLOGUE(mpn_sqr_diagonal)
+       ld      r8,0(r5)
+       daddiu  r6,r6,-2
+       dmultu  r8,r8
+       bltz    r6,$Lend1
+       nop
+       ld      r8,8(r5)
+       beq     r6,r0,$Lend2
+       nop
+
+$Loop: mflo    r10
+       mfhi    r9
+       daddiu  r6,r6,-1
+       sd      r10,0(r4)
+       sd      r9,8(r4)
+       dmultu  r8,r8
+       ld      r8,16(r5)
+       daddiu  r5,r5,8
+       bne     r6,r0,$Loop
+       daddiu  r4,r4,16
+
+$Lend2: mflo   r10
+       mfhi    r9
+       sd      r10,0(r4)
+       sd      r9,8(r4)
+       dmultu  r8,r8
+       mflo    r10
+       mfhi    r9
+       sd      r10,16(r4)
+       j       r31
+       sd      r9,24(r4)
+
+$Lend1: mflo   r10
+       mfhi    r9
+       sd      r10,0(r4)
+       j       r31
+       sd      r9,8(r4)
+EPILOGUE(mpn_sqr_diagonal)