Ускорение полинейного рекуррентного метода в подпространствах Крылова

На примере алгоритма LR1 полинейного рекуррентного метода [1, 2] рассматриваются два механизма его ускорения в подпространствах Крылова. В качестве ускоряющего метода используется алгоритм Bi-CGStab P ван дер Ворста. Показано, что традиционный подход: построение предобуславливателя на базе алгоритма LR1, не приводит к требуемому результату. В то время как прямое сочетание алгоритмов LR1 и Bi-CGStab P позволяет значительно повысить скорость сходимости решения.