tests/test-roundl-ieee.c

   1 /* Test of rounding to nearest, breaking ties away from zero.
   2    Copyright (C) 2010 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, write to the Free Software Foundation,
  16    Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.  */
  17
  18 #include <config.h>
  19
  20 #include <math.h>
  21
  22 #include "fpucw.h"
  23 #include "minus-zero.h"
  24 #include "macros.h"
  25
  26 int
  27 main ()
  28 {
  29   DECL_LONG_DOUBLE_ROUNDING
  30
  31   BEGIN_LONG_DOUBLE_ROUNDING ();
  32
  33   /* See IEEE 754, section 6.3:
  34        "the sign of the result of the round floating-point number to
  35         integral value operation is the sign of the operand. These rules
  36         shall apply even when operands or results are zero or infinite."  */
  37
  38   /* Zero.  */
  39   ASSERT (!signbit (roundl (0.0L)));
  40   ASSERT (!!signbit (roundl (minus_zerol)) == !!signbit (minus_zerol));
  41   /* Positive numbers.  */
  42   ASSERT (!signbit (roundl (0.3L)));
  43   ASSERT (!signbit (roundl (0.7L)));
  44   /* Negative numbers.  */
  45   ASSERT (!!signbit (roundl (-0.3L)) == !!signbit (minus_zerol));
  46   ASSERT (!!signbit (roundl (-0.7L)) == !!signbit (minus_zerol));
  47
  48   return 0;
  49 }