Optimistic mirror descent in saddle-point problems: Going the extra (gradient) mile

Owing to their connection with generative adversarial networks ( GANs ) , saddlepoint problems have recently attracted considerable interest in machine learning and beyond . By necessity , most theoretical guarantees revolve around convexconcave ( or even linear ) problems ; however , making theoretical inroads towards efficient GAN training depends crucially on moving beyond this classic framework . To make piecemeal progress along these lines , we analyze the behavior of mirror descent ( MD ) in a class of non-monotone problems whose solutions coincide with those of a naturally associated variational inequality – a property which we call coherence . We first show that ordinary , “ vanilla ” MD converges under a strict version of this condition , but not otherwise ; in particular , it may fail to converge even in bilinear models with a unique solution . We then show that this deficiency is mitigated by optimism : by taking an “ extra-gradient ” step , optimistic mirror descent ( OMD ) converges in all coherent problems . Our analysis generalizes and extends the results of Daskalakis et al . [ 2018 ] for optimistic gradient descent ( OGD ) in bilinear problems , and makes concrete headway for provable convergence beyond convex-concave games . We also provide stochastic analogues of these results , and we validate our analysis by numerical experiments in a wide array of GAN models ( including Gaussian mixture models , and the CelebA and CIFAR-10 datasets ) . 1 Introduction The surge of recent breakthroughs in deep learning has sparked significant interest in solving optimization problems that are universally considered hard . Accordingly , the need for an effective theory has two different sides : first , a deeper understanding would help demystify the reasons behind the success and/or failures of different training algorithms ; second , theoretical advances can inspire effective algorithmic tweaks leading to concrete performance gains . For instance , using tools from the theory of dynamical systems , Lee et al . [ 2016 , 2019 ] and Panageas & Piliouras [ 2017 ] showed that a wide variety of first-order methods ( including gradient descent and mirror descent ) almost always avoid saddle points . More generally , the optimization and machine learning communities alike have dedicated significant effort in understanding non-convex landscapes by searching for properties which could be leveraged for efficient training . As an example , the “ strict saddle ” property was shown to hold in a wide range of salient objective functions ranging from low-rank matrix factorization [ Bhojanapalli et al. , 2016 ; Ge et al. , 2017 ] and dictionary learning [ Sun et al. , 2017a , b ] , to principal component analysis [ Ge et al. , 2015 ] , and many other models . On the other hand , adversarial deep learning is nowhere near as well understood , especially in the case of generative adversarial networks ( GANs ) [ Goodfellow et al. , 2014 ] . Despite an immense amount of recent scrutiny , our theoretical understanding can not boast similar breakthroughs as in “ single-agent ” deep learning . Because of this , a considerable corpus of work has been devoted to exploring and enhancing the stability of GANs , including techniques as diverse as the use of Wasserstein metrics [ Arjovsky et al. , 2017 ] , critic gradient penalties [ Gulrajani et al. , 2017 ] , feature matching , minibatch discrimination , etc . [ Radford et al. , 2016 ; Salimans et al. , 2016 ] . Even before the advent of GANs , work on adaptive dynamics in general bilinear zero-sum games ( e.g . Rock-Paper-Scissors ) established that they lead to persistent , chaotic , recurrent ( i.e . cycle-like ) behavior [ Sato et al. , 2002 ; Piliouras & Shamma , 2014 ; Piliouras et al. , 2014 ] . Recently , simple specific instances of cycle-like behavior in bilinear games have been revisited mainly through the lens of GANs [ Mertikopoulos et al. , 2018 ; Daskalakis et al. , 2018 ; Mescheder et al. , 2018 ; Papadimitriou & Piliouras , 2018 ] . Two important recent results have established unified pictures about the behavior of continuous and discrete-time first order methods in bilinear games : First , Mertikopoulos et al . [ 2018 ] established that continuous-time descent methods in zero-sum games ( e.g. , gradient descent , follow-the-regularized-leader and the like ) are Poincaré recurrent , returning arbitrarily closely to their initial conditions infinitely many times . Second , Bailey & Piliouras [ 2018 ] examined the discrete-time analogues ( gradient descent , multiplicative weights and follow-the-regularized-leader ) showing that orbits spiral slowly outwards . These recurrent systems have formal connections to Hamiltonian dynamics and do not behave in a gradient-like fashion Balduzzi et al . [ 2018 ] ; Bailey & Piliouras [ 2019 ] . This is a critical failure of descent methods , but one which Daskalakis et al . [ 2018 ] showed can be overcome through “ optimism ” , interpreted in this context as an “ extra-gradient ” step that pushes the training process further along the incumbent gradient – as a result , optimistic gradient descent ( OGD ) succeeds in cases where vanilla gradient descent ( GD ) fails ( specifically , unconstrained bilinear saddle-point problems ) . A common theme in the above is that , to obtain a principled methodology for training GANs , it is beneficial to first establish improvements in a more restricted setting , and then test whether these gains carry over to more demanding learning environments . Following these theoretical breadcrumbs , we focus on a class of non-monotone problems whose solutions are related to those of a naturally associated variational inequality , a property which we call coherence . Then , hoping to overcome the shortcomings of ordinary descent methods by exploiting the problem ’ s geometry , we examine the convergence of MD in coherent problems . On the positive side , we show that if a problem is strictly coherent ( a condition satisfied by all strictly convex-concave problems ) , MD converges almost surely , even in stochastic problems ( Theorem 3.1 ) . However , under null coherence ( the “ saturated ” opposite to strict coherence ) , MD spirals outwards from the problem ’ s solutions and may cycle in perpetuity . The null coherence property covers all bilinear models , so this result encompasses fully the analysis of Bailey & Piliouras [ 2018 ] for GD and follow-the-regularized-leader ( FTRL ) in general bilinear zero-sum games within our coherence framework . Thus , in and by themselves , gradient/mirror descent methods do not suffice for training convoluted , adversarial deep learning models . To mitigate this deficiency , we consider the addition of an extra-gradient step which looks ahead and takes an additional step along a “ future ” gradient . This technique was first introduced by Korpelevich [ 1976 ] and subsequently gained great popularity as the basis of the mirror-prox algorithm of Nemirovski [ 2004 ] which achieves an optimal O ( 1/n ) convergence rate in Lipschitz monotone variational inequalities ( see also Nesterov , 2007 , for a primal-dual variant of the method and Juditsky et al. , 2011 , for an extension to stochastic variational inequalities and saddle-point problems ) . In the learning literature , the extra-gradient technique ( or , sometimes , a variant thereof ) is often referred to as optimistic mirror descent ( OMD ) [ Rakhlin & Sridharan , 2013 ] and its effectiveness in GAN training was recently examined by Daskalakis et al . [ 2018 ] and Yadav et al . [ 2018 ] ( the latter involving a damping mechanism for only one of the players ) . More recently , Gidel et al . [ 2018 ] considered a variant method which incorporates a mechanism that “ extrapolates from the past ” in order to circumvent the need for a second oracle call in the extra-gradient step . Specifically , Gidel et al . [ 2018 ] showed that the extra-gradient algorithm with gradient reuse converges a ) geometrically in strongly monotone , deterministic variational inequalities ; and b ) ergodically in general stochastic variational inequalities , achieving in that case an oracle complexity bound that is √ 13/7/2 ≈ 68 % of a bound previously established by Juditsky et al . [ 2011 ] for the mirror-prox algorithm . However , beyond convex-concave problems , averaging offers no tangible benefits because there is no way to relate the value of the ergodic average to the value of the iterates . As a result , moving closer to GAN training requires changing both the algorithm ’ s output as well as the accompanying analysis . With this as our guiding principle , we first show that the last iterate of OMD converges in all coherent problems , including null-coherent ones . As a special case , this generalizes and extends the results of Noor et al . [ 2011 ] for OGD in pseudo-monotone problems , and also settles in the affirmative an open question of Daskalakis et al . [ 2018 ] concerning the convergence of the last iterate of OGD in nonlinear problems . Going beyond deterministic problems , we also show that OMD converges with probability 1 even in stochastic saddle-point problems that are strictly coherent . These results complement the existing literature on the topic by showing that a cheap extra-gradient add-on can lead to significant performance gains when applied to state-of-the-art methods ( such as Adam ) . We validate this prediction for a wide array of standard GAN models in Section 5 . 2 Problem setup and preliminaries Saddle-point problems . Consider a saddle-point problem of the general form min x1∈X1 max x2∈X2 f ( x1 , x2 ) , ( SP ) where each feasible region Xi , i = 1 , 2 , is a compact convex subset of a finite-dimensional normed space Vi ≡ di , and f : X ≡ X1 × X2 → denotes the problem ’ s value function . From a gametheoretic standpoint , ( SP ) can be seen as a zero-sum game between two optimizing agents ( or players ) : Player 1 ( the minimizer ) seeks to incur the least possible loss , while Player 2 ( the maximizer ) seeks to obtain the highest possible reward – both determined by f ( x1 , x2 ) . To solve ( SP ) , we will focus on incremental processes that exploit the individual loss/reward gradients of f ( assumed throughout to be at least C1-smooth ) . Since the individual gradients of f will play a key role in our analysis , we will encode them in a single vector as g ( x ) = ( g1 ( x ) , g2 ( x ) ) = ( ∇x1 f ( x1 , x2 ) , −∇x2 f ( x1 , x2 ) ) ( 2.1 ) and , following standard conventions , we will treat g ( x ) as an element of Y ≡ V∗ , the dual of the ambient space V ≡ V1 × V2 , assumed to be endowed with the product norm ‖x‖2 = ‖x1‖2 + ‖x2‖2 . Variational inequalities and coherence . Most of the literature on saddle-point problems has focused on the monotone case , i.e. , when f is convex-concave . In such problems , it is well known that solutions of ( SP ) can be characterized equivalently as solutions of the Stampacchia variational inequality 〈g ( x∗ ) , x − x∗〉 ≥ 0 for all x ∈ X ( SVI ) or , in Minty form : 〈g ( x ) , x − x∗〉 ≥ 0 for all x ∈ X . ( MVI ) The equivalence between solutions of ( SP ) , ( SVI ) and ( MVI ) extends well beyond the realm of monotone problems : it trivially includes all bilinear problems ( f ( x1 , x2 ) = x > 1 Mx2 ) , pseudo-monotone objectives as in Noor et al . [ 2011 ] , etc . For a concrete example which is not even pseudo-monotone , consider the problem min x1∈ [ −1,1 ] max x2∈ [ −1,1 ] ( x41x 2 2 + x 2 1 + 1 ) ( x 2 1x 4 2 − x22 + 1 ) . ( 2.2 ) The only saddle-point of f is x∗ = ( 0 , 0 ) : it is easy to check that x∗ is also the unique solution of the corresponding variational inequality ( VI ) problems , despite the fact that f is not even pseudomonotone.1 This shows that this equivalence encompasses a wide range of phenomena that are innately incompatible with convexity/monotonicity , even in the lowest possible dimension ; for an in-depth discussion we refer the reader to Facchinei & Pang [ 2003 ] . Motivated by all this , we introduce below the following notion of coherence : Definition 2.1 . We say that ( SP ) is coherent if : 1 . Every solution of ( SVI ) also solves ( SP ) . 2 . There exists a solution p of ( SP ) that satisfies ( MVI ) . 3 . Every solution x∗ of ( SP ) satisfies ( MVI ) locally , i.e. , for all x sufficiently close to x∗ . In the above , if ( MVI ) holds as a strict inequality whenever x is not a solution thereof , ( SP ) will be called strictly coherent ; by contrast , if ( MVI ) holds as an equality for all x ∈ X , we will say that ( SP ) is null-coherent . The notion of coherence will play a central part in our considerations , so a few remarks are in order . To the best of our knowledge , its first antecedent is a gradient condition examined by Bottou [ 1998 ] in the context of nonlinear programming ; we borrow the term “ coherence ” from the more recent paper of Zhou et al . [ 2017b ] who used the term “ variational coherence ” for a stronger variant of the above definition . We should also note here that the set of solutions of a coherent problem does not need to be convex : for instance , if Player 1 controls ( x , y ) , and the objective function is f ( x , y ) = x2y2 ( i.e. , Player 2 has no impact in the game ) , the set of solutions is the non-convex set X ∗ = { ( x , y ) : x = 0 or y = 0 } . Moreover , regarding the distinction between coherence and strict coherence , we show in Appendix A that ( SP ) is strictly coherent when f is strictly convex-concave . At the other end of the spectrum , typical examples of problems that are null-coherent are bilinear objectives with an interior solution : for instance , f ( x1 , x2 ) = x1x2 with x1 , x2 ∈ [ −1 , 1 ] has 〈g ( x ) , x〉 = x1x2 − x2x1 = 0 for all x1 , x2 ∈ [ −1 , 1 ] , so it is null-coherent . Finally , neither strict , nor null coherence imply a unique solution to ( SP ) , a property which is particularly relevant for GANs ( the first example above is strictly coherent , but does not admit a unique solution ) . 3 Mirror descent The method . Motivated by its prolific success in convex programming , our starting point will be the well-known mirror descent ( MD ) method of Nemirovski & Yudin [ 1983 ] , suitably adapted to our saddle-point context . Several variants of the method exist , ranging from dual averaging [ Nesterov , 2009 ] to follow-the-regularized-leader ; for a survey , we refer the reader to Bubeck [ 2015 ] . The basic idea of mirror descent is to generate a new state variable x+ from some starting state x by taking a “ mirror step ” along a gradient-like vector y . To do this , let h : X → be a continuous and K-strongly convex distance-generating function ( DGF ) on X , i.e. , h ( tx + ( 1 − t ) x′ ) ≤ th ( x ) + ( 1 − t ) h ( x′ ) − 12 Kt ( 1 − t ) ‖x ′ − x‖2 , ( 3.1 ) for all x , x′ ∈ X and all t ∈ [ 0 , 1 ] . In terms of smoothness ( and in a slight abuse of notation ) , we also assume that the subdifferential of h admits a continuous selection , i.e. , a continuous function ∇h : dom ∂h→ Y such that ∇h ( x ) ∈ ∂h ( x ) for all x ∈ dom ∂h.2 Then , following Bregman [ 1967 ] , h generates a pseudo-distance on X via the relation D ( p , x ) = h ( p ) − h ( x ) − 〈∇h ( x ) , p − x〉 for all p ∈ X , x ∈ dom ∂h . ( 3.2 ) This pseudo-distance is known as the Bregman divergence . As we show in Appendix B , we have D ( p , x ) ≥ 12 K‖x − p‖2 , so the convergence of a sequence Xn to some target point p can be verified by showing that D ( p , Xn ) → 0 . On the other hand , D ( p , x ) typically fais to be symmetric and/or satisfy 1To see this , simply note that f ( x1 , x2 ) is multi-modal in x2 for certain values of x1 ( e.g. , when x1 = 1/2 ) . 2Recall here that the subdifferential of h at x ∈ X is defined as ∂h ( x ) ≡ { y ∈ Y : h ( x′ ) ≥ h ( x ) + 〈y , x′ − x〉 for all x ∈ V } , with the standard convention h ( x ) = ∞ for all x ∈ V \ X . Algorithm 1 : mirror descent ( MD ) for saddle-point problems Require : K-strongly convex regularizer h : X → , step-size sequence γn > 0 1 : choose X ∈ dom ∂h # initialization 2 : for n = 1 , 2 , . . . do 3 : oracle query at X returns g # gradient feedback 4 : set X ← PX ( −γng ) # new state 5 : end for 6 : return X the triangle inequality , so it is not a true distance function per se . Moreover , the level sets of D ( p , x ) may fail to form a neighborhood basis of p , so the convergence of Xn to p does not necessarily imply that D ( p , Xn ) → 0 ; we provide an example of this behavior in Appendix B . For technical reasons , it will be convenient to assume that such phenomena do not occur , i.e. , that D ( p , Xn ) → 0 whenever Xn → p. This mild regularity condition is known in the literature as “ Bregman reciprocity ” [ Chen & Teboulle , 1993 ; Alvarez et al. , 2004 ; Mertikopoulos & Staudigl , 2018 ; Bravo et al. , 2018 ] , and it will be our standing assumption in what follows ( note also that it holds trivially for both Examples 3.1 and 3.2 below ) . Now , as with standard Euclidean distances , the Bregman divergence generates an associated proxmapping defined as Px ( y ) = arg min x′∈X { 〈y , x − x′〉 + D ( x′ , x ) } for all x ∈ dom ∂h , y ∈ Y . ( 3.3 ) In analogy with the Euclidean case ( discussed below ) , the prox-mapping ( 3.3 ) produces a feasible point x+ = Px ( y ) by starting from x ∈ dom ∂h and taking a step along a dual ( gradient-like ) vector y ∈ Y . In this way , we obtain the mirror descent ( MD ) algorithm Xn+1 = PXn ( −γnĝn ) , ( MD ) where γn is a variable step-size sequence and ĝn is the calculated value of the gradient vector g ( Xn ) at the n-th stage of the algorithm ( for a pseudocode implementation , see Section 3 ) . For concreteness , two widely used examples of prox-mappings are as follows : Example 3.1 ( Euclidean projections ) . When X is endowed with the L2 norm ‖·‖2 , the archetypal prox-function is the ( square of the ) norm itself , i.e. , h ( x ) = 12‖x‖22 . In that case , D ( p , x ) = 1 2‖x − p‖2 and the induced prox-mapping is Px ( y ) = Π ( x + y ) , ( 3.4 ) with Π ( x ) = arg minx′∈X ‖x′ − x‖2 denoting the ordinary Euclidean projection onto X . Example 3.2 ( Entropic regularization ) . When X is a d-dimensional simplex , a widely used DGF is the ( negative ) Gibbs–Shannon entropy h ( x ) = ∑d j=1 x j log x j . This function is 1-strongly convex with respect to the L1 norm [ Shalev-Shwartz , 2011 ] and the associated pseudo-distance is the Kullback–Leibler divergence DKL ( p , x ) = ∑d j=1 p j log ( p j/x j ) ; in turn , this yields the prox-mapping Px ( y ) = ( x j exp ( y j ) ) dj=1∑d j=1 x j exp ( y j ) for all x ∈ X ◦ , y ∈ Y . ( 3.5 ) The update rule x← Px ( y ) is known in the literature as the multiplicative weights ( MW ) algorithm [ Arora et al. , 2012 ] , and is one of the centerpieces for learning in games [ Fudenberg & Levine , 1998 ; Freund & Schapire , 1999 ; Cohen et al. , 2017 ] , adversarial bandits [ Auer et al. , 1995 ] , etc . Regarding the gradient input sequence ĝn of ( MD ) , we assume that it is obtained by querying a first-order oracle which outputs an estimate of g ( Xn ) when called at Xn . This oracle could be either perfect , returning ĝn = g ( Xn ) for all n , or imperfect , providing noisy gradient estimations.3 By that token , we will make the following standard assumptions for the gradient feedback sequence ĝn [ Nesterov , 2007 ; Nemirovski et al. , 2009 ; Juditsky et al. , 2011 ] : a ) Unbiasedness : [ ĝn |Fn ] = g ( Xn ) . b ) Finite mean square error : [ ‖ĝn‖2∗ |Fn ] ≤ G2 for some finite G ≥ 0 . ( 3.6 ) 3The reason for this is that , depending on the application at hand , gradients might be difficult to compute directly e.g. , because they require huge amounts of data , the calculation of an unknown expectation , etc . In the above , ‖y‖∗ ≡ sup { 〈y , x〉 : x ∈ V , ‖x‖ ≤ 1 } denotes the dual norm on Y while Fn represents the history ( natural filtration ) of the generating sequence Xn up to stage n ( inclusive ) . Since ĝn is generated randomly from Xn at stage n , it is obviously not Fn-measurable , i.e. , ĝn = g ( Xn ) + Un+1 , where Un is an adapted martingale difference sequence with [ ‖Un+1‖2∗ |Fn ] ≤ σ2 for some finite σ ≥ 0 . Clearly , when σ = 0 , we recover the exact gradient feedback framework ĝn = g ( Xn ) . Convergence analysis . When ( SP ) is convex-concave , it is customary to take as the output of ( MD ) the so-called ergodic average X̄n = ∑n k=1 γkXk∑n k=1 γk , ( 3.7 ) or some other average of the sequence Xn where the objective is sampled . The reason for this is that convexity guarantees – via Jensen ’ s inequality and gradient monotonicity – that a regret-based analysis of ( MD ) can lead to explicit rates for the convergence of X̄n to the solution set of ( SP ) [ Nemirovski , 2004 ; Nesterov , 2007 ] . However , when the problem is not convex-concave , the standard proof techniques for establishing convergence of the method ’ s ergodic average no longer apply ; instead , we need to examine the convergence properties of the generating sequence Xn of ( MD ) directly . With all this in mind , our main result for ( MD ) may be stated is as follows : Theorem 3.1 . Suppose that ( MD ) is run with a gradient oracle satisfying ( 3.6 ) and a variable step-size sequence γn such that ∑∞ n=1 γn = ∞ . Then : a ) If f is strictly coherent and ∑∞ n=1 γ 2 n < ∞ , Xn converges ( a.s. ) to a solution of ( SP ) . b ) If f is null-coherent , the sequence [ D ( x∗ , Xn ) ] is non-decreasing for every solution x∗ of ( SP ) . This result establishes an important dichotomy between strict and null coherence : in strictly coherent problems , Xn is attracted to the solution set of ( SP ) ; in null-coherent problems , Xn drifts away and cycles without converging . In particular , this dichotomy leads to the following immediate corollaries : Corollary 3.2 . Suppose that f is strictly convex-concave . Then , with assumptions as above , Xn converges ( a.s. ) to the ( necessarily unique ) solution of ( SP ) . Corollary 3.3 . Suppose that f is bilinear and admits an interior saddle-point x∗ ∈ X ◦ . If X1 , x∗ and ( MD ) is run with exact gradient input ( σ = 0 ) , we have limn→∞ D ( x∗ , Xn ) > 0 . Since bilinear models include all finite two-player , zero-sum games , Corollary 3.3 also encapsulates the non-convergence results of Daskalakis et al . [ 2018 ] and Bailey & Piliouras [ 2018 ] for gradient descent and FTRL respectively ( for a more comprehensive formulation , see Proposition C.3 in Appendix C ) . The failure of ( MD ) to converge in this case is due to the fact that , witout a mitigating mechanism in place , a “ blind ” first-order step could overshoot and spiral outwards , even with a vanishing step-size . This becomes even more pronounced in GANs where it can lead to mode collapse and/or cycles between different modes ; the next two sections address precisely these issues . 4 Extra-gradient analysis The method . In convex-concave problems , taking an average of the algorithm ’ s generated samples as in ( 3.7 ) may resolve cycling phenomena by inducing an auxiliary sequence that gravitates towards the “ center of mass ” of the driving sequence Xn ( which orbits interior solutions ) . However , this technique can not be employed in problems that are not convex-concave because the structure of f can not be leveraged to establish convergence of the ergodic average of the process . In view of this , we replace averaging with an optimistic “ extra-gradient ” step which uses the obtained information to amortize the next prox step ( possibly by exiting the convex hull of generated states ) . The seed of this “ extra-gradient ” idea dates back to Korpelevich [ 1976 ] and Nemirovski [ 2004 ] , and has since found wide applications in optimization theory and beyond – for a survey , see Bubeck [ 2015 ] and references therein . In a nutshell , given a state x , the extra-gradient method first generates an intermediate , “ waiting ” state x̂ = Px ( −γg ( x ) ) by taking a prox step as usual . However , instead of continuing from x̂ , the method samples g ( x̂ ) and goes back to the original state x in order to generate a new state x+ = Px ( −γg ( x̂ ) ) . Based on this heuristic , we obtain the optimistic mirror descent ( OMD ) algorithm Xn+1/2 = PXn ( −γnĝn ) Xn+1 = PXn ( −γnĝn+1/2 ) ( OMD ) Algorithm 2 : optimistic mirror descent ( OMD ) for saddle-point problems Require : K-strongly convex regularizer h : X → , step-size sequence γn > 0 1 : choose X ∈ dom ∂h # initialization 2 : for n = 1 , 2 , . . . do 3 : oracle query at X returns g # gradient feedback 4 : set X+ ← PX ( −γng ) # waiting state 5 : oracle query at X+ returns g+ # gradient feedback 6 : set X ← PX ( −γng+ ) # new state 7 : end for 8 : return X where , in obvious notation , ĝn and ĝn+1/2 represent gradient oracle queries at the incumbent and intermediate states Xn and Xn+1/2 respectively . For a pseudocode implementation , see Algorithm 2 ; see also Rakhlin & Sridharan [ 2013 ] and Daskalakis et al . [ 2018 ] for a variant of the method with a “ momentum ” step , and Gidel et al . [ 2018 ] for a gradient reuse mechanism that replaces a second oracle call with a past gradient . Convergence analysis . In his original analysis , Nemirovski [ 2004 ] considered the ergodic average ( 3.7 ) of the algorithm ’ s iterates and established an O ( 1/n ) convergence rate in monotone problems . However , as we explained above , even though this kind of averaging is helpful in convex-concave problems , it does not provide any tangible benefits beyond this class : in more general problems , Xn appears to be the most natural solution candidate . Our first result below justifies this choice in the class of coherent problems : Theorem 4.1 . Suppose that ( SP ) is coherent and g is L-Lipschitz continuous . If ( OMD ) is run with exact gradient input ( σ = 0 ) and γn such that 0 < infn γn ≤ supn γn < K/L , the sequence Xn converges monotonically to a solution x∗ of ( SP ) , i.e. , D ( x∗ , Xn ) decreases monotonically to 0 . Corollary 4.2 . Suppose that f is bilinear . If ( OMD ) is run with assumptions as above , the sequence Xn converges monotonically to a solution of ( SP ) . Theorem 4.1 includes as a special case the analysis of Facchinei & Pang [ 2003 , Theorem 12.1.11 ] for optimistic gradient descent and the more recent results of Noor et al . [ 2011 ] for pseudo-monotone problems . Importantly , Theorem 4.1 shows that the extra-gradient step plays a crucial role in stabilizing ( MD ) : not only does ( OMD ) converge in problems where ( MD ) provably fails , but this convergence is , in fact , monotonic . In other words , at each iteration , ( OMD ) comes closer to a solution of ( SP ) , whereas ( MD ) may spiral outwards , ultimately converging to a limit cycle . This phenomenon is seen clearly in Fig . 1 , and also in the detailed analysis of Appendix C. Of course , except for very special cases , the monotonic convergence of Xn can not hold when the gradient input to ( OMD ) is imperfect : a single “ bad ” sample of ĝn would suffice to throw Xn off-track . In this case , we have : Theorem 4.3 . Suppose that ( SP ) is strictly coherent and ( OMD ) is run with a gradient oracle satisfying ( 3.6 ) and a variable step-size sequence γn such that ∑∞ n=1 γn = ∞ and ∑∞ n=1 γ 2 n < ∞ . Then , with probability 1 , Xn converges to a solution of ( SP ) . It is worth noting here that the step-size policy in Theorem 4.3 is different than that of Theorem 4.1 . This is due to a ) the lack of randomness ( which obviates the summability requirement ∑∞ n=1 γ 2 n < ∞ in Theorem 4.1 ) ; and b ) the lack of Lipschitz continuity assumption ( which , in the case of Theorem 4.1 guarantees monotonic decrease at each step , provided the step-size is not too big ) . Importantly , the maximum allowable step-size is also controlled by the strong convexity modulus of h , suggesting that the choice of distance-generating function can be fine-tuned further to allow for more aggressive step-size policies – a key benefit of mirror descent methods . 5 Experimental results Gaussian mixture models . For the experimental validation of our theoretical results , we began by evaluating the extra-gradient add-on in a highly multi-modal mixture of 16 Gaussians arranged in a 4 × 4 grid as in Metz et al . [ 2017 ] . The generator and discriminator have 6 fully connected layers with 384 neurons and Relu activations ( plus an additional layer for data space projection ) , and the generator generates 2-dimensional vectors . The output after { 4000 , 8000 , 12000 , 16000 , 20000 } iterations is shown in Fig . 2 . The networks were trained with RMSprop [ Tieleman & Hinton , 2012 ] and Adam [ Kingma & Ba , 2014 ] , and the results are compared to the corresponding extra-gradient variant ( for an explicit pseudocode representation in the case of Adam , see Daskalakis et al . [ 2018 ] and Appendix E ) . Learning rates and hyperparameters were chosen by an inspection of grid search results so as to enable a fair comparison between each method and its look-ahead version . Overall , the different optimization strategies without look-ahead exhibit mode collapse or oscillations throughout the training period ( we ran all models for at least 20000 iterations in order to evaluate the hopping behavior of the generator ) . In all cases , the extra-gradient add-on performs consistently better in learning the multi-modal distribution and greatly reduces occurrences of oscillatory behavior . Experiments with standard datasets . In our experiments with Gaussian mixture models ( GMMs ) , the most promising training method was Adam with an extra-gradient step ( a concrete pseudocode implementation is provided in Appendix E ) . Motivated by this , we trained a Wasserstein-GAN on the CelebA and CIFAR-10 datasets using Adam , both with and without an extra-gradient step . The architecture employed was a standard DCGAN ; hyperparameters and network architecture details may be found in Appendix E. Subsequently , to quantify the gains of the extra-gradient step , we employed the widely used inception score and Fréchet distance metrics , for which we report the results in Fig . 3 . Under both metrics , the extra-gradient add-on provides consistently higher scores after an initial warm-up period ( and is considerably more stable ) . For visualization purposes , we also present in Fig . 4 an ensemble of samples generated at the end of the training period . Overall , the generated samples provide accurate feature representation and low distortion ( especially in CelebA ) . 6 Conclusions Our results suggest that the implementation of an optimistic , extra-gradient step is a flexible add-on that can be easily attached to a wide variety of GAN training methods ( RMSProp , Adam , SGA , etc . ) , and provides noticeable gains in performance and stability . From a theoretical standpoint , the dichotomy between strict and null coherence provides a justification of why this is so : optimism eliminates cycles and , in so doing , stabilizes the method . We find this property particularly appealing because it paves the way to a local analysis with provable convergence guarantees in multi-modal settings , and beyond zero-sum games ; we intend to examine this question in future work . Acknowledgments P. Mertikopoulos was partially supported by the French National Research Agency ( ANR ) grant ORACLESS ( ANR–16–CE33–0004–01 ) , the CNRS PEPS program under grant MixedGAN , and the FMJH PGMO program under grant HEAVY.NET . G. Piliouras was partially supported by SUTD grant SRG ESD 2015 097 , MOE AcRF Tier 2 grant 2016-T2-1-170 , grant PIE-SGP-AI-2018-01 and NRF fellowship NRF-NRFF2018-07 . This work was partly funded by the deep learning 2.0 program at A * STAR . References Felipe Alvarez , Jérôme Bolte , and Olivier Brahic . Hessian Riemannian gradient flows in convex programming . SIAM Journal on Control and Optimization , 43 ( 2 ) :477–501 , 2004 . Martín Arjovsky , Soumith Chintala , and Léon Bottou . Wasserstein generative adversarial networks . In Proceedings of the 34th International Conference on Machine Learning , ICML 2017 , Sydney , NSW , Australia , 6-11 August 2017 , pp . 214–223 , 2017 . URL http : //proceedings.mlr.press/v70/arjovsky17a.html . Sanjeev Arora , Elad Hazan , and Satyen Kale . The multiplicative weights update method : A meta-algorithm and applications . Theory of Computing , 8 ( 1 ) :121–164 , 2012 . Peter Auer , Nicolò Cesa-Bianchi , Yoav Freund , and Robert E. Schapire . Gambling in a rigged casino : The adversarial multi-armed bandit problem . In Proceedings of the 36th Annual Symposium on Foundations of Computer Science , 1995 . James P Bailey and Georgios Piliouras . Multiplicative weights update in zero-sum games . In Proceedings of the 2018 ACM Conference on Economics and Computation , pp . 321–338 . ACM , 2018 . James P Bailey and Georgios Piliouras . Multi-agent learning in network zero-sum games is a Hamiltonian system . In Int . Conference on Autonomous Agents and Multiagent Systems ( AAMAS ) , 2019 . David Balduzzi , Sébastien Racanière , James Martens , Jakob N. Foerster , Karl Tuyls , and Thore Graepel . The mechanics of n-player differentiable games . In ICML , 2018 . Heinz H. Bauschke and Patrick L. Combettes . Convex Analysis and Monotone Operator Theory in Hilbert Spaces . Springer , New York , NY , USA , 2 edition , 2017 . Srinadh Bhojanapalli , Behnam Neyshabur , and Nati Srebro . Global optimality of local search for low rank matrix recovery . In Advances in Neural Information Processing Systems , pp . 3873–3881 , 2016 . Léon Bottou . Online learning and stochastic approximations . On-line learning in neural networks , 17 ( 9 ) :142 , 1998 . Mario Bravo , David S. Leslie , and Panayotis Mertikopoulos . Bandit learning in concave N-person games . In NIPS ’ 18 : Proceedings of the 32nd International Conference on Neural Information Processing Systems , 2018 . Lev M. Bregman . The relaxation method of finding the common point of convex sets and its application to the solution of problems in convex programming . USSR Computational Mathematics and Mathematical Physics , 7 ( 3 ) :200–217 , 1967 . Sébastien Bubeck . Convex optimization : Algorithms and complexity . Foundations and Trends in Machine Learning , 8 ( 3-4 ) :231–358 , 2015 . Gong Chen and Marc Teboulle . Convergence analysis of a proximal-like minimization algorithm using Bregman functions . SIAM Journal on Optimization , 3 ( 3 ) :538–543 , August 1993 . Johanne Cohen , Amélie Héliou , and Panayotis Mertikopoulos . Learning with bandit feedback in potential games . In NIPS ’ 17 : Proceedings of the 31st International Conference on Neural Information Processing Systems , 2017 . Constantinos Daskalakis , Andrew Ilyas , Vasilis Syrgkanis , and Haoyang Zeng . Training GANs with optimism . In ICLR ’ 18 : Proceedings of the 2018 International Conference on Learning Representations , 2018 . Francisco Facchinei and Jong-Shi Pang . Finite-Dimensional Variational Inequalities and Complementarity Problems . Springer Series in Operations Research . Springer , 2003 . Yoav Freund and Robert E. Schapire . Adaptive game playing using multiplicative weights . Games and Economic Behavior , 29:79–103 , 1999 . Drew Fudenberg and David K. Levine . The Theory of Learning in Games , volume 2 of Economic learning and social evolution . MIT Press , Cambridge , MA , 1998 . Rong Ge , Furong Huang , Chi Jin , and Yang Yuan . Escaping from saddle points online stochastic gradient for tensor decomposition . In Conference on Learning Theory , pp . 797–842 , 2015 . Rong Ge , Chi Jin , and Yi Zheng . No spurious local minima in nonconvex low rank problems : A unified geometric analysis . In ICML , 2017 . Gauthier Gidel , Hugo Berard , Pascal Vincent , and Simon Lacoste-Julien . A variational inequality perspective on generative adversarial networks . https : //arxiv.org/pdf/1802.10551.pdf , February 2018 . Ian Goodfellow , Jean Pouget-Abadie , Mehdi Mirza , Bing Xu , David Warde-Farley , Sherjil Ozair , Aaron Courville , and Yoshua Bengio . Generative adversarial nets . In Advances in neural information processing systems , pp . 2672–2680 , 2014 . Ishaan Gulrajani , Faruk Ahmed , Martin Arjovsky , Vincent Dumoulin , and Aaron Courville . Improved training of wasserstein gans . In Advances in Neural Information Processing Systems 30 ( NIPS 2017 ) , pp . 5769–5779 . Curran Associates , Inc. , December 2017 . URL https : //papers.nips.cc/paper/ 7159-improved-training-of-wasserstein-gans . arxiv : 1704.00028 . P. Hall and C. C. Heyde . Martingale Limit Theory and Its Application . Probability and Mathematical Statistics . Academic Press , New York , 1980 . Anatoli Juditsky , Arkadi Semen Nemirovski , and Claire Tauvel . Solving variational inequalities with stochastic mirror-prox algorithm . Stochastic Systems , 1 ( 1 ) :17–58 , 2011 . Diederik Kingma and Jimmy Ba . Adam : A method for stochastic optimization . 12 2014 . G. M. Korpelevich . The extragradient method for finding saddle points and other problems . Èkonom . i Mat . Metody , 12:747–756 , 1976 . Jason D Lee , Max Simchowitz , Michael I Jordan , and Benjamin Recht . Gradient descent only converges to minimizers . In Conference on Learning Theory , pp . 1246–1257 , 2016 . Jason D Lee , Ioannis Panageas , Georgios Piliouras , Max Simchowitz , Michael I Jordan , and Benjamin Recht . First-order methods almost always avoid strict saddle points . Mathematical Programming , 2019 . Panayotis Mertikopoulos and Mathias Staudigl . On the convergence of gradient-like flows with noisy gradient input . SIAM Journal on Optimization , 28 ( 1 ) :163–197 , January 2018 . Panayotis Mertikopoulos , Christos H. Papadimitriou , and Georgios Piliouras . Cycles in adversarial regularized learning . In SODA ’ 18 : Proceedings of the 29th annual ACM-SIAM Symposium on Discrete Algorithms , 2018 . Lars Mescheder , Andreas Geiger , and Sebastian Nowozin . Which training methods for GANs do actually converge ? In ICML , 2018 . Luke Metz , Ben Poole , David Pfau , and Jascha Sohl-Dickstein . Unrolled generative adversarial networks . ICLR Proceedings , 2017 . Arkadi Semen Nemirovski . Prox-method with rate of convergence O ( 1/t ) for variational inequalities with Lipschitz continuous monotone operators and smooth convex-concave saddle point problems . SIAM Journal on Optimization , 15 ( 1 ) :229–251 , 2004 . Arkadi Semen Nemirovski and David Berkovich Yudin . Problem Complexity and Method Efficiency in Optimization . Wiley , New York , NY , 1983 . Arkadi Semen Nemirovski , Anatoli Juditsky , Guanghui Lan , and Alexander Shapiro . Robust stochastic approximation approach to stochastic programming . SIAM Journal on Optimization , 19 ( 4 ) :1574–1609 , 2009 . Yurii Nesterov . Dual extrapolation and its applications to solving variational inequalities and related problems . Mathematical Programming , 109 ( 2 ) :319–344 , 2007 . Yurii Nesterov . Primal-dual subgradient methods for convex problems . Mathematical Programming , 120 ( 1 ) : 221–259 , 2009 . Muhammad Aslam Noor , Eisa Al-Said , Khalida Inayat Noor , and Yonghong Yao . Extragradient methods for solving nonconvex variational inequalities . Journal of Computational and Applied Mathematics , 235 ( 9 ) : 3104–3108 , March 2011 . Gerasimos Palaiopanos , Ioannis Panageas , and Georgios Piliouras . Multiplicative weights update with constant step-size in congestion games : Convergence , limit cycles and chaos . In NIPS ’ 17 : Proceedings of the 31st International Conference on Neural Information Processing Systems , 2017 . Ioannis Panageas and Georgios Piliouras . Gradient descent only converges to minimizers : Non-isolated critical points and invariant regions . In Innovations of Theoretical Computer Science ( ITCS ) , 2017 . Christos Papadimitriou and Georgios Piliouras . From nash equilibria to chain recurrent sets : An algorithmic solution concept for game theory . Entropy , 20 ( 10 ) , 2018 . ISSN 1099-4300. doi : 10.3390/e20100782 . URL http : //www.mdpi.com/1099-4300/20/10/782 . Georgios Piliouras and Jeff S Shamma . Optimization despite chaos : Convex relaxations to complex limit sets via poincaré recurrence . In Proceedings of the twenty-fifth annual ACM-SIAM symposium on Discrete algorithms , pp . 861–873 . SIAM , 2014 . Georgios Piliouras , Carlos Nieto-Granda , Henrik I. Christensen , and Jeff S. Shamma . Persistent patterns : Multi-agent learning beyond equilibrium and utility . In AAMAS , pp . 181–188 , 2014 . A. Radford , L. Metz , and S. Chintala . Unsupervised Representation Learning with Deep Convolutional Generative Adversarial Networks . In International Conference on Learning Representations ( ICLR ) , 2016 . Alexander Rakhlin and Karthik Sridharan . Optimization , learning , and games with predictable sequences . In NIPS ’ 13 : Proceedings of the 26th International Conference on Neural Information Processing Systems , 2013 . Ralph Tyrrell Rockafellar . Convex Analysis . Princeton University Press , Princeton , NJ , 1970 . Tim Salimans , Ian J. Goodfellow , Wojciech Zaremba , Vicki Cheung , Alec Radford , and Xi Chen . Improved techniques for training gans . In NIPS , pp . 2234–2242 , 2016 . Yuzuru Sato , Eizo Akiyama , and J. Doyne Farmer . Chaos in learning a simple two-person game . Proceedings of the National Academy of Sciences , 99 ( 7 ) :4748–4751 , 2002. doi : 10.1073/pnas.032086299 . URL http : //www.pnas.org/content/99/7/4748.abstract . Shai Shalev-Shwartz . Online learning and online convex optimization . Foundations and Trends in Machine Learning , 4 ( 2 ) :107–194 , 2011 . Ju Sun , Qing Qu , and John Wright . Complete dictionary recovery over the sphere i : Overview and the geometric picture . IEEE Transactions on Information Theory , 63 ( 2 ) :853–884 , 2017a . Ju Sun , Qing Qu , and John Wright . Complete dictionary recovery over the sphere ii : Recovery by riemannian trust-region method . IEEE Transactions on Information Theory , 63 ( 2 ) :885–914 , 2017b . T. Tieleman and G. Hinton . Lecture 6.5 - rmsprop , coursera : Neural networks for machine learning . 2012 . Abhay Yadav , Sohil Shah , Zheng Xu , David Jacobs , and Tom Goldstein . Stabilizing adversarial nets with predic- tion methods . In ICLR ’ 18 : Proceedings of the 2018 International Conference on Learning Representations , 2018 . Zhengyuan Zhou , Panayotis Mertikopoulos , Nicholas Bambos , Stephen Boyd , and Peter W. Glynn . On the convergence of mirror descent beyond stochastic convex programming . https : //arxiv.org/abs/1706 . 05681 , 2017a . Zhengyuan Zhou , Panayotis Mertikopoulos , Nicholas Bambos , Stephen Boyd , and Peter W. Glynn . Stochastic mirror descent for variationally coherent optimization problems . In NIPS ’ 17 : Proceedings of the 31st International Conference on Neural Information Processing Systems , 2017b . A Coherent saddle-point problems We begin our discussion with some basic results on coherence : Proposition A.1 . If f is convex-concave , ( SP ) is coherent . In addition , if f is strictly convex-concave , ( SP ) is strictly coherent . Proof . Let x∗ be a solution point of ( SP ) . Since f is convex-concave , first-order optimality gives 〈g1 ( x∗1 , x∗2 ) , x1 − x∗1〉 = 〈∇x1 f ( x∗1 , x∗2 ) , x1 − x∗1〉 ≥ 0 , ( A.1a ) and 〈g2 ( x∗1 , x∗2 ) , x2 − x∗2〉 = 〈−∇x2 f ( x∗1 , x∗2 ) , x2 − x∗2〉 ≥ 0 . ( A.1b ) Combining the two , we readily obtain the ( Stampacchia ) variational inequality 〈g ( x∗ ) , x − x∗〉 ≥ 0 for all x ∈ X . ( A.2 ) In addition to the above , the fact that f is convex-concave also implies that g ( x ) is monotone in the sense that 〈g ( x′ ) − g ( x ) , x′ − x〉 ≥ 0 ( A.3 ) for all x , x′ ∈ X [ Bauschke & Combettes , 2017 ] . Thus , setting x′ ← x∗ in ( A.3 ) and invoking ( SVI ) , we get 〈g ( x ) , x − x∗〉 ≥ 〈g ( x∗ ) , x − x∗〉 ≥ 0 , ( A.4 ) i.e. , ( MVI ) is satisfied . To establish the converse implication , focus for concreteness on the minimizer , and note that ( MVI ) implies that 〈g1 ( x ) , x1 − x∗1〉 ≥ 0 for all x1 ∈ X1 . ( A.5 ) Now , if we fix some x1 ∈ X1 and consider the function φ ( t ) = f ( x∗1 + t ( x1 − x∗1 ) , x∗2 ) , the inequality ( A.5 ) yields φ′ ( t ) = 〈g ( x∗1 + t ( x1 − x∗1 ) , x∗2 ) , x1 − x∗1〉 = 1 t 〈g ( x∗1 + t ( x1 − x∗1 ) , x∗2 ) , x∗1 + t ( x1 − x∗1 − x∗1〉 ≥ 0 , ( A.6 ) for all t ∈ [ 0 , 1 ] . This implies that φ is nondecreasing , so f ( x1 , x∗2 ) = φ ( 1 ) ≥ φ ( 0 ) = f ( x∗1 , x∗2 ) . The maximizing component follows similarly , showing that x∗ is a solution of ( SP ) and , in turn , establishing that ( SP ) is coherent . For the strict part of the claim , the same line of reasoning shows that if 〈g ( x ) , x − x∗〉 = 0 for some x that is not a saddle-point of f , the function φ ( t ) defined above must be constant on [ 0 , 1 ] , indicating in turn that f can not be strictly convex-concave , a contradiction . We proceed to show that the solution set of a coherent saddle-point problem is closed ( we will need this regularity result in the convergence analysis of Appendix C ) : Lemma A.2 . Let X ∗ denote the solution set of ( SP ) . If ( SP ) is coherent , X ∗ is closed . Proof . Let x∗n , n = 1 , 2 , . . . , be a sequence of solutions of ( SP ) converging to some limit point x ∗ ∈ X . To show that X ∗ is closed , it suffices to show that x∗ ∈ X . Indeed , given that ( SP ) is coherent , every solution thereof satisfies ( MVI ) , so we have 〈g ( x ) , x−x∗n〉 ≥ 0 for all x ∈ X . With x∗n → x∗ as n→ ∞ , it follows that 〈g ( x ) , x − x∗〉 = lim n→∞ 〈g ( x ) , x − x∗n〉 ≥ 0 for all x ∈ X , ( A.7 ) i.e. , x∗ satisfies ( MVI ) . By coherence , this implies that x∗ is a solution of ( SP ) , as claimed . B Properties of the Bregman divergence In this appendix , we provide some auxiliary results and estimates that are used throughout the convergence analysis of Appendix C. Some of the results we present here ( or close variants thereof ) are not new [ see e.g. , Nemirovski et al. , 2009 ; Juditsky et al. , 2011 ] . However , the hypotheses used to obtain them vary wildly in the literature , so we provide all the necessary details for completeness . To begin , recall that the Bregman divergence associated to a K-strongly convex distance-generating function h : X → is defined as D ( p , x ) = h ( p ) − h ( x ) − 〈∇h ( x ) , p − x〉 ( B.1 ) with ∇h ( x ) denoting a continuous selection of ∂h ( x ) . The induced prox-mapping is then given by Px ( y ) = arg min x′∈X { 〈y , x − x′〉 + D ( x′ , x ) } = arg max x′∈X { 〈y + ∇h ( x ) , x′〉 − h ( x′ ) } ( B.2 ) and is defined for all x ∈ dom ∂h , y ∈ Y ( recall here that Y ≡ V∗ denotes the dual of the ambient vector space V ) . In what follows , we will also make frequent use of the convex conjugate h∗ : Y → of h , defined as h∗ ( y ) = max x∈X { 〈y , x〉 − h ( x ) } . ( B.3 ) By standard results in convex analysis [ Rockafellar , 1970 , Chap . 26 ] , h∗ is differentiable on Y and its gradient satisfies the identity ∇h∗ ( y ) = arg max x∈X { 〈y , x〉 − h ( x ) } . ( B.4 ) For notational convenience , we will also write Q ( y ) = ∇h∗ ( y ) ( B.5 ) and we will refer to Q : Y → X as the mirror map generated by h. All these notions are related as follows : Lemma B.1 . Let h be a distance-generating function on X . Then , for all x ∈ dom ∂h , y ∈ Y , we have : a ) x = Q ( y ) ⇐⇒ y ∈ ∂h ( x ) . ( B.6a ) b ) x+ = Px ( y ) ⇐⇒ ∇h ( x ) + y ∈ ∂h ( x+ ) ⇐⇒ x+ = Q ( ∇h ( x ) + y ) . ( B.6b ) Finally , if x = Q ( y ) and p ∈ X , we have 〈∇h ( x ) , x − p〉 ≤ 〈y , x − p〉 . ( B.7 ) Remark . By ( B.6b ) , we have ∂h ( x+ ) , ∅ , i.e. , x+ ∈ dom ∂h . As a result , the update rule x← Px ( y ) is well-posed , i.e. , it can be iterated in perpetuity . Proof of Lemma B.1 . For ( B.6a ) , note that x solves ( B.3 ) if and only if y − ∂h ( x ) 3 0 , i.e. , if and only if y ∈ ∂h ( x ) . Similarly , comparing ( B.2 ) with ( B.3 ) , it follows that x+ solves ( B.2 ) if and only if ∇h ( x ) + y ∈ ∂h ( x+ ) , i.e. , if and only if x+ = Q ( ∇h ( x ) + y ) . For ( B.7 ) , by a simple continuity argument , it suffices to show that the inequality holds for interior p ∈ X ◦ . To establish this , let φ ( t ) = h ( x + t ( p − x ) ) − [ h ( x ) + 〈y , x + t ( p − x ) 〉 ] . ( B.8 ) Since h is strongly convex and y ∈ ∂h ( x ) by ( B.6a ) , it follows that φ ( t ) ≥ 0 with equality if and only if t = 0 . Since ψ ( t ) = 〈∇h ( x + t ( p − x ) ) − y , p − x〉 is a continuous selection of subgradients of φ and both φ and ψ are continuous on [ 0 , 1 ] , it follows that φ is continuously differentiable with φ′ = ψ on [ 0 , 1 ] . Hence , with φ convex and φ ( t ) ≥ 0 = φ ( 0 ) for all t ∈ [ 0 , 1 ] , we conclude that φ′ ( 0 ) = 〈∇h ( x ) − y , p − x〉 ≥ 0 , which proves our assertion . We continue with some basic bounds on the Bregman divergence before and after a prox step . The basic ingredient for these bounds is a generalization of the ( Euclidean ) law of cosines which is known in the literature as the “ three-point identity ” [ Chen & Teboulle , 1993 ] : Lemma B.2 . Let h be a distance-generating function on X . Then , for all p ∈ X and all x , x′ ∈ dom ∂h , we have D ( p , x′ ) = D ( p , x ) + D ( x , x′ ) + 〈∇h ( x′ ) − ∇h ( x ) , x − p〉 . ( B.9 ) Proof . By definition , we have : D ( p , x′ ) = h ( p ) − h ( x′ ) − 〈∇h ( x′ ) , p − x′〉 D ( p , x ) = h ( p ) − h ( x ) − 〈∇h ( x ) , p − x〉 D ( x , x′ ) = h ( x ) − h ( x′ ) − 〈∇h ( x′ ) , x − x′〉 . ( B.10 ) Our claim then follows by adding the last two lines and subtracting the first . With this identity at hand , we have the following series of upper and lower bounds : Proposition B.3 . Let h be a K-strongly convex distance-generating function on X , fix some p ∈ X , and let x+ = Px ( y ) for x ∈ dom ∂h , y ∈ Y . We then have : D ( p , x ) ≥ K 2 ‖x − p‖2 . ( B.11a ) D ( p , x+ ) ≤ D ( p , x ) − D ( x+ , x ) + 〈y , x+ − p〉 ( B.11b ) ≤ D ( p , x ) + 〈y , x − p〉 + 1 2K ‖y‖2∗ ( B.11c ) Proof of ( B.11a ) . By the strong convexity of h , we get h ( p ) ≥ h ( x ) + 〈∇h ( x ) , p − x〉 + K 2 ‖p − x‖2 ( B.12 ) so ( B.11a ) follows by gathering all terms involving h and recalling the definition of D ( p , x ) . Proof of ( B.11b ) and ( B.11c ) . By the three-point identity ( B.9 ) , we readily obtain D ( p , x ) = D ( p , x+ ) + D ( x+ , x ) + 〈∇h ( x ) − ∇h ( x+ ) , x+ − p〉 . ( B.13 ) In turn , this gives D ( p , x+ ) = D ( p , x ) − D ( x+ , x ) + 〈∇h ( x+ ) − ∇h ( x ) , x+ − p〉 ≤ D ( p , x ) − D ( x+ , x ) + 〈y , x+ − p〉 , ( B.14 ) where , in the last step , we used ( B.7 ) and the fact that x+ = Px ( y ) , so ∇h ( x ) + y ∈ ∂h ( x+ ) . The above is just ( B.11b ) , so the first part of our proof is complete . For ( B.11c ) , the bound ( B.14 ) gives D ( p , x+ ) ≤ D ( p , x ) + 〈y , x − p〉 + 〈y , x+ − x〉 − D ( x+ , x ) . ( B.15 ) Therefore , by Young ’ s inequality [ Rockafellar , 1970 ] , we get 〈y , x+ − x〉 ≤ K 2 ‖x+ − x‖2 + 1 2K ‖y‖2∗ , ( B.16 ) and hence D ( p , x+ ) ≤ D ( p , x ) + 〈y , x − p〉 + 1 2K ‖y‖2∗ + K 2 ‖x+ − x‖2 − D ( x+ , x ) ≤ D ( p , x ) + 〈y , x − p〉 + 1 2K ‖y‖2∗ , ( B.17 ) with the last step following from Lemma B.1 applied to x in place of p. The first part of Proposition B.3 shows that Xn converges to p if D ( p , Xn ) → 0 . However , as we mentioned in the main body of the paper , the converse may fail : in particular , we could have lim infn→∞ D ( p , Xn ) > 0 even if Xn → p. To see this , let X be the L2 ball of d and take h ( x ) = − √ 1 − ‖x‖22 . Then , a straightforward calculation gives D ( p , x ) = 1 − 〈p , x〉√ 1 − ‖x‖22 ( B.18 ) whenever ‖p‖2 = 1 . The corresponding level sets Lc ( p ) = { x ∈ d : D ( p , x ) = c } of D ( p , · ) are given by the equation 1 − 〈p , x〉 = c √ 1 − ‖x‖22 , ( B.19 ) which admits p as a solution for all c ≥ 0 ( so p belongs to the closure of Lc ( p ) even though D ( p , p ) = 0 by definition ) . As a result , under this distance-generating function , it is possible to have Xn → p even when lim infn→∞ D ( p , Xn ) > 0 ( simply take a sequence Xn that converges to p while remaining on the same level set of D ) . As we discussed in the main body of the paper , such pathologies are discarded by the Bregman reciprocity condition D ( p , Xn ) → 0 whenever Xn → p. ( B.20 ) This condition comes into play at the very last part of the proofs of Theorems 3.1 and 4.1 ; other than that , we will not need it in the rest of our analysis . Finally , for the analysis of the OMD algorithm , we will need to relate prox steps taken along different directions : Proposition B.4 . Let h be a K-strongly convex distance-generating function on X and fix some p ∈ X , x ∈ dom ∂h . Then : a ) For all y1 , y2 ∈ Y , we have : ‖Px ( y2 ) − Px ( y1 ) ‖ ≤ 1 K ‖y2 − y1‖∗ , ( B.21 ) i.e. , Px is ( 1/K ) -Lipschitz . b ) In addition , letting x+1 = Px ( y1 ) and x + 2 = Px ( y2 ) , we have : D ( p , x+2 ) ≤ D ( p , x ) + 〈y2 , x+1 − p〉 + [ 〈y2 , x+2 − x+1 〉 − D ( x+2 , x ) ] ( B.22a ) ≤ D ( p , x ) + 〈y2 , x+1 − p〉 + 1 2K ‖y2 − y1‖2∗ − K 2 ‖x+1 − x‖2 . ( B.22b ) Proof . We begin with the proof of the Lipschitz property of Px . Indeed , for all p ∈ X , ( B.7 ) gives 〈∇h ( x+1 ) − ∇h ( x ) − y1 , x+1 − p〉 ≤ 0 , ( B.23a ) and 〈∇h ( x+2 ) − ∇h ( x ) − y2 , x+2 − p〉 ≤ 0 . ( B.23b ) Therefore , setting p← x+2 in ( B.23a ) , p← x+1 in ( B.23b ) and rearranging , we obtain 〈∇h ( x+2 ) − ∇h ( x+1 ) , x+2 − x+1 〉 ≤ 〈y2 − y1 , x+2 − x+1 〉 . ( B.24 ) By the strong convexity of h , we also have K‖x+2 − x+1 ‖2 ≤ 〈∇h ( x+2 ) − ∇h ( x+1 ) , x+2 − x+1 〉 . ( B.25 ) Hence , combining ( B.24 ) and ( B.25 ) , we get K‖x+2 − x+1 ‖2 ≤ 〈y2 − y1 , x+2 − x+1 〉 ≤ ‖y2 − y1‖∗‖x+2 − x+1 ‖ , ( B.26 ) and our assertion follows . For the second part of our claim , the bound ( B.11b ) of Proposition B.3 applied to x+2 = Px ( y2 ) readily gives D ( p , x+2 ) ≤ D ( p , x ) − D ( x+2 , x ) + 〈y2 , x+2 − p〉 = D ( p , x ) + 〈y2 , x+1 − p〉 + [ 〈y2 , x+2 − x+1 〉 − D ( x+2 , x ) ] ( B.27 ) thus proving ( B.22a ) . To complete our proof , note that ( B.11b ) with p← x+2 gives D ( x+2 , x + 1 ) ≤ D ( x+2 , x ) + 〈y1 , x+1 − x+2 〉 − D ( x+1 , x ) , ( B.28 ) or , after rearranging , D ( x+2 , x ) ≥ D ( x+2 , x+1 ) + D ( x+1 , x ) + 〈y1 , x+2 − x+1 〉 . ( B.29 ) We thus obtain 〈y2 , x+2 − x+1 〉 − D ( x+2 , x ) ≤ 〈y2 − y1 , x+2 − x+1 〉 − D ( x+2 , x+1 ) − D ( x+1 , x ) ≤ ‖y2 − y1‖ 2 ∗ 2K + K 2 ‖x+2 − x+1 ‖2 − K 2 ‖x+2 − x+1 ‖2 − K 2 ‖x+1 − x‖2 ≤ 1 2K ‖y2 − y1‖2∗ − K 2 ‖x+1 − x‖2 , ( B.30 ) where we used Young ’ s inequality and ( B.11a ) in the second inequality . The bound ( B.22b ) then follows by substituting ( B.30 ) in ( B.27 ) . C Convergence analysis of mirror descent We begin by recalling the definition of the mirror descent algorithm . With notation as in the previous section , the algorithm is defined via the recursive scheme Xn+1 = PXn ( −γnĝn ) , ( MD ) where γn is a variable step-size sequence and ĝn is the calculated value of the gradient vector g ( Xn ) at the n-th stage of the algorithm . As we discussed in the main body of the paper , the gradient input sequence ĝn of ( MD ) is assumed to satisfy the standard oracle assumptions a ) Unbiasedness : [ ĝn |Fn ] = g ( Xn ) . b ) Finite mean square : [ ‖ĝn‖2∗ |Fn ] ≤ G2 for some finite G ≥ 0. where Fn represents the history ( natural filtration ) of the generating sequence Xn up to stage n ( inclusive ) . With this preliminaries at hand , our convergence proof for ( MD ) under strict coherence will hinge on the following results : Proposition C.1 . Suppose that ( SP ) is coherent and ( MD ) is run with a gradient oracle satisfying ( 3.6 ) and a variable step-size γn such that ∑∞ n=1 γ 2 n < ∞ . If x∗ ∈ X is a solution of ( SP ) , the Bregman divergence D ( x∗ , Xn ) converges ( a.s. ) to a random variable D ( x∗ ) with [ D ( x∗ ) ] < ∞ . Proposition C.2 . Suppose that ( SP ) is strictly coherent and ( MD ) is run with a gradient oracle satisfying ( 3.6 ) and a step-size γn such that ∑∞ n=1 γn = ∞ and ∑∞ n=1 γ 2 n < ∞ . Then , with probability 1 , there exists a ( possibly random ) solution x∗ of ( SP ) such that lim infn→∞ D ( x∗ , Xn ) = 0 . Proposition C.1 can be seen as a “ dichotomy ” result : it shows that the Bregman divergence is an asymptotic constant of motion , so ( MD ) either converges to a saddle-point x∗ ( if D ( x∗ ) = 0 ) or to some nonzero level set of the Bregman divergence ( with respect to x∗ ) . In this way , Proposition C.1 rules out more complicated chaotic or aperiodic behaviors that may arise in general – for instance , as in the analysis of Palaiopanos et al . [ 2017 ] for the long-run behavior of the multiplicative weights algorithm in two-player games . However , unless this limit value can be somehow predicted ( or estimated ) in advance , this result can not be easily applied . This is the main role of Proposition C.2 : it shows that ( MD ) admits a subsequence converging to a solution of ( SP ) so , by ( B.20 ) , the limit of D ( x∗ , Xn ) must be zero . Our first step is to prove Proposition C.2 . To do this , we first recall the following law of large numbers for L2 martingales : Theorem ( Hall & Heyde , 1980 , Theorem 2.18 ) . Let Yn = ∑n k=1 ζk be a martingale and Tn a nondecreasing sequence such that limn→∞ τn = ∞ . Then , lim n→∞ Yn/τn = 0 ( a.s. ) ( C.1 ) on the set ∑∞ n=1 τ −2 n [ ζ 2 n |Fn−1 ] < ∞ . With this in place , we have : Proof of Proposition C.2 . We begin with the technical observation that the solution set X ∗ of ( SP ) is closed – and hence , compact ( cf . Lemma A.2 in Appendix A ) . Clearly , if X ∗ = X , there is nothing to show ; hence , without loss of generality , we may assume in what follows that X ∗ , X . Assume now ad absurdum that , with positive probability , the sequence Xn generated by ( MD ) admits no limit points in X ∗ . Conditioning on this event , and given that X ∗ is compact , there exists a ( nonempty ) compact set C ⊂ X such that C ∩ X ∗ = ∅ and Xn ∈ C for all sufficiently large n. Moreover , letting p be as in Definition 2.1 , we have 〈g ( x ) , x − p〉 > 0 whenever x ∈ C. Therefore , by the continuity of g and the compactness of X ∗ and C , there exists some a > 0 such that 〈g ( x ) , x − p〉 ≥ a for all x ∈ C. ( C.2 ) To proceed , let Dn = D ( p , Xn ) . Then , by Proposition B.3 , we have Dn+1 = D ( p , PXn ( −γnĝn ) ) ≤ D ( p , Xn ) − γn〈ĝn , Xn − p〉 + γ2n 2K ‖ĝn‖2 = Dn − γn〈g ( Xn ) , Xn − p〉 − γn〈Un+1 , Xn − p〉 + γ2n 2K ‖ĝn‖2∗ ≤ Dn + γnξn+1 + γ2n 2K ‖ĝn‖2∗ , ( C.3 ) where , in the last line , we set Un+1 = ĝn − g ( Xn ) , ξn+1 = −〈Un+1 , Xn − p〉 , and we invoked the assumption that ( SP ) is coherent . Hence , telescoping ( C.3 ) yields the estimate Dn+1 ≤ D1 − n∑ k=1 γk〈g ( Xk ) , Xk − p〉 + n∑ k=1 γkξk+1 + n∑ k=1 γ2k 2K ‖ĝk‖2∗ . ( C.4 ) Subsequently , letting τn = ∑n k=1 γk and using ( C.2 ) , we obtain Dn+1 ≤ D1 − τn a − ∑nk=1 γkξk+1 τn − ( 2K ) −1 ∑n k=1 γ 2 k‖ĝk‖2∗ τn . ( C.5 ) By the unbiasedness hypothesis of ( 3.6 ) for Un , we have [ ξn+1 |Fn ] = 〈 [ Un+1 |Fn ] , Xn − p〉 = 0 ( recall that Xn is Fn-measurable by construction ) . Moreover , since Un is bounded in L2 and γn is ` 2 summable ( by assumption ) , it follows that ∞∑ n=1 γ2n [ ξ 2 n+1 |Fn ] ≤ ∞∑ n=1 γ2n‖Xn − p‖2 [ ‖Un+1‖2∗ |Fn ] ≤ diam ( X ) 2σ2 ∞∑ n=1 γ2n < ∞ . ( C.6 ) Therefore , by the law of large numbers for L2 martingales stated above [ Hall & Heyde , 1980 , Theorem 2.18 ] , we conclude that τ−1n ∑n k=1 γkξk+1 converges to 0 with probability 1 . Finally , for the last term of ( C.4 ) , let S n+1 = ∑n k=1 γ 2 k‖ĝk‖2∗ . Since ĝk is Fn-measurable for all k = 1 , 2 , . . . , n − 1 , we have [ S n+1 |Fn ] = n−1∑ k=1 γ2k‖ĝk‖2∗ + γ2n‖ĝn‖2∗ ∣∣∣∣∣∣∣Fn = S n + γ2n [ ‖ĝn‖2∗ |Fn ] ≥ S n , ( C.7 ) i.e. , S n is a submartingale with respect to Fn . Furthermore , by the law of total expectation , we also have [ S n+1 ] = [ [ S n+1 |Fn ] ] ≤ G2 n∑ k=1 γ2n ≤ G2 ∞∑ k=1 γ2n < ∞ , ( C.8 ) so S n is bounded in L1 . Hence , by Doob ’ s submartingale convergence theorem [ Hall & Heyde , 1980 , Theorem 2.5 ] , we conclude that S n converges to some ( almost surely finite ) random variable S∞ with [ S∞ ] < ∞ , implying in turn that limn→∞ S n+1/τn = 0 ( a.s. ) . Applying all of the above , the estimate ( C.4 ) gives Dn+1 ≤ D1 − aτn/2 for sufficiently large n , so D ( p , Xn ) → −∞ , a contradiction . Going back to our original assumption , this shows that , at least one of the limit points of Xn must lie in X ∗ ( a.s. ) , as claimed . We now turn to the proof of Proposition C.1 : Proof of Proposition C.1 . Let x∗ ∈ X ∗ be a limit point of Xn , as guaranteed by Proposition C.2 , and let Dn = D ( x∗ , Xn ) . Then , by Proposition B.3 , we have Dn+1 = D ( x∗ , PXn ( −γnĝn ) ) ≤ D ( x∗ , Xn ) − γn〈ĝn , Xn − x∗〉 + γ2n 2K ‖ĝn‖2 = Dn − γn〈g ( Xn ) , Xn − x∗〉 − γn〈Un+1 , Xn − x∗〉 + γ2n 2K ‖ĝn‖2∗ ( C.9 ) and hence , for large enough n : Dn+1 ≤ Dn + γnξn+1 + γ2n 2K ‖ĝn‖2∗ , ( C.10 ) where we used the ansatz that 〈g ( Xn ) , Xn − x∗〉 ≤ 0 for sufficiently large n ( to be proved below ) , and , as in the proof of Proposition C.2 , we set Un+1 = ĝn − g ( Xn ) , ξn+1 = −〈Un+1 , Xn − x∗〉 . Thus , conditioning on Fn and taking expectations , we get [ Dn+1 |Fn ] ≤ Dn + [ ξn+1 |Fn ] + γ2n 2K [ ‖ĝn‖2∗ |Fn ] ≤ Dn + G2 2K γ2n , ( C.11 ) where we used the oracle assumptions ( 3.6 ) and the fact that Xn is Fn-measurable ( by definition ) . Now , letting Rn = Dn + ( 2K ) −1G2 ∑∞ k=n γ 2 k , the estimate ( C.10 ) gives [ Rn+1 |Fn ] = [ Dn+1 |Fn ] + G2 2K ∞∑ k=n+1 γ2k ≤ Dn + G2 2K ∞∑ k=n γ2k = Rn , ( C.12 ) i.e. , Rn is an Fn-adapted supermartingale . Since ∑∞ n=1 γ 2 n < ∞ , it follows that [ Rn ] = [ [ Rn |Fn−1 ] ] ≤ [ Rn−1 ] ≤ · · · ≤ [ R1 ] ≤ [ D1 ] + G2 2K ∞∑ n=1 γ2n < ∞ , ( C.13 ) i.e. , Rn is uniformly bounded in L1 . Thus , by Doob ’ s convergence theorem for supermartingales [ Hall & Heyde , 1980 , Theorem 2.5 ] , it follows that Rn converges ( a.s. ) to some finite random variable R∞ with [ R∞ ] < ∞ . In turn , by inverting the definition of Rn , this shows that Dn converges ( a.s. ) to some random variable D ( x∗ ) with [ D ( x∗ ) ] < ∞ , as claimed . It remains to be shown that 〈g ( Xn ) , Xn − x∗〉 ≥ 0 for sufficiently large n. By Definition 2.1 , this amounts to showing that , for all large enough n , Xn lies in a neighborhood U of x∗ such that ( MVI ) holds . Since x∗ has been chosen so that lim inf D ( x∗ , Xn ) = 0 , it follows that , for all ε > 0 , there exists some n0 such that ∑∞ n=n0 γ 2 n < ε and Xn0 ∈ U . Hence , arguing in the same way as in the proof of Theorem 5.2 of Zhou et al . [ 2017a ] , we conclude that ( Xn ∈ U for all n ≥ n0 ) = 1 , implying in turn that 〈g ( Xn ) , Xn − x∗〉 ≥ 0 for all n ≥ n0 . This proves our last claim and concludes our proof . With all this at hand , we are finally in a position to prove our main result for ( MD ) : Proof of Theorem 3.1 ( a ) . Proposition C.2 shows that , with probability 1 , there exists a ( possibly random ) solution x∗ of ( SP ) such that lim infn→∞‖Xn − x∗‖ = 0 and , hence , lim infn→∞ D ( x∗ , Xn ) = 0 ( by Bregman reciprocity ) . Since limn→∞ D ( x∗ , Xn ) exists with probability 1 ( by Proposition C.1 ) , it follows that limn→∞ D ( x∗ , Xn ) = lim infn→∞ D ( x∗ , Xn ) = 0 , i.e. , Xn converges to x∗ . We proceed with the negative result hinted at in the main body of the paper , namely the failure of ( MD ) to converge under null coherence : Proof of Theorem 3.1 ( b ) . The evolution of the Bregman divergence under ( MD ) satisfies the identity D ( x∗ , Xn+1 ) = D ( x∗ , Xn ) + D ( Xn , Xn+1 ) + γn〈ĝn , Xn − x∗〉 = D ( x∗ , Xn ) + D ( Xn , Xn+1 ) + 〈Un+1 , Xn − x∗〉 ( C.14 ) where , in the last line , we used the null coherence assumption 〈g ( x ) , x − x∗〉 = 0 for all x ∈ X . Since D ( Xn , Xn+1 ) ≥ 0 , taking expecations above shows that D ( x∗ , Xn ) is nondecreasing , as claimed . With Theorem 3.1 at hand , the proof of Corollary 3.2 is an immediate consequence of the fact that strictly convex-concave problems satisfy strict coherence ( Proposition A.1 ) . As for Corollary 3.3 , we provide below a more general result for two-player , zero-sum finite games . To state it , let Ai = { 1 , . . . , Ai } , i = 1 , 2 , be two finite sets of pure strategies , and let Xi = ∆ ( Ai ) denote the set of mixed strategies of player i . A finite , two-player zero-sum game is then defined by a matrix M ∈ A1×A2 so that the loss of Player 1 and the reward of Player 2 in the mixed strategy profile x = ( x1 , x2 ) ∈ X are concurrently given by f ( x1 , x2 ) = x > 1 Mx2 ( C.15 ) Then , writing Γ ≡ Γ ( A1 , A2 , M ) for the resulting game , we have : Proposition C.3 . Let Γ be a two-player zero-sum game with an interior Nash equilibrium x∗ . If X1 , x∗ and ( MD ) is run with exact gradient input ( σ2 = 0 ) , we have limn→∞ D ( x∗ , Xn ) > 0 . If , in addition , ∑∞ n=1 γ 2 n < ∞ , limn→∞ D ( x∗ , Xn ) is finite . Remark . Note that non-convergence does not require any summability assumptions on γn . In words , Proposition C.3 states that ( MD ) does not converge in finite zero-sum games with a unique interior equilibrium and exact gradient input : instead , Xn cycles at positive Bregman distance from the game ’ s Nash equilibrium . Heuristically , the reason for this behavior is that , for small γ → 0 , the incremental step Vγ ( x ) = Px ( −γg ( x ) ) − x of ( MD ) is essentially tangent to the level set of D ( x∗ , · ) that passes through x.4 For finite γ > 0 , things are even worse because Vγ ( x ) points noticeably away from x , i.e. , towards higher level sets of D. As a result , the “ best-case scenario ” for ( MD ) is to orbit x∗ ( when γ → 0 ) ; in practice , for finite γ , the algorithm takes small outward steps throughout its runtime , eventually converging to some limit cycle farther away from x∗ . We make this intuition precise below ( for a schematic illustration , see also Fig . 1 above ) : Proof of Proposition C.3 . Write v1 ( x ) = −Mx2 and v2 ( x ) = x > 1 M for the players ’ payoff vectors under the mixed strategy profile x = ( x1 , x2 ) . By construction , we have g ( x ) = − ( v1 ( x ) , v2 ( x ) ) . Furthermore , since x∗ is an interior equilibrium of f , elementary game-theoretic considerations show that v1 ( x∗ ) and v2 ( x∗ ) are both proportional to the constant vector of ones . We thus get 〈g ( x ) , x − x∗〉 = 〈v1 ( x ) , x1 − x∗1〉 + 〈v2 ( x ) , x2 − x∗2〉 = −x > 1 Mx2 + ( x∗1 ) > Mx2 + x > 1 Mx2 − x > 1 Mx∗2 = 0 , ( C.16 ) where , in the last line , we used the fact that x∗ is interior . This shows that f satisfies null coherence , so our claim follows from Theorem 3.1 ( b ) . For our second claim , arguing as above and using ( B.11c ) , we get D ( x∗ , Xn+1 ) ≤ D ( x∗ , Xn ) + γn〈g ( Xn ) , Xn − x∗〉 + γ2n 2K ‖g ( Xn ) ‖2∗ ≤ D ( x∗ , Xn ) + γ2nG 2 2K ( C.17 ) with G = maxx1∈X1 , x2∈X2‖ ( −Mx2 , x > 1 M ) ‖∗ . Telescoping this last bound yields sup n D ( x∗ , Xn ) ≤ D ( x∗ , X1 ) + ∞∑ k=1 γ2nG 2 2K < ∞ , ( C.18 ) 4This observation was also the starting point of Mertikopoulos et al . [ 2018 ] who showed that FTRL in continuous time exhibits a similar cycling behavior in zero-sum games with an interior equilibrium . so D ( x∗ , Xn ) is also bounded from above . Therefore , with D ( x∗ , Xn ) nondecreasing , bounded from above and D ( x∗ , X1 ) > 0 , it follows that limn→∞ D ( x∗ , Xn ) > 0 , as claimed . D Convergence analysis of optimistic mirror descent We now turn to the optimistic mirror descent ( OMD ) algorithm , as defined by the recursion Xn+1/2 = PXn ( −γnĝn ) Xn+1 = PXn ( −γnĝn+1/2 ) ( OMD ) with X1 initialized arbitrarily in dom ∂h , and ĝn , ĝn+1/2 representing gradient oracle queries at the incumbent and intermediate states Xn and Xn+1/2 respectively . The heavy lifting for our analysis is provided by Proposition B.4 , which leads to the following crucial lemma : Lemma D.1 . Suppose that ( SP ) is coherent and g is L-Lipschitz continuous . With notation as above and exact gradient input ( σ = 0 ) , we have D ( p , Xn+1 ) ≤ D ( p , Xn ) − 1 2 ( K − γ 2 nL 2 K ) ‖Xn+1/2 − Xn‖2 , ( D.1 ) with p as in Definition 2.1 . Proof . Substituting x← Xn , y1 ← −γng ( Xn ) , and y2 ← −γng ( Xn+1/2 ) in Proposition B.4 , we obtain the estimate : D ( p , Xn+1 ) ≤ D ( p , Xn ) − γn〈g ( Xn+1/2 ) , Xn+1/2 − p〉 + γ2n 2K ‖g ( Xn+1/2 ) − g ( Xn ) ‖2∗ − K 2 ‖Xn+1/2 − Xn‖2 ≤ D ( p , Xn ) + γ2nL 2 2K ‖Xn+1/2 − Xn‖2 − K 2 ‖Xn+1/2 − Xn‖2 , ( D.2 ) where we used the fact that g is L-Lipschitz and that p is a solution of ( SP ) such that ( MVI ) holds for all x ∈ X . We are now finally in a position to prove Theorem 4.1 ( reproduced below for convenience ) : Theorem . Suppose that ( SP ) is coherent and g is L-Lipschitz continuous . If ( OMD ) is run with exact gradient input and a step-size sequence γn such that 0 < limn→∞ γn ≤ supn γn < K/L , ( D.3 ) the sequence Xn converges monotonically to a solution x∗ of ( SP ) , i.e. , D ( x∗ , Xn ) is non-increasing and converges to 0 . Proof . Let p be a solution of ( SP ) such that ( MVI ) holds for all x ∈ X ( that such a solution exists is a consequence of Definition 2.1 ) . Then , by the stated assumptions for γn , Lemma D.1 yields D ( p , Xn+1 ) ≤ D ( p , Xn ) − 1 2 K ( 1 − α2 ) ‖Xn+1/2 − Xn‖2 , ( D.4 ) where α ∈ ( 0 , 1 ) is such that γ2n < αK/L for all n ( that such an α exists is a consequence of the assumption that supn γn < K/L ) . Now , telescoping ( D.1 ) , we obtain D ( p , Xn+1 ) ≤ D ( p , X1 ) − 1 2 n∑ k=1 K − γ2k L2K ‖Xk+1/2 − Xk‖2 , ( D.5 ) and hence : n∑ k=1 1 − γ2k L2K2 ‖Xk+1/2 − Xk‖2 ≤ 2K D ( p , X1 ) . ( D.6 ) With supn γn < K/L , the above estimate readily yields ∑∞ n=1‖Xn+1/2−Xn‖2 < ∞ , which in turn implies that ‖Xn+1/2 − Xn‖ → 0 as n→ ∞ . By the compactness of X , we further infer that Xn admits an accumulation point x∗ , i.e. , there exists a subsequence nk such that Xnk → x∗ as k → ∞ . Since ‖Xnk+1/2 − Xnk‖ → 0 , this also implies that Xnk+1/2 converges to x ∗ as k → ∞ . Further , by passing to a subsequence if necessary , we may also assume without loss of generality that γnk converges to some limit value γ > 0 . Then , by the Lipschitz continuity of the prox-mapping ( cf . Proposition B.4 ) , we readily obtain x∗ = lim k→∞ Xnk+1/2 = limk→∞ PXnk ( −γnkg ( Xnk ) ) = Px∗ ( −γg ( x ∗ ) ) , ( D.7 ) i.e. , x∗ is a solution of ( SP ) . Since ( MVI ) holds locally around x∗ ( by Definition 2.1 ) , the same reasoning as above shows that D ( x∗ , Xn+1 ) ≤ D ( x∗ , Xn ) − 1 2 K ( 1 − α2 ) ‖Xn+1/2 − Xn‖2 ≤ D ( x∗ , Xn ) , ( D.8 ) for all sufficiently large n. This shows that D ( x∗ , Xn ) is non-decreasing ; since lim infn→∞ D ( x∗ , Xn ) = 0 ( by Bregman reciprocity ) , we ultimately conclude that limn→∞ D ( x∗ , Xn ) = 0 , i.e. , Xn → x∗ . Our last result concerns the convergence of ( OMD ) in strictly coherent problems with a stochastic gradient oracle : Proof of Theorem 4.3 . Our argument hinges on the inequality D ( x∗ , Xn+1 ) ≤ D ( x∗ , Xn ) − γn〈ĝn+1/2 , Xn+1/2 − x∗〉 + γ2n/ ( 2K ) ‖ĝn+1/2 − ĝn‖2∗ ( D.9 ) which is obtained from the two-point estimate ( B.22b ) by substituting x ← x∗ , x1 ← Xn , y1 ← ĝn , x+1 ← Xn+1/2 = PXn ( −γnĝn ) , y2 ← ĝn+1/2 , and x+2 ← Xn = PXn ( −γnĝn+1/2 ) . Then , working as in the proof of Proposition C.1 , we obtain the following estimate for the sequence Dn = D ( x∗ , Xn ) : Dn+1 ≤ Dn − γn〈g ( Xn+1/2 ) , Xn+1/2 − x∗〉 − γn〈U+n+1 , Xn − x∗〉 + γ2n 2K ‖ĝn+1/2 − ĝn‖2∗ ≤ Dn + γnξ+n+1 + γ2n K [ ‖ĝn‖2∗ + ‖ĝ2n+1/2‖∗ ] , ( D.10 ) where U+n+1 = ĝn+1/2 − g ( Xn+1/2 ) denotes the martingale part of ĝn+1/2 and we have set ξ+n+1 = 〈U+n+1 , Xn+1/2 − x∗〉 . Since [ ‖gn‖2∗ | Xn , . . . , X1 ] and [ ‖gn+1/2‖2∗ | Xn+1/2 , . . . , X1 ] are both bounded by G2 , we get the bound [ Dn+1 |Fn ] ≤ Dn + G K γ2n . ( D.11 ) Following the same steps as in the proof of Proposition C.1 , it then follows that Dn converges to some limit value D∞ . To proceed , telescoping ( D.10 ) also yields Dn+1 ≤ D1 − n∑ k=1 γk〈g ( Xk+1/2 ) , Xk+1/2 − x∗〉 + n∑ k=1 γkξ + k+1 + n∑ k=1 γ2k 2K ‖ĝk+1/2 − ĝk‖2∗ . ( D.12 ) Each term in the above bound can be controlled in the same way as the corresponding terms in ( C.4 ) . Thus , repeating the steps in the proof of Proposition C.2 , it follows that there exists a subsequence of Xn+1/2 ( and hence also of Xn ) which converges to x∗ . Our claim then follows by combining the two intermediate results above in the same way as in the proof of Theorem 3.1 ( a ) ; to avoid needless repetition , we omit the details . E Experimental results E.1 Adam with extra-gradient step For most of our experiments , the method that seemed to generate the best results was Adam and its optimistic version [ Daskalakis et al. , 2018 ] ; for a pseudocode iplementation , see Algorithm 3 below . We also noticed empirically that it was more efficient to use two different sets of moment estimates ( mt , vt ) and ( m′t , v ′ t ) for the first and the second gradient steps . We used this algorithm for our experiments with both GMMs and the CelebA/CIFAR-10 datasets . Algorithm 3 : Adam with extra-gradient add-on ( optimistic Adam ) Compute stochastic gradient : ∇θ , t Update biased estimate of 1st momentum : mt = β1mt−1 + ( 1 − β1 ) ∇θ , t Update biased estimate of 2nd momentum : vt = β2vt−1 + ( 1 − β2 ) ∇2θ , t Compute bias corrected 1st moment : m̂t = mt1−βt1 Compute bias corrected 2nd moment : v̂t = vt1−βt2 Perform : θ ′ t = θt−1 − η m̂t√v̂t+ Compute stochastic gradient : ∇θ′ , t Update biased estimate of 1st momentum : m′t = β1m ′ t−1 + ( 1 − β1 ) ∇θ′ , t Update biased estimate of 2nd momentum : v′t = β2v ′ t−1 + ( 1 − β1 ) ∇2θ′ , t Compute bias corrected 1st moment : m̂′t = m′t 1−βt1 Compute bias corrected 2nd moment : v̂′t = v′t 1−βt1 Perform : θt = θt−1 − η′ m̂ ′ t√ v̂′t + Return θt E.2 Experiments with standards datasets In this section we present the results of our image experiments using OMD training techniques . Inception and FID scores obtained by our model during training were reported in Fig . 3 : as can be seen there , the extra-gradient add-on improves the performance of GAN training and efficiently stabilizes the model ; without the extra-gradient step , performance tends to drop noticeably after approximately 100k steps . For ease of comparison , we provide below a collection of samples generated by Adam and optimistic Adam in the CelebA and CIFAR-10 datasets . Especially in the case of CelebA , the generated samples are consistently more representative and faithful to the target data distribution . E.2.1 Network Architecture and hyperparameters For the reproducibility of our experiments , we provide Table 1 and Table 2 the network architectures and the hyperparameters of the GANs that we used . The architecture employed is a standard DCGAN architecture with a 5-layer generator with batchnorm , and an 8-layer discriminator . The generated samples were 32×32×3 RGB images .

Summary: This paper is trying to find a saddle-point of a Lagrangian using mirror descent. Mirror descent based methods use Bregman divergence to encode the convexity and smoothness of objective function beyond the euclidean structure. The main contribution of this paper is adding an extra gradient step to the standard MD, i.e., step 5 in Algorithm 2 as well as stochastic versions. Numerical experiments support their results.