1/**
 2 * This file has no copyright assigned and is placed in the Public Domain.
 3 * This file is part of the mingw-w64 runtime package.
 4 * No warranty is given; refer to the file DISCLAIMER.PD within this package.
 5 */
 6#include "cephes_mconf.h"
 7
 8#if __SIZEOF_LONG_DOUBLE__ == __SIZEOF_DOUBLE__
 9#include <math.h>
10
11long double coshl(long double x)
12{
13  return cosh(x);
14}
15#else
16
17#ifndef _SET_ERRNO
18#define _SET_ERRNO(x)
19#endif
20
21long double coshl(long double x)
22{
23  long double y;
24  int x_class = fpclassify (x);
25  if (x_class == FP_NAN)
26    {
27      errno = EDOM;
28      return x;
29    }
30  else if (x_class == FP_INFINITE)
31    {
32       errno = ERANGE;
33       return INFINITY;
34    }
35  x = fabsl (x);
36  if (x > (MAXLOGL + LOGE2L))
37    {
38      errno = ERANGE;
39#ifdef INFINITIES
40      return (INFINITYL);
41#else
42      return (MAXNUML);
43#endif
44    }
45  if (x >= (MAXLOGL - LOGE2L))
46    {
47      y = expl(0.5L * x);
48      y = (0.5L * y) * y;
49      return y;
50    }
51  y = expl(x);
52  y = 0.5L * (y + 1.0L / y);
53  return y;
54}
55#endif