Multigrid (MG) puede usarse para resolver un sistema lineal construyendo una conjetura inicial x 0 y repitiendo lo siguiente para i = 0 , 1 .. hasta la convergencia:
- Calcule el residuo
- Aplique un ciclo de cuadrículas múltiples para obtener una aproximación , donde A e i = r i .
- Actualizar
El ciclo de múltiples cuadrículas es una secuencia de operaciones de suavizado, interpolación, restricción y resolución exacta de grillas gruesas aplicadas a para producir Δ x i . Esto es típicamente un ciclo V o un ciclo W. Esta es una operación lineal, por lo que escribimos Δ x i = B r i .
Uno puede interpretar este proceso como una iteración precondicionada de Richardson. Es decir, actualizamos .
La iteración de Richardson es un método prototípico del subespacio de Krylov, que sugiere el uso de ciclos de múltiples cuadrículas para preacondicionar otros métodos del subespacio de Krylov. Esto a veces se denomina multirredes "acelerado" con un método de Krylov, o alternativamente se puede ver como una opción de un preacondicionador para un método de Krylov.
Otra forma de extender el algoritmo anterior es emplear Full Multigrid (FMG). Vea esta respuesta para una descripción concisa.
¿En qué situaciones es preferible MG acelerado por Krylov a MG o FMG?