(minor) might be nicer just to assigning 1 to the diagonal

maybe something like
```
cov = np.empty((dim, dim)), np.float64)
cov.fill(self.dparams)
cov[np.diag_indices(dim, dim)] = 1
```
minus syntax errors