Mar*_*rkD 8 c++ iteration matrix linear-algebra solver
我目前正在研究一个基于C++的大型稀疏线性代数问题库(是的,我知道很多这样的库存在,但我主要是为了学习迭代求解器,稀疏存储容器等等).
我正处于我在其他编程项目中使用求解器的地步,并希望针对不属于我自己的问题测试解算器.首先,我希望测试对称的稀疏系统是正定的.我找到了这样的系统矩阵的几个来源,例如:
话虽这么说,我还没有找到任何包含整个系统 - 系统矩阵和RHS的良好测试矩阵的来源.为了检查结果,这将是很好的.关于我可以在哪里找到这样的完整系统的任何提示,或者我可以做什么来为我可以上网的系统矩阵生成"好的"RHS?我目前只是填充一个随机值或所有值的矩阵,但怀疑这不一定是最好的方法.
我还没有使用它,我正准备使用它,但 GiNAC 似乎是我为 C++ 找到的最好的东西。它是Maple背后用于CAS的库,我不知道它的性能如何。