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#include <errno.h>
  8
  9#if __SIZEOF_LONG_DOUBLE__ == __SIZEOF_DOUBLE__
 10#include <math.h>
 11
 12long double sinhl(long double x)
 13{
 14  return sinh(x);
 15}
 16#else
 17
 18#ifdef UNK
 19static uLD P[] = {
 20  { { 1.7550769032975377032681E-6L } },
 21  { { 4.1680702175874268714539E-4L } },
 22  { { 3.0993532520425419002409E-2L } },
 23  { { 9.9999999999999999998002E-1L } }
 24};
 25static long double Q[] = {
 26  { { 1.7453965448620151484660E-8L } },
 27  { { -5.9116673682651952419571E-6L } },
 28  { { 1.0599252315677389339530E-3L } },
 29  { { -1.1403880487744749056675E-1L } },
 30  { { 6.0000000000000000000200E0L } }
 31};
 32#endif
 33
 34#ifdef IBMPC
 35static const uLD P[] = {
 36  { { 0xec6a,0xd942,0xfbb3,0xeb8f,0x3feb, 0, 0, 0 } },
 37  { { 0x365e,0xb30a,0xe437,0xda86,0x3ff3, 0, 0, 0 } },
 38  { { 0x8890,0x01f6,0x2612,0xfde6,0x3ff9, 0, 0, 0 } },
 39  { { 0x0000,0x0000,0x0000,0x8000,0x3fff, 0, 0, 0 } }
 40};
 41static const uLD Q[] = {
 42  { { 0x4edd,0x4c21,0xad09,0x95ed,0x3fe5, 0, 0, 0 } },
 43  { { 0x4376,0x9b70,0xd605,0xc65c,0xbfed, 0, 0, 0 } },
 44  { { 0xc8ad,0x5d21,0x3069,0x8aed,0x3ff5, 0, 0, 0 } },
 45  { { 0x9c32,0x6374,0x2d4b,0xe98d,0xbffb, 0, 0, 0 } },
 46  { { 0x0000,0x0000,0x0000,0xc000,0x4001, 0, 0, 0 } }
 47};
 48#endif
 49
 50#ifdef MIEEE
 51static uLD P[] = {
 52  { { 0x3feb0000,0xeb8ffbb3,0xd942ec6a, 0 } },
 53  { { 0x3ff30000,0xda86e437,0xb30a365e, 0 } },
 54  { { 0x3ff90000,0xfde62612,0x01f68890, 0 } },
 55  { { 0x3fff0000,0x80000000,0x00000000, 0 } }
 56};
 57static uLD Q[] = {
 58  { { 0x3fe50000,0x95edad09,0x4c214edd, 0 } },
 59  { { 0xbfed0000,0xc65cd605,0x9b704376, 0 } },
 60  { { 0x3ff50000,0x8aed3069,0x5d21c8ad, 0 } },
 61  { { 0xbffb0000,0xe98d2d4b,0x63749c32, 0 } },
 62  { { 0x40010000,0xc0000000,0x00000000, 0 } }
 63};
 64#endif
 65
 66long double sinhl(long double x)
 67{
 68  long double a;
 69  int x_class = fpclassify (x);
 70
 71  if (x_class == FP_NAN)
 72    {
 73      errno = EDOM;
 74      return x;
 75    }
 76  if (x_class == FP_ZERO)
 77    return x;
 78  if (x_class == FP_INFINITE ||
 79      (fabsl (x) > (MAXLOGL + LOGE2L)))
 80  {
 81    errno = ERANGE;
 82#ifdef INFINITIES
 83    return (signbit (x) ? -INFINITYL : INFINITYL);
 84#else
 85    return (signbit (x) ? -MAXNUML : MAXNUML);
 86#endif
 87  }
 88  a = fabsl (x);
 89  if (a > 1.0L)
 90  {
 91    if (a >= (MAXLOGL - LOGE2L))
 92    {
 93      a = expl(0.5L*a);
 94      a = (0.5L * a) * a;
 95      if (x < 0.0L)
 96	a = -a;
 97      return (a);
 98    }
 99    a = expl(a);
100    a = 0.5L*a - (0.5L/a);
101    if (x < 0.0L)
102      a = -a;
103    return (a);
104  }
105
106  a *= a;
107  return (x + x * a * (polevll(a,P,3)/polevll(a,Q,4)));
108}
109#endif