sysdeps/ieee754/ldbl-128ibm/w_expl.c


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24

#include <math.h>
#include <math_private.h>
#include <math_ldbl_opt.h>

static const long double o_thres = 709.78271289338399678773454114191496482L;
static const long double u_thres = -744.44007192138126231410729844608163411L;

long double __expl(long double x)	/* wrapper exp  */
{
  long double z;
  z = __ieee754_expl(x);
  if (_LIB_VERSION == _IEEE_)
    return z;
  if (isfinite(x))
    {
      if (x >= o_thres)
	return __kernel_standard_l(x,x,206); /* exp overflow  */
      else if (x <= u_thres)
	return __kernel_standard_l(x,x,207); /* exp underflow  */
    }
  return z;
}
hidden_def (__expl)
long_double_symbol (libm, __expl, expl);