Follow-up Comment #1, bug #43015 (project gsl):
The problem is most likely due to these lines of eigen/test.c:
1270 gsl_eigen_nonsymmv_sort(w->eval, NULL, GSL_EIGEN_SORT_ABS_ASC);
1271 gsl_eigen_nonsymmv_sort(w->evalv, NULL, GSL_EIGEN_SORT_ABS_ASC);
1272 test_eigenvalues_complex(w->evalv, w->eval, "gen", desc);
With the extra optimization, the magnitudes of the eigenvalues could be
slightly different leading to a different sort result. Perhaps the solution is
to sort first by real part, then by imaginary part
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?43015>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/
