寻找迭代线性求解器的测试矩阵/系统

Mar*_*rkD 8 c++ iteration matrix linear-algebra solver

我目前正在研究一个基于C++的大型稀疏线性代数问题库(是的,我知道很多这样的库存在,但我主要是为了学习迭代求解器,稀疏存储容器等等).

我正处于我在其他编程项目中使用求解器的地步,并希望针对不属于我自己的问题测试解算器.首先,我希望测试对称的稀疏系统是正定的.我找到了这样的系统矩阵的几个来源,例如:

Matrix Market UF稀疏矩阵集合

话虽这么说,我还没有找到任何包含整个系统 - 系统矩阵和RHS的良好测试矩阵的来源.为了检查结果,这将是很好的.关于我可以在哪里找到这样的完整系统的任何提示,或者我可以做什么来为我可以上网的系统矩阵生成"好的"RHS?我目前只是填充一个随机值或所有值的矩阵,但怀疑这不一定是最好的方法.

nor*_*lli 0

我还没有使用它,我正准备使用它,但 GiNAC 似乎是我为 C++ 找到的最好的东西。它是Maple背后用于CAS的库,我不知道它的性能如何。

http://www.ginac.de/