Fast Algorithm for Non-Stationary Gaussian Process Prediction