Exact Gaussian Processes for Massive Datasets via Non-Stationary Sparsity-Discovering Kernels