-
Notifications
You must be signed in to change notification settings - Fork 11
/
arith.h
39 lines (33 loc) · 1022 Bytes
/
arith.h
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
#ifndef F2C_ARITH_H
#define F2C_ARITH_H
#include <config.h>
#include <math.h>
#ifdef _MSC_VER
#define isnan _isnan
#define isinf(x) (!_finite(x))
#endif
#ifndef isnan
# define isnan(x) \
(sizeof (x) == sizeof (long double) ? isnan_ld (x) \
: sizeof (x) == sizeof (double) ? isnan_d (x) \
: isnan_f (x))
static inline int isnan_f (float x) { return x != x; }
static inline int isnan_d (double x) { return x != x; }
static inline int isnan_ld (long double x) { return x != x; }
#endif
#ifndef isinf
# define isinf(x) \
(sizeof (x) == sizeof (long double) ? isinf_ld (x) \
: sizeof (x) == sizeof (double) ? isinf_d (x) \
: isinf_f (x))
static inline int isinf_f (float x)
{ return !isnan (x) && isnan (x - x); }
static inline int isinf_d (double x)
{ return !isnan (x) && isnan (x - x); }
static inline int isinf_ld (long double x)
{ return !isnan (x) && isnan (x - x); }
#endif
#ifndef signbit
#define signbit(x) (((x) < 0)? 1 : 0)
#endif
#endif