1 #ifndef STAN_MATH_PRIM_MAT_FUN_TRACE_QUAD_FORM_HPP 2 #define STAN_MATH_PRIM_MAT_FUN_TRACE_QUAD_FORM_HPP 13 template <
int RA,
int CA,
int RB,
int CB>
15 const Eigen::Matrix<double, RB, CB> &B) {
19 return (B.transpose() * A * B).
trace();
fvar< T > trace_quad_form(const Eigen::Matrix< fvar< T >, RA, CA > &A, const Eigen::Matrix< fvar< T >, RB, CB > &B)
void check_square(const char *function, const char *name, const matrix_cl &y)
Check if the matrix_cl is square.
void check_multiplicable(const char *function, const char *name1, const T1 &y1, const char *name2, const T2 &y2)
Check if the matrices can be multiplied.
T trace(const Eigen::Matrix< T, Eigen::Dynamic, Eigen::Dynamic > &m)
Returns the trace of the specified matrix.