OneScience-Group
/

AlphaFold3

#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);
}

	#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);
	}