File: llt.cpp

package info (click to toggle)
eigen3 3.4.0-5
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 20,492 kB
  • sloc: cpp: 151,572; ansic: 75,396; fortran: 24,137; sh: 971; python: 244; javascript: 205; makefile: 42
file content (15 lines) | stat: -rw-r--r-- 298 bytes parent folder | download | duplicates (8) 
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);
}