gsl_eigen_nonsymmv_workspace has no member n_evals
issue: Diagonalizing double data_3 [] = { 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 } ; I get eigenvalue = 0 +0i eigenvector = 1 +0i 0 +0i 0 +0i eigenvalue = 0 +0i eigenvector = 0 +0i 1 +0i 0 +0i eigenvalue = 0 +0i eigenvector = -1 +0i 0 +0i 3.00625e-292 +0i which is wrong. The last eigenvector is not correct because this matrix is not diagonalizable. I need to identify such matrices. =========================================================================== Patrick DUPRÉ | | email: pdu...@gmx.com Laboratoire interdisciplinaire Carnot de Bourgogne 9 Avenue Alain Savary, BP 47870, 21078 DIJON Cedex FRANCE Tel: +33 (0)380395988 =========================================================================== > Sent: Tuesday, January 19, 2021 at 6:56 PM > From: "Patrick Alken" <patrick.al...@colorado.edu> > To: help-gsl@gnu.org > Subject: Re: eigensystem > > What do you mean by handle it? According to the documentation, if the > function cannot compute all eigenvalues, an error code is returned. In > the case of gsl_eigen_nonsymm, the number of converged eigenvalues is > stored in w->n_evals. > > Patrick > > On 1/19/21 10:33 AM, Patrick Dupre wrote: > > Hello, > > > > Is there a way to handle the possible error of gsl_eigen_nonsymmv ? > > > > For example, when the matrix is not diagonalizable. > > > > Thanks > > > > =========================================================================== > > Patrick DUPRÉ | | email: pdu...@gmx.com > > Laboratoire interdisciplinaire Carnot de Bourgogne > > 9 Avenue Alain Savary, BP 47870, 21078 DIJON Cedex FRANCE > > Tel: +33 (0)380395988 > > =========================================================================== > > > > > > >