divide.hpp

Go to the documentation of this file.

#ifndef STAN_MATH_PRIM_FUN_DIVIDE_HPP
#define STAN_MATH_PRIM_FUN_DIVIDE_HPP
 
#include <stan/math/prim/fun/Eigen.hpp>
#include <stan/math/prim/meta.hpp>
#include <stan/math/prim/err.hpp>
#include <stan/math/prim/fun/as_array_or_scalar.hpp>
#include <cstddef>
#include <cstdlib>
#include <type_traits>
 
namespace stan {
namespace math {
 
template <typename Scal1, typename Scal2,
          require_all_stan_scalar_t<Scal1, Scal2>* = nullptr>
inline return_type_t<Scal1, Scal2> divide(const Scal1& x, const Scal2& y) {
  return x / y;
}
 
inline int divide(int x, int y) {
  if (unlikely(y == 0)) {
    throw_domain_error("divide", "denominator is", y, "");
  }
  return x / y;
}
 
template <typename T1, typename T2, require_any_eigen_t<T1, T2>* = nullptr,
          require_all_not_st_var<T1, T2>* = nullptr>
inline auto divide(T1&& m, T2&& c) {
  return make_holder(
      [](auto&& m_, auto&& c_) {
        return (as_array_or_scalar(std::forward<decltype(m_)>(m_))
                / as_array_or_scalar(std::forward<decltype(c_)>(c_)))
            .matrix();
      },
      std::forward<T1>(m), std::forward<T2>(c));
}
 
}  // namespace math
}  // namespace stan
 
#endif