[alasdairj]

Thanks! I implemented this for the pseudo-inverse calculation, and also in the bias in the final calculation of h. What surprised me was that the Ein got slightly worse rather than better, which is not what I expected. Or do I have an implementation issue?