gl/float+.h

   1 /* Supplemental information about the floating-point formats.
   2    Copyright (C) 2007, 2009, 2010 Free Software Foundation, Inc.
   3    Written by Bruno Haible <bruno@clisp.org>, 2007.
   4
   5    This program is free software; you can redistribute it and/or modify
   6    it under the terms of the GNU General Public License as published by
   7    the Free Software Foundation; either version 3, or (at your option)
   8    any later version.
   9
  10    This program is distributed in the hope that it will be useful,
  11    but WITHOUT ANY WARRANTY; without even the implied warranty of
  12    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  13    GNU General Public License for more details.
  14
  15    You should have received a copy of the GNU General Public License
  16    along with this program; if not, write to the Free Software Foundation,
  17    Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.  */
  18
  19 #ifndef _FLOATPLUS_H
  20 #define _FLOATPLUS_H
  21
  22 #include <float.h>
  23 #include <limits.h>
  24
  25 /* Number of bits in the mantissa of a floating-point number, including the
  26    "hidden bit".  */
  27 #if FLT_RADIX == 2
  28 # define FLT_MANT_BIT FLT_MANT_DIG
  29 # define DBL_MANT_BIT DBL_MANT_DIG
  30 # define LDBL_MANT_BIT LDBL_MANT_DIG
  31 #elif FLT_RADIX == 4
  32 # define FLT_MANT_BIT (FLT_MANT_DIG * 2)
  33 # define DBL_MANT_BIT (DBL_MANT_DIG * 2)
  34 # define LDBL_MANT_BIT (LDBL_MANT_DIG * 2)
  35 #elif FLT_RADIX == 16
  36 # define FLT_MANT_BIT (FLT_MANT_DIG * 4)
  37 # define DBL_MANT_BIT (DBL_MANT_DIG * 4)
  38 # define LDBL_MANT_BIT (LDBL_MANT_DIG * 4)
  39 #endif
  40
  41 /* Bit mask that can be used to mask the exponent, as an unsigned number.  */
  42 #define FLT_EXP_MASK ((FLT_MAX_EXP - FLT_MIN_EXP) | 7)
  43 #define DBL_EXP_MASK ((DBL_MAX_EXP - DBL_MIN_EXP) | 7)
  44 #define LDBL_EXP_MASK ((LDBL_MAX_EXP - LDBL_MIN_EXP) | 7)
  45
  46 /* Number of bits used for the exponent of a floating-point number, including
  47    the exponent's sign.  */
  48 #define FLT_EXP_BIT \
  49   (FLT_EXP_MASK < 0x100 ? 8 : \
  50    FLT_EXP_MASK < 0x200 ? 9 : \
  51    FLT_EXP_MASK < 0x400 ? 10 : \
  52    FLT_EXP_MASK < 0x800 ? 11 : \
  53    FLT_EXP_MASK < 0x1000 ? 12 : \
  54    FLT_EXP_MASK < 0x2000 ? 13 : \
  55    FLT_EXP_MASK < 0x4000 ? 14 : \
  56    FLT_EXP_MASK < 0x8000 ? 15 : \
  57    FLT_EXP_MASK < 0x10000 ? 16 : \
  58    FLT_EXP_MASK < 0x20000 ? 17 : \
  59    FLT_EXP_MASK < 0x40000 ? 18 : \
  60    FLT_EXP_MASK < 0x80000 ? 19 : \
  61    FLT_EXP_MASK < 0x100000 ? 20 : \
  62    FLT_EXP_MASK < 0x200000 ? 21 : \
  63    FLT_EXP_MASK < 0x400000 ? 22 : \
  64    FLT_EXP_MASK < 0x800000 ? 23 : \
  65    FLT_EXP_MASK < 0x1000000 ? 24 : \
  66    FLT_EXP_MASK < 0x2000000 ? 25 : \
  67    FLT_EXP_MASK < 0x4000000 ? 26 : \
  68    FLT_EXP_MASK < 0x8000000 ? 27 : \
  69    FLT_EXP_MASK < 0x10000000 ? 28 : \
  70    FLT_EXP_MASK < 0x20000000 ? 29 : \
  71    FLT_EXP_MASK < 0x40000000 ? 30 : \
  72    FLT_EXP_MASK <= 0x7fffffff ? 31 : \
  73    32)
  74 #define DBL_EXP_BIT \
  75   (DBL_EXP_MASK < 0x100 ? 8 : \
  76    DBL_EXP_MASK < 0x200 ? 9 : \
  77    DBL_EXP_MASK < 0x400 ? 10 : \
  78    DBL_EXP_MASK < 0x800 ? 11 : \
  79    DBL_EXP_MASK < 0x1000 ? 12 : \
  80    DBL_EXP_MASK < 0x2000 ? 13 : \
  81    DBL_EXP_MASK < 0x4000 ? 14 : \
  82    DBL_EXP_MASK < 0x8000 ? 15 : \
  83    DBL_EXP_MASK < 0x10000 ? 16 : \
  84    DBL_EXP_MASK < 0x20000 ? 17 : \
  85    DBL_EXP_MASK < 0x40000 ? 18 : \
  86    DBL_EXP_MASK < 0x80000 ? 19 : \
  87    DBL_EXP_MASK < 0x100000 ? 20 : \
  88    DBL_EXP_MASK < 0x200000 ? 21 : \
  89    DBL_EXP_MASK < 0x400000 ? 22 : \
  90    DBL_EXP_MASK < 0x800000 ? 23 : \
  91    DBL_EXP_MASK < 0x1000000 ? 24 : \
  92    DBL_EXP_MASK < 0x2000000 ? 25 : \
  93    DBL_EXP_MASK < 0x4000000 ? 26 : \
  94    DBL_EXP_MASK < 0x8000000 ? 27 : \
  95    DBL_EXP_MASK < 0x10000000 ? 28 : \
  96    DBL_EXP_MASK < 0x20000000 ? 29 : \
  97    DBL_EXP_MASK < 0x40000000 ? 30 : \
  98    DBL_EXP_MASK <= 0x7fffffff ? 31 : \
  99    32)
 100 #define LDBL_EXP_BIT \
 101   (LDBL_EXP_MASK < 0x100 ? 8 : \
 102    LDBL_EXP_MASK < 0x200 ? 9 : \
 103    LDBL_EXP_MASK < 0x400 ? 10 : \
 104    LDBL_EXP_MASK < 0x800 ? 11 : \
 105    LDBL_EXP_MASK < 0x1000 ? 12 : \
 106    LDBL_EXP_MASK < 0x2000 ? 13 : \
 107    LDBL_EXP_MASK < 0x4000 ? 14 : \
 108    LDBL_EXP_MASK < 0x8000 ? 15 : \
 109    LDBL_EXP_MASK < 0x10000 ? 16 : \
 110    LDBL_EXP_MASK < 0x20000 ? 17 : \
 111    LDBL_EXP_MASK < 0x40000 ? 18 : \
 112    LDBL_EXP_MASK < 0x80000 ? 19 : \
 113    LDBL_EXP_MASK < 0x100000 ? 20 : \
 114    LDBL_EXP_MASK < 0x200000 ? 21 : \
 115    LDBL_EXP_MASK < 0x400000 ? 22 : \
 116    LDBL_EXP_MASK < 0x800000 ? 23 : \
 117    LDBL_EXP_MASK < 0x1000000 ? 24 : \
 118    LDBL_EXP_MASK < 0x2000000 ? 25 : \
 119    LDBL_EXP_MASK < 0x4000000 ? 26 : \
 120    LDBL_EXP_MASK < 0x8000000 ? 27 : \
 121    LDBL_EXP_MASK < 0x10000000 ? 28 : \
 122    LDBL_EXP_MASK < 0x20000000 ? 29 : \
 123    LDBL_EXP_MASK < 0x40000000 ? 30 : \
 124    LDBL_EXP_MASK <= 0x7fffffff ? 31 : \
 125    32)
 126
 127 /* Number of bits used for a floating-point number: the mantissa (not
 128    counting the "hidden bit", since it may or may not be explicit), the
 129    exponent, and the sign.  */
 130 #define FLT_TOTAL_BIT ((FLT_MANT_BIT - 1) + FLT_EXP_BIT + 1)
 131 #define DBL_TOTAL_BIT ((DBL_MANT_BIT - 1) + DBL_EXP_BIT + 1)
 132 #define LDBL_TOTAL_BIT ((LDBL_MANT_BIT - 1) + LDBL_EXP_BIT + 1)
 133
 134 /* Number of bytes used for a floating-point number.
 135    This can be smaller than the 'sizeof'.  For example, on i386 systems,
 136    'long double' most often have LDBL_MANT_BIT = 64, LDBL_EXP_BIT = 16, hence
 137    LDBL_TOTAL_BIT = 80 bits, i.e. 10 bytes of consecutive memory, but
 138    sizeof (long double) = 12 or = 16.  */
 139 #define SIZEOF_FLT ((FLT_TOTAL_BIT + CHAR_BIT - 1) / CHAR_BIT)
 140 #define SIZEOF_DBL ((DBL_TOTAL_BIT + CHAR_BIT - 1) / CHAR_BIT)
 141 #define SIZEOF_LDBL ((LDBL_TOTAL_BIT + CHAR_BIT - 1) / CHAR_BIT)
 142
 143 /* Verify that SIZEOF_FLT <= sizeof (float) etc.  */
 144 typedef int verify_sizeof_flt[2 * (SIZEOF_FLT <= sizeof (float)) - 1];
 145 typedef int verify_sizeof_dbl[2 * (SIZEOF_DBL <= sizeof (double)) - 1];
 146 typedef int verify_sizeof_ldbl[2 * (SIZEOF_LDBL <= sizeof (long double)) - 1];
 147
 148 #endif /* _FLOATPLUS_H */