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

◆ eigenvalues() [4/4]

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

Definition at line 442 of file eigen.cpp.

442 {
443 auto old_prec = Real::get_default_prec();
444 Real::set_default_prec(A->ring().get_precision());
445
446 MatrixXmpRRR AXmp(A->numRows(), A->numColumns());
447 fill_to_MatrixXmp(*A, AXmp);
448
449 Eigen::EigenSolver<MatrixXmpRRR> es(AXmp,false);
450 fill_from_MatrixXmp(es.eigenvalues(), *eigenvals);
451
452 Real::set_default_prec(old_prec);
453 return true;
454}
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-RRR.hpp:82
Eigen::Matrix< Real, Eigen::Dynamic, Eigen::Dynamic > MatrixXmpRRR
Definition eigen.cpp:25
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::ARingRRR::get_precision(), DMat< ACoeffRing >::numColumns(), DMat< ACoeffRing >::numRows(), and DMat< ACoeffRing >::ring().

Referenced by MatrixOps::eigenvalues(), MatrixOps::eigenvalues(), MatrixOps::eigenvalues(), and MatrixOps::eigenvalues().