_C_STD_BEGIN
_MRTIMP2_NCEEPURE long double __CLRCALL_PURE_OR_CDECL _LCosh(long double x, long double y)
{ /* compute y * cosh(x), |y| <= 1 */
switch (_LDtest(&x))
{ /* test for special codes */
case _NANCODE:
case _INFCODE:
return (x);
case 0:
return (y);
default: /* finite */
if (y == 0.0L)
return (y);
if (x < 0.0)
x = -x;
if (x < _LXbig)
{ /* worth adding in exp(-x) */
_LExp(&x, 1.0L, -1);
return (y * (x + 0.25L / x));
}
switch (_LExp(&x, y, -1))
{ /* report over/underflow */
case 0:
_Feraise(_FE_UNDERFLOW);
break;
case _INFCODE:
_Feraise(_FE_OVERFLOW);
}
return (x);
}
}
_MRTIMP2_NCEEPURE long double __CLRCALL_PURE_OR_CDECL _LSinh(long double x, long double y)
{ /* compute y*sinh(x), |y| <= 1 */
short neg;
switch (_LDtest(&x))
{ /* test for special codes */
case _NANCODE:
return (x);
case _INFCODE:
return (y != 0.0L ? x : LSIGN(x) ? -y : y);
case 0:
return (x * y);
default: /* finite */
if (y == 0.0L)
return (x < 0.0L ? -y : y);
if (x < 0.0L)
x = -x, neg = 1;
else
neg = 0;
if (x < _LRteps._Long_double)
x *= y; /* x tiny */
else if (x < 1.0L)
{
long double w = x * x;
x += x * w * _LPoly(w, p, NP - 1);
x *= y;
}
else if (x < _LXbig)
{ /* worth adding in exp(-x) */
_LExp(&x, 1.0L, -1);
x = y * (x - 0.25L / x);
}
else
switch (_LExp(&x, y, -1))
{ /* report over/underflow */
case 0:
_Feraise(_FE_UNDERFLOW);
break;
case _INFCODE:
_Feraise(_FE_OVERFLOW);
}
return (neg ? -x : x);
}
}
_CRTIMP2 long double __cdecl _LSinh(long double x, long double y)
{ /* compute y*sinh(x), |y| <= 1 */
switch (_LDtest(&x))
{ /* test for special codes */
case NAN:
errno = EDOM;
return (x);
case INF:
if (y == 0)
return (0);
errno = ERANGE;
return (LSIGN(x) ? -_LInf._L : _LInf._L);
case 0:
return (0);
default: /* finite */
{ /* compute sinh(finite) */
short neg;
if (x < 0)
x = -x, neg = 1;
else
neg = 0;
if (x < _LRteps._L)
x *= y; /* x tiny */
else if (x < 1)
{
long double w = x * x;
x += x * w * _LPoly(w, p, NP - 1);
x *= y;
}
else if (x < _LXbig)
{ /* worth adding in exp(-x) */
_LExp(&x, 1, -1);
x = y * (x - 0.25 / x);
}
else if (0 <= _LExp(&x, y, -1))
errno = ERANGE; /* x large */
return (neg ? -x : x);
}
}
}