Subversion Repositories Kolibri OS

Rev

Go to most recent revision | Blame | Last modification | View Log | Download | RSS feed

  1. #include <math.h>
  2. #include <errno.h>
  3. #include "fastmath.h"
  4.  
  5. /* atanh (x) = 0.5 * log ((1.0 + x)/(1.0 - x)) */
  6. float atanhf (float x)
  7. {
  8.   float z;
  9.   if isnan (x)
  10.     return x;
  11.   z = fabsf (x);
  12.   if (z == 1.0)
  13.     {
  14.       errno  = ERANGE;
  15.       return (x > 0 ? INFINITY : -INFINITY);
  16.     }
  17.   if ( z > 1.0)
  18.     {
  19.       errno = EDOM;
  20.       return nanf("");
  21.     }
  22.   /* Rearrange formula to avoid precision loss for small x.
  23.  
  24.   atanh(x) = 0.5 * log ((1.0 + x)/(1.0 - x))
  25.            = 0.5 * log1p ((1.0 + x)/(1.0 - x) - 1.0)
  26.            = 0.5 * log1p ((1.0 + x - 1.0 + x) /(1.0 - x))
  27.            = 0.5 * log1p ((2.0 * x ) / (1.0 - x))  */
  28.   z = 0.5 * __fast_log1p ((z + z) / (1.0 - z));
  29.   return x >= 0 ? z : -z;
  30. }
  31.