cregit-Linux how code gets into the kernel

Release 4.16 lib/mpi/longlong.h

Directory: lib/mpi
/* longlong.h -- definitions for mixed size 32/64 bit arithmetic.
 * Note: I added some stuff for use with gnupg
 *
 * Copyright (C) 1991, 1992, 1993, 1994, 1996, 1998,
 *      2000, 2001, 2002, 2003 Free Software Foundation, Inc.
 *
 * This file is free software; you can redistribute it and/or modify
 * it under the terms of the GNU Library General Public License as published by
 * the Free Software Foundation; either version 2 of the License, 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 Library General Public
 * License for more details.
 *
 * You should have received a copy of the GNU Library General Public License
 * along with this file; see the file COPYING.LIB.  If not, write to
 * the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
 * MA 02111-1307, USA. */

#include <linux/count_zeros.h>

/* You have to define the following before including this file:
 *
 * UWtype -- An unsigned type, default type for operations (typically a "word")
 * UHWtype -- An unsigned type, at least half the size of UWtype.
 * UDWtype -- An unsigned type, at least twice as large a UWtype
 * W_TYPE_SIZE -- size in bits of UWtype
 *
 * SItype, USItype -- Signed and unsigned 32 bit types.
 * DItype, UDItype -- Signed and unsigned 64 bit types.
 *
 * On a 32 bit machine UWtype should typically be USItype;
 * on a 64 bit machine, UWtype should typically be UDItype.
*/


#define __BITS4 (W_TYPE_SIZE / 4)

#define __ll_B ((UWtype) 1 << (W_TYPE_SIZE / 2))

#define __ll_lowpart(t) ((UWtype) (t) & (__ll_B - 1))

#define __ll_highpart(t) ((UWtype) (t) >> (W_TYPE_SIZE / 2))

/* This is used to make sure no undesirable sharing between different libraries
        that use this file takes place.  */
#ifndef __MPN

#define __MPN(x) __##x
#endif

/* Define auxiliary asm macros.
 *
 * 1) umul_ppmm(high_prod, low_prod, multipler, multiplicand) multiplies two
 * UWtype integers MULTIPLER and MULTIPLICAND, and generates a two UWtype
 * word product in HIGH_PROD and LOW_PROD.
 *
 * 2) __umulsidi3(a,b) multiplies two UWtype integers A and B, and returns a
 * UDWtype product.  This is just a variant of umul_ppmm.

 * 3) udiv_qrnnd(quotient, remainder, high_numerator, low_numerator,
 * denominator) divides a UDWtype, composed by the UWtype integers
 * HIGH_NUMERATOR and LOW_NUMERATOR, by DENOMINATOR and places the quotient
 * in QUOTIENT and the remainder in REMAINDER.  HIGH_NUMERATOR must be less
 * than DENOMINATOR for correct operation.  If, in addition, the most
 * significant bit of DENOMINATOR must be 1, then the pre-processor symbol
 * UDIV_NEEDS_NORMALIZATION is defined to 1.
 * 4) sdiv_qrnnd(quotient, remainder, high_numerator, low_numerator,
 * denominator).  Like udiv_qrnnd but the numbers are signed.  The quotient
 * is rounded towards 0.
 *
 * 5) count_leading_zeros(count, x) counts the number of zero-bits from the
 * msb to the first non-zero bit in the UWtype X.  This is the number of
 * steps X needs to be shifted left to set the msb.  Undefined for X == 0,
 * unless the symbol COUNT_LEADING_ZEROS_0 is defined to some value.
 *
 * 6) count_trailing_zeros(count, x) like count_leading_zeros, but counts
 * from the least significant end.
 *
 * 7) add_ssaaaa(high_sum, low_sum, high_addend_1, low_addend_1,
 * high_addend_2, low_addend_2) adds two UWtype integers, composed by
 * HIGH_ADDEND_1 and LOW_ADDEND_1, and HIGH_ADDEND_2 and LOW_ADDEND_2
 * respectively.  The result is placed in HIGH_SUM and LOW_SUM.  Overflow
 * (i.e. carry out) is not stored anywhere, and is lost.
 *
 * 8) sub_ddmmss(high_difference, low_difference, high_minuend, low_minuend,
 * high_subtrahend, low_subtrahend) subtracts two two-word UWtype integers,
 * composed by HIGH_MINUEND_1 and LOW_MINUEND_1, and HIGH_SUBTRAHEND_2 and
 * LOW_SUBTRAHEND_2 respectively.  The result is placed in HIGH_DIFFERENCE
 * and LOW_DIFFERENCE.  Overflow (i.e. carry out) is not stored anywhere,
 * and is lost.
 *
 * If any of these macros are left undefined for a particular CPU,
 * C macros are used.  */

/* The CPUs come in alphabetical order below.
 *
 * Please add support for more CPUs here, or improve the current support
 * for the CPUs below!  */

#if defined(__GNUC__) && !defined(NO_ASM)

/* We sometimes need to clobber "cc" with gcc2, but that would not be
        understood by gcc1.     Use cpp to avoid major code duplication.  */
#if __GNUC__ < 2

#define __CLOBBER_CC

#define __AND_CLOBBER_CC
#else /* __GNUC__ >= 2 */

#define __CLOBBER_CC : "cc"

#define __AND_CLOBBER_CC , "cc"
#endif /* __GNUC__ < 2 */

/***************************************
        **************  A29K  *****************
        ***************************************/
#if (defined(__a29k__) || defined(_AM29K)) && W_TYPE_SIZE == 32

#define add_ssaaaa(sh, sl, ah, al, bh, bl) \
	__asm__ ("add %1,%4,%5\n" \
                "addc %0,%2,%3" \
        : "=r" ((USItype)(sh)), \
                "=&r" ((USItype)(sl)) \
        : "%r" ((USItype)(ah)), \
                "rI" ((USItype)(bh)), \
                "%r" ((USItype)(al)), \
                "rI" ((USItype)(bl)))

#define sub_ddmmss(sh, sl, ah, al, bh, bl) \
	__asm__ ("sub %1,%4,%5\n" \
                "subc %0,%2,%3" \
        : "=r" ((USItype)(sh)), \
                "=&r" ((USItype)(sl)) \
        : "r" ((USItype)(ah)), \
                "rI" ((USItype)(bh)), \
                "r" ((USItype)(al)), \
                "rI" ((USItype)(bl)))

#define umul_ppmm(xh, xl, m0, m1) \
do { \
                USItype __m0 = (m0), __m1 = (m1); \
                __asm__ ("multiplu %0,%1,%2" \
                : "=r" ((USItype)(xl)) \
                : "r" (__m0), \
                        "r" (__m1)); \
                __asm__ ("multmu %0,%1,%2" \
                : "=r" ((USItype)(xh)) \
                : "r" (__m0), \
                        "r" (__m1)); \
} while (0)

#define udiv_qrnnd(q, r, n1, n0, d) \
	__asm__ ("dividu %0,%3,%4" \
        : "=r" ((USItype)(q)), \
                "=q" ((USItype)(r)) \
        : "1" ((USItype)(n1)), \
                "r" ((USItype)(n0)), \
                "r" ((USItype)(d)))
#endif /* __a29k__ */

#if defined(__alpha) && W_TYPE_SIZE == 64

#define umul_ppmm(ph, pl, m0, m1)			\
do {                                                    \
        UDItype __m0 = (m0), __m1 = (m1);               \
        (ph) = __builtin_alpha_umulh(__m0, __m1);       \
        (pl) = __m0 * __m1;                             \
} while (0)

