Estou interessado nos casos em que o gradiente conjugado funciona muito melhor que o método GMRES. Em geral, o CG é a escolha preferível em muitos casos de SPD (simétrico-positivo-definido) porque requer menos armazenamento e o limite teórico da taxa de convergência para o CG é o dobro do GMRES....
Em métodos como gmres ou bicgstab, pode ser atraente usar outro método de krylov como pré-condicionador. Afinal, eles são fáceis de implementar de maneira livre de matriz e em um ambiente paralelo. Por exemplo, um exemplo pode usar algumas (digamos ~ 5) iterações da bigcstab não condicionada como...
Atualmente, estou trabalhando na solução de sistemas simétricos muito grandes (mas não positivos definidos), gerados por alguns algoritmos. Essas matrizes têm uma boa esparsidade de bloco que pode ser usada para solução paralela. Mas não consigo decidir se devo usar uma abordagem direta (como...