from alexandr.shadchin at gmail, with some tweaks
ARPACK is a collection of Fortran77 subroutines designed to solve large scale eigenvalue problems.