1 #ifndef STAN_MATH_PRIM_MAT_FUN_QUAD_FORM_HPP 2 #define STAN_MATH_PRIM_MAT_FUN_QUAD_FORM_HPP 13 template <
int RA,
int CA,
int RB,
int CB,
typename T>
14 inline Eigen::Matrix<T, CB, CB>
quad_form(
const Eigen::Matrix<T, RA, CA>& A,
15 const Eigen::Matrix<T, RB, CB>& B) {
18 return B.transpose() * A * B;
21 template <
int RA,
int CA,
int RB,
typename T>
22 inline T
quad_form(
const Eigen::Matrix<T, RA, CA>& A,
23 const Eigen::Matrix<T, RB, 1>& 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.
Eigen::Matrix< T, CB, CB > quad_form(const Eigen::Matrix< T, RA, CA > &A, const Eigen::Matrix< T, RB, CB > &B)
Compute B^T A B.