Subversion Repositories Kolibri OS

Rev

Blame | Last modification | View Log | Download | RSS feed

  1. /*                                                      coshl.c
  2.  *
  3.  *      Hyperbolic cosine, long double precision
  4.  *
  5.  *
  6.  *
  7.  * SYNOPSIS:
  8.  *
  9.  * long double x, y, coshl();
  10.  *
  11.  * y = coshl( x );
  12.  *
  13.  *
  14.  *
  15.  * DESCRIPTION:
  16.  *
  17.  * Returns hyperbolic cosine of argument in the range MINLOGL to
  18.  * MAXLOGL.
  19.  *
  20.  * cosh(x)  =  ( exp(x) + exp(-x) )/2.
  21.  *
  22.  *
  23.  *
  24.  * ACCURACY:
  25.  *
  26.  *                      Relative error:
  27.  * arithmetic   domain     # trials      peak         rms
  28.  *    IEEE     +-10000      30000       1.1e-19     2.8e-20
  29.  *
  30.  *
  31.  * ERROR MESSAGES:
  32.  *
  33.  *   message         condition              value returned
  34.  * cosh overflow    |x| > MAXLOGL+LOGE2L      INFINITYL
  35.  *
  36.  *
  37.  */
  39.  
  40. /*
  41. Cephes Math Library Release 2.7:  May, 1998
  42. Copyright 1985, 1991, 1998 by Stephen L. Moshier
  43. */
  44.  
  45. /*
  46. Modified for mingw
  47. 2002-07-22 Danny Smith <dannysmith@users.sourceforge.net>
  48. */
  49.  
  50. #ifdef __MINGW32__
  51. #include "cephes_mconf.h"
  52. #else
  53. #include "mconf.h"
  54. #endif
  55.  
  56. #ifndef _SET_ERRNO
  57. #define _SET_ERRNO(x)
  58. #endif
  59.  
  60.  
  61. #ifndef __MINGW32__
  62. extern long double MAXLOGL, MAXNUML, LOGE2L;
  63. #ifdef ANSIPROT
  64. extern long double expl ( long double );
  65. extern int isnanl ( long double );
  66. #else
  67. long double expl(), isnanl();
  68. #endif
  69. #ifdef INFINITIES
  70. extern long double INFINITYL;
  71. #endif
  72. #ifdef NANS
  73. extern long double NANL;
  74. #endif
  75. #endif /* __MINGW32__ */
  76.  
  77. long double coshl(x)
  78. long double x;
  79. {
  80. long double y;
  81.  
  82. #ifdef NANS
  83. if( isnanl(x) )
  84.         {
  85.         _SET_ERRNO(EDOM);
  86.         return(x);
  87.         }
  88. #endif
  89. if( x < 0 )
  90.         x = -x;
  91. if( x > (MAXLOGL + LOGE2L) )
  92.         {
  93.         mtherr( "coshl", OVERFLOW );
  94.         _SET_ERRNO(ERANGE);
  95. #ifdef INFINITIES
  96.         return( INFINITYL );
  97. #else
  98.         return( MAXNUML );
  99. #endif
  100.         }      
  101. if( x >= (MAXLOGL - LOGE2L) )
  102.         {
  103.         y = expl(0.5L * x);
  104.         y = (0.5L * y) * y;
  105.         return(y);
  106.         }
  107. y = expl(x);
  108. y = 0.5L * (y + 1.0L / y);
  109. return( y );
  110. }
  111.