1 /*- 2 * Copyright (c) 1992, 1993 3 * The Regents of the University of California. All rights reserved. 4 * 5 * This software was developed by the Computer Systems Engineering group 6 * at Lawrence Berkeley Laboratory under DARPA contract BG 91-66 and 7 * contributed to Berkeley. 8 * 9 * Redistribution and use in source and binary forms, with or without 10 * modification, are permitted provided that the following conditions 11 * are met: 12 * 1. Redistributions of source code must retain the above copyright 13 * notice, this list of conditions and the following disclaimer. 14 * 2. Redistributions in binary form must reproduce the above copyright 15 * notice, this list of conditions and the following disclaimer in the 16 * documentation and/or other materials provided with the distribution. 17 * 3. Neither the name of the University nor the names of its contributors 18 * may be used to endorse or promote products derived from this software 19 * without specific prior written permission. 20 * 21 * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND 22 * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE 23 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE 24 * ARE DISCLAIMED. IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE 25 * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL 26 * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS 27 * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) 28 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT 29 * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY 30 * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF 31 * SUCH DAMAGE. 32 * 33 * $OpenBSD: quad.h,v 1.9 2019/11/10 22:23:28 guenther Exp $ 34 */ 35 36 /* 37 * Quad arithmetic. 38 * 39 * This library makes the following assumptions: 40 * 41 * - The type long long (aka quad_t) exists. 42 * 43 * - A quad variable is exactly twice as long as `int'. 44 * 45 * - The machine's arithmetic is two's complement. 46 * 47 * This library can provide 128-bit arithmetic on a machine with 128-bit 48 * quads and 64-bit ints, for instance, or 96-bit arithmetic on machines 49 * with 48-bit ints. 50 */ 51 52 #include <sys/types.h> 53 #include <limits.h> 54 55 /* 56 * Depending on the desired operation, we view a `long long' (aka quad_t) in 57 * one or more of the following formats. 58 */ 59 union uu { 60 quad_t q; /* as a (signed) quad */ 61 u_quad_t uq; /* as an unsigned quad */ 62 int sl[2]; /* as two signed ints */ 63 u_int ul[2]; /* as two unsigned ints */ 64 }; 65 66 /* 67 * Define high and low parts of a quad_t. 68 */ 69 #define H _QUAD_HIGHWORD 70 #define L _QUAD_LOWWORD 71 72 /* 73 * Total number of bits in a quad_t and in the pieces that make it up. 74 * These are used for shifting, and also below for halfword extraction 75 * and assembly. 76 */ 77 #define QUAD_BITS (sizeof(quad_t) * CHAR_BIT) 78 #define INT_BITS (sizeof(int) * CHAR_BIT) 79 #define HALF_BITS (sizeof(int) * CHAR_BIT / 2) 80 81 /* 82 * Extract high and low shortwords from longword, and move low shortword of 83 * longword to upper half of long, i.e., produce the upper longword of 84 * ((quad_t)(x) << (number_of_bits_in_int/2)). (`x' must actually be u_int.) 85 * 86 * These are used in the multiply code, to split a longword into upper 87 * and lower halves, and to reassemble a product as a quad_t, shifted left 88 * (sizeof(int)*CHAR_BIT/2). 89 */ 90 #define HHALF(x) ((u_int)(x) >> HALF_BITS) 91 #define LHALF(x) ((u_int)(x) & (((int)1 << HALF_BITS) - 1)) 92 #define LHUP(x) ((u_int)(x) << HALF_BITS) 93 94 typedef unsigned int qshift_t; 95 96 quad_t __adddi3(quad_t, quad_t) __dso_protected; 97 quad_t __anddi3(quad_t, quad_t) __dso_protected; 98 quad_t __ashldi3(quad_t, qshift_t) __dso_protected; 99 quad_t __ashrdi3(quad_t, qshift_t) __dso_protected; 100 int __cmpdi2(quad_t, quad_t) __dso_protected; 101 quad_t __divdi3(quad_t, quad_t) __dso_protected; 102 quad_t __fixdfdi(double) __dso_protected; 103 quad_t __fixsfdi(float) __dso_protected; 104 u_quad_t __fixunsdfdi(double) __dso_protected; 105 u_quad_t __fixunssfdi(float) __dso_protected; 106 double __floatdidf(quad_t) __dso_protected; 107 float __floatdisf(quad_t) __dso_protected; 108 double __floatunsdidf(u_quad_t) __dso_protected; 109 quad_t __iordi3(quad_t, quad_t) __dso_protected; 110 quad_t __lshldi3(quad_t, qshift_t) __dso_protected; 111 quad_t __lshrdi3(quad_t, qshift_t) __dso_protected; 112 quad_t __moddi3(quad_t, quad_t) __dso_protected; 113 quad_t __muldi3(quad_t, quad_t) __dso_protected; 114 quad_t __negdi2(quad_t) __dso_protected; 115 quad_t __one_cmpldi2(quad_t) __dso_protected; 116 u_quad_t __qdivrem(u_quad_t, u_quad_t, u_quad_t *) __dso_protected; 117 quad_t __subdi3(quad_t, quad_t) __dso_protected; 118 int __ucmpdi2(u_quad_t, u_quad_t) __dso_protected; 119 u_quad_t __udivdi3(u_quad_t, u_quad_t ) __dso_protected; 120 u_quad_t __umoddi3(u_quad_t, u_quad_t ) __dso_protected; 121 quad_t __xordi3(quad_t, quad_t) __dso_protected; 122