1

在我看来,在 Sage 中,创建密集矩阵和稀疏矩阵的唯一区别在于传递给构造函数的标志(“sparse=True”)。特别是,这意味着如果我想创建一个大小为 200,000 x 200,000 但只有 500,000 个非零条目的稀疏矩阵,我首先需要创建一个大小为 40,000,000,000 的列表,这是荒谬的。有没有更好的办法?

此外,由于我只对求解由矩阵定义的线性方程组感兴趣,因此我可以对使用“黑匣子”表示感到满意(即给出一个函数 g 使得 g(v)=Av 其中 A 是矩阵)。这在底层的 LinBox 库中是支持的,但我不知道对应的 Sage 接口是什么。

4

2 回答 2

3

使用字典指定稀疏矩阵

matrix(5000,5000,{(1,2): 1})
于 2012-02-11T20:51:19.307 回答
1

有关矩阵构造函数的更多示例,请参阅Sage 参考手册。

于 2012-05-17T16:21:36.217 回答