1 /* SPDX-License-Identifier: GPL-2.0-or-later W !! 1 /* libgcc1 routines for 68000 w/o floating-point hardware. 2 #include <linux/linkage.h> !! 2 Copyright (C) 1994, 1996, 1997, 1998 Free Software Foundation, Inc. 3 #include <asm/asmmacro.h> << 4 #include <asm/core.h> << 5 << 6 ENTRY(__udivsi3) << 7 << 8 abi_entry_default << 9 #if XCHAL_HAVE_DIV32 << 10 quou a2, a2, a3 << 11 #else << 12 bltui a3, 2, .Lle_one /* check if th << 13 << 14 mov a6, a2 /* keep divide << 15 do_nsau a5, a6, a2, a7 /* dividend_sh << 16 do_nsau a4, a3, a2, a7 /* divisor_shi << 17 bgeu a5, a4, .Lspecial << 18 << 19 sub a4, a4, a5 /* count = div << 20 ssl a4 << 21 sll a3, a3 /* divisor <<= << 22 movi a2, 0 /* quotient = << 23 << 24 /* test-subtract-and-shift loop; one q << 25 #if XCHAL_HAVE_LOOPS << 26 loopnez a4, .Lloopend << 27 #endif /* XCHAL_HAVE_LOOPS */ << 28 .Lloop: << 29 bltu a6, a3, .Lzerobit << 30 sub a6, a6, a3 << 31 addi a2, a2, 1 << 32 .Lzerobit: << 33 slli a2, a2, 1 << 34 srli a3, a3, 1 << 35 #if !XCHAL_HAVE_LOOPS << 36 addi a4, a4, -1 << 37 bnez a4, .Lloop << 38 #endif /* !XCHAL_HAVE_LOOPS */ << 39 .Lloopend: << 40 << 41 bltu a6, a3, .Lreturn << 42 addi a2, a2, 1 /* increment q << 43 .Lreturn: << 44 abi_ret_default << 45 << 46 .Lle_one: << 47 beqz a3, .Lerror /* if divisor << 48 abi_ret_default << 49 << 50 .Lspecial: << 51 /* return dividend >= divisor */ << 52 bltu a6, a3, .Lreturn0 << 53 movi a2, 1 << 54 abi_ret_default << 55 << 56 .Lerror: << 57 /* Divide by zero: Use an illegal inst << 58 The subsequent "DIV0" string can be << 59 handler to identify the real cause << 60 ill << 61 .ascii "DIV0" << 62 << 63 .Lreturn0: << 64 movi a2, 0 << 65 #endif /* XCHAL_HAVE_DIV32 */ << 66 abi_ret_default << 67 3 68 ENDPROC(__udivsi3) !! 4 This file is part of GNU CC. 69 EXPORT_SYMBOL(__udivsi3) !! 5 >> 6 GNU CC is free software; you can redistribute it and/or modify it >> 7 under the terms of the GNU General Public License as published by the >> 8 Free Software Foundation; either version 2, or (at your option) any >> 9 later version. >> 10 >> 11 In addition to the permissions in the GNU General Public License, the >> 12 Free Software Foundation gives you unlimited permission to link the >> 13 compiled version of this file with other programs, and to distribute >> 14 those programs without any restriction coming from the use of this >> 15 file. (The General Public License restrictions do apply in other >> 16 respects; for example, they cover modification of the file, and >> 17 distribution when not linked into another program.) >> 18 >> 19 This file is distributed in the hope that it will be useful, but >> 20 WITHOUT ANY WARRANTY; without even the implied warranty of >> 21 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU >> 22 General Public License for more details. */ >> 23 >> 24 /* As a special exception, if you link this library with files >> 25 compiled with GCC to produce an executable, this does not cause >> 26 the resulting executable to be covered by the GNU General Public License. >> 27 This exception does not however invalidate any other reasons why >> 28 the executable file might be covered by the GNU General Public License. */ >> 29 >> 30 /* Use this one for any 680x0; assumes no floating point hardware. >> 31 The trailing " '" appearing on some lines is for ANSI preprocessors. Yuk. >> 32 Some of this code comes from MINIX, via the folks at ericsson. >> 33 D. V. Henkel-Wallace (gumby@cygnus.com) Fete Bastille, 1992 >> 34 */ >> 35 #include <linux/export.h> >> 36 /* These are predefined by new versions of GNU cpp. */ >> 37 >> 38 #ifndef __USER_LABEL_PREFIX__ >> 39 #define __USER_LABEL_PREFIX__ _ >> 40 #endif >> 41 >> 42 #ifndef __REGISTER_PREFIX__ >> 43 #define __REGISTER_PREFIX__ >> 44 #endif >> 45 >> 46 #ifndef __IMMEDIATE_PREFIX__ >> 47 #define __IMMEDIATE_PREFIX__ # >> 48 #endif >> 49 >> 50 /* ANSI concatenation macros. */ >> 51 >> 52 #define CONCAT1(a, b) CONCAT2(a, b) >> 53 #define CONCAT2(a, b) a ## b >> 54 >> 55 /* Use the right prefix for global labels. */ >> 56 >> 57 #define SYM(x) CONCAT1 (__USER_LABEL_PREFIX__, x) >> 58 >> 59 /* Use the right prefix for registers. */ >> 60 >> 61 #define REG(x) CONCAT1 (__REGISTER_PREFIX__, x) >> 62 >> 63 /* Use the right prefix for immediate values. */ >> 64 >> 65 #define IMM(x) CONCAT1 (__IMMEDIATE_PREFIX__, x) >> 66 >> 67 #define d0 REG (d0) >> 68 #define d1 REG (d1) >> 69 #define d2 REG (d2) >> 70 #define d3 REG (d3) >> 71 #define d4 REG (d4) >> 72 #define d5 REG (d5) >> 73 #define d6 REG (d6) >> 74 #define d7 REG (d7) >> 75 #define a0 REG (a0) >> 76 #define a1 REG (a1) >> 77 #define a2 REG (a2) >> 78 #define a3 REG (a3) >> 79 #define a4 REG (a4) >> 80 #define a5 REG (a5) >> 81 #define a6 REG (a6) >> 82 #define fp REG (fp) >> 83 #define sp REG (sp) >> 84 >> 85 .text >> 86 .proc >> 87 .globl SYM (__udivsi3) >> 88 SYM (__udivsi3): >> 89 #if !(defined(__mcf5200__) || defined(__mcoldfire__)) >> 90 movel d2, sp@- >> 91 movel sp@(12), d1 /* d1 = divisor */ >> 92 movel sp@(8), d0 /* d0 = dividend */ >> 93 >> 94 cmpl IMM (0x10000), d1 /* divisor >= 2 ^ 16 ? */ >> 95 jcc L3 /* then try next algorithm */ >> 96 movel d0, d2 >> 97 clrw d2 >> 98 swap d2 >> 99 divu d1, d2 /* high quotient in lower word */ >> 100 movew d2, d0 /* save high quotient */ >> 101 swap d0 >> 102 movew sp@(10), d2 /* get low dividend + high rest */ >> 103 divu d1, d2 /* low quotient */ >> 104 movew d2, d0 >> 105 jra L6 >> 106 >> 107 L3: movel d1, d2 /* use d2 as divisor backup */ >> 108 L4: lsrl IMM (1), d1 /* shift divisor */ >> 109 lsrl IMM (1), d0 /* shift dividend */ >> 110 cmpl IMM (0x10000), d1 /* still divisor >= 2 ^ 16 ? */ >> 111 jcc L4 >> 112 divu d1, d0 /* now we have 16 bit divisor */ >> 113 andl IMM (0xffff), d0 /* mask out divisor, ignore remainder */ >> 114 >> 115 /* Multiply the 16 bit tentative quotient with the 32 bit divisor. Because of >> 116 the operand ranges, this might give a 33 bit product. If this product is >> 117 greater than the dividend, the tentative quotient was too large. */ >> 118 movel d2, d1 >> 119 mulu d0, d1 /* low part, 32 bits */ >> 120 swap d2 >> 121 mulu d0, d2 /* high part, at most 17 bits */ >> 122 swap d2 /* align high part with low part */ >> 123 tstw d2 /* high part 17 bits? */ >> 124 jne L5 /* if 17 bits, quotient was too large */ >> 125 addl d2, d1 /* add parts */ >> 126 jcs L5 /* if sum is 33 bits, quotient was too large */ >> 127 cmpl sp@(8), d1 /* compare the sum with the dividend */ >> 128 jls L6 /* if sum > dividend, quotient was too large */ >> 129 L5: subql IMM (1), d0 /* adjust quotient */ >> 130 >> 131 L6: movel sp@+, d2 >> 132 rts >> 133 >> 134 #else /* __mcf5200__ || __mcoldfire__ */ >> 135 >> 136 /* Coldfire implementation of non-restoring division algorithm from >> 137 Hennessy & Patterson, Appendix A. */ >> 138 link a6,IMM (-12) >> 139 moveml d2-d4,sp@ >> 140 movel a6@(8),d0 >> 141 movel a6@(12),d1 >> 142 clrl d2 | clear p >> 143 moveq IMM (31),d4 >> 144 L1: addl d0,d0 | shift reg pair (p,a) one bit left >> 145 addxl d2,d2 >> 146 movl d2,d3 | subtract b from p, store in tmp. >> 147 subl d1,d3 >> 148 jcs L2 | if no carry, >> 149 bset IMM (0),d0 | set the low order bit of a to 1, >> 150 movl d3,d2 | and store tmp in p. >> 151 L2: subql IMM (1),d4 >> 152 jcc L1 >> 153 moveml sp@,d2-d4 | restore data registers >> 154 unlk a6 | and return >> 155 rts >> 156 #endif /* __mcf5200__ || __mcoldfire__ */ >> 157 EXPORT_SYMBOL(__udivsi3)
Linux® is a registered trademark of Linus Torvalds in the United States and other countries.
TOMOYO® is a registered trademark of NTT DATA CORPORATION.