blob: d55b7d80340a7a07c1090da6f82bb97b515b4842 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
|
#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);
}
|
