Algorithm for the smallest (algebraic) eigenvalues of a symmetric (sparse) matrix Hi,
I'm looking for a way to get the negative eigenspace of a large (sparse) symmetric matrix. This matrix is basically a discretized version of the operator $-\Delta + V$, $V$ negative, on some domain $[-L,L]$ with Dirichlet BC, so its spectrum consists of a few negative eigenvalues (which I want to find), and a lot of positive ones (whose distribution is roughly known). 
The way I currently do it is to use the shift-invert mode of ARPACK (so Lanczos), with a negative shift and 'LM' mode (lowest magnitude). This requires me to choose a good shift: too large a shift might miss negative eigenvalues, and too small a shift leads to slow convergence. The 'LA' mode (lowest algebraic) is just not an option, it's too slow/imprecise.
Is there any better method out there?
 A: If you are interested in software for doing this problem quickly, you might consider the SLEPc package which includes various sparse and parallel eigenvalue algorithms.  It is quite difficult to get set up but once you do you will be able to experiment with different methods easily.
A: The chebyshev-davidson method is very likely the method of choice for this type of problems. Since you are interested in the negative eigenvalues, and the eigenvalues are known to be inside  [-L, L], you can input the bounds [0, L] to the chebyshev filters at each davidson iteration step, so that eigenvalues in [0, L] will be mapped into [-1, 1] for damping, while the eigenvalues you are interested in will be magnified (since they are mapped outside [-1, 1]) and converged quickly. The method is an accelerated version of Lanczos, with an advantage that it requires much smaller memory than krylov methods (memory efficiency is important when the number of eigenvalues to be computed is large).  
A: Arpack with which = 'SR'?
