我将不得不解决一个大型线性系统。我现在正在寻找一个“无矩阵”的求解器（所以我只需要指定一个矩阵向量乘积，而不是矩阵）。据我了解（我不是专家。）LAPACK 不提供此选项。

你能推荐任何图书馆吗？还是我实现这种无矩阵的方法根本不合理？

备注：我计划使用 C++，矩阵来自 PDE 系统的 FEM，但不是对称的。矩阵是稀疏的，但由于我计划自己提供矩阵向量乘积，并且生成的向量将是密集的，因此稀疏不应该在这里发挥作用——据我所知。

问候，迈克尔