Iterated total variation regularization with finite element methods for reconstruction the source term in elliptic systems