¿Cómo afectan las modificaciones de bajo rango a la convergencia del método Krylov?

14

Decir que tengo un sistema lineal Ax=b , que converge rápidamente usando un método de Krylov adecuado (tal como CG o GMRES) para todos . Si es una matriz con un rango bajo , ¿el mismo método de Krylov en el sistema también convergerá rápidamente (idealmente con un número adicional de iteraciones que depende aproximadamente de )?bBr(A+B)x=br

Un ejemplo de dicho sistema sería una elasticidad y flexión de la membrana bien acondicionadas más términos de presión de aire no acondicionados con una estructura exterior densa del producto.

Tenga en cuenta que la pregunta es la misma con o sin preacondicionamiento, ya que es una modificación de rango de .P(A+B)Q=PAQ+PBQrPAQ

Geoffrey Irving
fuente

Respuestas:

7

Si su subespacio de Krylov se basa en las potencias de , la convergencia se retrasará en varias iteraciones como máximo en el rango de la corrección. Si se basa en potencias de , como máximo dos veces este número.A T AAATA

No he visto tal declaración en la literatura. Pero para ver la validez en el primer caso, es suficiente para mostrar que el ésimo espacio Krylov de la matriz donde tiene columnas está contenido en el espacio correspondiente sin correcciones de bajo rango, pero con una correspondientemente mayor índice . Esto es sencillo de verificar.A + U S V T U , V r k + rkA+USVTU,Vrk+r

Arnold Neumaier
fuente
¿Puedes explicar qué quieres decir con "basado en los poderes de "? El solucionador de Krylov solo recibe información sobre A + B , no A directamente. AA+BA
Geoffrey Irving
No importa: presumiblemente te refieres a los poderes de la matriz en cuestión, por lo que en este caso. A+B
Geoffrey Irving
Si. El método tiene una matriz como un parámetro, y esta matriz se denota generalmente por . A
Arnold Neumaier el
Tal vez para mayor interés podría reescribir su ecuación (o la solución) con algunos requisitos en para que podría ayudar si es nilpotente o de pequeña norma. También se reconoce la dependencia de la solución del problema. BB A - 1 Bx=(E+k=1(A1B)k)A1bBA1Bundisturbed
Bastian Ebeling