glibc/sysdeps/alpha/udiv_qrnnd.S

160 lines
2.8 KiB
ArmAsm
Raw Normal View History

1995-02-17 20:13:09 +00:00
# Alpha 21064 __udiv_qrnnd
1996-03-01 18:43:45 +00:00
# Copyright (C) 1992, 1994, 1995 Free Software Foundation, Inc.
1995-02-17 20:13:09 +00:00
# This file is part of the GNU MP Library.
# The GNU MP Library is free software; you can redistribute it and/or modify
2001-07-07 10:13:33 +00:00
# it under the terms of the GNU Lesser General Public License as published by
# the Free Software Foundation; either version 2.1 of the License, or (at your
1995-02-17 20:13:09 +00:00
# option) any later version.
# The GNU MP Library is distributed in the hope that it will be useful, but
# WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
2001-07-07 10:13:33 +00:00
# or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
1995-02-17 20:13:09 +00:00
# License for more details.
2001-07-07 10:13:33 +00:00
# You should have received a copy of the GNU Lesser General Public License
2012-03-09 23:56:38 +00:00
# along with the GNU MP Library. If not, see <http://www.gnu.org/licenses/>.
1995-02-17 20:13:09 +00:00
#include <sysdep.h>
1995-02-17 20:13:09 +00:00
.set noreorder
.set noat
.text
LEAF(__udiv_qrnnd, 0)
1996-11-06 04:24:11 +00:00
#ifdef PROF
ldgp gp, 0(pv)
lda AT, _mcount
jsr AT, (AT), _mcount
.prologue 1
#else
1995-02-17 20:13:09 +00:00
.prologue 0
1996-11-06 04:24:11 +00:00
#endif
1995-02-17 20:13:09 +00:00
#define cnt $2
#define tmp $3
#define rem_ptr $16
#define n1 $17
#define n0 $18
#define d $19
#define qb $20
ldiq cnt,16
1996-11-06 04:24:11 +00:00
blt d,$largedivisor
1995-02-17 20:13:09 +00:00
1996-11-06 04:24:11 +00:00
$loop1: cmplt n0,0,tmp
1995-02-17 20:13:09 +00:00
addq n1,n1,n1
bis n1,tmp,n1
addq n0,n0,n0
cmpule d,n1,qb
subq n1,d,tmp
cmovne qb,tmp,n1
bis n0,qb,n0
cmplt n0,0,tmp
addq n1,n1,n1
bis n1,tmp,n1
addq n0,n0,n0
cmpule d,n1,qb
subq n1,d,tmp
cmovne qb,tmp,n1
bis n0,qb,n0
cmplt n0,0,tmp
addq n1,n1,n1
bis n1,tmp,n1
addq n0,n0,n0
cmpule d,n1,qb
subq n1,d,tmp
cmovne qb,tmp,n1
bis n0,qb,n0
cmplt n0,0,tmp
addq n1,n1,n1
bis n1,tmp,n1
addq n0,n0,n0
cmpule d,n1,qb
subq n1,d,tmp
cmovne qb,tmp,n1
bis n0,qb,n0
subq cnt,1,cnt
1996-11-06 04:24:11 +00:00
bgt cnt,$loop1
1995-02-17 20:13:09 +00:00
stq n1,0(rem_ptr)
bis $31,n0,$0
ret $31,($26),1
1996-11-06 04:24:11 +00:00
$largedivisor:
1995-02-17 20:13:09 +00:00
and n0,1,$4
srl n0,1,n0
sll n1,63,tmp
or tmp,n0,n0
srl n1,1,n1
and d,1,$6
srl d,1,$5
addq $5,$6,$5
1996-11-06 04:24:11 +00:00
$loop2: cmplt n0,0,tmp
1995-02-17 20:13:09 +00:00
addq n1,n1,n1
bis n1,tmp,n1
addq n0,n0,n0
cmpule $5,n1,qb
subq n1,$5,tmp
cmovne qb,tmp,n1
bis n0,qb,n0
cmplt n0,0,tmp
addq n1,n1,n1
bis n1,tmp,n1
addq n0,n0,n0
cmpule $5,n1,qb
subq n1,$5,tmp
cmovne qb,tmp,n1
bis n0,qb,n0
cmplt n0,0,tmp
addq n1,n1,n1
bis n1,tmp,n1
addq n0,n0,n0
cmpule $5,n1,qb
subq n1,$5,tmp
cmovne qb,tmp,n1
bis n0,qb,n0
cmplt n0,0,tmp
addq n1,n1,n1
bis n1,tmp,n1
addq n0,n0,n0
cmpule $5,n1,qb
subq n1,$5,tmp
cmovne qb,tmp,n1
bis n0,qb,n0
subq cnt,1,cnt
1996-11-06 04:24:11 +00:00
bgt cnt,$loop2
1995-02-17 20:13:09 +00:00
addq n1,n1,n1
addq $4,n1,n1
1996-11-06 04:24:11 +00:00
bne $6,$Odd
1995-02-17 20:13:09 +00:00
stq n1,0(rem_ptr)
bis $31,n0,$0
ret $31,($26),1
1996-11-06 04:24:11 +00:00
$Odd:
1995-10-14 02:29:43 +00:00
/* q' in n0. r' in n1 */
1995-02-17 20:13:09 +00:00
addq n1,n0,n1
1996-11-06 04:24:11 +00:00
1995-02-17 20:13:09 +00:00
cmpult n1,n0,tmp # tmp := carry from addq
1996-11-06 04:24:11 +00:00
subq n1,d,AT
addq n0,tmp,n0
cmovne tmp,AT,n1
cmpult n1,d,tmp
addq n0,1,AT
cmoveq tmp,AT,n0
subq n1,d,AT
cmoveq tmp,AT,n1
1995-02-17 20:13:09 +00:00
stq n1,0(rem_ptr)
bis $31,n0,$0
ret $31,($26),1
.end __udiv_qrnnd