This subroutine uses the Lanczos algorithm to compute the part of the spectrum of a large symmetric matrix A that lies in a specified interval, that is, it computes eigenvalues without regard to multiplicities. The user need only provide in the form of code which computes for any given vectors and . Auxiliary calls allow corresponding eigenvectors to be found.
Precision: At least 8-byte arithmetic is recommended.