We study and develop (stochastic) primal-dual block-coordinate descent methods for convex problems based on the method due to Chambolle and Pock. Our methods have known convergence rates for the iterates and the ergodic gap of O(1/N 2 ) if each block is strongly convex, O(1/N ) if no convexity is present, and more generally a mixed rate O(1/N 2 ) + O(1/N ) for strongly convex blocks if only some blocks are strongly convex. Additional novelties of our methods include blockwise-adapted step lengths and acceleration as well as the ability to update both the primal and dual variables randomly in blocks under a very light compatibility condition. In other words, these variants of our methods are doubly-stochastic. We test the proposed methods on various image processing problems, where we employ pixelwise-adapted acceleration.BLOCK-PROXIMAL METHODS WITH SPATIALLY ADAPTED ACCELERATION 17 conditions on the various step length and testing parameters. These conditions need to be verified through the development of explicit parameter update rules. We do this in Section 4 along with proving the claimed convergence rates (Theorem 4.5 and its corollaries). We also present there the final detailed versions of our proposed algorithms: Algorithm 1 (doubly stochastic) and Algorithm 2 (simplified). We finish with numerical experiments in Section 5.2. Background and overall structure of the algorithm. To make the notation definite, we write L(X; Y ) for the space of bounded linear operators between Hilbert spaces X and Y . The identity operator we denote by I. For T, S ∈ L(X; X), we use T ≥ S to indicate that T − S is positive semi-definite; in particular T ≥ 0 means that T is positive semi-definite. Also for possibly non-self-adjoint T , we introduce the inner product and norm-like notationsx, z T := T x, z andThe iteration index is off-by-one for σ ,i+1 and ψ ,i+1 for reasons of the historical development of the Chambolle-Pock method, when it was not written as a preconditioned proximal point method.