blob: 101e5186cc7c4b8042e5c04ddfa9d8c3a38bb88e (
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
|
# define TGSOURCE "log1p.c"
#include "_tgmath.h"
#include <math.h>
#include <fenv.h>
TYPE TGFN(log1p)(TYPE x)
{
switch (fpclassify(x)) {
case FP_ZERO: return x;
case FP_INFINITE: if (!signbit(x)) { return x; } break;
default: break;
}
if (x == -1) {
feraiseexcept(FE_DIVBYZERO);
return - INFINITY;
}
if (x < -1) {
feraiseexcept(FE_INVALID);
return NAN;
}
return x;
}
/*
STDC(199901)
XOPEN(400)
LINK(m)
*/
|
