Efficient computation of the exponential operator for large, sparse, symmetric matrices