Partial Refactorization in Sparse Matrix Solution: A New Possibility for Faster Nonlinear Finite Element Analysis

Algorithm 4

Simple 𝑗𝑘𝑖 form of 𝐋𝐃𝐋^T refactorization.

do row 𝑗=1, 𝑛

 if (CHANGE(𝑗) = 1) then

  do 𝑘 = all appropriate row

   RowTask (𝑘,𝑗)

  end do

  RowTask (𝑗,𝑗)

 end if

end do