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