python - 为什么我在 Sage 中遇到 "large"矩阵的问题？

Question

我正在尝试使用 sage 中的 22x22 矩阵进行一些计算。这似乎并不是那么糟糕的事情，特别是因为矩阵非常稀疏。但是，当我尝试对矩阵执行任何操作时，我会收到“IndexError:字符串索引超出范围”或无意义的计算。给出了什么？

即使非常简单，如果我尝试以下代码:

M = matrix(ZZ,20); M

这应该显示一个 20x20 的零矩阵。然而，它却给出了相同的 IndexError，其中

M = matrix(ZZ,19); M

不会，并且给出的正是您所期望的。

有趣的是，输入我的矩阵后，任何显示它的尝试都会给出 IndexError。然而，尽管我对它们的正确性持怀疑态度，但我仍然可以做一些事情来给出“答案”。例如，我可以做类似的事情

M.parent()

这似乎是有道理的。然而，

M.determinant()

输出一个数字，但我大约 99% 确定它给出的数字与我输入的矩阵的行列式几乎没有关系。

那怎么办？ Sage 是否无法处理大于 19x19 的矩阵？

编辑:这是在 Mac OS X 10.9.4 上，我的 sage 版本是 5.10。这是在重新启动 sage 之后，这给了我同样的错误。然而，看来我应该更新 sage 看看它是否能解决这个问题......

Answer 1

这是一个已知且已修复的错误，请参阅票证 #14785和 #14579在 Sage's trac ，我们了解到该错误是在 Python 中，请参阅issue 17526在 Python's bug tracker ,它已在上游解决，并且自版本 5.11.beta3 以来在 Sage 中运行良好。

我同意 John Palmieri 鼓励更新的观点。 Sage 变得越来越好，所以它总是值得使用最新版本。目前 Sage 6.2 已发布，Sage 6.3.rc1 已发布已经发布并且运行良好，所以 Sage 6.3 应该很快就会发布。

关于行列式，如果计算错误我会感到非常惊讶，但是你为什么不在另一个计算机代数系统中计算它来检查呢？如果确实如此不正确，请报告错误。在网络上搜索“wims 决定因素”以获取在线信息行列式计算器。

编辑(2014-08-10): Sage 6.3 is out!

当我编辑时，我会链接到 wims matrix tool ,只是为了完整性。我还承认我还没有检查 Wims 和 Sage 呼吁谁来计算这样的行列式；它最终可能是相同的软件，所以不会有太多一张支票。很高兴您以其他方式检查了它。