cov_matrix_constrain_lkj.hpp

Go to the documentation of this file.

#ifndef STAN_MATH_PRIM_CONSTRAINT_COV_MATRIX_CONSTRAIN_LKJ_HPP
#define STAN_MATH_PRIM_CONSTRAINT_COV_MATRIX_CONSTRAIN_LKJ_HPP
 
#include <stan/math/prim/fun/Eigen.hpp>
#include <stan/math/prim/constraint/corr_constrain.hpp>
#include <stan/math/prim/constraint/positive_constrain.hpp>
#include <stan/math/prim/fun/read_cov_matrix.hpp>
#include <stan/math/prim/fun/to_ref.hpp>
 
namespace stan {
namespace math {
 
template <typename T, require_eigen_vector_t<T>* = nullptr>
inline Eigen::Matrix<value_type_t<T>, Eigen::Dynamic, Eigen::Dynamic>
cov_matrix_constrain_lkj(const T& x, size_t k) {
  size_t k_choose_2 = (k * (k - 1)) / 2;
  const auto& x_ref = to_ref(x);
  return read_cov_matrix(corr_constrain(x_ref.head(k_choose_2)),
                         positive_constrain(x_ref.tail(k)));
}
 
template <typename T, require_eigen_vector_t<T>* = nullptr>
inline Eigen::Matrix<value_type_t<T>, Eigen::Dynamic, Eigen::Dynamic>
cov_matrix_constrain_lkj(const T& x, size_t k, return_type_t<T>& lp) {
  size_t k_choose_2 = (k * (k - 1)) / 2;
  const auto& x_ref = x;
  return read_cov_matrix(corr_constrain(x_ref.head(k_choose_2)),
                         positive_constrain(x_ref.tail(k)), lp);
}
 
template <bool Jacobian, typename T, require_not_std_vector_t<T>* = nullptr>
inline auto cov_matrix_constrain_lkj(const T& x, size_t k,
                                     return_type_t<T>& lp) {
  if (Jacobian) {
    return cov_matrix_constrain_lkj(x, k, lp);
  } else {
    return cov_matrix_constrain_lkj(x, k);
  }
}
 
template <bool Jacobian, typename T, require_std_vector_t<T>* = nullptr>
inline auto cov_matrix_constrain_lkj(const T& x, size_t k,
                                     return_type_t<T>& lp) {
  return apply_vector_unary<T>::apply(x, [&lp, k](auto&& v) {
    return cov_matrix_constrain_lkj<Jacobian>(v, k, lp);
  });
}
 
}  // namespace math
}  // namespace stan
 
#endif