#define UMUL_TIME 46
#ifndef LONGLONG_STANDALONE

#define udiv_qrnnd(q, r, n1, n0, d) \
do { UDItype __r; \
        (q) = __udiv_qrnnd(&__r, (n1), (n0), (d)); \
        (r) = __r; \
} while (0)
extern UDItype __udiv_qrnnd(UDItype *, UDItype, UDItype, UDItype);

#define UDIV_TIME 220
#endif /* LONGLONG_STANDALONE */
#endif /* __alpha */

/***************************************
        **************  ARM  ******************
        ***************************************/
#if defined(__arm__) && W_TYPE_SIZE == 32

#define add_ssaaaa(sh, sl, ah, al, bh, bl) \
	__asm__ ("adds %1, %4, %5\n" \
                "adc  %0, %2, %3" \
        : "=r" (sh), \
                "=&r" (sl) \
        : "%r" ((USItype)(ah)), \
                "rI" ((USItype)(bh)), \
                "%r" ((USItype)(al)), \
                "rI" ((USItype)(bl)))

#define sub_ddmmss(sh, sl, ah, al, bh, bl) \
	__asm__ ("subs %1, %4, %5\n" \
                "sbc  %0, %2, %3" \
        : "=r" (sh), \
                "=&r" (sl) \
        : "r" ((USItype)(ah)), \
                "rI" ((USItype)(bh)), \
                "r" ((USItype)(al)), \
                "rI" ((USItype)(bl)))
#if defined __ARM_ARCH_2__ || defined __ARM_ARCH_3__

#define umul_ppmm(xh, xl, a, b) \
	__asm__ ("@ Inlined umul_ppmm\n" \
                "mov    %|r0, %2, lsr #16               @ AAAA\n" \
                "mov    %|r2, %3, lsr #16               @ BBBB\n" \
                "bic    %|r1, %2, %|r0, lsl #16         @ aaaa\n" \
                "bic    %0, %3, %|r2, lsl #16           @ bbbb\n" \
                "mul    %1, %|r1, %|r2                  @ aaaa * BBBB\n" \
                "mul    %|r2, %|r0, %|r2                @ AAAA * BBBB\n" \
                "mul    %|r1, %0, %|r1                  @ aaaa * bbbb\n" \
                "mul    %0, %|r0, %0                    @ AAAA * bbbb\n" \
                "adds   %|r0, %1, %0                    @ central sum\n" \
                "addcs  %|r2, %|r2, #65536\n" \
                "adds   %1, %|r1, %|r0, lsl #16\n" \
                "adc    %0, %|r2, %|r0, lsr #16" \
        : "=&r" (xh), \
                "=r" (xl) \
        : "r" ((USItype)(a)), \
                "r" ((USItype)(b)) \
        : "r0", "r1", "r2")
#else

#define umul_ppmm(xh, xl, a, b) \
	__asm__ ("@ Inlined umul_ppmm\n" \
                "umull %1, %0, %2, %3" \
        : "=&r" (xh), \
                "=&r" (xl) \
        : "r" ((USItype)(a)), \
                "r" ((USItype)(b)) \
        : "r0", "r1")
#endif

#define UMUL_TIME 20

#define UDIV_TIME 100
#endif /* __arm__ */

/***************************************
        **************  CLIPPER  **************
        ***************************************/
#if defined(__clipper__) && W_TYPE_SIZE == 32

#define umul_ppmm(w1, w0, u, v) \
	({union {UDItype __ll; \
                struct {USItype __l, __h; } __i; \
        } __xx; \
        __asm__ ("mulwux %2,%0" \
        : "=r" (__xx.__ll) \
        : "%0" ((USItype)(u)), \
                "r" ((USItype)(v))); \
        (w1) = __xx.__i.__h; (w0) = __xx.__i.__l; })

#define smul_ppmm(w1, w0, u, v) \
	({union {DItype __ll; \
                struct {SItype __l, __h; } __i; \
        } __xx; \
        __asm__ ("mulwx %2,%0" \
        : "=r" (__xx.__ll) \
        : "%0" ((SItype)(u)), \
                "r" ((SItype)(v))); \
        (w1) = __xx.__i.__h; (w0) = __xx.__i.__l; })

#define __umulsidi3(u, v) \
	({UDItype __w; \
        __asm__ ("mulwux %2,%0" \
        : "=r" (__w) \
        : "%0" ((USItype)(u)), \
                "r" ((USItype)(v))); \
        __w; })
#endif /* __clipper__ */

/***************************************
        **************  GMICRO  ***************
        ***************************************/
#if defined(__gmicro__) && W_TYPE_SIZE == 32

#define add_ssaaaa(sh, sl, ah, al, bh, bl) \
	__asm__ ("add.w %5,%1\n" \
                "addx %3,%0" \
        : "=g" ((USItype)(sh)), \
                "=&g" ((USItype)(sl)) \
        : "%0" ((USItype)(ah)), \
                "g" ((USItype)(bh)), \
                "%1" ((USItype)(al)), \
                "g" ((USItype)(bl)))

#define sub_ddmmss(sh, sl, ah, al, bh, bl) \
	__asm__ ("sub.w %5,%1\n" \
                "subx %3,%0" \
        : "=g" ((USItype)(sh)), \
                "=&g" ((USItype)(sl)) \
        : "0" ((USItype)(ah)), \
                "g" ((USItype)(bh)), \
                "1" ((USItype)(al)), \
                "g" ((USItype)(bl)))

#define umul_ppmm(ph, pl, m0, m1) \
	__asm__ ("mulx %3,%0,%1" \
        : "=g" ((USItype)(ph)), \
                "=r" ((USItype)(pl)) \
        : "%0" ((USItype)(m0)), \
                "g" ((USItype)(m1)))

#define udiv_qrnnd(q, r, nh, nl, d) \
	__asm__ ("divx %4,%0,%1" \
        : "=g" ((USItype)(q)), \
                "=r" ((USItype)(r)) \
        : "1" ((USItype)(nh)), \
                "0" ((USItype)(nl)), \
                "g" ((USItype)(d)))
#endif

/***************************************
        **************  HPPA  *****************
        ***************************************/
#if defined(__hppa) && W_TYPE_SIZE == 32

#define add_ssaaaa(sh, sl, ah, al, bh, bl) \
	__asm__ ("add %4,%5,%1\n" \
                   "addc %2,%3,%0" \
        : "=r" ((USItype)(sh)), \
             "=&r" ((USItype)(sl)) \
        : "%rM" ((USItype)(ah)), \
             "rM" ((USItype)(bh)), \
             "%rM" ((USItype)(al)), \
             "rM" ((USItype)(bl)))

#define sub_ddmmss(sh, sl, ah, al, bh, bl) \
	__asm__ ("sub %4,%5,%1\n" \
           "subb %2,%3,%0" \
        : "=r" ((USItype)(sh)), \
             "=&r" ((USItype)(sl)) \
        : "rM" ((USItype)(ah)), \
             "rM" ((USItype)(bh)), \
             "rM" ((USItype)(al)), \
             "rM" ((USItype)(bl)))
#if 0 && defined(_PA_RISC1_1)
/* xmpyu uses floating point register which is not allowed in Linux kernel. */
#define umul_ppmm(wh, wl, u, v) \
do { \
        union {UDItype __ll; \
        struct {USItype __h, __l; } __i; \
        } __xx; \
        __asm__ ("xmpyu %1,%2,%0" \
        : "=*f" (__xx.__ll) \
        : "*f" ((USItype)(u)), \
               "*f" ((USItype)(v))); \
        (wh) = __xx.__i.__h; \
        (wl) = __xx.__i.__l; \
} while (0)
#define UMUL_TIME 8
#define UDIV_TIME 60
#else

