python – 创建对角稀疏矩阵的有效方法

前端之家收集整理的这篇文章主要介绍了python – 创建对角稀疏矩阵的有效方法前端之家小编觉得挺不错的,现在分享给大家,也给大家做个参考。
我使用Numpy在 Python中有以下代码
p = np.diag(1.0 / np.array(x))

如何在不创建p的情况下对其进行变换以使稀疏矩阵p2具有与p相同的值?

解决方法

使用 scipy.sparse.spdiags(这很多,因此最初可能会令人困惑),scipy.sparse.dia_matrix和/或 scipy.sparse.lil_diags.(取决于 format你想要稀疏矩阵……)

例如.使用spdiags:

import numpy as np
import scipy as sp
import scipy.sparse

x = np.arange(10)

# "0" here indicates the main diagonal...
# "y" will be a dia_matrix type of sparse array,by default
y = sp.sparse.spdiags(x,x.size,x.size)

猜你在找的Python相关文章