Practical evaluation of condensed phase quantum correlation functions: A Feynman–Kleinert variational linearized path integral method