158f0484fSRodney W. Grimes /*- 28a16b7a1SPedro F. Giffuni * SPDX-License-Identifier: BSD-3-Clause 38a16b7a1SPedro F. Giffuni * 458f0484fSRodney W. Grimes * Copyright (c) 1992, 1993 558f0484fSRodney W. Grimes * The Regents of the University of California. All rights reserved. 658f0484fSRodney W. Grimes * 758f0484fSRodney W. Grimes * This software was developed by the Computer Systems Engineering group 858f0484fSRodney W. Grimes * at Lawrence Berkeley Laboratory under DARPA contract BG 91-66 and 958f0484fSRodney W. Grimes * contributed to Berkeley. 1058f0484fSRodney W. Grimes * 1158f0484fSRodney W. Grimes * Redistribution and use in source and binary forms, with or without 1258f0484fSRodney W. Grimes * modification, are permitted provided that the following conditions 1358f0484fSRodney W. Grimes * are met: 1458f0484fSRodney W. Grimes * 1. Redistributions of source code must retain the above copyright 1558f0484fSRodney W. Grimes * notice, this list of conditions and the following disclaimer. 1658f0484fSRodney W. Grimes * 2. Redistributions in binary form must reproduce the above copyright 1758f0484fSRodney W. Grimes * notice, this list of conditions and the following disclaimer in the 1858f0484fSRodney W. Grimes * documentation and/or other materials provided with the distribution. 19fbbd9655SWarner Losh * 3. Neither the name of the University nor the names of its contributors 2058f0484fSRodney W. Grimes * may be used to endorse or promote products derived from this software 2158f0484fSRodney W. Grimes * without specific prior written permission. 2258f0484fSRodney W. Grimes * 2358f0484fSRodney W. Grimes * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND 2458f0484fSRodney W. Grimes * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE 2558f0484fSRodney W. Grimes * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE 2658f0484fSRodney W. Grimes * ARE DISCLAIMED. IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE 2758f0484fSRodney W. Grimes * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL 2858f0484fSRodney W. Grimes * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS 2958f0484fSRodney W. Grimes * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) 3058f0484fSRodney W. Grimes * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT 3158f0484fSRodney W. Grimes * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY 3258f0484fSRodney W. Grimes * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF 3358f0484fSRodney W. Grimes * SUCH DAMAGE. 3458f0484fSRodney W. Grimes */ 3558f0484fSRodney W. Grimes 3658f0484fSRodney W. Grimes /* 3758f0484fSRodney W. Grimes * Quad arithmetic. 3858f0484fSRodney W. Grimes * 3958f0484fSRodney W. Grimes * This library makes the following assumptions: 4058f0484fSRodney W. Grimes * 4158f0484fSRodney W. Grimes * - The type long long (aka quad_t) exists. 4258f0484fSRodney W. Grimes * 4358f0484fSRodney W. Grimes * - A quad variable is exactly twice as long as `long'. 4458f0484fSRodney W. Grimes * 4558f0484fSRodney W. Grimes * - The machine's arithmetic is two's complement. 4658f0484fSRodney W. Grimes * 4758f0484fSRodney W. Grimes * This library can provide 128-bit arithmetic on a machine with 128-bit 4858f0484fSRodney W. Grimes * quads and 64-bit longs, for instance, or 96-bit arithmetic on machines 4958f0484fSRodney W. Grimes * with 48-bit longs. 5058f0484fSRodney W. Grimes */ 5158f0484fSRodney W. Grimes 5258f0484fSRodney W. Grimes #include <sys/types.h> 5358f0484fSRodney W. Grimes #include <limits.h> 5458f0484fSRodney W. Grimes 5558f0484fSRodney W. Grimes /* 5658f0484fSRodney W. Grimes * Depending on the desired operation, we view a `long long' (aka quad_t) in 5758f0484fSRodney W. Grimes * one or more of the following formats. 5858f0484fSRodney W. Grimes */ 5958f0484fSRodney W. Grimes union uu { 6058f0484fSRodney W. Grimes quad_t q; /* as a (signed) quad */ 6158f0484fSRodney W. Grimes quad_t uq; /* as an unsigned quad */ 6258f0484fSRodney W. Grimes long sl[2]; /* as two signed longs */ 6358f0484fSRodney W. Grimes u_long ul[2]; /* as two unsigned longs */ 6458f0484fSRodney W. Grimes }; 6558f0484fSRodney W. Grimes 6658f0484fSRodney W. Grimes /* 6758f0484fSRodney W. Grimes * Define high and low longwords. 6858f0484fSRodney W. Grimes */ 6958f0484fSRodney W. Grimes #define H _QUAD_HIGHWORD 7058f0484fSRodney W. Grimes #define L _QUAD_LOWWORD 7158f0484fSRodney W. Grimes 7258f0484fSRodney W. Grimes /* 7358f0484fSRodney W. Grimes * Total number of bits in a quad_t and in the pieces that make it up. 7458f0484fSRodney W. Grimes * These are used for shifting, and also below for halfword extraction 7558f0484fSRodney W. Grimes * and assembly. 7658f0484fSRodney W. Grimes */ 7758f0484fSRodney W. Grimes #define QUAD_BITS (sizeof(quad_t) * CHAR_BIT) 7858f0484fSRodney W. Grimes #define LONG_BITS (sizeof(long) * CHAR_BIT) 7958f0484fSRodney W. Grimes #define HALF_BITS (sizeof(long) * CHAR_BIT / 2) 8058f0484fSRodney W. Grimes 8158f0484fSRodney W. Grimes /* 8258f0484fSRodney W. Grimes * Extract high and low shortwords from longword, and move low shortword of 8358f0484fSRodney W. Grimes * longword to upper half of long, i.e., produce the upper longword of 8458f0484fSRodney W. Grimes * ((quad_t)(x) << (number_of_bits_in_long/2)). (`x' must actually be u_long.) 8558f0484fSRodney W. Grimes * 8658f0484fSRodney W. Grimes * These are used in the multiply code, to split a longword into upper 8758f0484fSRodney W. Grimes * and lower halves, and to reassemble a product as a quad_t, shifted left 8858f0484fSRodney W. Grimes * (sizeof(long)*CHAR_BIT/2). 8958f0484fSRodney W. Grimes */ 9058f0484fSRodney W. Grimes #define HHALF(x) ((x) >> HALF_BITS) 91*feb1d550SDimitry Andric #define LHALF(x) ((x) & ((1L << HALF_BITS) - 1)) 9258f0484fSRodney W. Grimes #define LHUP(x) ((x) << HALF_BITS) 9358f0484fSRodney W. Grimes 94a3b5200eSStefan Farfeleder int __cmpdi2(quad_t a, quad_t b); 95c05ac53bSDavid E. O'Brien quad_t __divdi3(quad_t a, quad_t b); 96c05ac53bSDavid E. O'Brien quad_t __moddi3(quad_t a, quad_t b); 97c05ac53bSDavid E. O'Brien u_quad_t __qdivrem(u_quad_t u, u_quad_t v, u_quad_t *rem); 98a3b5200eSStefan Farfeleder int __ucmpdi2(u_quad_t a, u_quad_t b); 99c05ac53bSDavid E. O'Brien u_quad_t __udivdi3(u_quad_t a, u_quad_t b); 100c05ac53bSDavid E. O'Brien u_quad_t __umoddi3(u_quad_t a, u_quad_t b); 10158f0484fSRodney W. Grimes 10258f0484fSRodney W. Grimes typedef unsigned int qshift_t; 103