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.