Implement Moore-Penrose pseudo-inverse

[lambday]

for compute function where the covariance matrix is singular.