Check for definiteness in generalised eigensolvers. Use I in eigs. Change default nev in eigs because different restrictions apply to symmetric and non-symmetric problems. Use factorize instead of lufact for inverse operations.
Merge pull request #6964 from JuliaLang/anj/geneig

Check for definiteness in generalised eigensolvers.