#define UMUL_TIME 40

#define UDIV_TIME 80
#endif
#if 0 /* #ifndef LONGLONG_STANDALONE */
#define udiv_qrnnd(q, r, n1, n0, d) \
do { USItype __r; \
        (q) = __udiv_qrnnd(&__r, (n1), (n0), (d)); \
        (r) = __r; \
} while (0)
extern USItype __udiv_qrnnd();
#endif /* LONGLONG_STANDALONE */
#endif /* hppa */

/***************************************
        **************  I370  *****************
        ***************************************/
#if (defined(__i370__) || defined(__mvs__)) && W_TYPE_SIZE == 32

#define umul_ppmm(xh, xl, m0, m1) \
do { \
        union {UDItype __ll; \
           struct {USItype __h, __l; } __i; \
        } __xx; \
        USItype __m0 = (m0), __m1 = (m1); \
        __asm__ ("mr %0,%3" \
        : "=r" (__xx.__i.__h), \
               "=r" (__xx.__i.__l) \
        : "%1" (__m0), \
               "r" (__m1)); \
        (xh) = __xx.__i.__h; (xl) = __xx.__i.__l; \
        (xh) += ((((SItype) __m0 >> 31) & __m1) \
             + (((SItype) __m1 >> 31) & __m0)); \
} while (0)

#define smul_ppmm(xh, xl, m0, m1) \
do { \
        union {DItype __ll; \
           struct {USItype __h, __l; } __i; \
        } __xx; \
        __asm__ ("mr %0,%3" \
        : "=r" (__xx.__i.__h), \
               "=r" (__xx.__i.__l) \
        : "%1" (m0), \
               "r" (m1)); \
        (xh) = __xx.__i.__h; (xl) = __xx.__i.__l; \
} while (0)

#define sdiv_qrnnd(q, r, n1, n0, d) \
do { \
        union {DItype __ll; \
           struct {USItype __h, __l; } __i; \
        } __xx; \
        __xx.__i.__h = n1; __xx.__i.__l = n0; \
        __asm__ ("dr %0,%2" \
        : "=r" (__xx.__ll) \
        : "0" (__xx.__ll), "r" (d)); \
        (q) = __xx.__i.__l; (r) = __xx.__i.__h; \
} while (0)
#endif

/***************************************
        **************  I386  *****************
        ***************************************/
#undef __i386__
#if (defined(__i386__) || defined(__i486__)) && W_TYPE_SIZE == 32

#define add_ssaaaa(sh, sl, ah, al, bh, bl) \
	__asm__ ("addl %5,%1\n" \
           "adcl %3,%0" \
        : "=r" ((USItype)(sh)), \
             "=&r" ((USItype)(sl)) \
        : "%0" ((USItype)(ah)), \
             "g" ((USItype)(bh)), \
             "%1" ((USItype)(al)), \
             "g" ((USItype)(bl)))

#define sub_ddmmss(sh, sl, ah, al, bh, bl) \
	__asm__ ("subl %5,%1\n" \
           "sbbl %3,%0" \
        : "=r" ((USItype)(sh)), \
             "=&r" ((USItype)(sl)) \
        : "0" ((USItype)(ah)), \
             "g" ((USItype)(bh)), \
             "1" ((USItype)(al)), \
             "g" ((USItype)(bl)))

#define umul_ppmm(w1, w0, u, v) \
	__asm__ ("mull %3" \
        : "=a" ((USItype)(w0)), \
             "=d" ((USItype)(w1)) \
        : "%0" ((USItype)(u)), \
             "rm" ((USItype)(v)))

#define udiv_qrnnd(q, r, n1, n0, d) \
	__asm__ ("divl %4" \
        : "=a" ((USItype)(q)), \
             "=d" ((USItype)(r)) \
        : "0" ((USItype)(n0)), \
             "1" ((USItype)(n1)), \
             "rm" ((USItype)(d)))
#ifndef UMUL_TIME

#define UMUL_TIME 40
#endif
#ifndef UDIV_TIME

#define UDIV_TIME 40
#endif
#endif /* 80x86 */

/***************************************
        **************  I860  *****************
        ***************************************/
#if defined(__i860__) && W_TYPE_SIZE == 32

#define rshift_rhlc(r, h, l, c) \
	__asm__ ("shr %3,r0,r0\n" \
        "shrd %1,%2,%0" \
           "=r" (r) : "r" (h), "r" (l), "rn" (c))
#endif /* i860 */

/***************************************
        **************  I960  *****************
        ***************************************/
#if defined(__i960__) && W_TYPE_SIZE == 32

#define add_ssaaaa(sh, sl, ah, al, bh, bl) \
	__asm__ ("cmpo 1,0\n" \
        "addc %5,%4,%1\n" \
        "addc %3,%2,%0" \
        : "=r" ((USItype)(sh)), \
             "=&r" ((USItype)(sl)) \
        : "%dI" ((USItype)(ah)), \
             "dI" ((USItype)(bh)), \
             "%dI" ((USItype)(al)), \
             "dI" ((USItype)(bl)))

#define sub_ddmmss(sh, sl, ah, al, bh, bl) \
	__asm__ ("cmpo 0,0\n" \
        "subc %5,%4,%1\n" \
        "subc %3,%2,%0" \
        : "=r" ((USItype)(sh)), \
             "=&r" ((USItype)(sl)) \
        : "dI" ((USItype)(ah)), \
             "dI" ((USItype)(bh)), \
             "dI" ((USItype)(al)), \
             "dI" ((USItype)(bl)))

#define umul_ppmm(w1, w0, u, v) \
	({union {UDItype __ll; \
           struct {USItype __l, __h; } __i; \
        } __xx; \
        __asm__ ("emul        %2,%1,%0" \
        : "=d" (__xx.__ll) \
        : "%dI" ((USItype)(u)), \
             "dI" ((USItype)(v))); \
        (w1) = __xx.__i.__h; (w0) = __xx.__i.__l; })

#define __umulsidi3(u, v) \
	({UDItype __w; \
        __asm__ ("emul      %2,%1,%0" \
        : "=d" (__w) \
        : "%dI" ((USItype)(u)), \
               "dI" ((USItype)(v))); \
        __w; })

#define udiv_qrnnd(q, r, nh, nl, d) \
do { \
        union {UDItype __ll; \
           struct {USItype __l, __h; } __i; \
        } __nn; \
        __nn.__i.__h = (nh); __nn.__i.__l = (nl); \
        __asm__ ("ediv %d,%n,%0" \
        : "=d" (__rq.__ll) \
        : "dI" (__nn.__ll), \
             "dI" ((USItype)(d))); \
        (r) = __rq.__i.__l; (q) = __rq.__i.__h; \
} while (0)
#if defined(__i960mx)		/* what is the proper symbol to test??? */

#define rshift_rhlc(r, h, l, c) \
do { \
        union {UDItype __ll; \
           struct {USItype __l, __h; } __i; \
        } __nn; \
        __nn.__i.__h = (h); __nn.__i.__l = (l); \
        __asm__ ("shre %2,%1,%0" \
        : "=d" (r) : "dI" (__nn.__ll), "dI" (c)); \
}
#endif /* i960mx */
#endif /* i960 */

/***************************************
        **************  68000   ****************
        ***************************************/
