Algoritmo de subespacios de Krylov para matrices hermíticas (simétricas) que genera una matriz tridiagonal mediante recurrencias cortas de tres términos, permitiendo extraer de forma eficiente eigenvalores extremos e información espectral aproximada de grandes operadores dispersos.

Principio

Aplicar recurrencias Lanczos de tres términos sucesivas para construir una base ortonormal del subespacio de Krylov mientras se genera una proyección tridiagonal simétrica cuyas eigenvalores aproximan las del operador hermítico original; aprovechar las recurrencias cortas para bajo coste por paso pero vigilar la pérdida de ortogonalidad numérica.

Demostración

Calcular los pocos eigenvalores mayores de una gran matriz de rigidez simétrica y dispersa ejecutando iteraciones de Lanczos con reortogonalización completa selectiva o reortogonalización parcial para evitar eigenvalores fantasma causados por la precisión finita.

Aplicación incorrecta

Ejecutar Lanczos indefinidamente sin reortogonalización en precisión finita, lo que conduce a pérdida de ortogonalidad y eigenvalores espurios repetidos (fantasmas) o multiplicidades incorrectas en el espectro calculado.

Consecuencia

Proporciona una representación tridiagonal compacta a partir de la cual se pueden obtener de forma económica algunos pares propios extremos, facilitando aproximación espectral, construcción de precondicionadores y reducción de modelos para grandes problemas hermíticos.

Inversión

Uso de métodos de ortogonalización completa como Arnoldi con Gram-Schmidt completo para problemas no simétricos, o solvers directos densos de eigenvalores que no explotan la esparcidad y requieren mucha más memoria y cálculo.

Límite

Se aplica principalmente a operadores hermíticos (realmente simétricos) o hermíticos complejos; el comportamiento en precisión finita requiere estrategias de reortogonalización; no es directamente aplicable a matrices no simétricas sin modificación (entonces es preferible Arnoldi).

Tensión semántica

Tensión con Arnoldi/otros métodos de Krylov: Lanczos usa recurrencias cortas de tres términos óptimas para problemas hermíticos pero es vulnerable a la pérdida de ortogonalidad numérica; Arnoldi maneja casos no hermíticos con recurrencias más largas y mayor estabilidad a mayor coste.

Síntesis

Método eficiente de aproximación espectral para operadores hermíticos que cambia recurrencias cortas y de bajo coste por vulnerabilidad a la pérdida de ortogonalidad en precisión finita, produciendo una proyección tridiagonal cuyos eigenvalores aproximan el espectro dominante del operador.