123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303 |
- /* libgcc routines for the MCore.
- Copyright (C) 1993-2022 Free Software Foundation, Inc.
- This file is part of GCC.
- GCC 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 3, or (at your option) any
- later version.
- 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.
- Under Section 7 of GPL version 3, you are granted additional
- permissions described in the GCC Runtime Library Exception, version
- 3.1, as published by the Free Software Foundation.
- You should have received a copy of the GNU General Public License and
- a copy of the GCC Runtime Library Exception along with this program;
- see the files COPYING3 and COPYING.RUNTIME respectively. If not, see
- <http://www.gnu.org/licenses/>. */
- #define CONCAT1(a, b) CONCAT2(a, b)
- #define CONCAT2(a, b) a ## b
- /* Use the right prefix for global labels. */
- #define SYM(x) CONCAT1 (__, x)
- #ifdef __ELF__
- #define TYPE(x) .type SYM (x),@function
- #define SIZE(x) .size SYM (x), . - SYM (x)
- #else
- #define TYPE(x)
- #define SIZE(x)
- #endif
- .macro FUNC_START name
- .text
- .globl SYM (\name)
- TYPE (\name)
- SYM (\name):
- .endm
- .macro FUNC_END name
- SIZE (\name)
- .endm
- #ifdef L_udivsi3
- FUNC_START udiv32
- FUNC_START udivsi32
- movi r1,0 // r1-r2 form 64 bit dividend
- movi r4,1 // r4 is quotient (1 for a sentinel)
- cmpnei r3,0 // look for 0 divisor
- bt 9f
- trap 3 // divide by 0
- 9:
- // control iterations; skip across high order 0 bits in dividend
- mov r7,r2
- cmpnei r7,0
- bt 8f
- movi r2,0 // 0 dividend
- jmp r15 // quick return
- 8:
- ff1 r7 // figure distance to skip
- lsl r4,r7 // move the sentinel along (with 0's behind)
- lsl r2,r7 // and the low 32 bits of numerator
- // appears to be wrong...
- // tested out incorrectly in our OS work...
- // mov r7,r3 // looking at divisor
- // ff1 r7 // I can move 32-r7 more bits to left.
- // addi r7,1 // ok, one short of that...
- // mov r1,r2
- // lsr r1,r7 // bits that came from low order...
- // rsubi r7,31 // r7 == "32-n" == LEFT distance
- // addi r7,1 // this is (32-n)
- // lsl r4,r7 // fixes the high 32 (quotient)
- // lsl r2,r7
- // cmpnei r4,0
- // bf 4f // the sentinel went away...
- // run the remaining bits
- 1: lslc r2,1 // 1 bit left shift of r1-r2
- addc r1,r1
- cmphs r1,r3 // upper 32 of dividend >= divisor?
- bf 2f
- sub r1,r3 // if yes, subtract divisor
- 2: addc r4,r4 // shift by 1 and count subtracts
- bf 1b // if sentinel falls out of quotient, stop
- 4: mov r2,r4 // return quotient
- mov r3,r1 // and piggyback the remainder
- jmp r15
- FUNC_END udiv32
- FUNC_END udivsi32
- #endif
- #ifdef L_umodsi3
- FUNC_START urem32
- FUNC_START umodsi3
- movi r1,0 // r1-r2 form 64 bit dividend
- movi r4,1 // r4 is quotient (1 for a sentinel)
- cmpnei r3,0 // look for 0 divisor
- bt 9f
- trap 3 // divide by 0
- 9:
- // control iterations; skip across high order 0 bits in dividend
- mov r7,r2
- cmpnei r7,0
- bt 8f
- movi r2,0 // 0 dividend
- jmp r15 // quick return
- 8:
- ff1 r7 // figure distance to skip
- lsl r4,r7 // move the sentinel along (with 0's behind)
- lsl r2,r7 // and the low 32 bits of numerator
- 1: lslc r2,1 // 1 bit left shift of r1-r2
- addc r1,r1
- cmphs r1,r3 // upper 32 of dividend >= divisor?
- bf 2f
- sub r1,r3 // if yes, subtract divisor
- 2: addc r4,r4 // shift by 1 and count subtracts
- bf 1b // if sentinel falls out of quotient, stop
- mov r2,r1 // return remainder
- jmp r15
- FUNC_END urem32
- FUNC_END umodsi3
- #endif
- #ifdef L_divsi3
- FUNC_START div32
- FUNC_START divsi3
- mov r5,r2 // calc sign of quotient
- xor r5,r3
- abs r2 // do unsigned divide
- abs r3
- movi r1,0 // r1-r2 form 64 bit dividend
- movi r4,1 // r4 is quotient (1 for a sentinel)
- cmpnei r3,0 // look for 0 divisor
- bt 9f
- trap 3 // divide by 0
- 9:
- // control iterations; skip across high order 0 bits in dividend
- mov r7,r2
- cmpnei r7,0
- bt 8f
- movi r2,0 // 0 dividend
- jmp r15 // quick return
- 8:
- ff1 r7 // figure distance to skip
- lsl r4,r7 // move the sentinel along (with 0's behind)
- lsl r2,r7 // and the low 32 bits of numerator
- // tested out incorrectly in our OS work...
- // mov r7,r3 // looking at divisor
- // ff1 r7 // I can move 32-r7 more bits to left.
- // addi r7,1 // ok, one short of that...
- // mov r1,r2
- // lsr r1,r7 // bits that came from low order...
- // rsubi r7,31 // r7 == "32-n" == LEFT distance
- // addi r7,1 // this is (32-n)
- // lsl r4,r7 // fixes the high 32 (quotient)
- // lsl r2,r7
- // cmpnei r4,0
- // bf 4f // the sentinel went away...
- // run the remaining bits
- 1: lslc r2,1 // 1 bit left shift of r1-r2
- addc r1,r1
- cmphs r1,r3 // upper 32 of dividend >= divisor?
- bf 2f
- sub r1,r3 // if yes, subtract divisor
- 2: addc r4,r4 // shift by 1 and count subtracts
- bf 1b // if sentinel falls out of quotient, stop
- 4: mov r2,r4 // return quotient
- mov r3,r1 // piggyback the remainder
- btsti r5,31 // after adjusting for sign
- bf 3f
- rsubi r2,0
- rsubi r3,0
- 3: jmp r15
- FUNC_END div32
- FUNC_END divsi3
- #endif
- #ifdef L_modsi3
- FUNC_START rem32
- FUNC_START modsi3
- mov r5,r2 // calc sign of remainder
- abs r2 // do unsigned divide
- abs r3
- movi r1,0 // r1-r2 form 64 bit dividend
- movi r4,1 // r4 is quotient (1 for a sentinel)
- cmpnei r3,0 // look for 0 divisor
- bt 9f
- trap 3 // divide by 0
- 9:
- // control iterations; skip across high order 0 bits in dividend
- mov r7,r2
- cmpnei r7,0
- bt 8f
- movi r2,0 // 0 dividend
- jmp r15 // quick return
- 8:
- ff1 r7 // figure distance to skip
- lsl r4,r7 // move the sentinel along (with 0's behind)
- lsl r2,r7 // and the low 32 bits of numerator
- 1: lslc r2,1 // 1 bit left shift of r1-r2
- addc r1,r1
- cmphs r1,r3 // upper 32 of dividend >= divisor?
- bf 2f
- sub r1,r3 // if yes, subtract divisor
- 2: addc r4,r4 // shift by 1 and count subtracts
- bf 1b // if sentinel falls out of quotient, stop
- mov r2,r1 // return remainder
- btsti r5,31 // after adjusting for sign
- bf 3f
- rsubi r2,0
- 3: jmp r15
- FUNC_END rem32
- FUNC_END modsi3
- #endif
- /* GCC expects that {__eq,__ne,__gt,__ge,__le,__lt}{df2,sf2}
- will behave as __cmpdf2. So, we stub the implementations to
- jump on to __cmpdf2 and __cmpsf2.
-
- All of these shortcircuit the return path so that __cmp{sd}f2
- will go directly back to the caller. */
- .macro COMPARE_DF_JUMP name
- .import SYM (cmpdf2)
- FUNC_START \name
- jmpi SYM (cmpdf2)
- FUNC_END \name
- .endm
-
- #ifdef L_eqdf2
- COMPARE_DF_JUMP eqdf2
- #endif /* L_eqdf2 */
- #ifdef L_nedf2
- COMPARE_DF_JUMP nedf2
- #endif /* L_nedf2 */
- #ifdef L_gtdf2
- COMPARE_DF_JUMP gtdf2
- #endif /* L_gtdf2 */
- #ifdef L_gedf2
- COMPARE_DF_JUMP gedf2
- #endif /* L_gedf2 */
- #ifdef L_ltdf2
- COMPARE_DF_JUMP ltdf2
- #endif /* L_ltdf2 */
-
- #ifdef L_ledf2
- COMPARE_DF_JUMP ledf2
- #endif /* L_ledf2 */
- /* SINGLE PRECISION FLOATING POINT STUBS */
- .macro COMPARE_SF_JUMP name
- .import SYM (cmpsf2)
- FUNC_START \name
- jmpi SYM (cmpsf2)
- FUNC_END \name
- .endm
-
- #ifdef L_eqsf2
- COMPARE_SF_JUMP eqsf2
- #endif /* L_eqsf2 */
-
- #ifdef L_nesf2
- COMPARE_SF_JUMP nesf2
- #endif /* L_nesf2 */
-
- #ifdef L_gtsf2
- COMPARE_SF_JUMP gtsf2
- #endif /* L_gtsf2 */
-
- #ifdef L_gesf2
- COMPARE_SF_JUMP __gesf2
- #endif /* L_gesf2 */
-
- #ifdef L_ltsf2
- COMPARE_SF_JUMP __ltsf2
- #endif /* L_ltsf2 */
-
- #ifdef L_lesf2
- COMPARE_SF_JUMP lesf2
- #endif /* L_lesf2 */
|