#if (defined(__mc68000__) || defined(__mc68020__) || defined(__NeXT__) || defined(mc68020)) && W_TYPE_SIZE == 32

#define add_ssaaaa(sh, sl, ah, al, bh, bl) \
	__asm__ ("add%.l %5,%1\n" \
           "addx%.l %3,%0" \
        : "=d" ((USItype)(sh)), \
             "=&d" ((USItype)(sl)) \
        : "%0" ((USItype)(ah)), \
             "d" ((USItype)(bh)), \
             "%1" ((USItype)(al)), \
             "g" ((USItype)(bl)))

#define sub_ddmmss(sh, sl, ah, al, bh, bl) \
	__asm__ ("sub%.l %5,%1\n" \
           "subx%.l %3,%0" \
        : "=d" ((USItype)(sh)), \
             "=&d" ((USItype)(sl)) \
        : "0" ((USItype)(ah)), \
             "d" ((USItype)(bh)), \
             "1" ((USItype)(al)), \
             "g" ((USItype)(bl)))
#if (defined(__mc68020__) || defined(__NeXT__) || defined(mc68020))

#define umul_ppmm(w1, w0, u, v) \
	__asm__ ("mulu%.l %3,%1:%0" \
        : "=d" ((USItype)(w0)), \
             "=d" ((USItype)(w1)) \
        : "%0" ((USItype)(u)), \
             "dmi" ((USItype)(v)))

#define UMUL_TIME 45

#define udiv_qrnnd(q, r, n1, n0, d) \
	__asm__ ("divu%.l %4,%1:%0" \
        : "=d" ((USItype)(q)), \
             "=d" ((USItype)(r)) \
        : "0" ((USItype)(n0)), \
             "1" ((USItype)(n1)), \
             "dmi" ((USItype)(d)))

#define UDIV_TIME 90

#define sdiv_qrnnd(q, r, n1, n0, d) \
	__asm__ ("divs%.l %4,%1:%0" \
        : "=d" ((USItype)(q)), \
             "=d" ((USItype)(r)) \
        : "0" ((USItype)(n0)), \
             "1" ((USItype)(n1)), \
             "dmi" ((USItype)(d)))
#else /* not mc68020 */

#define umul_ppmm(xh, xl, a, b) \
do { USItype __umul_tmp1, __umul_tmp2; \
        __asm__ ("| Inlined umul_ppmm\n" \
        "move%.l %5,%3\n" \
        "move%.l %2,%0\n" \
        "move%.w %3,%1\n" \
        "swap   %3\n" \
        "swap   %0\n" \
        "mulu   %2,%1\n" \
        "mulu   %3,%0\n" \
        "mulu   %2,%3\n" \
        "swap   %2\n" \
        "mulu   %5,%2\n" \
        "add%.l %3,%2\n" \
        "jcc    1f\n" \
        "add%.l %#0x10000,%0\n" \
        "1:     move%.l %2,%3\n" \
        "clr%.w %2\n" \
        "swap   %2\n" \
        "swap   %3\n" \
        "clr%.w %3\n" \
        "add%.l %3,%1\n" \
        "addx%.l %2,%0\n" \
        "| End inlined umul_ppmm" \
        : "=&d" ((USItype)(xh)), "=&d" ((USItype)(xl)), \
                "=d" (__umul_tmp1), "=&d" (__umul_tmp2) \
        : "%2" ((USItype)(a)), "d" ((USItype)(b))); \
} while (0)

#define UMUL_TIME 100

#define UDIV_TIME 400
#endif /* not mc68020 */
#endif /* mc68000 */

/***************************************
        **************  88000   ****************
        ***************************************/
#if defined(__m88000__) && W_TYPE_SIZE == 32

#define add_ssaaaa(sh, sl, ah, al, bh, bl) \
	__asm__ ("addu.co %1,%r4,%r5\n" \
           "addu.ci %0,%r2,%r3" \
        : "=r" ((USItype)(sh)), \
             "=&r" ((USItype)(sl)) \
        : "%rJ" ((USItype)(ah)), \
             "rJ" ((USItype)(bh)), \
             "%rJ" ((USItype)(al)), \
             "rJ" ((USItype)(bl)))

#define sub_ddmmss(sh, sl, ah, al, bh, bl) \
	__asm__ ("subu.co %1,%r4,%r5\n" \
           "subu.ci %0,%r2,%r3" \
        : "=r" ((USItype)(sh)), \
             "=&r" ((USItype)(sl)) \
        : "rJ" ((USItype)(ah)), \
             "rJ" ((USItype)(bh)), \
             "rJ" ((USItype)(al)), \
             "rJ" ((USItype)(bl)))
#if defined(__m88110__)

#define umul_ppmm(wh, wl, u, v) \
do { \
        union {UDItype __ll; \
           struct {USItype __h, __l; } __i; \
        } __x; \
        __asm__ ("mulu.d %0,%1,%2" : "=r" (__x.__ll) : "r" (u), "r" (v)); \
        (wh) = __x.__i.__h; \
        (wl) = __x.__i.__l; \
} while (0)

#define udiv_qrnnd(q, r, n1, n0, d) \
	({union {UDItype __ll; \
           struct {USItype __h, __l; } __i; \
        } __x, __q; \
        __x.__i.__h = (n1); __x.__i.__l = (n0); \
        __asm__ ("divu.d %0,%1,%2" \
        : "=r" (__q.__ll) : "r" (__x.__ll), "r" (d)); \
        (r) = (n0) - __q.__l * (d); (q) = __q.__l; })

#define UMUL_TIME 5

#define UDIV_TIME 25
#else

#define UMUL_TIME 17

#define UDIV_TIME 150
#endif /* __m88110__ */
#endif /* __m88000__ */

/***************************************
        **************  MIPS  *****************
        ***************************************/
#if defined(__mips__) && W_TYPE_SIZE == 32
#if (__GNUC__ >= 5) || (__GNUC__ >= 4 && __GNUC_MINOR__ >= 4)

#define umul_ppmm(w1, w0, u, v)			\
do {                                            \
        UDItype __ll = (UDItype)(u) * (v);      \
        w1 = __ll >> 32;                        \
        w0 = __ll;                              \
} while (0)
#elif __GNUC__ > 2 || __GNUC_MINOR__ >= 7

#define umul_ppmm(w1, w0, u, v) \
	__asm__ ("multu %2,%3" \
        : "=l" ((USItype)(w0)), \
             "=h" ((USItype)(w1)) \
        : "d" ((USItype)(u)), \
             "d" ((USItype)(v)))
#else

#define umul_ppmm(w1, w0, u, v) \
	__asm__ ("multu %2,%3\n" \
           "mflo %0\n" \
           "mfhi %1" \
        : "=d" ((USItype)(w0)), \
             "=d" ((USItype)(w1)) \
        : "d" ((USItype)(u)), \
             "d" ((USItype)(v)))
#endif

#define UMUL_TIME 10

#define UDIV_TIME 100
#endif /* __mips__ */

/***************************************
        **************  MIPS/64  **************
        ***************************************/
#if (defined(__mips) && __mips >= 3) && W_TYPE_SIZE == 64
#if defined(__mips_isa_rev) && __mips_isa_rev >= 6
/*
 * GCC ends up emitting a __multi3 intrinsic call for MIPS64r6 with the plain C
 * code below, so we special case MIPS64r6 until the compiler can do better.
 */

