Abstract

In this paper we propose optimisation methods for variational regularisation problems based on discretising the inverse scale space flow with discrete gradient methods. Inverse scale space flow generalises gradient flows by incorporating a generalised Bregman distance as the underlying metric. Its discrete-time counterparts, Bregman iterations and linearised Bregman iterations are popular regularisation schemes for inverse problems that incorporate a priori information without loss of contrast. Discrete gradient methods are tools from geometric numerical integration for preserving energy dissipation of dissipative differential systems. The resultant Bregman discrete gradient methods are unconditionally dissipative and achieve rapid convergence rates by exploiting structures of the problem such as sparsity. Building on previous work on discrete gradients for non-smooth, non-convex optimisation, we prove convergence guarantees for these methods in a Clarke subdifferential framework. Numerical results for convex and non-convex examples are presented.

Highlights

  • This article is dedicated to Mila Nikolova whose keen interest and inspiring discussions have encouraged our research on geometric integration for optimisation

  • We study the Itoh–Abe discrete gradient method applied to the inverse scale space (ISS) flow

  • We review papers on discrete gradient methods for optimisation based on gradient flows

Read more

Summary

Introduction

This article is dedicated to Mila Nikolova whose keen interest and inspiring discussions have encouraged our research on geometric integration for optimisation. We propose and study optimisation schemes by using tools from geometric numerical integration to solve the inverse scale space (ISS) flow. The ISS flow is a dissipative system, and its dissipative structure is determined by the function J This allows one to solve (1.1) while incorporating a priori information into the optimisation scheme, with the benefits of converging to superior solutions, and doing so faster. We propose to discretise the inverse scale space flow with discrete gradient methods These are methods from geometric numerical integration that preserve the aforementioned geometric structures in a general setting. The theoretical convergence rates of the discrete gradient methods match those of explicit gradient descent and coordinate descent [21] The drawback of these methods is that the updates are in general implicit.

Related Literature
Structure and Contributions
Notation and Preliminaries
Convex Analysis
Non-Convex Subdifferential Analysis
The Discrete Gradient Method for the Inverse Scale Space Flow
Inverse Scale Space Flow and Bregman Methods
The Bregman Discrete Gradient Method
Theoretical Results
Well-Posedness and Preliminary Results
Amended Scheme
Main Convergence Theorem
Examples of Bregman Itoh–Abe Discrete Gradient Schemes
Sparse SOR Method
Equivalence of Iterative Methods for Linear Systems
Sparse SOR
Numerical Examples
Student-t Regularised Image Denoising
Conclusion
Else if
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call