asinh.hpp

Go to the documentation of this file.

#ifndef STAN_MATH_REV_FUN_ASINH_HPP
#define STAN_MATH_REV_FUN_ASINH_HPP
 
#include <stan/math/prim/core.hpp>
#include <stan/math/prim/meta.hpp>
#include <stan/math/prim/fun/asinh.hpp>
#include <stan/math/prim/fun/isinf.hpp>
#include <stan/math/prim/fun/is_inf.hpp>
#include <stan/math/rev/core.hpp>
#include <stan/math/rev/meta.hpp>
#include <stan/math/rev/fun/value_of_rec.hpp>
#include <stan/math/rev/fun/abs.hpp>
#include <stan/math/rev/fun/arg.hpp>
#include <stan/math/rev/fun/cosh.hpp>
#include <stan/math/rev/fun/is_inf.hpp>
#include <stan/math/rev/fun/is_nan.hpp>
#include <stan/math/rev/fun/log.hpp>
#include <stan/math/rev/fun/polar.hpp>
#include <stan/math/rev/fun/sqrt.hpp>
#include <cmath>
#include <complex>
 
namespace stan {
namespace math {
 
inline var asinh(const var& x) {
  return make_callback_var(std::asinh(x.val()), [x](const auto& vi) mutable {
    x.adj() += vi.adj() / std::sqrt(x.val() * x.val() + 1.0);
  });
}
 
template <typename VarMat, require_var_matrix_t<VarMat>* = nullptr>
inline auto asinh(const VarMat& x) {
  return make_callback_var(
      x.val().unaryExpr([](const auto x) { return asinh(x); }),
      [x](const auto& vi) mutable {
        x.adj().array()
            += vi.adj().array() / (x.val().array().square() + 1.0).sqrt();
      });
}
 
inline std::complex<var> asinh(const std::complex<var>& z) {
  return stan::math::internal::complex_asinh(z);
}
 
}  // namespace math
}  // namespace stan
#endif