1 #ifndef STAN_MATH_PRIM_MAT_FUN_CHOLESKY_FACTOR_FREE_HPP
2 #define STAN_MATH_PRIM_MAT_FUN_CHOLESKY_FACTOR_FREE_HPP
23 Eigen::Matrix<T, Eigen::Dynamic, 1>
25 <T, Eigen::Dynamic, Eigen::Dynamic>& y) {
29 "y is not a Cholesky factor");
32 Eigen::Matrix<T, Eigen::Dynamic, 1> x((N * (N + 1)) / 2 + (M - N) * N);
35 for (
int m = 0; m < N; ++m) {
36 for (
int n = 0; n < m; ++n)
39 x(pos++) =
log(y(m, m));
42 for (
int m = N; m < M; ++m)
43 for (
int n = 0; n < N; ++n)
fvar< T > log(const fvar< T > &x)
Eigen::Matrix< T, Eigen::Dynamic, 1 > cholesky_factor_free(const Eigen::Matrix< T, Eigen::Dynamic, Eigen::Dynamic > &y)
Return the unconstrained vector of parameters correspdonding to the specified Cholesky factor...
void domain_error(const char *function, const char *name, const T &y, const char *msg1, const char *msg2)
Throw a domain error with a consistently formatted message.
bool check_cholesky_factor(const char *function, const char *name, const Eigen::Matrix< T_y, Dynamic, Dynamic > &y)
Return true if the specified matrix is a valid Cholesky factor.