#define umul_ppmm(w1, w0, u, v)						\
do {                                                                    \
        __asm__ ("dmulu %0,%1,%2"                                       \
                 : "=d" ((UDItype)(w0))                                 \
                 : "d" ((UDItype)(u)),                                  \
                   "d" ((UDItype)(v)));                                 \
        __asm__ ("dmuhu %0,%1,%2"                                       \
                 : "=d" ((UDItype)(w1))                                 \
                 : "d" ((UDItype)(u)),                                  \
                   "d" ((UDItype)(v)));                                 \
} while (0)
#elif (__GNUC__ >= 5) || (__GNUC__ >= 4 && __GNUC_MINOR__ >= 4)

#define umul_ppmm(w1, w0, u, v) \
do {                                                                    \
        typedef unsigned int __ll_UTItype __attribute__((mode(TI)));    \
        __ll_UTItype __ll = (__ll_UTItype)(u) * (v);                    \
        w1 = __ll >> 64;                                                \
        w0 = __ll;                                                      \
} while (0)
#elif __GNUC__ > 2 || __GNUC_MINOR__ >= 7

#define umul_ppmm(w1, w0, u, v) \
	__asm__ ("dmultu %2,%3" \
        : "=l" ((UDItype)(w0)), \
             "=h" ((UDItype)(w1)) \
        : "d" ((UDItype)(u)), \
             "d" ((UDItype)(v)))
#else

#define umul_ppmm(w1, w0, u, v) \
	__asm__ ("dmultu %2,%3\n" \
           "mflo %0\n" \
           "mfhi %1" \
        : "=d" ((UDItype)(w0)), \
             "=d" ((UDItype)(w1)) \
        : "d" ((UDItype)(u)), \
             "d" ((UDItype)(v)))
#endif

#define UMUL_TIME 20

#define UDIV_TIME 140
#endif /* __mips__ */

/***************************************
        **************  32000   ****************
        ***************************************/
#if defined(__ns32000__) && W_TYPE_SIZE == 32

#define umul_ppmm(w1, w0, u, v) \
	({union {UDItype __ll; \
           struct {USItype __l, __h; } __i; \
        } __xx; \
        __asm__ ("meid %2,%0" \
        : "=g" (__xx.__ll) \
        : "%0" ((USItype)(u)), \
             "g" ((USItype)(v))); \
        (w1) = __xx.__i.__h; (w0) = __xx.__i.__l; })

#define __umulsidi3(u, v) \
	({UDItype __w; \
        __asm__ ("meid %2,%0" \
        : "=g" (__w) \
        : "%0" ((USItype)(u)), \
               "g" ((USItype)(v))); \
        __w; })

#define udiv_qrnnd(q, r, n1, n0, d) \
	({union {UDItype __ll; \
           struct {USItype __l, __h; } __i; \
        } __xx; \
        __xx.__i.__h = (n1); __xx.__i.__l = (n0); \
        __asm__ ("deid %2,%0" \
        : "=g" (__xx.__ll) \
        : "0" (__xx.__ll), \
             "g" ((USItype)(d))); \
        (r) = __xx.__i.__l; (q) = __xx.__i.__h; })
#endif /* __ns32000__ */

/***************************************
        **************  PPC  ******************
        ***************************************/
#if (defined(_ARCH_PPC) || defined(_IBMR2)) && W_TYPE_SIZE == 32

#define add_ssaaaa(sh, sl, ah, al, bh, bl) \
do { \
        if (__builtin_constant_p(bh) && (bh) == 0) \
                __asm__ ("{a%I4|add%I4c} %1,%3,%4\n\t{aze|addze} %0,%2" \
                : "=r" ((USItype)(sh)), \
                "=&r" ((USItype)(sl)) \
                : "%r" ((USItype)(ah)), \
                "%r" ((USItype)(al)), \
                "rI" ((USItype)(bl))); \
        else if (__builtin_constant_p(bh) && (bh) == ~(USItype) 0) \
                __asm__ ("{a%I4|add%I4c} %1,%3,%4\n\t{ame|addme} %0,%2" \
                : "=r" ((USItype)(sh)), \
                "=&r" ((USItype)(sl)) \
                : "%r" ((USItype)(ah)), \
                "%r" ((USItype)(al)), \
                "rI" ((USItype)(bl))); \
        else \
                __asm__ ("{a%I5|add%I5c} %1,%4,%5\n\t{ae|adde} %0,%2,%3" \
                : "=r" ((USItype)(sh)), \
                "=&r" ((USItype)(sl)) \
                : "%r" ((USItype)(ah)), \
                "r" ((USItype)(bh)), \
                "%r" ((USItype)(al)), \
                "rI" ((USItype)(bl))); \
} while (0)

#define sub_ddmmss(sh, sl, ah, al, bh, bl) \
do { \
        if (__builtin_constant_p(ah) && (ah) == 0) \
                __asm__ ("{sf%I3|subf%I3c} %1,%4,%3\n\t{sfze|subfze} %0,%2" \
                : "=r" ((USItype)(sh)), \
                "=&r" ((USItype)(sl)) \
                : "r" ((USItype)(bh)), \
                "rI" ((USItype)(al)), \
                "r" ((USItype)(bl))); \
        else if (__builtin_constant_p(ah) && (ah) == ~(USItype) 0) \
                __asm__ ("{sf%I3|subf%I3c} %1,%4,%3\n\t{sfme|subfme} %0,%2" \
                : "=r" ((USItype)(sh)), \
                "=&r" ((USItype)(sl)) \
                : "r" ((USItype)(bh)), \
                "rI" ((USItype)(al)), \
                "r" ((USItype)(bl))); \
        else if (__builtin_constant_p(bh) && (bh) == 0) \
                __asm__ ("{sf%I3|subf%I3c} %1,%4,%3\n\t{ame|addme} %0,%2" \
                : "=r" ((USItype)(sh)), \
                "=&r" ((USItype)(sl)) \
                : "r" ((USItype)(ah)), \
                "rI" ((USItype)(al)), \
                "r" ((USItype)(bl))); \
        else if (__builtin_constant_p(bh) && (bh) == ~(USItype) 0) \
                __asm__ ("{sf%I3|subf%I3c} %1,%4,%3\n\t{aze|addze} %0,%2" \
                : "=r" ((USItype)(sh)), \
                "=&r" ((USItype)(sl)) \
                : "r" ((USItype)(ah)), \
                "rI" ((USItype)(al)), \
                "r" ((USItype)(bl))); \
        else \
                __asm__ ("{sf%I4|subf%I4c} %1,%5,%4\n\t{sfe|subfe} %0,%3,%2" \
                : "=r" ((USItype)(sh)), \
                "=&r" ((USItype)(sl)) \
                : "r" ((USItype)(ah)), \
                "r" ((USItype)(bh)), \
                "rI" ((USItype)(al)), \
                "r" ((USItype)(bl))); \
} while (0)
#if defined(_ARCH_PPC)

#define umul_ppmm(ph, pl, m0, m1) \
do { \
        USItype __m0 = (m0), __m1 = (m1); \
        __asm__ ("mulhwu %0,%1,%2" \
        : "=r" ((USItype) ph) \
        : "%r" (__m0), \
        "r" (__m1)); \
        (pl) = __m0 * __m1; \
} while (0)

#define UMUL_TIME 15

#define smul_ppmm(ph, pl, m0, m1) \
do { \
        SItype __m0 = (m0), __m1 = (m1); \
        __asm__ ("mulhw %0,%1,%2" \
        : "=r" ((SItype) ph) \
        : "%r" (__m0), \
        "r" (__m1)); \
        (pl) = __m0 * __m1; \
} while (0)

