/* Test of rounding towards positive infinity.
Copyright (C) 2007-2008 Free Software Foundation, Inc.
This program is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation; either version 3 of the License, or
(at your option) any later version.
This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
You should have received a copy of the GNU General Public License
along with this program. If not, see . */
/* Written by Bruno Haible , 2007. */
#include
#include
#include
#include
#include "fpucw.h"
#include "isnanl-nolibm.h"
#define ASSERT(expr) \
do \
{ \
if (!(expr)) \
{ \
fprintf (stderr, "%s:%d: assertion failed\n", __FILE__, __LINE__); \
fflush (stderr); \
abort (); \
} \
} \
while (0)
/* HP cc on HP-UX 10.20 has a bug with the constant expression -0.0L.
So we use -zero instead. */
long double zero = 0.0L;
int
main ()
{
DECL_LONG_DOUBLE_ROUNDING
BEGIN_LONG_DOUBLE_ROUNDING ();
/* Zero. */
ASSERT (ceill (0.0L) == 0.0L);
ASSERT (ceill (-zero) == 0.0L);
/* Positive numbers. */
ASSERT (ceill (0.3L) == 1.0L);
ASSERT (ceill (0.7L) == 1.0L);
ASSERT (ceill (1.0L) == 1.0L);
ASSERT (ceill (1.001L) == 2.0L);
ASSERT (ceill (1.5L) == 2.0L);
ASSERT (ceill (1.999L) == 2.0L);
ASSERT (ceill (2.0L) == 2.0L);
ASSERT (ceill (65535.999L) == 65536.0L);
ASSERT (ceill (65536.0L) == 65536.0L);
ASSERT (ceill (2.341e31L) == 2.341e31L);
/* Negative numbers. */
ASSERT (ceill (-0.3L) == 0.0L);
ASSERT (ceill (-0.7L) == 0.0L);
ASSERT (ceill (-1.0L) == -1.0L);
ASSERT (ceill (-1.5L) == -1.0L);
ASSERT (ceill (-1.999L) == -1.0L);
ASSERT (ceill (-2.0L) == -2.0L);
ASSERT (ceill (-65535.999L) == -65535.0L);
ASSERT (ceill (-65536.0L) == -65536.0L);
ASSERT (ceill (-2.341e31L) == -2.341e31L);
/* Infinite numbers. */
ASSERT (ceill (1.0L / 0.0L) == 1.0L / 0.0L);
ASSERT (ceill (-1.0L / 0.0L) == -1.0L / 0.0L);
/* NaNs. */
ASSERT (isnanl (ceill (0.0L / 0.0L)));
return 0;
}