Parallel Computation of Sturm-Liouville Spectral Density Functions