1  /* Test for ldbl-128ibm fmodl etc. handling of equal values.
       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 <fenv.h>
      20  #include <float.h>
      21  #include <math.h>
      22  #include <stdio.h>
      23  
      24  /* FUNC is defined to be the name of the function to test.  */
      25  #define STRX(x) #x
      26  #define STR(x) STRX (x)
      27  #define SFUNC STR (FUNC)
      28  
      29  union u
      30  {
      31    long double ld;
      32    double d[2];
      33  };
      34  
      35  volatile union u p1 = { .d = { DBL_MIN, 0.0 } };
      36  volatile union u p2 = { .d = { DBL_MIN, -0.0 } };
      37  volatile union u m1 = { .d = { -DBL_MIN, 0.0 } };
      38  volatile union u m2 = { .d = { -DBL_MIN, -0.0 } };
      39  
      40  static int
      41  test_func (const char *s, long double x, long double y, long double expected)
      42  {
      43    volatile long double r;
      44    r = FUNC (x, y);
      45    if (r != expected || copysignl (1.0, r) != copysignl (1.0, expected))
      46      {
      47        printf ("FAIL: " SFUNC " (%s)\n", s);
      48        return 1;
      49      }
      50    else
      51      {
      52        printf ("PASS: " SFUNC " (%s)\n", s);
      53        return 0;
      54      }
      55  }
      56  
      57  #define TEST_FUNC(a, b, e) test_func (#a ", " #b, a, b, e)
      58  
      59  static int
      60  do_test (void)
      61  {
      62    int result = 0;
      63    SETUP;
      64    result |= TEST_FUNC (p1.ld, p1.ld, 0.0L);
      65    result |= TEST_FUNC (p1.ld, p2.ld, 0.0L);
      66    result |= TEST_FUNC (p1.ld, m1.ld, 0.0L);
      67    result |= TEST_FUNC (p1.ld, m2.ld, 0.0L);
      68    result |= TEST_FUNC (p2.ld, p1.ld, 0.0L);
      69    result |= TEST_FUNC (p2.ld, p2.ld, 0.0L);
      70    result |= TEST_FUNC (p2.ld, m1.ld, 0.0L);
      71    result |= TEST_FUNC (p2.ld, m2.ld, 0.0L);
      72    result |= TEST_FUNC (m1.ld, p1.ld, -0.0L);
      73    result |= TEST_FUNC (m1.ld, p2.ld, -0.0L);
      74    result |= TEST_FUNC (m1.ld, m1.ld, -0.0L);
      75    result |= TEST_FUNC (m1.ld, m2.ld, -0.0L);
      76    result |= TEST_FUNC (m2.ld, p1.ld, -0.0L);
      77    result |= TEST_FUNC (m2.ld, p2.ld, -0.0L);
      78    result |= TEST_FUNC (m2.ld, m1.ld, -0.0L);
      79    result |= TEST_FUNC (m2.ld, m2.ld, -0.0L);
      80    return result;
      81  }
      82  
      83  #define TEST_FUNCTION do_test ()
      84  #include "../../../test-skeleton.c"