Before revision:  The resulting linear algebraic systems are shown to be symmetric and therefore can be very efficiently solved by means of a standard matrix-free conjugate gradient algorithm.
After revision:  The resulting linear algebraic systems are shown to be symmetric and therefore can be solved by means of a standard matrix-free conjugate gradient algorithm.