def pinv(x): #return np.linalg.pinv(x, max(x.shape) * np.spacing(np.linalg.norm(x))) #return scipy.linalg.pinv(x, max(x.shape) * np.spacing(np.linalg.norm(x))) return scipy.linalg.pinv(x, 1e-6)