Me pregunto si el algoritmo de Thomas es la forma más rápida (¿probablemente?) De resolver un sistema tridiagonal disperso simétrico que domina diagonalmente en términos de complejidad algorítmica (no busca paquetes de implementación como LAPACK, etc.). Sé que tanto el algoritmo de Thomas como la multigrid son complejidad, pero ¿tal vez el factor constante para multigrid sea menor? No me parece que la multigrid pueda ser más rápida, pero no soy positivo.
Nota: Estoy considerando el caso donde las matrices son muy grandes. Se aceptan métodos directos o iterativos.