1 /* PSPP - a program for statistical analysis.
2 Copyright (C) 1997-9, 2000 Free Software Foundation, Inc.
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.
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.
14 You should have received a copy of the GNU General Public License
15 along with this program. If not, see <http://www.gnu.org/licenses/>. */
23 #define EPSILON (10 * DBL_EPSILON)
25 /* HUGE_VAL is traditionally defined as positive infinity, or
26 alternatively, DBL_MAX. */
27 #if !HAVE_ISINF && !defined isinf
28 #define isinf(X) (fabs (X) == HUGE_VAL)
31 /* A Not a Number is not equal to itself. */
32 #if !HAVE_ISNAN && !defined isnan
33 #define isnan(X) ((X) != (X))
36 /* Finite numbers are not infinities or NaNs. */
37 #if !HAVE_FINITE && !defined finite
38 #define finite(X) (!isinf (X) && !isnan (X))
40 #include <ieeefp.h> /* Declares finite() under Solaris. */
43 /* Divides nonnegative X by positive Y, rounding up. */
44 #define DIV_RND_UP(X, Y) (((X) + ((Y) - 1)) / (Y))
46 /* Returns nonnegative difference between {nonnegative X} and {the
47 least multiple of positive Y greater than or equal to X}. */
48 #define REM_RND_UP(X, Y) ((X) % (Y) ? (Y) - (X) % (Y) : 0)
50 /* Rounds X up to the next multiple of Y. */
51 #define ROUND_UP(X, Y) (((X) + ((Y) - 1)) / (Y) * (Y))
53 /* Rounds X down to the previous multiple of Y. */
54 #define ROUND_DOWN(X, Y) ((X) / (Y) * (Y))
56 int intlog10 (unsigned);
58 /* Returns the square of X. */
65 /* Returns the cube of X. */
72 /* Returns the fourth power of X. */
81 #endif /* math/misc.h */