summaryrefslogtreecommitdiff
path: root/support/dktools/dk4math.c
blob: 55cddeb8bb54518c719d88aa9bfe3e6a6d948e9f (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
/*
Copyright (C) 2017-2020, Dirk Krause
SPDX-License-Identifier: BSD-3-Clause
*/

/*
	WARNING: This file was generated by the dkct program (see
	http://dktools.sourceforge.net/ for details).
	Changes you make here will be lost if dkct is run again!
	You should modify the original source and run dkct on it.
	Original source: dk4math.ctr
*/

/**	@file dk4math.c The dk4math module.
*/


#line 243 "dk4math.ctr"

#include "dk4conf.h"
#include "dk4math.h"



#if (!(DK4_HAVE_RINT))
double
dk4ma_rint(double x)
{
	return (floor(x+0.5));
}
#endif



#if	(!(DK4_HAVE_ISGREATER))
int
dk4ma_isgreater(double a, double b)
{
	int		back = 0;
#if DK4_HAVE_ISFINITE
	if ((isfinite(a)) && (isfinite(b))) {
#else
#if	DK4_HAVE__FINITE
	if ((_finite(a)) && (_finite(b))) {
#endif
#endif
		if (a > b) {
			back = 1;
		}
#if DK4_HAVE_ISFINITE
	}
#else
#if	DK4_HAVE__FINITE
	}
#endif
#endif
	return back;
}
#endif



#if	(!(DK4_HAVE_ISGREATEREQUAL))
int
dk4ma_isgreaterequal(double a, double b)
{
	int		back = 0;
#if DK4_HAVE_ISFINITE
	if ((isfinite(a)) && (isfinite(b))) {
#else
#if	DK4_HAVE__FINITE
	if ((_finite(a)) && (_finite(b))) {
#endif
#endif
		if (a >= b) {
			back = 1;
		}
#if DK4_HAVE_ISFINITE
	}
#else
#if	DK4_HAVE__FINITE
	}
#endif
#endif
	return back;
}
#endif



#if (!(DK4_HAVE_ISLESS))
int
dk4ma_isless(double a, double b)
{
	int		back = 0;
#if DK4_HAVE_ISFINITE
	if ((isfinite(a)) && (isfinite(b))) {
#else
#if	DK4_HAVE__FINITE
	if ((_finite(a)) && (_finite(b))) {
#endif
#endif
		if (a < b) {
			back = 1;
		}
#if DK4_HAVE_ISFINITE
	}
#else
#if	DK4_HAVE__FINITE
	}
#endif
#endif
	return back;
}
#endif



#if	(!(DK4_HAVE_ISLESSEQUAL))
int
dk4ma_islessequal(double a, double b)
{
	int		back = 0;
#if DK4_HAVE_ISFINITE
	if ((isfinite(a)) && (isfinite(b))) {
#else
#if	DK4_HAVE__FINITE
	if ((_finite(a)) && (_finite(b))) {
#endif
#endif
		if (a <= b) {
			back = 1;
		}
#if DK4_HAVE_ISFINITE
	}
#else
#if	DK4_HAVE__FINITE
	}
#endif
#endif
	return back;
}
#endif



int
dk4ma_is_finite(double d) 
{
  int		 back	= 0;
#if _WIN32
  back = ((0 != _finite(d)) ? 1 : 0);
#else
  back = ((0 != isfinite(d)) ? 1 : 0);
#endif
  return back;
}



int
dk4ma_is_normal(double d)
{
  int		 back = 0;
#if _WIN32
  switch (_fpclass(d)) {
    case _FPCLASS_NN : case _FPCLASS_NZ :
    case _FPCLASS_PZ : case _FPCLASS_PN :
    {
      back = 1;
    } break;
  }
#else
  back = ((0 != isnormal(d)) ? 1 : 0);
#endif
  return back;
}



/* vim: set ai sw=4 ts=4 : */