跳过numpy.linalg.eigh非收敛错误

时间:2016-02-15 09:36:55

标签: python numpy scipy linear-algebra

我有一个for - 循环,使用函数生成随机Hermitian矩阵,对角化并存储特征值。但是,对于某些实例,对角化无法收敛并且代码停止。

我想要的是,如果对角化失败,它应该跳到下一个实例。作为伪代码:

vals_list = []
for i in range(1000):
    M = random_matrix(....) #generates a 2D numpy array
    vals,vecs = np.linalg.eigh(M)
    vals_list.append(vals)

如何修改它以在对角化不收敛时跳过实例?

0 个答案:

没有答案