llt.cpp 298 B

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