Path integral solution handled by Fast Gauss Transform