Solving non-linear elliptic difference equations by extendable sparse factorization procedures