least_squares() [2/4]

bool EigenM2::least_squares	(	const LMatrixCCC *	A,
		const LMatrixCCC *	B,
		LMatrixCCC *	X )

Definition at line 579 of file eigen.cpp.

{
  auto old_prec = Real::get_default_prec(); 
  Real::set_default_prec(A->ring().get_precision());
 
  MatrixXmpCCC AXmp(A->numRows(), A->numColumns());
  fill_to_MatrixXmp(*A, AXmp);
  MatrixXmpCCC BXmp(B->numRows(), B->numColumns());
  fill_to_MatrixXmp(*B, BXmp);
 
  Eigen::BDCSVD<MatrixXmpCCC> svd(AXmp, Eigen::ComputeThinU | Eigen::ComputeThinV);
  
  fill_from_MatrixXmp(svd.solve(BXmp), *X);
 
  Real::set_default_prec(old_prec);
  return true;
}

References fill_from_MatrixXmp(), fill_to_MatrixXmp(), M2::ARingCCC::get_precision(), DMat< ACoeffRing >::numColumns(), DMat< ACoeffRing >::numRows(), and DMat< ACoeffRing >::ring().