LU decomposition

From CFD-Wiki

Revision as of 00:47, 27 September 2005 (view source) Zxaar (Talk | contribs) ← Older edit		Revision as of 06:22, 3 October 2005 (view source) Zxaar (Talk | contribs) Newer edit →
Line 14:		Line 14:
	----		----
		+
		+
		+	----
		+	<i> Return to [[Numerical methods | Numerical Methods]] </i>

---

Revision as of 06:22, 3 October 2005

LU Solvers

For the system of equations Ax=b.
The solvers based on factorization are widely popular. The factorization of a non singular matrix A in two matrices L and U, called lower and upper matrices, leads to a direct procedure for the evaluation of the inverse of the matrix. Thus making it possible to calculate the solution vector with given source vector. It now remains a two step process: a forward substitution followed by a backward substitution.

Algorithm

Calculate LU factors of A

Ax = (LU)x = L(Ux) = Ly = b

Solve Ly = b by forward substitution

y = L^-1b

Solve Ux=y by backward substitution

x = U^-1y

---

---

Return to Numerical Methods