// -*-C++-*-

#ifndef MATHFUNCS_ASINH_H
#define MATHFUNCS_ASINH_H

#include "mathfuncs_base.h"

#include <cmath>

namespace vecmathlib {

template <typename realvec_t>
realvec_t mathfuncs<realvec_t>::vml_asinh(realvec_t x) {
  // Reduce range
  realvec_t r = fabs(x);
  r = log(r + sqrt(r * r + RV(1.0)));
  r = copysign(r, x);
  return r;
}

template <typename realvec_t>
realvec_t mathfuncs<realvec_t>::vml_acosh(realvec_t x) {
  return log(x + sqrt(x * x - RV(1.0)));
}

template <typename realvec_t>
realvec_t mathfuncs<realvec_t>::vml_atanh(realvec_t x) {
  // Reduce range
  realvec_t r = fabs(x);
  r = RV(0.5) * log((RV(1.0) + r) / (RV(1.0) - r));
  r = copysign(r, x);
  return r;
}

}; // namespace vecmathlib

#endif // #ifndef MATHFUNCS_ASINH_H