问题

是否有实现块 Krylov 子空间方法的库？（我无法从简单的 Google 搜索中找到任何内容。）

背景

现在，我正在使用一个解决多个系统形式的代码

在哪里$A$是$n$经过$n$,$i = 1, \ldots, m$，一般来说，$n \gg m$($n$至少是 10 倍以上$m$，并且可能比$m$）。此外，$A$足够大，迭代方法更可取。

目前，代码解决了$Ax_{i} = b_{i}$使用给定的预条件子对每个单独的右侧使用 Krylov 子空间方法，这种方法似乎工作得很好。（代码还会重新计算预处理器，这是不必要的。）

为了加快代码速度，似乎值得考虑解决

相反，在哪里$B = [b_{1}\, b_{2}\, \ldots\, b_{m}]$是一个$n$经过$m$矩阵收集所有$b_{i}$， 和$X$也是一个$n$经过$m$矩阵收集$x_{i}$.

我想尝试解决$AX = B$使用块 Krylov 子空间方法（使用与以前相同的预处理器）来查看它是否优于当前的求解方法$Ax_{i} = b_{i}$并循环$i$. 解决方案既可能也不太可能$Ax_{i} = b_{i}$然后可以用作系统的猜测$Ax_{j} = b_{j}$为了$i \neq j$; 我不希望各种系统之间存在关系。