$\mathbf{v}_{j+1}=\ifrac{\mathbf{L}^{-1}\left(\mathbf{L}^{-1}\right)^{\rm T}%
\mathbf{u}}{\beta_{j+1}},$