tests/test-roundl-ieee.c

   1 /* Test of rounding to nearest, breaking ties away from zero.
   2    Copyright (C) 2010-2020 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, or (at your option)
   7    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 #include <config.h>
  18
  19 #include <math.h>
  20
  21 #include "fpucw.h"
  22 #include "isnanl-nolibm.h"
  23 #include "minus-zero.h"
  24 #include "infinity.h"
  25 #include "nan.h"
  26 #include "macros.h"
  27
  28 int
  29 main ()
  30 {
  31   DECL_LONG_DOUBLE_ROUNDING
  32
  33   BEGIN_LONG_DOUBLE_ROUNDING ();
  34
  35   /* See IEEE 754, section 6.3:
  36        "the sign of the result of the round floating-point number to
  37         integral value operation is the sign of the operand. These rules
  38         shall apply even when operands or results are zero or infinite."  */
  39
  40   /* Zero.  */
  41   ASSERT (!signbit (roundl (0.0L)));
  42   ASSERT (!!signbit (roundl (minus_zerol)) == !!signbit (minus_zerol));
  43   /* Positive numbers.  */
  44   ASSERT (!signbit (roundl (0.3L)));
  45   ASSERT (!signbit (roundl (0.7L)));
  46   /* Negative numbers.  */
  47   ASSERT (!!signbit (roundl (-0.3L)) == !!signbit (minus_zerol));
  48   ASSERT (!!signbit (roundl (-0.7L)) == !!signbit (minus_zerol));
  49
  50   /* [MX] shaded specification in POSIX.  */
  51
  52   /* NaN.  */
  53   ASSERT (isnanl (roundl (NaNl ())));
  54   /* Infinity.  */
  55   ASSERT (roundl (Infinityl ()) == Infinityl ());
  56   ASSERT (roundl (- Infinityl ()) == - Infinityl ());
  57
  58   return 0;
  59 }