Geometric iterative methods, including progressive approximation and geometric interpolation are efficient for fitting a given data set. With the development of big technology, number points has become massive, least-squares (LSPIA) is generally applied to fit mass data. Combining Schulz method calculating Moore–Penrose generalized inverse matrix with traditional LSPIA method, this paper presen...