Macaulay2 Engine
Loading...
Searching...
No Matches

◆ eigenvalues() [2/4]

bool EigenM2::eigenvalues ( const LMatrixCCC * A,
LMatrixCCC * eigenvals )

Definition at line 456 of file eigen.cpp.

456 {
457 auto old_prec = Real::get_default_prec();
458 Real::set_default_prec(A->ring().get_precision());
459
460 MatrixXmpCCC AXmp(A->numRows(), A->numColumns());
461 fill_to_MatrixXmp(*A, AXmp);
462
463 Eigen::ComplexEigenSolver<MatrixXmpCCC> ces(AXmp, false);
464 fill_from_MatrixXmp(ces.eigenvalues(), *eigenvals);
465
466 Real::set_default_prec(old_prec);
467 return true;
468}
size_t numRows() const
Definition dmat.hpp:144
const ACoeffRing & ring() const
Definition dmat.hpp:143
size_t numColumns() const
Definition dmat.hpp:145
unsigned long get_precision() const
Definition aring-CCC.hpp:88
Eigen::Matrix< Complex, Eigen::Dynamic, Eigen::Dynamic > MatrixXmpCCC
Definition eigen.cpp:26
void fill_to_MatrixXmp(const LMatrixRRR &orig, MatrixXmpRRR &result)
Definition eigen.cpp:313
void fill_from_MatrixXmp(const MatrixXmpRRR &orig, LMatrixRRR &result)
Definition eigen.cpp:327

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