add alternative computing methods for complex sigma
maybe a better way to create a complex matrix
shift-and-invert mode for sparse matrices
correct typo
improved documentation
improved documentation