I have been using the routine DGGEV to solve a generalised eigenvalue problem for the values E satisfying A.x = E B.x
The routine DGGEV works fine and gives me the right results for the system, however when I try to make the problem complex and use ZGGEV I get results that seem to bear no relation to those generated by DGGEV. The input has only changed in that some arguments have been assigned as complex16, where previously they were double precision, (and that ZGGEV drops an argument) however the values of these arguments contain no complex part at present and I would have thought that the outcome should consequently be the same.
In any case, any help would be much appreciated!