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 !libpspp_misc_h
18 #define libpspp_misc_h 1
19
20 #include <stddef.h>
21 #include <float.h>
22 #include <math.h>
23
24 #define EPSILON (10 * DBL_EPSILON)
25
26 /* Divides nonnegative X by positive Y, rounding up. */
27 #define DIV_RND_UP(X, Y) (((X) + ((Y) - 1)) / (Y))
28
29 /* Returns nonnegative difference between {nonnegative X} and {the
30 least multiple of positive Y greater than or equal to X}. */
31 #define REM_RND_UP(X, Y) ((X) % (Y) ? (Y) - (X) % (Y) : 0)
32
33 /* Rounds X up to the next multiple of Y. */
34 #define ROUND_UP(X, Y) (((X) + ((Y) - 1)) / (Y) * (Y))
35
36 /* Rounds X down to the previous multiple of Y. */
37 #define ROUND_DOWN(X, Y) ((X) / (Y) * (Y))
38
39 int intlog10 (unsigned);
40
41 /* Returns the square of X. */
42 static inline double
pow2(double x)43 pow2 (double x)
44 {
45 return x * x;
46 }
47
48 /* Returns the cube of X. */
49 static inline double
pow3(double x)50 pow3 (double x)
51 {
52 return x * x * x;
53 }
54
55 /* Returns the fourth power of X. */
56 static inline double
pow4(double x)57 pow4 (double x)
58 {
59 double y = x * x;
60 y *= y;
61 return y;
62 }
63
64 /* Set *DEST to the lower of *DEST and SRC */
65 static inline void
minimize(double * dest,double src)66 minimize (double *dest, double src)
67 {
68 if (src < *dest)
69 *dest = src;
70 }
71
72
73 /* Set *DEST to the greater of *DEST and SRC */
74 static inline void
maximize(double * dest,double src)75 maximize (double *dest, double src)
76 {
77 if (src > *dest)
78 *dest = src;
79 }
80
81
82 /* Set *DEST to the lower of *DEST and SRC */
83 static inline void
minimize_int(int * dest,int src)84 minimize_int (int *dest, int src)
85 {
86 if (src < *dest)
87 *dest = src;
88 }
89
90
91 /* Set *DEST to the greater of *DEST and SRC */
92 static inline void
maximize_int(int * dest,int src)93 maximize_int (int *dest, int src)
94 {
95 if (src > *dest)
96 *dest = src;
97 }
98
99 int c_dtoastr (char *buf, size_t bufsize, int flags, int width, double x);
100
101
102 #endif /* libpspp/misc.h */
103