我正在研究一些 c++ 代码来解决（共轭梯度、最小二乘共轭梯度、LSQR、..）大型稀疏超定线性方程组。我的矩阵和它们运行的​​系统有一个扭曲，这就是我自己编写求解器的原因。

出于测试、调试和基准测试的目的，我正在寻找一些具有众所周知的解决方案的测试矩阵向量问题。我希望能够为不同的矩阵大小生成这样的系统N x M。（通常我处理大小为非零元素的矩阵10^7 x 10^5。10^9）

您知道任何可以轻松生成的矩阵向量系统吗？