search

Macaulay2 / M2

The primary source code repository for Macaulay2, a system for computing in commutative algebra, algebraic geometry and related fields.
https://macaulay2.com
345 stars 230 forks source link

Generalized eigenvalue problem #976

Open antonleykin opened 5 years ago

antonleykin commented 5 years ago

@JoseMath and @ltrbuse are interested in having a method to find generalized eigen- values/vectors.

Either wait until eigen can be used in M2 or wrap lapack routines: dggev zggev

mahrud commented 4 years ago

Is it beneficial to replace all lapack routines with eigen counterparts?

DanGrayson commented 4 years ago

Probably, since it's siad to be faster, by virtue of using vector operations on certain architectures. But I don't know.

We might also be able to use https://www.gnu.org/software/gsl/.

mikestillman commented 4 years ago

We should definitely do so time comparisons between eigen and lapack.

antonleykin commented 4 years ago

Perhaps we should ask @ofloveandhate for an opinion.

ofloveandhate commented 4 years ago

I think you can easily move entirely to Eigen -- because you can configure Eigen to use BLAS and Lapack routines . This would make doing timings easy, because it would only be setting whether Eigen uses Lapack or not, and the rest of the code would be identical.

Does this help?