我正在寻找简单的GMRES实现。不幸的是,我有几个要求:
到目前为止,netlib 上的“模板”代码似乎很合适。但是,它不处理复杂的数据(尽管我认为将一些声明从 real 更改为 complex 可以解决这个问题),而且它显然更好地实现了操作留给 BLAS(我认为例如这个代码的Givens 旋转可能会更好被drot()调用)。另一种选择 PETSc 似乎有点矫枉过正。
drot()
我期待着您的建议。谢谢!