rmgarnett/gpml_extensions — pull request #1 by dougalsutherland (2014-06-04 17:05:58+00:00)
Title: fixed_discrete_covariance: get diagonal indices faster

Profiling repeated calls to model(), way too much time was being spent on `diagonal = diag(K)`. This way doesn't copy the diagonal.