#define SMUL_TIME 14

#define UDIV_TIME 120
#else

#define umul_ppmm(xh, xl, m0, m1) \
do { \
        USItype __m0 = (m0), __m1 = (m1); \
        __asm__ ("mul %0,%2,%3" \
        : "=r" ((USItype)(xh)), \
        "=q" ((USItype)(xl)) \
        : "r" (__m0), \
        "r" (__m1)); \
        (xh) += ((((SItype) __m0 >> 31) & __m1) \
        + (((SItype) __m1 >> 31) & __m0)); \
} while (0)

#define UMUL_TIME 8

#define smul_ppmm(xh, xl, m0, m1) \
	__asm__ ("mul %0,%2,%3" \
        : "=r" ((SItype)(xh)), \
        "=q" ((SItype)(xl)) \
        : "r" (m0), \
        "r" (m1))

#define SMUL_TIME 4

#define sdiv_qrnnd(q, r, nh, nl, d) \
	__asm__ ("div %0,%2,%4" \
        : "=r" ((SItype)(q)), "=q" ((SItype)(r)) \
        : "r" ((SItype)(nh)), "1" ((SItype)(nl)), "r" ((SItype)(d)))

#define UDIV_TIME 100
#endif
#endif /* Power architecture variants.  */

/***************************************
        **************  PYR  ******************
        ***************************************/
#if defined(__pyr__) && W_TYPE_SIZE == 32

#define add_ssaaaa(sh, sl, ah, al, bh, bl) \
	__asm__ ("addw        %5,%1\n" \
        "addwc  %3,%0" \
        : "=r" ((USItype)(sh)), \
        "=&r" ((USItype)(sl)) \
        : "%0" ((USItype)(ah)), \
        "g" ((USItype)(bh)), \
        "%1" ((USItype)(al)), \
        "g" ((USItype)(bl)))

#define sub_ddmmss(sh, sl, ah, al, bh, bl) \
	__asm__ ("subw        %5,%1\n" \
        "subwb  %3,%0" \
        : "=r" ((USItype)(sh)), \
        "=&r" ((USItype)(sl)) \
        : "0" ((USItype)(ah)), \
        "g" ((USItype)(bh)), \
        "1" ((USItype)(al)), \
        "g" ((USItype)(bl)))
	/* This insn works on Pyramids with AP, XP, or MI CPUs, but not with SP.  */

#define umul_ppmm(w1, w0, u, v) \
	({union {UDItype __ll; \
        struct {USItype __h, __l; } __i; \
        } __xx; \
        __asm__ ("movw %1,%R0\n" \
        "uemul %2,%0" \
        : "=&r" (__xx.__ll) \
        : "g" ((USItype) (u)), \
        "g" ((USItype)(v))); \
        (w1) = __xx.__i.__h; (w0) = __xx.__i.__l; })
#endif /* __pyr__ */

/***************************************
        **************  RT/ROMP  **************
        ***************************************/
#if defined(__ibm032__) /* RT/ROMP */	&& W_TYPE_SIZE == 32

#define add_ssaaaa(sh, sl, ah, al, bh, bl) \
	__asm__ ("a %1,%5\n" \
        "ae %0,%3" \
        : "=r" ((USItype)(sh)), \
        "=&r" ((USItype)(sl)) \
        : "%0" ((USItype)(ah)), \
        "r" ((USItype)(bh)), \
        "%1" ((USItype)(al)), \
        "r" ((USItype)(bl)))

#define sub_ddmmss(sh, sl, ah, al, bh, bl) \
	__asm__ ("s %1,%5\n" \
        "se %0,%3" \
        : "=r" ((USItype)(sh)), \
        "=&r" ((USItype)(sl)) \
        : "0" ((USItype)(ah)), \
        "r" ((USItype)(bh)), \
        "1" ((USItype)(al)), \
        "r" ((USItype)(bl)))

#define umul_ppmm(ph, pl, m0, m1) \
do { \
        USItype __m0 = (m0), __m1 = (m1); \
        __asm__ ( \
        "s       r2,r2\n" \
        "mts    r10,%2\n" \
        "m      r2,%3\n" \
        "m      r2,%3\n" \
        "m      r2,%3\n" \
        "m      r2,%3\n" \
        "m      r2,%3\n" \
        "m      r2,%3\n" \
        "m      r2,%3\n" \
        "m      r2,%3\n" \
        "m      r2,%3\n" \
        "m      r2,%3\n" \
        "m      r2,%3\n" \
        "m      r2,%3\n" \
        "m      r2,%3\n" \
        "m      r2,%3\n" \
        "m      r2,%3\n" \
        "m      r2,%3\n" \
        "cas    %0,r2,r0\n" \
        "mfs    r10,%1" \
        : "=r" ((USItype)(ph)), \
        "=r" ((USItype)(pl)) \
        : "%r" (__m0), \
        "r" (__m1) \
        : "r2"); \
        (ph) += ((((SItype) __m0 >> 31) & __m1) \
        + (((SItype) __m1 >> 31) & __m0)); \
} while (0)

#define UMUL_TIME 20

#define UDIV_TIME 200
#endif /* RT/ROMP */

/***************************************
        **************  SH2  ******************
        ***************************************/
#if (defined(__sh2__) || defined(__sh3__) || defined(__SH4__)) \
	&& W_TYPE_SIZE == 32

#define umul_ppmm(w1, w0, u, v) \
	__asm__ ( \
        "dmulu.l %2,%3\n" \
        "sts    macl,%1\n" \
        "sts    mach,%0" \
        : "=r" ((USItype)(w1)), \
        "=r" ((USItype)(w0)) \
        : "r" ((USItype)(u)), \
        "r" ((USItype)(v)) \
        : "macl", "mach")

#define UMUL_TIME 5
#endif

/***************************************
        **************  SPARC   ****************
        ***************************************/
#if defined(__sparc__) && W_TYPE_SIZE == 32

#define add_ssaaaa(sh, sl, ah, al, bh, bl) \
	__asm__ ("addcc %r4,%5,%1\n" \
        "addx %r2,%3,%0" \
        : "=r" ((USItype)(sh)), \
        "=&r" ((USItype)(sl)) \
        : "%rJ" ((USItype)(ah)), \
        "rI" ((USItype)(bh)), \
        "%rJ" ((USItype)(al)), \
        "rI" ((USItype)(bl)) \
        __CLOBBER_CC)

#define sub_ddmmss(sh, sl, ah, al, bh, bl) \
	__asm__ ("subcc %r4,%5,%1\n" \
        "subx %r2,%3,%0" \
        : "=r" ((USItype)(sh)), \
        "=&r" ((USItype)(sl)) \
        : "rJ" ((USItype)(ah)), \
        "rI" ((USItype)(bh)), \
        "rJ" ((USItype)(al)), \
        "rI" ((USItype)(bl)) \
        __CLOBBER_CC)
#if defined(__sparc_v8__)
/* Don't match immediate range because, 1) it is not often useful,
        2) the 'I' flag thinks of the range as a 13 bit signed interval,
        while we want to match a 13 bit interval, sign extended to 32 bits,
        but INTERPRETED AS UNSIGNED.  */

#define umul_ppmm(w1, w0, u, v) \
	__asm__ ("umul %2,%3,%1;rd %%y,%0" \
        : "=r" ((USItype)(w1)), \
        "=r" ((USItype)(w0)) \
        : "r" ((USItype)(u)), \
        "r" ((USItype)(v)))

