#include "gemm_common.h" | |
#include <Eigen/Cholesky> | |
EIGEN_DONT_INLINE void llt(const Mat &A, const Mat &B, Mat &C) { | |
C = A; | |
C.diagonal().array() += 1000; | |
Eigen::internal::llt_inplace<Mat::Scalar, Lower>::blocked(C); | |
} | |
int main(int argc, char **argv) { return main_gemm(argc, argv, llt); } |