GE Solver with shifts. Works for GEigsMode = Cholesky or RegularInverse. See the Spectra Documentation (SymGEigsSolver) for more information