Método del gradiente biconjugado estabilizado

Método del gradiente biconjugado estabilizado

En álgebra lineal numérica, el método del gradiente biconjugado estabilizado, generalmente abreviado como BiCGSTAB (del inglés «biconjugate gradient stabilized method»), es un método iterativo propuesto por H. A. van der Vorst para la resolución numérica de los sistemas de ecuaciones lineales no simétricos. Es una variante del método del gradiente biconjugado (BiCG) y ofrece convergencia más rápida y suave que el original BiCG así como otras variantes como el método del gradiente conjugado cuadrado (CGS). Es un método del subespacio de Krylov.

Contenido

Pasos algorítmicos

BiCGSTAB sin precondicionamiento

Para resolver el sistema \boldsymbol{Ax}=\boldsymbol{b}, el BiCGSTAB comienza con una aproximación inicial \boldsymbol{x}_0 y procede como sigue:

  1. \boldsymbol{r}_0=\boldsymbol{b}-\boldsymbol{Ax}
  2. Elige un vector arbitrario \boldsymbol{\hat{r}}_0 tal que (\boldsymbol{\hat{r}}_0,\boldsymbol{r}_0) \neq0, por ejemplo, \boldsymbol{\hat{r}}_0=\boldsymbol{r}_0
  3. \rho_0=\alpha=\omega_0=1\;
  4. \boldsymbol{v}_0=\boldsymbol{p}_0=\boldsymbol{0}
  5. Para i=1,2,3,\ldots\;
    1. \rho_i=(\boldsymbol{\hat{r}}_0,\boldsymbol{r}_{i-1})
    2. \beta=(\rho_i/\rho_{i-1})(\alpha/\omega_{i-1})\;
    3. \boldsymbol{p}_i=\boldsymbol{r}_{i-1}+\beta(\boldsymbol{p}_{i-1}-\omega_{i-1}\boldsymbol{v}_{i-1})
    4. \boldsymbol{v}_i=\boldsymbol{Ap}_i
    5. \alpha=\rho_i/(\boldsymbol{\hat{r}}_0,\boldsymbol{v}_i)
    6. \boldsymbol{s}=\boldsymbol{r}_{i-1}-\alpha\boldsymbol{v}_i
    7. \boldsymbol{t}=\boldsymbol{As}
    8. \omega_i=(\boldsymbol{t},\boldsymbol{s})/(\boldsymbol{s},\boldsymbol{s})
    9. \boldsymbol{x}_i=\boldsymbol{x}_{i-1}+\alpha\boldsymbol{p}_i+\omega_i\boldsymbol{s}
    10. Termina si \boldsymbol{x}_i es lo suficientemente preciso
    11. \boldsymbol{r}_i=\boldsymbol{s}-\omega_i\boldsymbol{t}

BiCGSTAB precondicionado

Generalmente se utiliza los precondicionadores para acelerar la convergencia de los métodos iterativos. Para resolver el sistema \boldsymbol{Ax}=\boldsymbol{b} con un precondicionador \boldsymbol{K}=\boldsymbol{K}_1\boldsymbol{K}_2\approx\boldsymbol{A}, el BiCGSTAB precondicionado comienza con una aproximación inicial \boldsymbol{x}_0 y procede como sigue:

  1. \boldsymbol{r}_0=\boldsymbol{b}-\boldsymbol{Ax}
  2. Elige un vector arbitrario \boldsymbol{\hat{r}}_0 tal que (\boldsymbol{\hat{r}}_0,\boldsymbol{r}_0) \neq0, por ejemplo, \boldsymbol{\hat{r}}_0=\boldsymbol{r}_0
  3. \rho_0=\alpha=\omega_0=1\;
  4. \boldsymbol{v}_0=\boldsymbol{p}_0=\boldsymbol{0}
  5. Para i=1,2,3,\ldots\;
    1. \rho_i=(\boldsymbol{\hat{r}}_0,\boldsymbol{r}_{i-1})
    2. \beta=(\rho_i/\rho_{i-1})(\alpha/\omega_{i-1})\;
    3. \boldsymbol{p}_i=\boldsymbol{r}_{i-1}+\beta(\boldsymbol{p}_{i-1}-\omega_{i-1}\boldsymbol{v}_{i-1})
    4. \boldsymbol{y}=\boldsymbol{K}^{-1}\boldsymbol{p}_i
    5. \boldsymbol{v}_i=\boldsymbol{Ay}
    6. \alpha=\rho_i/(\boldsymbol{\hat{r}}_0,\boldsymbol{v}_i)
    7. \boldsymbol{s}=\boldsymbol{r}_i-\alpha\boldsymbol{v}_i
    8. \boldsymbol{z}=\boldsymbol{As}
    9. \boldsymbol{t}=\boldsymbol{K}^{-1}\boldsymbol{z}
    10. \omega_i=(\boldsymbol{K}_1^{-1}\boldsymbol{t},\boldsymbol{K}_1^{-1}\boldsymbol{s})/(\boldsymbol{K}_1^{-1}\boldsymbol{s},\boldsymbol{K}_1^{-1}\boldsymbol{s})
    11. \boldsymbol{x}_i=\boldsymbol{x}_{i-1}+\alpha\boldsymbol{y}+\omega_i\boldsymbol{z}
    12. Termina si \boldsymbol{x}_i es lo suficientemente preciso
    13. \boldsymbol{r}_i=\boldsymbol{s}-\omega_i\boldsymbol{t}

Esta formulación es equivalente a aplicar el BiCGSTAB sin precondicionamiento al sistema explícitamente precondicionado

\boldsymbol{\tilde{A}\tilde{x}}=\boldsymbol{\tilde{b}}

con \boldsymbol{\tilde{A}}=\boldsymbol{K}_1^{-1}\boldsymbol{AK}_2^{-1}, \boldsymbol{\tilde{x}}=\boldsymbol{K}_2\boldsymbol{x}, \boldsymbol{\tilde{b}}=\boldsymbol{K}_1^{-1}\boldsymbol{b}. En otras palabras, precondicionamiento a ambas la izquierda y la derecha es posible en esta formulación.

Generalización

BiCGSTAB puede ser visto como una combinación de BiCG y GMRES en que cada paso de BiCG se sigue por un paso de GMRES(1) (GMRES reiniciado en cada paso) para reparar el comportamiento irregular de convergencia de CGS, de lo cual BiCGSTAB fue desarrollado como una mejora. No obstante, debido al uso de los polinomios del residuo mínimo de grado uno, la dicha reparación puede no ser eficaz si la matriz \boldsymbol{A} tiene pares propios complejos grandes. En tales casos, es probable que BiCGSTAB se estanca como lo confirman los experimentos numéricos.

Se puede esparar que los polinomios del residuo mínimo de grado más alto pueda mejor manejar esta situación. Esto da lugar a los métodos que incluyen BiCGSTAB2[1] La plantilla {{ref}} está obsoleta, véase el nuevo sistema de referencias. y el más general BiCGSTAB(l)[2] La plantilla {{ref}} está obsoleta, véase el nuevo sistema de referencias.. En BiCGSTAB(l), un paso de GMRES(l) sigue cada l pasos de BiCG. BiCGSTAB2 es equivalente a BiCGSTAB(l) con l = 2.

Véase también

Referencias


Wikimedia foundation. 2010.

Игры ⚽ Поможем сделать НИР

Compartir el artículo y extractos

Link directo
Do a right-click on the link above
and select “Copy Link”