迭代Krylov子空间求解器的性能

Question

我正在使用几个预处理的迭代Krylov求解器（gmres，bicg，bicgstab，cgs，qmr）来成功解决中等大小（10000x10000）的线性系统，但是有一个问题困扰着我。

当我使用相同的krylov解算器重复解决同一系统时，结果是相同的，但每次重复的CPU时间不同，最多1秒。