We present a new method for calculating the density of states of large-scal
e molecular systems. The new algorithm does not require calculation of all
eigenvalues of the force constant matrix. The histogram for the eigenvalue
distribution is constructed by performing a sequence of sparse LDLT matrix
factorizations. Moreover, the entire calculation can be fully parallelized
on a multi-processor computer system.