On the repeated solution of matrix equations in imaging techniques

In various imaging procedures or in linear optimization it is necessary to solve large systems of linear algebraic equations with a few equations being changed many times. Once a complete solution of an n × n system is obtained by applying a classical Gaussian elimination method [1]-[4], for instance, the updated solution of the system, when its last equation is changed, requires 2n2 arithmetic operations for large systems.