#define UMUL_TIME 5
#ifndef SUPERSPARC		/* SuperSPARC's udiv only handles 53 bit dividends */

#define udiv_qrnnd(q, r, n1, n0, d) \
do { \
        USItype __q; \
        __asm__ ("mov %1,%%y;nop;nop;nop;udiv %2,%3,%0" \
        : "=r" ((USItype)(__q)) \
        : "r" ((USItype)(n1)), \
        "r" ((USItype)(n0)), \
        "r" ((USItype)(d))); \
        (r) = (n0) - __q * (d); \
        (q) = __q; \
} while (0)

#define UDIV_TIME 25
#endif /* SUPERSPARC */
#else /* ! __sparc_v8__ */
#if defined(__sparclite__)
/* This has hardware multiply but not divide.  It also has two additional
        instructions scan (ffs from high bit) and divscc.  */

#define umul_ppmm(w1, w0, u, v) \
	__asm__ ("umul %2,%3,%1;rd %%y,%0" \
        : "=r" ((USItype)(w1)), \
        "=r" ((USItype)(w0)) \
        : "r" ((USItype)(u)), \
        "r" ((USItype)(v)))

#define UMUL_TIME 5

#define udiv_qrnnd(q, r, n1, n0, d) \
	__asm__ ("! Inlined udiv_qrnnd\n" \
        "wr     %%g0,%2,%%y     ! Not a delayed write for sparclite\n" \
        "tst    %%g0\n" \
        "divscc %3,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%%g1\n" \
        "divscc %%g1,%4,%0\n" \
        "rd     %%y,%1\n" \
        "bl,a 1f\n" \
        "add    %1,%4,%1\n" \
        "1:     ! End of inline udiv_qrnnd" \
        : "=r" ((USItype)(q)), \
        "=r" ((USItype)(r)) \
        : "r" ((USItype)(n1)), \
        "r" ((USItype)(n0)), \
        "rI" ((USItype)(d)) \
        : "%g1" __AND_CLOBBER_CC)

#define UDIV_TIME 37
#endif /* __sparclite__ */
#endif /* __sparc_v8__ */
	/* Default to sparc v7 versions of umul_ppmm and udiv_qrnnd.  */
#ifndef umul_ppmm

#define umul_ppmm(w1, w0, u, v) \
	__asm__ ("! Inlined umul_ppmm\n" \
        "wr     %%g0,%2,%%y     ! SPARC has 0-3 delay insn after a wr\n" \
        "sra    %3,31,%%g2      ! Don't move this insn\n" \
        "and    %2,%%g2,%%g2    ! Don't move this insn\n" \
        "andcc  %%g0,0,%%g1     ! Don't move this insn\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,%3,%%g1\n" \
        "mulscc %%g1,0,%%g1\n" \
        "add    %%g1,%%g2,%0\n" \
        "rd     %%y,%1" \
        : "=r" ((USItype)(w1)), \
        "=r" ((USItype)(w0)) \
        : "%rI" ((USItype)(u)), \
        "r" ((USItype)(v)) \
        : "%g1", "%g2" __AND_CLOBBER_CC)

#define UMUL_TIME 39		
/* 39 instructions */
/* It's quite necessary to add this much assembler for the sparc.
   The default udiv_qrnnd (in C) is more than 10 times slower!  */

#define udiv_qrnnd(q, r, n1, n0, d) \
  __asm__ ("! Inlined udiv_qrnnd\n\t"                                     \
           "mov 32,%%g1\n\t"                                            \
           "subcc       %1,%2,%%g0\n\t"                                 \
           "1:  bcs     5f\n\t"                                         \
           "addxcc %0,%0,%0     ! shift n1n0 and a q-bit in lsb\n\t"    \
           "sub %1,%2,%1        ! this kills msb of n\n\t"              \
           "addx        %1,%1,%1        ! so this can't give carry\n\t" \
           "subcc       %%g1,1,%%g1\n\t"                                \
           "2:  bne     1b\n\t"                                         \
           "subcc       %1,%2,%%g0\n\t"                                 \
           "bcs 3f\n\t"                                                 \
           "addxcc %0,%0,%0     ! shift n1n0 and a q-bit in lsb\n\t"    \
           "b           3f\n\t"                                         \
           "sub %1,%2,%1        ! this kills msb of n\n\t"              \
           "4:  sub     %1,%2,%1\n\t"                                   \
           "5:  addxcc  %1,%1,%1\n\t"                                   \
           "bcc 2b\n\t"                                                 \
           "subcc       %%g1,1,%%g1\n\t"                                \
           "! Got carry from n.  Subtract next step to cancel this carry.\n\t" \
           "bne 4b\n\t"                                                 \
           "addcc       %0,%0,%0        ! shift n1n0 and a 0-bit in lsb\n\t" \
           "sub %1,%2,%1\n\t"                                           \
           "3:  xnor    %0,0,%0\n\t"                                    \
           "! End of inline udiv_qrnnd\n"                               \
           : "=&r" ((USItype)(q)),                                      \
             "=&r" ((USItype)(r))                                       \
           : "r" ((USItype)(d)),                                        \
             "1" ((USItype)(n1)),                                       \
             "0" ((USItype)(n0)) : "%g1", "cc")

#define UDIV_TIME (3+7*32)      
/* 7 instructions/iteration. 32 iterations.  */
#endif
#endif /* __sparc__ */

/***************************************
        **************  VAX  ******************
        ***************************************/
#if defined(__vax__) && W_TYPE_SIZE == 32

#define add_ssaaaa(sh, sl, ah, al, bh, bl) \
	__asm__ ("addl2 %5,%1\n" \
        "adwc %3,%0" \
        : "=g" ((USItype)(sh)), \
        "=&g" ((USItype)(sl)) \
        : "%0" ((USItype)(ah)), \
        "g" ((USItype)(bh)), \
        "%1" ((USItype)(al)), \
        "g" ((USItype)(bl)))

#define sub_ddmmss(sh, sl, ah, al, bh, bl) \
	__asm__ ("subl2 %5,%1\n" \
        "sbwc %3,%0" \
        : "=g" ((USItype)(sh)), \
        "=&g" ((USItype)(sl)) \
        : "0" ((USItype)(ah)), \
        "g" ((USItype)(bh)), \
        "1" ((USItype)(al)), \
        "g" ((USItype)(bl)))

#define umul_ppmm(xh, xl, m0, m1) \
do { \
        union {UDItype __ll; \
        struct {USItype __l, __h; } __i; \
        } __xx; \
        USItype __m0 = (m0), __m1 = (m1); \
        __asm__ ("emul %1,%2,$0,%0" \
        : "=g" (__xx.__ll) \
        : "g" (__m0), \
        "g" (__m1)); \
        (xh) = __xx.__i.__h; (xl) = __xx.__i.__l; \
        (xh) += ((((SItype) __m0 >> 31) & __m1) \
        + (((SItype) __m1 >> 31) & __m0)); \
} while (0)

#define sdiv_qrnnd(q, r, n1, n0, d) \
do { \
        union {DItype __ll; \
        struct {SItype __l, __h; } __i; \
        } __xx; \
        __xx.__i.__h = n1; __xx.__i.__l = n0; \
        __asm__ ("ediv %3,%2,%0,%1" \
        : "=g" (q), "=g" (r) \
        : "g" (__xx.__ll), "g" (d)); \
} while (0)
#endif /* __vax__ */

/***************************************
        **************  Z8000   ****************
        ***************************************/
#if defined(__z8000__) && W_TYPE_SIZE == 16

