If I use this model to solve the eigenvalue problem with the qz solver, I get the independent eigenvalues of the blocks, namely, the eigenvalues of (K1,M1) and (K2,M2)? The coupling blocks seem to be not taken into consideration. What could be the cause of this?
What should be the path to use the QZ nonsymmetric eigenvalue solver?
Following the documentation, I guess I should do,
+ balance to improve for performance: ggbal + then reduced to hessenberg format: gghrd + use the solver: hgeqz + recover: ggbak
Is this correct?
If this is correct, then there is a problem because the problem gives me eigenvalues of the blocks only as explained.
Usually it is much easier to help if the test is provided. From the first glance, the sequence of calls you refer to looks OK. But at the same time, describing the matrices you draw them as block upper triangular and block lower triangular. But for the QZ solver one of matrices should be in the upper Hessenberg form, the second - in the upper triangular form. This might be the root cause of the issue you observe.