All,

I have the following problem with dgges. For version 3.1.1 and sooner, I get a reasonable result, for version 3.2 and 3.2.1 I get info=n+2.

I am separating eigenvalues in the unit circle from one outside the unique circle.

The two D and E matrices have relatively well separated null spaces, the minimum angle is acos(0.97). However, if i calculate condition numbers of E-lambda*D of lambda=[-1:0.01:1], they are quite bad.

The matrices are attached with a small c++ program which calls dgges and sorts eigenvalues for a better comparison. There is also a Makefile, which links with different version of lapack.

I use the reference blas.

My question is if it is a bug in Lapack introduced between 3.2. and 3.1.1 or the matrix is just very bad and in version 3.1.1 I was just lucky to get a reasonable solution.

Many thanks for any help.

Ondra Kamenik