Kriging (Gaussian Process Regression) is a powerful technique for creating surrogate models used in engineering design optimization when dealing with computationally expensive simulations. Despite its effectiveness, traditional Kriging implementations face computational challenges when scaling to large datasets, particularly for repeated inference operations required by optimization and uncertainty quantification applications. We present a high-performance PyTorch implementation of Kriging that significantly accelerates both training and inference while maintaining predictive accuracy. Our implementation leverages PyTorch's optimized tensor operations, automatic differentiation capabilities, and efficient memory management to achieve up to 20x speedup during inference with large datasets compared to traditional NumPy implementations. Benchmark results across datasets ranging from 100-2000+ training points demonstrate comparable interpolation accuracy with substantial runtime improvements, particularly for larger prediction set sizes. The results indicate that numerical stability is largely maintained between the NumPy implementation and PyTorch. This implementation is particularly valuable for optimization, uncertainty quantification, adaptive sampling, and other applications requiring repeated surrogate model evaluations. Our work demonstrates that modern deep learning frameworks can significantly benefit classical statistical learning algorithms beyond neural networks.