summaryrefslogtreecommitdiff
path: root/src/math/ldexp.c
blob: c56fe40508cd36f5975f0242c274309e2524573d (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
# define TGSOURCE "ldexp.c"
#include <errno.h>
#include <math.h>
#include "_tgmath.h"

/** multiply by a power of 2 **/

TYPE TGFN(ldexp)(TYPE x, int exp)
{
	switch (fpclassify(x)) {
	case FP_ZERO:           return x;
	case FP_INFINITE:       return x;
	default:                break;
	}

	if (exp == 0) {
		return x;
	}

	if (0) {
		errno = ERANGE; /* The result cannot be represented */
		/* RETURN_FAILURE(CONSTANT(HUGE_VAL), A range error occurred); */
		return TGHUGE;
	}

	/* RETURN_SUCCESS(ARGUMENT(x) * POW(2,ARGUMENT(exp))); */
	return x * TGFN(pow)(2, exp);
}

/***
multiply the floating-point number specifed by ARGUMENT(x)
by 2 raised to the power ARGUMENT(exp).
***/

/*
IMPLEMENTATION(The value returned on a domain error, CONSTANT(HUGE_VAL))
LINK(m)
STDC(1)
*/