Evaluate diagonal only on predict std

[mhangaard]

Predicting mean and covariance on a test set, X_test, scale N**2, because we are constructing K(X_test, X_test). In case X_test is large, we would prefer to not calculate the covariance matrix, but just the standard deviation based on the diagonal (See gpfunctions.uncertainty)