Skip to content

ENH: Improve polynomial equations and analytical eigendecompositions

Kutalmış Berçin requested to merge feature-analytic-eigen into develop


Analytical eigendecomposition routines show irregular numerical instabilities particularly for cases wherein off-diagonal elements of a tensor involve noise. The issue revealed to be caused by the polynomialEqn containers, particularly cubicEqn.

In addition, the analytical eigendecomposition routines seem to be mathematically inconsistent since tensors are not allowed to return complex types, which is nevertheless the norm for asymmetric matrices.

Accordingly, this set of commits aims to improve the numerical stability of analytical eigendecompositions, and to provide a group of verification tests to prevent future breaks in the workflow.

Resolved bugs

#1311 (closed) #1312 (closed) #1527 #1575 #1596



polynomialEqns and analytical eigendecomposition routines were considerably isolated from the rest of the code.


Edited by Kutalmış Berçin

Merge request reports