#22 Initial work on replacing scalar variances with matrices

Also swap out all uses of np.eye with ssp.eye as it’s 10-times faster.