Vectorized computation of correlation functions from phase space trajectories generated by molecular dynamics calculations