math/test-iseqsig-excess-precision.c

   1 /* Test iseqsig with excess precision.
   2    Copyright (C) 2016-2023 Free Software Foundation, Inc.
   3    This file is part of the GNU C Library.
   4
   5    The GNU C Library is free software; you can redistribute it and/or
   6    modify it under the terms of the GNU Lesser General Public
   7    License as published by the Free Software Foundation; either
   8    version 2.1 of the License, or (at your option) any later version.
   9
  10    The GNU C Library 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 GNU
  13    Lesser General Public License for more details.
  14
  15    You should have received a copy of the GNU Lesser General Public
  16    License along with the GNU C Library; if not, see
  17    <https://www.gnu.org/licenses/>.  */
  18
  19 #include <float.h>
  20 #include <math.h>
  21 #include <stdio.h>
  22
  23 static int
  24 do_test (void)
  25 {
  26   int result = 0;
  27
  28   if (FLT_EVAL_METHOD == 1 || FLT_EVAL_METHOD == 2 || FLT_EVAL_METHOD > 32)
  29     {
  30       /* Excess precision for float.  */
  31       if (iseqsig (1.0f, 1.0f + (float) DBL_EPSILON))
  32         {
  33           puts ("iseqsig removes excess precision float -> double");
  34           result = 1;
  35         }
  36       else
  37         puts ("iseqsig preserves excess precision float -> double");
  38       if (iseqsig (__builtin_inff (), FLT_MAX * FLT_MAX))
  39         {
  40           puts ("iseqsig removes excess range float -> double");
  41           result = 1;
  42         }
  43       else
  44         puts ("iseqsig preserves excess range float -> double");
  45     }
  46
  47   if (FLT_EVAL_METHOD == 2 || FLT_EVAL_METHOD > 64)
  48     {
  49       /* Excess precision for float and double.  */
  50       if (iseqsig (1.0f, 1.0f + (float) LDBL_EPSILON))
  51         {
  52           puts ("iseqsig removes excess precision float -> long double");
  53           result = 1;
  54         }
  55       else
  56         puts ("iseqsig preserves excess precision float -> long double");
  57       if (iseqsig (1.0, 1.0 + (double) LDBL_EPSILON))
  58         {
  59           puts ("iseqsig removes excess precision double -> long double");
  60           result = 1;
  61         }
  62       else
  63         puts ("iseqsig preserves excess precision double -> long double");
  64       if (LDBL_MAX_EXP >= 2 * DBL_MAX_EXP)
  65         {
  66           if (iseqsig (__builtin_inf (), DBL_MAX * DBL_MAX))
  67             {
  68               puts ("iseqsig removes excess range double -> long double");
  69               result = 1;
  70             }
  71             else
  72               puts ("iseqsig preserves excess range double -> long double");
  73         }
  74     }
  75
  76   return result;
  77 }
  78
  79 #define TEST_FUNCTION do_test ()
  80 #include "../test-skeleton.c"