fixed_discrete_covariance: diagonal indices faster
Merge pull request #1 from dougalsutherland/master

fixed_discrete_covariance: get diagonal indices faster