#define add_ssaaaa(sh, sl, ah, al, bh, bl) \
	__asm__ ("add %H1,%H5\n\tadc  %H0,%H3" \
        : "=r" ((unsigned int)(sh)), \
        "=&r" ((unsigned int)(sl)) \
        : "%0" ((unsigned int)(ah)), \
        "r" ((unsigned int)(bh)), \
        "%1" ((unsigned int)(al)), \
        "rQR" ((unsigned int)(bl)))

#define sub_ddmmss(sh, sl, ah, al, bh, bl) \
	__asm__ ("sub %H1,%H5\n\tsbc  %H0,%H3" \
        : "=r" ((unsigned int)(sh)), \
        "=&r" ((unsigned int)(sl)) \
        : "0" ((unsigned int)(ah)), \
        "r" ((unsigned int)(bh)), \
        "1" ((unsigned int)(al)), \
        "rQR" ((unsigned int)(bl)))

#define umul_ppmm(xh, xl, m0, m1) \
do { \
        union {long int __ll; \
        struct {unsigned int __h, __l; } __i; \
        } __xx; \
        unsigned int __m0 = (m0), __m1 = (m1); \
        __asm__ ("mult      %S0,%H3" \
        : "=r" (__xx.__i.__h), \
        "=r" (__xx.__i.__l) \
        : "%1" (__m0), \
        "rQR" (__m1)); \
        (xh) = __xx.__i.__h; (xl) = __xx.__i.__l; \
        (xh) += ((((signed int) __m0 >> 15) & __m1) \
        + (((signed int) __m1 >> 15) & __m0)); \
} while (0)
#endif /* __z8000__ */

#endif /* __GNUC__ */

/***************************************
        ***********  Generic Versions   ********
        ***************************************/
#if !defined(umul_ppmm) && defined(__umulsidi3)

#define umul_ppmm(ph, pl, m0, m1) \
{ \
        UDWtype __ll = __umulsidi3(m0, m1); \
        ph = (UWtype) (__ll >> W_TYPE_SIZE); \
        pl = (UWtype) __ll; \
}
#endif

#if !defined(__umulsidi3)

#define __umulsidi3(u, v) \
	({UWtype __hi, __lo; \
        umul_ppmm(__hi, __lo, u, v); \
        ((UDWtype) __hi << W_TYPE_SIZE) | __lo; })
#endif

	/* If this machine has no inline assembler, use C macros.  */

#if !defined(add_ssaaaa)

#define add_ssaaaa(sh, sl, ah, al, bh, bl) \
do { \
        UWtype __x; \
        __x = (al) + (bl); \
        (sh) = (ah) + (bh) + (__x < (al)); \
        (sl) = __x; \
} while (0)
#endif

#if !defined(sub_ddmmss)

#define sub_ddmmss(sh, sl, ah, al, bh, bl) \
do { \
        UWtype __x; \
        __x = (al) - (bl); \
        (sh) = (ah) - (bh) - (__x > (al)); \
        (sl) = __x; \
} while (0)
#endif

#if !defined(umul_ppmm)

#define umul_ppmm(w1, w0, u, v) \
do { \
        UWtype __x0, __x1, __x2, __x3; \
        UHWtype __ul, __vl, __uh, __vh; \
        UWtype __u = (u), __v = (v); \
        \
        __ul = __ll_lowpart(__u); \
        __uh = __ll_highpart(__u); \
        __vl = __ll_lowpart(__v); \
        __vh = __ll_highpart(__v); \
        \
        __x0 = (UWtype) __ul * __vl; \
        __x1 = (UWtype) __ul * __vh; \
        __x2 = (UWtype) __uh * __vl; \
        __x3 = (UWtype) __uh * __vh; \
        \
        __x1 += __ll_highpart(__x0);/* this can't give carry */ \
        __x1 += __x2;		/* but this indeed can */ \
        if (__x1 < __x2)		/* did we get it? */ \
        __x3 += __ll_B;		/* yes, add it in the proper pos. */ \
        \
        (w1) = __x3 + __ll_highpart(__x1); \
        (w0) = (__ll_lowpart(__x1) << W_TYPE_SIZE/2) + __ll_lowpart(__x0); \
} while (0)
#endif

#if !defined(umul_ppmm)

#define smul_ppmm(w1, w0, u, v) \
do { \
        UWtype __w1; \
        UWtype __m0 = (u), __m1 = (v); \
        umul_ppmm(__w1, w0, __m0, __m1); \
        (w1) = __w1 - (-(__m0 >> (W_TYPE_SIZE - 1)) & __m1) \
        - (-(__m1 >> (W_TYPE_SIZE - 1)) & __m0); \
} while (0)
#endif

	/* Define this unconditionally, so it can be used for debugging.  */

#define __udiv_qrnnd_c(q, r, n1, n0, d) \
do { \
        UWtype __d1, __d0, __q1, __q0, __r1, __r0, __m; \
        __d1 = __ll_highpart(d); \
        __d0 = __ll_lowpart(d); \
        \
        __r1 = (n1) % __d1; \
        __q1 = (n1) / __d1; \
        __m = (UWtype) __q1 * __d0; \
        __r1 = __r1 * __ll_B | __ll_highpart(n0); \
        if (__r1 < __m) { \
                __q1--, __r1 += (d); \
                if (__r1 >= (d)) /* i.e. we didn't get carry when adding to __r1 */ \
                if (__r1 < __m) \
                        __q1--, __r1 += (d); \
        } \
        __r1 -= __m; \
        \
        __r0 = __r1 % __d1; \
        __q0 = __r1 / __d1; \
        __m = (UWtype) __q0 * __d0; \
        __r0 = __r0 * __ll_B | __ll_lowpart(n0); \
        if (__r0 < __m) { \
                __q0--, __r0 += (d); \
                if (__r0 >= (d)) \
                        if (__r0 < __m) \
                                __q0--, __r0 += (d); \
        } \
        __r0 -= __m; \
        \
        (q) = (UWtype) __q1 * __ll_B | __q0; \
        (r) = __r0; \
} while (0)

/* If the processor has no udiv_qrnnd but sdiv_qrnnd, go through
        __udiv_w_sdiv (defined in libgcc or elsewhere).  */
#if !defined(udiv_qrnnd) && defined(sdiv_qrnnd)

#define udiv_qrnnd(q, r, nh, nl, d) \
do { \
        UWtype __r; \
        (q) = __MPN(udiv_w_sdiv) (&__r, nh, nl, d); \
        (r) = __r; \
} while (0)
#endif

	/* If udiv_qrnnd was not defined for this processor, use __udiv_qrnnd_c.  */
#if !defined(udiv_qrnnd)

#define UDIV_NEEDS_NORMALIZATION 1

#define udiv_qrnnd __udiv_qrnnd_c
#endif

#ifndef UDIV_NEEDS_NORMALIZATION

#define UDIV_NEEDS_NORMALIZATION 0
#endif

Overall Contributors

PersonTokensPropCommitsCommitProp
Dmitry Kasatkin225695.07%110.00%
Manuel Lauss461.94%110.00%
James Hogan251.05%110.00%
Jaedon Shin160.67%110.00%
Richard Henderson110.46%110.00%
Helge Deller80.34%110.00%
Stefan Agner40.17%110.00%
David S. Miller40.17%110.00%
David Howells20.08%110.00%
Christoph Hellwig10.04%110.00%
Total2373100.00%10100.00%
Directory: lib/mpi
Information contained on this website is for historical information purposes only and does not indicate or represent copyright ownership.
Created with cregit.