我已经用Lis,Intel mkl和试过了PETSc。在任何地方您都需要将实际矩阵传递给GMRES或其他类似的迭代求解器，而我只想传递将我的矩阵乘以向量的函数（我的矩阵具有特殊结构并且完全形成它是低效的）。

我自己编写可以接受我的自定义矩阵向量产品的简单求解器没有问题，但我认为这种问题可能已经在许多库中得到解决和实现。