1 #include "gemm_common.h" 2 #include <Eigen/Cholesky> 3 4 EIGEN_DONT_INLINE llt(const Mat & A,const Mat & B,Mat & C)5void llt(const Mat &A, const Mat &B, Mat &C) 6 { 7 C = A; 8 C.diagonal().array() += 1000; 9 Eigen::internal::llt_inplace<Mat::Scalar, Lower>::blocked(C); 10 } 11 main(int argc,char ** argv)12int main(int argc, char **argv) 13 { 14 return main_gemm(argc, argv, llt); 15 } 16