Large-Scale Low-Rank Gaussian Process Prediction with Support Points