dnl SPARC v9 64-bit 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/.
include(`../config.m4')
C cycles/limb
C UltraSPARC 1&2: 22
C UltraSPARC 3: 36
C This was generated by the Sun C compiler. It runs at 22 cycles/limb on the
C UltraSPARC-1/2, three cycles slower than theoretically possible for optimal
C code using the same algorithm. For 1-3 limbs, a special loop was generated,
C which causes performance problems in particular for 2 and 3 limbs.
C Ultimately, this should be replaced by hand-written code in the same software
C pipeline style as e.g., addmul_1.asm.
ASM_START()
REGISTER(%g2,#scratch)
REGISTER(%g3,#scratch)
PROLOGUE(mpn_sqr_diagonal)
save %sp, -240, %sp
sethi %hi(0x1ffc00), %o0
sethi %hi(0x3ffc00), %o1
add %o0, 1023, %o7
cmp %i2, 4
add %o1, 1023, %o4
or %g0, %i1, %g1
or %g0, %i0, %o0
bl,pn %xcc, .Lsmall
or %g0, 0, %g2
ldx [%i1], %o1
add %i1, 24, %g1
or %g0, 3, %g2
srlx %o1, 42, %g3
stx %g3, [%sp+2279]
and %o1, %o7, %o2
stx %o2, [%sp+2263]
srlx %o1, 21, %o1
ldd [%sp+2279], %f0
and %o1, %o7, %o1
stx %o1, [%sp+2271]
ldx [%i1+8], %o2
fxtod %f0, %f12
srlx %o2, 21, %o1
and %o2, %o7, %g3
ldd [%sp+2263], %f2
fmuld %f12, %f12, %f10
srlx %o2, 42, %o2
ldd [%sp+2271], %f0
and %o1, %o7, %o1
fxtod %f2, %f8
stx %o2, [%sp+2279]
stx %o1, [%sp+2271]
fxtod %f0, %f0
stx %g3, [%sp+2263]
fdtox %f10, %f14
fmuld %f12, %f8, %f6
ldx [%i1+16], %o2
std %f14, [%sp+2255]
fmuld %f0, %f0, %f2
fmuld %f8, %f8, %f10
srlx %o2, 42, %o1
faddd %f6, %f6, %f6
fmuld %f12, %f0, %f12
fmuld %f0, %f8, %f8
ldd [%sp+2279], %f0
ldd [%sp+2263], %f4
fdtox %f10, %f10
std %f10, [%sp+2239]
faddd %f2, %f6, %f6
ldd [%sp+2271], %f2
fdtox %f12, %f12
std %f12, [%sp+2247]
fdtox %f8, %f8
std %f8, [%sp+2231]
fdtox %f6, %f6
std %f6, [%sp+2223]
.Loop: srlx %o2, 21, %g3
stx %o1, [%sp+2279]
add %g2, 1, %g2
and %g3, %o7, %o1
ldx [%sp+2255], %g4
cmp %g2, %i2
stx %o1, [%sp+2271]
add %g1, 8, %g1
add %o0, 16, %o0
ldx [%sp+2239], %o1
fxtod %f0, %f10
fxtod %f4, %f14
ldx [%sp+2231], %i0
ldx [%sp+2223], %g5
ldx [%sp+2247], %g3
and %o2, %o7, %o2
fxtod %f2, %f8
fmuld %f10, %f10, %f0
stx %o2, [%sp+2263]
fmuld %f10, %f14, %f6
ldx [%g1-8], %o2
fmuld %f10, %f8, %f12
fdtox %f0, %f2
ldd [%sp+2279], %f0
fmuld %f8, %f8, %f4
faddd %f6, %f6, %f6
fmuld %f14, %f14, %f10
std %f2, [%sp+2255]
sllx %g4, 20, %g4
ldd [%sp+2271], %f2
fmuld %f8, %f14, %f8
sllx %i0, 22, %i1
fdtox %f12, %f12
std %f12, [%sp+2247]
sllx %g5, 42, %i0
add %o1, %i1, %o1
faddd %f4, %f6, %f6
ldd [%sp+2263], %f4
add %o1, %i0, %o1
add %g3, %g4, %g3
fdtox %f10, %f10
std %f10, [%sp+2239]
srlx %o1, 42, %g4
and %g5, %o4, %i0
fdtox %f8, %f8
std %f8, [%sp+2231]
srlx %g5, 22, %g5
sub %g4, %i0, %g4
fdtox %f6, %f6
std %f6, [%sp+2223]
srlx %g4, 63, %g4
add %g3, %g5, %g3
add %g3, %g4, %g3
stx %o1, [%o0-16]
srlx %o2, 42, %o1
bl,pt %xcc, .Loop
stx %g3, [%o0-8]
stx %o1, [%sp+2279]
srlx %o2, 21, %o1
fxtod %f0, %f16
ldx [%sp+2223], %g3
fxtod %f4, %f6
and %o2, %o7, %o3
stx %o3, [%sp+2263]
fxtod %f2, %f4
and %o1, %o7, %o1
ldx [%sp+2231], %o2
sllx %g3, 42, %g4
fmuld %f16, %f16, %f14
stx %o1, [%sp+2271]
fmuld %f16, %f6, %f8
add %o0, 48, %o0
ldx [%sp+2239], %o1
sllx %o2, 22, %o2
fmuld %f4, %f4, %f10
ldx [%sp+2255], %o3
fdtox %f14, %f14
fmuld %f4, %f6, %f2
std %f14, [%sp+2255]
faddd %f8, %f8, %f12
add %o1, %o2, %o2
fmuld %f16, %f4, %f4
ldd [%sp+2279], %f0
sllx %o3, 20, %g5
add %o2, %g4, %o2
fmuld %f6, %f6, %f6
srlx %o2, 42, %o3
and %g3, %o4, %g4
srlx %g3, 22, %g3
faddd %f10, %f12, %f16
ldd [%sp+2271], %f12
ldd [%sp+2263], %f8
fxtod %f0, %f0
sub %o3, %g4, %o3
ldx [%sp+2247], %o1
srlx %o3, 63, %o3
fdtox %f2, %f10
fxtod %f8, %f8
std %f10, [%sp+2231]
fdtox %f6, %f6
std %f6, [%sp+2239]
add %o1, %g5, %o1
fmuld %f0, %f0, %f2
fdtox %f16, %f16
std %f16, [%sp+2223]
add %o1, %g3, %o1
fdtox %f4, %f4
std %f4, [%sp+2247]
fmuld %f0, %f8, %f10
fxtod %f12, %f12
add %o1, %o3, %o1
stx %o2, [%o0-48]
fmuld %f8, %f8, %f6
stx %o1, [%o0-40]
fdtox %f2, %f2
ldx [%sp+2231], %o2
faddd %f10, %f10, %f10
ldx [%sp+2223], %g3
fmuld %f12, %f12, %f4
fdtox %f6, %f6
ldx [%sp+2239], %o1
sllx %o2, 22, %o2
fmuld %f12, %f8, %f8
sllx %g3, 42, %g5
ldx [%sp+2255], %o3
fmuld %f0, %f12, %f0
add %o1, %o2, %o2
faddd %f4, %f10, %f4
ldx [%sp+2247], %o1
add %o2, %g5, %o2
and %g3, %o4, %g4
fdtox %f8, %f8
sllx %o3, 20, %g5
std %f8, [%sp+2231]
fdtox %f0, %f0
srlx %o2, 42, %o3
add %o1, %g5, %o1
fdtox %f4, %f4
srlx %g3, 22, %g3
sub %o3, %g4, %o3
std %f6, [%sp+2239]
std %f4, [%sp+2223]
srlx %o3, 63, %o3
add %o1, %g3, %o1
std %f2, [%sp+2255]
add %o1, %o3, %o1
std %f0, [%sp+2247]
stx %o2, [%o0-32]
stx %o1, [%o0-24]
ldx [%sp+2231], %o2
ldx [%sp+2223], %o3
ldx [%sp+2239], %o1
sllx %o2, 22, %o2
sllx %o3, 42, %g5
ldx [%sp+2255], %g4
and %o3, %o4, %g3
add %o1, %o2, %o2
ldx [%sp+2247], %o1
add %o2, %g5, %o2
stx %o2, [%o0-16]
sllx %g4, 20, %g4
srlx %o2, 42, %o2
add %o1, %g4, %o1
srlx %o3, 22, %o3
sub %o2, %g3, %o2
srlx %o2, 63, %o2
add %o1, %o3, %o1
add %o1, %o2, %o1
stx %o1, [%o0-8]
ret
restore %g0, %g0, %g0
.Lsmall:
ldx [%g1], %o2
.Loop0:
and %o2, %o7, %o1
stx %o1, [%sp+2263]
add %g2, 1, %g2
srlx %o2, 21, %o1
add %g1, 8, %g1
srlx %o2, 42, %o2
stx %o2, [%sp+2279]
and %o1, %o7, %o1
ldd [%sp+2263], %f0
cmp %g2, %i2
stx %o1, [%sp+2271]
fxtod %f0, %f6
ldd [%sp+2279], %f0
ldd [%sp+2271], %f4
fxtod %f0, %f2
fmuld %f6, %f6, %f0
fxtod %f4, %f10
fmuld %f2, %f6, %f4
fdtox %f0, %f0
std %f0, [%sp+2239]
fmuld %f10, %f6, %f8
fmuld %f10, %f10, %f0
faddd %f4, %f4, %f6
fmuld %f2, %f2, %f4
fdtox %f8, %f8
std %f8, [%sp+2231]
fmuld %f2, %f10, %f2
faddd %f0, %f6, %f0
fdtox %f4, %f4
std %f4, [%sp+2255]
fdtox %f2, %f2
std %f2, [%sp+2247]
fdtox %f0, %f0
std %f0, [%sp+2223]
ldx [%sp+2239], %o1
ldx [%sp+2255], %g4
ldx [%sp+2231], %o2
sllx %g4, 20, %g4
ldx [%sp+2223], %o3
sllx %o2, 22, %o2
sllx %o3, 42, %g5
add %o1, %o2, %o2
ldx [%sp+2247], %o1
add %o2, %g5, %o2
stx %o2, [%o0]
and %o3, %o4, %g3
srlx %o2, 42, %o2
add %o1, %g4, %o1
srlx %o3, 22, %o3
sub %o2, %g3, %o2
srlx %o2, 63, %o2
add %o1, %o3, %o1
add %o1, %o2, %o1
stx %o1, [%o0+8]
add %o0, 16, %o0
bl,a,pt %xcc, .Loop0
ldx [%g1], %o2
ret
restore %g0, %g0, %g0
EPILOGUE(mpn_sqr_diagonal)
syntax highlighted by Code2HTML, v. 0.9.1