八度错误：SparseMatrix :: solve数值分解失败

Question

我有一个很大的方程组要解决。 系数存储在尺寸为320001 x 320001元素的稀疏矩阵CM中，其中18536032为非零。 结果向量B为320001个元素长。

执行时

I=CM\\B

八度错误：SparseMatrix :: solve数值分解失败

我收到上述错误消息。 简要查看源代码并没有给我任何线索。

有谁知道是什么导致了该错误？

顺便说一句：用较小的矩阵（例如180001x180001）解决相同的问题时，程序可以正常运行。

约翰内斯

Answer 1

Octave使用UMFPACK库来解决稀疏线性系统。 检查源表明错误消息是由于错误状态为负值引起的。 错误代码列表可在用户指南中找到 。 其中之一与内存不足有关：

UMFPACK错误的内存不足，（-1）：没有足够的内存。 ANSI C malloc或realloc例程失败。