CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Wiki > LU decomposition

LU decomposition

From CFD-Wiki

(Difference between revisions)
Jump to: navigation, search

Revision as of 00:47, 27 September 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

My wiki