Pivoting e stabilizzazione

Il pivoting - ossia lo scambio preliminare tra righe - è importante per la stabilizzazione dell'algoritmo:

Si consideri:

$\displaystyle Ax = b \Leftrightarrow PAx = Pb $

Poichè $ PA = LU $:

$\displaystyle \begin{split}Ax = b & \Leftrightarrow PAx = Pb \\ & = LUx = Pb \\ & \Leftrightarrow Ux = \underbrace{PL^{-1}b}_{\beta} \end{split} $

È allora possibile risolvere il sistema trasformandolo in due sistemi triangolari nel seguente modo:

$\displaystyle \begin{cases}Ly = Pb \\ Ux = y \end{cases} $



Matteo Lisotto, Tobia Tesan - CC-BY 2.0