robust_inv_vTv function

Possibly use ridge-regularized inversion of crossprod(v)