Interpolating moving least-squares methods for fitting potential energy surfaces: using classical trajectories to explore configuration space.