1 #ifndef STAN_MATH_PRIM_MAT_FUN_TRACE_QUAD_FORM_HPP
2 #define STAN_MATH_PRIM_MAT_FUN_TRACE_QUAD_FORM_HPP
4 #include <boost/utility/enable_if.hpp>
5 #include <boost/type_traits.hpp>
15 template<
int RA,
int CA,
int RB,
int CB>
18 const Eigen::Matrix<double, RB, CB> &B) {
24 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)
bool check_multiplicable(const char *function, const char *name1, const T1 &y1, const char *name2, const T2 &y2)
Return true if the matrices can be multiplied.
T trace(const Eigen::Matrix< T, Eigen::Dynamic, Eigen::Dynamic > &m)
Returns the trace of the specified matrix.
bool check_square(const char *function, const char *name, const Eigen::Matrix< T_y, Eigen::Dynamic, Eigen::Dynamic > &y)
Return true if the specified matrix is square.