You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
16 lines
298 B
16 lines
298 B
2 years ago
|
#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);
|
||
|
|
}
|