主題
Search

廣義最小殘差法


廣義最小殘差 (GMRES) 法(Saad 和 Schultz 1986)是最小殘差法 (MINRES) 的擴充套件,後者僅適用於對稱系統,而前者適用於非對稱系統。與 MINRES 類似,它生成一個正交向量序列,但在缺少對稱性的情況下,這不再能用短遞推關係完成;相反,正交序列中所有先前計算的向量都必須保留。因此,使用了該方法的“重啟”版本。

共軛梯度法中,殘差構成該空間的正交基

 span{r^((0)),Ar^((0)),A^2r^((0))}.

在 GMRES 中,這個基被顯式地構造出來

讀者可能會認識到這是一個改進的 Gram-Schmidt 正交化。應用於 Krylov 序列 {A^kr^((0))} 的這種正交化被稱為“Arnoldi 方法”(Arnoldi 1951)。內積係數 (w^((i)),v^((k)))|w^((i))| 儲存在一個上 Hessenberg 矩陣中。

GMRES 迭代被構造為

 x^((i))=x^((0))+y_1v^((1))+...+y_iv^((i)),

其中係數 y_k 的選擇是為了最小化殘差範數 |b-Ax^((i))|。GMRES 演算法的特性是,即使迭代尚未形成,也可以計算此殘差範數。因此,形成迭代的昂貴操作可以推遲到殘差範數被認為足夠小的時候。

在此方案中,UPDATE(x^~,i) 替換了以下計算

廣義最小殘差法旨在求解非對稱線性系統(Saad 和 Schultz 1986)。GMRES 最流行的形式是基於改進的 Gram-Schmidt 正交化程式,並使用重啟來控制儲存需求。

如果不使用重啟,GMRES(像任何正交化 Krylov 子空間方法一樣)將在不超過 n 步內收斂(假設精確算術)。當然,當 n 很大時,這沒有實際價值;此外,在沒有重啟的情況下,儲存和計算需求是令人望而卻步的。實際上,成功應用 GMRES(m) 的關鍵要素圍繞何時重啟的決定;也就是 m 的選擇。不幸的是,存在一些例子,其中該方法停滯不前,並且收斂僅發生在第 n 步。對於此類系統,任何小於 nm 選擇都無法收斂。

Saad 和 Schultz (1986) 已經證明了幾個有用的結果。特別是,他們表明,如果係數矩陣 A 是實數且接近正定,則可以選擇 m 的“合理”值。下面討論 m 選擇的含義。

Saad 和 Schultz (1986) 提出了一種常見的 GMRES 實現,它依賴於使用改進的 Gram-Schmidt 正交化。Householder 變換雖然成本相對較高但很穩定,也被提出。Householder 方法導致與內積和向量更新相關的工作量增加三倍(不是與矩陣向量乘積);然而,收斂性可能會更好,特別是對於病態系統(Walker 1988)。從並行性的角度來看,Gram-Schmidt 正交化可能是首選,犧牲一些穩定性以獲得更好的並行化特性(Demmel et al. 1993)。在這裡,我們採用改進的 Gram-Schmidt 方法。

GMRES 的主要缺點是,每次迭代所需的工作量和儲存量隨著迭代次數線性增加。除非足夠幸運地獲得極快的收斂速度,否則成本將迅速變得令人望而卻步。克服此限制的常用方法是重啟迭代。在選定的迭代次數 m 後,累積的資料被清除,中間結果用作下一次 m 迭代的初始資料。重複此過程,直到實現收斂。困難在於為 m 選擇合適的值。如果 m 太小,則 GMRES(m) 可能收斂緩慢,或者完全無法收斂。大於必要的 m 值會涉及過多的工作(並使用更多的儲存空間)。不幸的是,沒有明確的規則來指導 m 的選擇;選擇何時重啟是一項經驗問題。

有關向量和共享記憶體計算機的 GMRES 的討論,請參閱 Dongarra et al. (1991)。有關更通用的架構,請參閱 Demmel et al. (1993)。


另請參閱

雙共軛梯度法, Chebyshev 迭代, 正規方程上的共軛梯度法 共軛梯度法, 共軛梯度平方方法, 靈活的廣義最小殘差法, 線性方程組, 最小殘差法, 非平穩迭代法, 預處理器, 擬最小殘差法 平穩迭代法, 對稱 LQ 方法, 無轉置擬最小殘差法

此條目由 Noel Black 和 Shirley Moore 貢獻,改編自 Barrett et al. (1994) (作者連結)

使用 探索

參考文獻

Arnoldi, W. "The Principle of Minimized Iterations in the Solution of the Matrix Eigenvalue Problem." Quart. Appl. Math. 9, 17-29, 1951.Barrett, R.; Berry, M.; Chan, T. F.; Demmel, J.; Donato, J.; Dongarra, J.; Eijkhout, V.; Pozo, R.; Romine, C.; and van der Vorst, H. Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods, 2nd ed. Philadelphia, PA: SIAM, 1994. http://www.netlib.org/linalg/html_templates/Templates.html.Demmel, J.; Heath, M.; and van der Vorst, H. "Parallel Numerical Linear Algebra." In Acta Numerica, Vol. 2. Cambridge, England: Cambridge University Press, 1993.Dongarra, J.; Duff, I.; Sorensen, D.; and van der Vorst, H. Solving Linear Systems on Vector and Shared Memory Computers. Philadelphia, PA: SIAM, 1991.Saad, Y. and Schultz, M. "GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems." SIAM J. Sci. Statist. Comput. 7, 856-869, 1986.Walker, H. F. "Implementation of the GMRES Method using Householder Transformations." SIAM J. Sci. Statist. Comput. 9, 152-163, 1988.

在 上引用

廣義最小殘差法

引用為

Black, NoelMoore, Shirley. "廣義最小殘差法。" 來自 Web 資源,由 Eric W. Weisstein 建立。 https://mathworld.tw/GeneralizedMinimalResidualMethod.html

主題分類