158 lines
4.9 KiB
ArmAsm
158 lines
4.9 KiB
ArmAsm
/* libgcc1 routines for 68000 w/o floating-point hardware.
|
|
Copyright (C) 1994, 1996, 1997, 1998 Free Software Foundation, Inc.
|
|
|
|
This file is part of GNU CC.
|
|
|
|
GNU CC is free software; you can redistribute it and/or modify it
|
|
under the terms of the GNU General Public License as published by the
|
|
Free Software Foundation; either version 2, or (at your option) any
|
|
later version.
|
|
|
|
In addition to the permissions in the GNU General Public License, the
|
|
Free Software Foundation gives you unlimited permission to link the
|
|
compiled version of this file with other programs, and to distribute
|
|
those programs without any restriction coming from the use of this
|
|
file. (The General Public License restrictions do apply in other
|
|
respects; for example, they cover modification of the file, and
|
|
distribution when not linked into another program.)
|
|
|
|
This file is distributed in the hope that it will be useful, but
|
|
WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
|
|
General Public License for more details. */
|
|
|
|
/* As a special exception, if you link this library with files
|
|
compiled with GCC to produce an executable, this does not cause
|
|
the resulting executable to be covered by the GNU General Public License.
|
|
This exception does not however invalidate any other reasons why
|
|
the executable file might be covered by the GNU General Public License. */
|
|
|
|
/* Use this one for any 680x0; assumes no floating point hardware.
|
|
The trailing " '" appearing on some lines is for ANSI preprocessors. Yuk.
|
|
Some of this code comes from MINIX, via the folks at ericsson.
|
|
D. V. Henkel-Wallace (gumby@cygnus.com) Fete Bastille, 1992
|
|
*/
|
|
#include <asm/export.h>
|
|
/* These are predefined by new versions of GNU cpp. */
|
|
|
|
#ifndef __USER_LABEL_PREFIX__
|
|
#define __USER_LABEL_PREFIX__ _
|
|
#endif
|
|
|
|
#ifndef __REGISTER_PREFIX__
|
|
#define __REGISTER_PREFIX__
|
|
#endif
|
|
|
|
#ifndef __IMMEDIATE_PREFIX__
|
|
#define __IMMEDIATE_PREFIX__ #
|
|
#endif
|
|
|
|
/* ANSI concatenation macros. */
|
|
|
|
#define CONCAT1(a, b) CONCAT2(a, b)
|
|
#define CONCAT2(a, b) a ## b
|
|
|
|
/* Use the right prefix for global labels. */
|
|
|
|
#define SYM(x) CONCAT1 (__USER_LABEL_PREFIX__, x)
|
|
|
|
/* Use the right prefix for registers. */
|
|
|
|
#define REG(x) CONCAT1 (__REGISTER_PREFIX__, x)
|
|
|
|
/* Use the right prefix for immediate values. */
|
|
|
|
#define IMM(x) CONCAT1 (__IMMEDIATE_PREFIX__, x)
|
|
|
|
#define d0 REG (d0)
|
|
#define d1 REG (d1)
|
|
#define d2 REG (d2)
|
|
#define d3 REG (d3)
|
|
#define d4 REG (d4)
|
|
#define d5 REG (d5)
|
|
#define d6 REG (d6)
|
|
#define d7 REG (d7)
|
|
#define a0 REG (a0)
|
|
#define a1 REG (a1)
|
|
#define a2 REG (a2)
|
|
#define a3 REG (a3)
|
|
#define a4 REG (a4)
|
|
#define a5 REG (a5)
|
|
#define a6 REG (a6)
|
|
#define fp REG (fp)
|
|
#define sp REG (sp)
|
|
|
|
.text
|
|
.proc
|
|
.globl SYM (__udivsi3)
|
|
SYM (__udivsi3):
|
|
#if !(defined(__mcf5200__) || defined(__mcoldfire__))
|
|
movel d2, sp@-
|
|
movel sp@(12), d1 /* d1 = divisor */
|
|
movel sp@(8), d0 /* d0 = dividend */
|
|
|
|
cmpl IMM (0x10000), d1 /* divisor >= 2 ^ 16 ? */
|
|
jcc L3 /* then try next algorithm */
|
|
movel d0, d2
|
|
clrw d2
|
|
swap d2
|
|
divu d1, d2 /* high quotient in lower word */
|
|
movew d2, d0 /* save high quotient */
|
|
swap d0
|
|
movew sp@(10), d2 /* get low dividend + high rest */
|
|
divu d1, d2 /* low quotient */
|
|
movew d2, d0
|
|
jra L6
|
|
|
|
L3: movel d1, d2 /* use d2 as divisor backup */
|
|
L4: lsrl IMM (1), d1 /* shift divisor */
|
|
lsrl IMM (1), d0 /* shift dividend */
|
|
cmpl IMM (0x10000), d1 /* still divisor >= 2 ^ 16 ? */
|
|
jcc L4
|
|
divu d1, d0 /* now we have 16 bit divisor */
|
|
andl IMM (0xffff), d0 /* mask out divisor, ignore remainder */
|
|
|
|
/* Multiply the 16 bit tentative quotient with the 32 bit divisor. Because of
|
|
the operand ranges, this might give a 33 bit product. If this product is
|
|
greater than the dividend, the tentative quotient was too large. */
|
|
movel d2, d1
|
|
mulu d0, d1 /* low part, 32 bits */
|
|
swap d2
|
|
mulu d0, d2 /* high part, at most 17 bits */
|
|
swap d2 /* align high part with low part */
|
|
tstw d2 /* high part 17 bits? */
|
|
jne L5 /* if 17 bits, quotient was too large */
|
|
addl d2, d1 /* add parts */
|
|
jcs L5 /* if sum is 33 bits, quotient was too large */
|
|
cmpl sp@(8), d1 /* compare the sum with the dividend */
|
|
jls L6 /* if sum > dividend, quotient was too large */
|
|
L5: subql IMM (1), d0 /* adjust quotient */
|
|
|
|
L6: movel sp@+, d2
|
|
rts
|
|
|
|
#else /* __mcf5200__ || __mcoldfire__ */
|
|
|
|
/* Coldfire implementation of non-restoring division algorithm from
|
|
Hennessy & Patterson, Appendix A. */
|
|
link a6,IMM (-12)
|
|
moveml d2-d4,sp@
|
|
movel a6@(8),d0
|
|
movel a6@(12),d1
|
|
clrl d2 | clear p
|
|
moveq IMM (31),d4
|
|
L1: addl d0,d0 | shift reg pair (p,a) one bit left
|
|
addxl d2,d2
|
|
movl d2,d3 | subtract b from p, store in tmp.
|
|
subl d1,d3
|
|
jcs L2 | if no carry,
|
|
bset IMM (0),d0 | set the low order bit of a to 1,
|
|
movl d3,d2 | and store tmp in p.
|
|
L2: subql IMM (1),d4
|
|
jcc L1
|
|
moveml sp@,d2-d4 | restore data registers
|
|
unlk a6 | and return
|
|
rts
|
|
#endif /* __mcf5200__ || __mcoldfire__ */
|
|
EXPORT_SYMBOL(__udivsi3)
|