Product Kernel Interpolation for Scalable Gaussian Processes