Generalization Error Minimization in SubSampling for Gaussian Processes // Implements the GEMSS algorithm for sequential subdata selection in large-scale Gaussian process modeling (Chang, Hua, and Wu, 2026) doi:10.1080/00401706.2026.2670596. The method selects data points by a criterion consisting of predictive and space-filling parts, enabling efficient surrogate modeling for massive datasets.