A fast, optimal spatial-prediction method for massive datasets