def _sigma(matrix, mu):
n, p = matrix.shape
if p > n:
tmp = np.dot(matrix, matrix.T)
else:
tmp = np.dot(matrix.T, matrix)
return (la.norm(tmp, 2)/n) + mu
微信
公众号
扫码关注公众号