我有兴趣在有限域(可能是接近 1000n
或所以)。从文献来看,块 Lanczos 方法可能是最合适的。Ax = b
GF(q)
q
我有 Linbox,它应该有这样的方法,但无法让 BlockLanczos 求解器在那里工作,一份报告说自 2003 年以来这已经被打破。该SparseElimination
方法确实有效,但似乎这不会很好因为矩阵的填充很大n
。
那么,有什么方法可以解决这些问题呢?
我有兴趣在有限域(可能是接近 1000n
或所以)。从文献来看,块 Lanczos 方法可能是最合适的。Ax = b
GF(q)
q
我有 Linbox,它应该有这样的方法,但无法让 BlockLanczos 求解器在那里工作,一份报告说自 2003 年以来这已经被打破。该SparseElimination
方法确实有效,但似乎这不会很好因为矩阵的填充很大n
。
那么,有什么方法可以解决这些问题呢?
Julia 支持有限域。我的教授对他的方法有一个简短的说明。它在第 37 行。LU decomp 和其他命令是内置的,并且是从 GF 类型派生的。