"Fossies" - the Fresh Open Source Software Archive

Member "gnuastro-0.8/bootstrapped/tests/float+.h" (28 Dec 2018, 5606 Bytes) of package /linux/privat/gnuastro-0.8.tar.lz:


As a special service "Fossies" has tried to format the requested source page into HTML format using (guessed) C and C++ source code syntax highlighting (style: standard) with prefixed line numbers and code folding option. Alternatively you can here view or download the uninterpreted source code file. See also the last Fossies "Diffs" side-by-side code changes report for "float+.h": 0.5_vs_0.6.

    1 /* Supplemental information about the floating-point formats.
    2    Copyright (C) 2007, 2009-2018 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, see <https://www.gnu.org/licenses/>.  */
   17 
   18 #ifndef _FLOATPLUS_H
   19 #define _FLOATPLUS_H
   20 
   21 #include <float.h>
   22 #include <limits.h>
   23 
   24 /* Number of bits in the mantissa of a floating-point number, including the
   25    "hidden bit".  */
   26 #if FLT_RADIX == 2
   27 # define FLT_MANT_BIT FLT_MANT_DIG
   28 # define DBL_MANT_BIT DBL_MANT_DIG
   29 # define LDBL_MANT_BIT LDBL_MANT_DIG
   30 #elif FLT_RADIX == 4
   31 # define FLT_MANT_BIT (FLT_MANT_DIG * 2)
   32 # define DBL_MANT_BIT (DBL_MANT_DIG * 2)
   33 # define LDBL_MANT_BIT (LDBL_MANT_DIG * 2)
   34 #elif FLT_RADIX == 16
   35 # define FLT_MANT_BIT (FLT_MANT_DIG * 4)
   36 # define DBL_MANT_BIT (DBL_MANT_DIG * 4)
   37 # define LDBL_MANT_BIT (LDBL_MANT_DIG * 4)
   38 #endif
   39 
   40 /* Bit mask that can be used to mask the exponent, as an unsigned number.  */
   41 #define FLT_EXP_MASK ((FLT_MAX_EXP - FLT_MIN_EXP) | 7)
   42 #define DBL_EXP_MASK ((DBL_MAX_EXP - DBL_MIN_EXP) | 7)
   43 #define LDBL_EXP_MASK ((LDBL_MAX_EXP - LDBL_MIN_EXP) | 7)
   44 
   45 /* Number of bits used for the exponent of a floating-point number, including
   46    the exponent's sign.  */
   47 #define FLT_EXP_BIT \
   48   (FLT_EXP_MASK < 0x100 ? 8 : \
   49    FLT_EXP_MASK < 0x200 ? 9 : \
   50    FLT_EXP_MASK < 0x400 ? 10 : \
   51    FLT_EXP_MASK < 0x800 ? 11 : \
   52    FLT_EXP_MASK < 0x1000 ? 12 : \
   53    FLT_EXP_MASK < 0x2000 ? 13 : \
   54    FLT_EXP_MASK < 0x4000 ? 14 : \
   55    FLT_EXP_MASK < 0x8000 ? 15 : \
   56    FLT_EXP_MASK < 0x10000 ? 16 : \
   57    FLT_EXP_MASK < 0x20000 ? 17 : \
   58    FLT_EXP_MASK < 0x40000 ? 18 : \
   59    FLT_EXP_MASK < 0x80000 ? 19 : \
   60    FLT_EXP_MASK < 0x100000 ? 20 : \
   61    FLT_EXP_MASK < 0x200000 ? 21 : \
   62    FLT_EXP_MASK < 0x400000 ? 22 : \
   63    FLT_EXP_MASK < 0x800000 ? 23 : \
   64    FLT_EXP_MASK < 0x1000000 ? 24 : \
   65    FLT_EXP_MASK < 0x2000000 ? 25 : \
   66    FLT_EXP_MASK < 0x4000000 ? 26 : \
   67    FLT_EXP_MASK < 0x8000000 ? 27 : \
   68    FLT_EXP_MASK < 0x10000000 ? 28 : \
   69    FLT_EXP_MASK < 0x20000000 ? 29 : \
   70    FLT_EXP_MASK < 0x40000000 ? 30 : \
   71    FLT_EXP_MASK <= 0x7fffffff ? 31 : \
   72    32)
   73 #define DBL_EXP_BIT \
   74   (DBL_EXP_MASK < 0x100 ? 8 : \
   75    DBL_EXP_MASK < 0x200 ? 9 : \
   76    DBL_EXP_MASK < 0x400 ? 10 : \
   77    DBL_EXP_MASK < 0x800 ? 11 : \
   78    DBL_EXP_MASK < 0x1000 ? 12 : \
   79    DBL_EXP_MASK < 0x2000 ? 13 : \
   80    DBL_EXP_MASK < 0x4000 ? 14 : \
   81    DBL_EXP_MASK < 0x8000 ? 15 : \
   82    DBL_EXP_MASK < 0x10000 ? 16 : \
   83    DBL_EXP_MASK < 0x20000 ? 17 : \
   84    DBL_EXP_MASK < 0x40000 ? 18 : \
   85    DBL_EXP_MASK < 0x80000 ? 19 : \
   86    DBL_EXP_MASK < 0x100000 ? 20 : \
   87    DBL_EXP_MASK < 0x200000 ? 21 : \
   88    DBL_EXP_MASK < 0x400000 ? 22 : \
   89    DBL_EXP_MASK < 0x800000 ? 23 : \
   90    DBL_EXP_MASK < 0x1000000 ? 24 : \
   91    DBL_EXP_MASK < 0x2000000 ? 25 : \
   92    DBL_EXP_MASK < 0x4000000 ? 26 : \
   93    DBL_EXP_MASK < 0x8000000 ? 27 : \
   94    DBL_EXP_MASK < 0x10000000 ? 28 : \
   95    DBL_EXP_MASK < 0x20000000 ? 29 : \
   96    DBL_EXP_MASK < 0x40000000 ? 30 : \
   97    DBL_EXP_MASK <= 0x7fffffff ? 31 : \
   98    32)
   99 #define LDBL_EXP_BIT \
  100   (LDBL_EXP_MASK < 0x100 ? 8 : \
  101    LDBL_EXP_MASK < 0x200 ? 9 : \
  102    LDBL_EXP_MASK < 0x400 ? 10 : \
  103    LDBL_EXP_MASK < 0x800 ? 11 : \
  104    LDBL_EXP_MASK < 0x1000 ? 12 : \
  105    LDBL_EXP_MASK < 0x2000 ? 13 : \
  106    LDBL_EXP_MASK < 0x4000 ? 14 : \
  107    LDBL_EXP_MASK < 0x8000 ? 15 : \
  108    LDBL_EXP_MASK < 0x10000 ? 16 : \
  109    LDBL_EXP_MASK < 0x20000 ? 17 : \
  110    LDBL_EXP_MASK < 0x40000 ? 18 : \
  111    LDBL_EXP_MASK < 0x80000 ? 19 : \
  112    LDBL_EXP_MASK < 0x100000 ? 20 : \
  113    LDBL_EXP_MASK < 0x200000 ? 21 : \
  114    LDBL_EXP_MASK < 0x400000 ? 22 : \
  115    LDBL_EXP_MASK < 0x800000 ? 23 : \
  116    LDBL_EXP_MASK < 0x1000000 ? 24 : \
  117    LDBL_EXP_MASK < 0x2000000 ? 25 : \
  118    LDBL_EXP_MASK < 0x4000000 ? 26 : \
  119    LDBL_EXP_MASK < 0x8000000 ? 27 : \
  120    LDBL_EXP_MASK < 0x10000000 ? 28 : \
  121    LDBL_EXP_MASK < 0x20000000 ? 29 : \
  122    LDBL_EXP_MASK < 0x40000000 ? 30 : \
  123    LDBL_EXP_MASK <= 0x7fffffff ? 31 : \
  124    32)
  125 
  126 /* Number of bits used for a floating-point number: the mantissa (not
  127    counting the "hidden bit", since it may or may not be explicit), the
  128    exponent, and the sign.  */
  129 #define FLT_TOTAL_BIT ((FLT_MANT_BIT - 1) + FLT_EXP_BIT + 1)
  130 #define DBL_TOTAL_BIT ((DBL_MANT_BIT - 1) + DBL_EXP_BIT + 1)
  131 #define LDBL_TOTAL_BIT ((LDBL_MANT_BIT - 1) + LDBL_EXP_BIT + 1)
  132 
  133 /* Number of bytes used for a floating-point number.
  134    This can be smaller than the 'sizeof'.  For example, on i386 systems,
  135    'long double' most often have LDBL_MANT_BIT = 64, LDBL_EXP_BIT = 16, hence
  136    LDBL_TOTAL_BIT = 80 bits, i.e. 10 bytes of consecutive memory, but
  137    sizeof (long double) = 12 or = 16.  */
  138 #define SIZEOF_FLT ((FLT_TOTAL_BIT + CHAR_BIT - 1) / CHAR_BIT)
  139 #define SIZEOF_DBL ((DBL_TOTAL_BIT + CHAR_BIT - 1) / CHAR_BIT)
  140 #define SIZEOF_LDBL ((LDBL_TOTAL_BIT + CHAR_BIT - 1) / CHAR_BIT)
  141 
  142 /* Verify that SIZEOF_FLT <= sizeof (float) etc.  */
  143 typedef int verify_sizeof_flt[SIZEOF_FLT <= sizeof (float) ? 1 : -1];
  144 typedef int verify_sizeof_dbl[SIZEOF_DBL <= sizeof (double) ? 1 : - 1];
  145 typedef int verify_sizeof_ldbl[SIZEOF_LDBL <= sizeof (long double) ? 1 : - 1];
  146 
  147 #endif /* _FLOATPLUS_H */