1  /* Test of <math.h> substitute.
       2     Copyright (C) 2007-2022 Free Software Foundation, Inc.
       3  
       4     This program is free software: you can redistribute it and/or modify
       5     it under the terms of the GNU General Public License as published by
       6     the Free Software Foundation, either version 3 of the License, or
       7     (at your option) any later version.
       8  
       9     This program is distributed in the hope that it will be useful,
      10     but WITHOUT ANY WARRANTY; without even the implied warranty of
      11     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
      12     GNU General Public License for more details.
      13  
      14     You should have received a copy of the GNU General Public License
      15     along with this program.  If not, see <https://www.gnu.org/licenses/>.  */
      16  
      17  /* Written by Bruno Haible <bruno@clisp.org>, 2007.  */
      18  
      19  #include <config.h>
      20  
      21  #include <math.h>
      22  
      23  #ifndef NAN
      24  # error NAN should be defined
      25  choke me
      26  #endif
      27  
      28  #ifndef HUGE_VALF
      29  # error HUGE_VALF should be defined
      30  choke me
      31  #endif
      32  
      33  #ifndef HUGE_VAL
      34  # error HUGE_VAL should be defined
      35  choke me
      36  #endif
      37  
      38  #ifndef HUGE_VALL
      39  # error HUGE_VALL should be defined
      40  choke me
      41  #endif
      42  
      43  #ifndef FP_ILOGB0
      44  # error FP_ILOGB0 should be defined
      45  choke me
      46  #endif
      47  
      48  #ifndef FP_ILOGBNAN
      49  # error FP_ILOGBNAN should be defined
      50  choke me
      51  #endif
      52  
      53  #include <limits.h>
      54  
      55  #include "macros.h"
      56  
      57  #if 0
      58  /* Check that NAN expands into a constant expression.  */
      59  static float n = NAN;
      60  #endif
      61  
      62  /* Compare two numbers with ==.
      63     This is a separate function because IRIX 6.5 "cc -O" miscompiles an
      64     'x == x' test.  */
      65  static int
      66  numeric_equalf (float x, float y)
      67  {
      68    return x == y;
      69  }
      70  static int
      71  numeric_equald (double x, double y)
      72  {
      73    return x == y;
      74  }
      75  static int
      76  numeric_equall (long double x, long double y)
      77  {
      78    return x == y;
      79  }
      80  
      81  int
      82  main (void)
      83  {
      84    double d = NAN;
      85    double zero = 0.0;
      86    ASSERT (!numeric_equald (d, d));
      87  
      88    d = HUGE_VAL;
      89    ASSERT (numeric_equald (d, 1.0 / zero));
      90  
      91    ASSERT (numeric_equalf (HUGE_VALF, HUGE_VALF + HUGE_VALF));
      92  
      93    ASSERT (numeric_equald (HUGE_VAL, HUGE_VAL + HUGE_VAL));
      94  
      95    ASSERT (numeric_equall (HUGE_VALL, HUGE_VALL + HUGE_VALL));
      96  
      97    /* Check the value of FP_ILOGB0.  */
      98    ASSERT (FP_ILOGB0 == INT_MIN || FP_ILOGB0 == - INT_MAX);
      99  
     100    /* Check the value of FP_ILOGBNAN.  */
     101    ASSERT (FP_ILOGBNAN == INT_MIN || FP_ILOGBNAN == INT_MAX);
     102  
     103    return 0;
     104  }