I have a generalized eigenvalue problem Ax = lambda*Bx, where A and B are real multiprecision matrices. I used eig(A,B) and expected that non-real eigenvalues came out in perfectly complex conjugate pairs (because I want to group them in pairs), but this is not the case. I'm not sure this should be considered a bug, or rather a question about whether there is a way to enforce that EIG returns perfectly conjugate pairs when A and B are real mp matrices?
Here is a simple example (tested in MATLAB 2015A and MCT version 18.104.22.16843):
H = [ 2.285927555469414e-03 -9.683344516608315e-03 0 0 ;
Customer support service by UserEcho