Discrete gradient methods for solving variational image regularisation models. (28th June 2017)