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

Algorithm 3

Spread the effect of modification in 𝑗𝑘𝑖 factorization.

do 𝑗=1, 𝑛

 if (CHANGE(𝑗) = 1) then

  do 𝑖=𝑗+1, 𝑛

   if (𝑢_𝑗𝑖 is nonzero) then CHANGE(𝑖) = 1

  end do

 end if

end do