src/libpspp/misc.h

   1 /* PSPP - a program for statistical analysis.
   2    Copyright (C) 1997-9, 2000 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 of the License, or
   7    (at your option) 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 <http://www.gnu.org/licenses/>. */
  16
  17 #if !math_misc_h
  18 #define math_misc_h 1
  19
  20 #include <float.h>
  21 #include <math.h>
  22
  23 #define EPSILON (10 * DBL_EPSILON)
  24
  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)
  29 #endif
  30
  31 /* A Not a Number is not equal to itself. */
  32 #if !HAVE_ISNAN && !defined isnan
  33 #define isnan(X) ((X) != (X))
  34 #endif
  35
  36 /* Finite numbers are not infinities or NaNs. */
  37 #if !HAVE_FINITE && !defined finite
  38 #define finite(X) (!isinf (X) && !isnan (X))
  39 #elif HAVE_IEEEFP_H
  40 #include <ieeefp.h>             /* Declares finite() under Solaris. */
  41 #endif
  42
  43 /* Divides nonnegative X by positive Y, rounding up. */
  44 #define DIV_RND_UP(X, Y) (((X) + ((Y) - 1)) / (Y))
  45
  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)
  49
  50 /* Rounds X up to the next multiple of Y. */
  51 #define ROUND_UP(X, Y) (((X) + ((Y) - 1)) / (Y) * (Y))
  52
  53 /* Rounds X down to the previous multiple of Y. */
  54 #define ROUND_DOWN(X, Y) ((X) / (Y) * (Y))
  55
  56 int intlog10 (unsigned);
  57
  58 /* Returns the square of X. */
  59 static inline double
  60 pow2 (double x)
  61 {
  62   return x * x;
  63 }
  64
  65 /* Returns the cube of X. */
  66 static inline double
  67 pow3 (double x)
  68 {
  69   return x * x * x;
  70 }
  71
  72 /* Returns the fourth power of X. */
  73 static inline double
  74 pow4 (double x)
  75 {
  76   double y = x * x;
  77   y *= y;
  78   return y;
  79 }
  80
  81 #endif /* math/misc.h */