An earlier paper proved the convergence of a variable stepsize Bregman operator splitting algorithm (BOSVS) for minimizing $\phi(Bu)+H(u)$ where $H$ and $\phi$ are convex functions, and $\phi$ is possibly nonsmooth. The algorithm was shown to be relatively efficient when applied to partially parallel magnetic resonance image reconstruction problems. In this paper, the convergence rate of BOSVS is analyzed. When $H(u) = \|Au-f\|^2$ where $A$ is a matrix, it is shown that for an ergodic approximation ${\bf u}_k$ obtained by averaging $k$ BOSVS iterates, the error in the objective value $\phi (B\m{u}_k) + H(\m{u}_k)$ is $\C{O}(1/k)$. When the optimization problem has a unique solution $u^*$, we obtain the estimate $\|\m{u}_k - u^*\| = \C{O}(1/\sqrt{k})$. The theoretical analysis is compared to observed convergence rates for partially parallel magnetic resonance image reconstruction problems where $A$ is a large dense ill conditioned matrix.
Citation
University of Florida, September 23, 2013
Article
View An O(1/k) Convergence Rate for the Variable Stepsize Bregman Operator Splitting Algorithm