Title: STABLE OPPONENT SHAPING IN DIFFERENTIABLE GAMES

Abstract: A growing number of learning methods are actually differentiable games whose players optimise multiple, interdependent objectives in parallel -from GANs and intrinsic curiosity to multi-agent RL. Opponent shaping is a powerful approach to improve learning dynamics in these games, accounting for player influence on others' updates. Learning with Opponent-Learning Awareness (LOLA) is a recent algorithm that exploits this response and leads to cooperation in settings like the Iterated Prisoner's Dilemma. Although experimentally successful, we show that LOLA agents can exhibit 'arrogant' behaviour directly at odds with convergence. In fact, remarkably few algorithms have theoretical guarantees applying across all (n-player, non-convex) games. In this paper we present Stable Opponent Shaping (SOS), a new method that interpolates between LOLA and a stable variant named LookAhead. We prove that LookAhead converges locally to equilibria and avoids strict saddles in all differentiable games. SOS inherits these essential guarantees, while also shaping the learning of opponents and consistently either matching or outperforming LOLA experimentally.

1. INTRODUCTION:

Problem Setting. While machine learning has traditionally focused on optimising single objectives, generative adversarial nets (GANs) (Goodfellow et al., 2014) have showcased the potential of architectures dealing with multiple interacting goals. They have since then proliferated substantially, including intrinsic curiosity (Pathak et al., 2017) , imaginative agents (Racanière et al., 2017) , synthetic gradients (Jaderberg et al., 2017) , hierarchical reinforcement learning (RL) (Wayne & Abbott, 2014; Vezhnevets et al., 2017) and multi-agent RL in general (Busoniu et al., 2008) . These can effectively be viewed as differentiable games played by cooperating and competing agents -which may simply be different internal components of a single system, like the generator and discriminator in GANs. The difficulty is that each loss depends on all parameters, including those of other agents. While gradient descent on single functions has been widely successful, converging to local minima under rather mild conditions (Lee et al., 2017) , its simultaneous generalisation can fail even in simple two-player, two-parameter zero-sum games. No algorithm has yet been shown to converge, even locally, in all differentiable games. Related Work. Convergence has widely been studied in convex n-player games, see especially Rosen (1965) ; Facchinei & Kanzow (2007) . However, the recent success of non-convex games exemplified by GANs calls for a better understanding of this general class where comparatively little is known. Mertikopoulos & Zhou (2018) recently prove local convergence of no-regreat learning to variationally stable equilibria, though under a number of regularity assumptions.

Conversely, a number of algorithms have been successful in the non-convex setting for restricted classes of games. These include policy prediction in two-player two-action bimatrix games (Zhang & Lesser, 2010) ; WoLF in two-player two-action games (Bowling & Veloso, 2001) ; AWESOME in repeated games (Conitzer & Sandholm, 2007) ; Optimistic Mirror Descent in two-player bilinear zero-sum games (Daskalakis et al., 2018) and Consensus Optimisation (CO) in two-player zerosum games (Mescheder et al., 2017) . An important body of work including Heusel et al. (2017) ; Nagarajan & Kolter (2017) has also appeared for the specific case of GANs.

Working towards bridging this gap, some of the authors recently proposed Symplectic Gradient Adjustment (SGA), see Balduzzi et al. (2018) . This algorithm is provably 'attracted' to stable fixed points while 'repelled' from unstable ones in all differentiable games (n-player, non-convex). Nonetheless, these results are weaker than strict convergence guarantees. Moreover, SGA agents may act against their own self-interest by prioritising stability over individual loss. SGA was also discovered independently by Gemp & Mahadevan (2018) , drawing on variational inequalities.

In a different direction, Learning with Opponent-Learning Awareness (LOLA) (Foerster et al., 2018) modifies the learning objective by predicting and differentiating through opponent learning steps. This is intuitively appealing and experimentally successful, encouraging cooperation in settings like the Iterated Prisoner's Dilemma (IPD) where more stable algorithms like SGA defect. However, LOLA has no guarantees of converging or even preserving fixed points of the game.

Contribution. We begin by constructing the first explicit tandem game where LOLA agents adopt 'arrogant' behaviour and converge to non-fixed points. We pinpoint the cause of failure and show that a natural variant named LookAhead (LA), discovered before LOLA by Zhang & Lesser (2010) , successfully preserves fixed points. We then prove that LookAhead locally converges and avoids strict saddles in all differentiable games, filling a theoretical gap in multi-agent learning. This is enabled through a unified approach based on fixed-point iterations and dynamical systems. These techniques apply equally well to algorithms like CO and SGA, though this is not our present focus.

While LookAhead is theoretically robust, the shaping component endowing LOLA with a capacity to exploit opponent dynamics is lost. We solve this dilemma with an algorithm named Stable Opponent Shaping (SOS), trading between stability and exploitation by interpolating between LookAhead and LOLA. Using an intuitive and theoretically grounded criterion for this interpolation parameter, SOS inherits both strong convergence guarantees from LA and opponent shaping from LOLA.

On the experimental side, we show that SOS plays tit-for-tat in the IPD on par with LOLA, while all other methods mostly defect. We display the practical consequences of our theoretical guarantees in the tandem game, where SOS always outperforms LOLA. Finally we implement a more involved GAN setup, testing for mode collapse and mode hopping when learning Gaussian mixture distributions. SOS successfully spreads mass across all Gaussians, at least matching dedicated algorithms like CO, while LA is significantly slower and simultaneous gradient descent fails entirely.

2. BACKGROUND:



2.1. DIFFERENTIABLE GAMES:

We frame the problem of multi-agent learning as a game. Adapted from Balduzzi et al. (2018) , the following definition insists only on differentiability for gradient-based methods to apply. This concept is strictly more general than stochastic games, whose parameters are usually restricted to action-state transition probabilities or functional approximations thereof. Definition 1. A differentiable game is a set of n players with parameters θ = (θ 1 , . . . , θ n ) ∈ R d and twice continuously differentiable losses L i : R d → R, where θ i ∈ R di for each i and i d i = d.

Crucially, note that each loss is a function of all parameters. From the viewpoint of player i, parameters can be written as θ = (θ i , θ -i ) where θ -i contains all other players' parameters. We do not make the common assumption that each L i is convex as a function of θ i alone, for any fixed opponent parameters θ -i , nor do we restrict θ to the probability simplex -though this restriction can be recovered via projection or sigmoid functions σ : R → [0, 1]. If n = 1, the 'game' is simply to minimise a given loss function. In this case one can reach local minima by (possibly stochastic) gradient descent (GD). For arbitrary n, the standard solution concept is that of Nash equilibria. Definition 2. A point θ ∈ R d is a (local) Nash equilibrium if for each i, there are neighbourhoods U i of θi such that L i (θ i , θ-i ) ≥ L i ( θ) for all θ i ∈ U i . In other words, each player's strategy is a local best response to current opponent strategies.

We write:

∇ i L k = ∇ θ i L k and ∇ ij L k = ∇ θ j ∇ θ i L k

for any i, j, k. Define the simultaneous gradient of the game as the concatenation of each player's gradient,

ξ = ∇ 1 L 1 , . . . , ∇ n L n ∈ R d .

The ith component of ξ is the direction of greatest increase in L i with respect to θ i . If each agent minimises their loss independently from others, they perform GD on their component ∇ i L i with learning rate α i . Hence, the parameter update for all agents is given by θ ← θ -α ξ, where α = (α 1 , . . . , α n ) and is element-wise multiplication. This is also called naive learning (NL), reducing to θ ← θ -αξ if agents have the same learning rate. This is assumed for notational simplicity, though irrelevant to our results. The following example shows that NL can fail to converge.

Example 1. Consider L 1/2 = ±xy, where players control the x and y parameters respectively. The origin is a (global and unique) Nash equilibrium. The simultaneous gradient is ξ = (y, -x) and cycles around the origin. Explicitly, a gradient step from (x, y) yields (x, y) ← (x, y) -α(y, -x) = (x -αy, y + αx) which has distance from the origin (1 + α 2 )(x 2 + y 2 ) > (x 2 + y 2 ) for any α > 0 and (x, y) = 0. It follows that agents diverge away from the origin for any α > 0. The cause of failure is that ξ is not the gradient of a single function, implying that each agent's loss is inherently dependent on others. This results in a contradiction between the non-stationarity of each agent, and the optimisation of each loss independently from others. Failure of convergence in this simple two-player zero-sum game shows that gradient descent does not generalise well to differentiable games. We consider an alternative solution concept to Nash equilibria before introducing LOLA.

2.2. STABLE FIXED POINTS:

Consider the game given by L 1 = L 2 = xy where players control the x and y parameters respectively. The optimal solution is (x, y) → ±(∞, -∞), since then L 1 = L 2 → -∞. However the origin is a global Nash equilibrium, while also a saddle point of xy. It is highly undesirable to converge to the origin in this game, since infinitely better losses can be reached in the anti-diagonal direction. In this light, Nash equilibria cannot be the right solution concept to aim for in multi-agent learning. To define stable fixed points, first introduce the 'Hessian' of the game as the block matrix

H = ∇ξ =     ∇ 11 L 1 • • • ∇ 1n L 1 . . . . . . . . . ∇ n1 L n • • • ∇ nn L n     ∈ R d×d .

This can equivalently be viewed as the Jacobian of the vector field ξ. Importantly, note that H is not symmetric in general unless n = 1, in which case we recover the usual Hessian

H = ∇ 2 L. Definition 3. A point θ is a fixed point if ξ( θ) = 0. It is stable if H( θ) 0, unstable if H( θ) ≺ 0

and a strict saddle if H( θ) has an eigenvalue with negative real part.

The name 'fixed point' is coherent with GD, since ξ( θ) = 0 implies a fixed update θ ← θ -αξ( θ) = θ. Though Nash equilibria were shown to be inadequate above, it is not obvious that stable fixed points (SFPs) are a better solution concept. In Appendix A we provide intuition for why SFPs are both closer to local minima in the context of multi-loss optimisation, and more tractable for convergence proofs. Moreover, this definition is an improved variant on that in Balduzzi et al. (2018) , assuming positive semi-definiteness only at θ instead of holding in a neighbourhood. This makes the class of SFPs as large as possible, while sufficient for all our theoretical results.

Assuming invertibility of H( θ) at SFPs is crucial to all convergence results in this paper. The same assumption is present in related work including Mescheder et al. (2017) , and cannot be avoided. Even for single losses, a fixed point with singular Hessian can be a local minimum, maximum, or saddle point. Invertibility is thus necessary to ensure that SFPs really are 'local minima'. This is omitted from now on. Finally note that unstable fixed points are a subset of strict saddles, making Theorem 6 both stronger and more general than results for SGA by Balduzzi et al. (2018) .

2.3. LEARNING WITH OPPONENT-LEARNING AWARENESS (LOLA):

Accounting for nonstationarity, Learning with Opponent-Learning Awareness (LOLA) modifies the learning objective by predicting and differentiating through opponent learning steps (Foerster et al., 2018) . For simplicity, if n = 2 then agent 1 optimises L 1 (θ 1 , θ 2 + ∆θ 2 ) with respect to θ 1 , where ∆θ 2 is the predicted learning step for agent 2. Foerster et al. (2018) assume that opponents are naive learners, namely ∆θ 2 = -α 2 ∇ 2 L 2 . After first-order Taylor expansion, the loss is approximately given by L 1 + ∇ 2 L 1 • ∆θ 2 . By minimising this quantity, agent 1 learns parameters that align the opponent learning step ∆θ 2 with the direction of greatest decrease in L 1 , exploiting opponent dynamics to further reduce one's losses. Differentiating with respect to θ 1 , the adjustment is

∇ 1 L 1 + ∇ 21 L 1 ∆θ 2 + ∇ 1 ∆θ 2 ∇ 2 L 1 .

By explicitly differentiating through ∆θ 2 in the rightmost term, LOLA agents actively shape opponent learning. This has proven effective in reaching cooperative equilibria in multi-agent learning, finding success in a number of games including tit-for-tat in the IPD. The middle term above was originally dropped by the authors because "LOLA focuses on this shaping of the learning direction of the opponent". We choose not to eliminate this term, as also inherent in LOLA-DiCE (Foerster et al., 2018) . Preserving both terms will in fact be key to developing stable opponent shaping. 

= (I -αH o )ξ -α χ .

While experimentally successful, LOLA fails to preserve fixed points θ of the game since

(I -αH o )ξ( θ) -α χ ( θ) = -α χ ( θ) = 0

in general. Even if θ is a Nash equilibrium, the update θ ← θ -αLOLA = θ can push them away despite parameters being optimal. This may worsen the losses for all agents, as in the game below.

Example 2 (Tandem). Imagine a tandem controlled by agents facing opposite directions, who feed x and y force into their pedals respectively. Negative numbers correspond to pedalling backwards.

Moving coherently requires x ≈ -y, embodied by a quadratic loss (x+y) 2 . However it is easier for agents to pedal forwards, translated by linear losses -2x and -2y. The game is thus given by L 1 (x, y) = (x + y) 2 -2x and L 2 (x, y) = (x + y) 2 -2y. These sub-goals are incompatible, so agents cannot simply accelerate forwards. The SFPs are given by {x + y = 1}.

Computing χ (x, 1 -x) = (4, 4) = 0, none of these are preserved by LOLA. Instead, we show in Appendix C that LOLA can only converge to sub-optimal scenarios with worse losses for both agents, for any α. Intuitively, the root of failure is that LOLA agents try to shape opponent learning and enforce compliance by accelerating forwards, assuming a dynamic response from their opponent. The other agent does the same, so they become 'arrogant' and suffer by pushing strongly in opposite directions.

3. METHOD:



3.1. LOOKAHEAD:

The shaping term χ prevents LOLA from preserving fixed points. Consider removing this component entirely, giving (I -αH o )ξ. This variant preserves fixed points, but what does it mean from the perspective of each agent? Note that LOLA optimises L 1 (θ 1 , θ 2 + ∆θ 2 ) with respect to θ 1 , while ∆θ 2 is a function of θ 1 . In other words, we assume that our opponent's learning step depends on our current optimisation with respect to θ 1 . This is inaccurate, since opponents cannot see our updated parameters until the next step. Instead, assume we optimise L 1 (θ 1 , θ2 + ∆θ 2 ( θ1 , θ2 )) where θ1 , θ2 are the current parameters. After Taylor expansion, the gradient with respect to θ 1 is given by

∇ 1 L 1 + ∇ 21 L 1 ∆θ 2

since ∆θ 2 ( θ1 , θ2 ) does not depend on θ 1 . In vectorial form, we recover the variant (I -αH o )ξ since the shaping term corresponds precisely to differentiating through ∆θ 2 . We name this LookAhead, which was discovered before LOLA by Zhang & Lesser (2010) though not explicitly named. Using the stop-gradient operator ⊥1 , this can be reformulated as optimising L 1 (θ 1 , θ 2 + ⊥∆θ 2 ) where ⊥ prevents gradient flowing from ∆θ 2 upon differentiation.

The main result of Zhang & Lesser (2010) is that LookAhead converges to Nash equilibria in the small class of two-player, two-action bimatrix games. We will prove local convergence to SFP and non-convergence to strict saddles in all differentiable games. On the other hand, by discarding the problematic shaping term, we also eliminated LOLA's capacity to exploit opponent dynamics and encourage cooperation. This will be witnessed in the IPD, where LookAhead agents mostly defect.

3.2. STABLE OPPONENT SHAPING (SOS):

We propose Stable Opponent Shaping (SOS), an algorithm preserving both advantages at once. Define the partial stop-gradient operator ⊥ p := p⊥ + (1 -p)I, where I is the identity and p stands for partial. A p-LOLA agent optimises the modified objective

L 1 (θ 1 , θ 2 + ⊥ 1-p ∆θ 2 , . . . , θ n + ⊥ 1-p ∆θ n

) , collapsing to LookAhead at p = 0 and LOLA at p = 1. The resulting gradient is given by ξ p := p-LOLA = (I -αH o )ξ -pα χ with ξ 0 = LA. We obtain an algorithm trading between shaping and stability as a function of p. Note however that preservation of fixed points only holds if p is infinitesimal, in which case p-LOLA is almost identical to LookAhead -losing the very purpose of interpolation. Instead we propose a two-part criterion for p at each learning step, through which all guarantees descend.

First choose p such that ξ p points in the same direction as LookAhead. This will not be enough to prove convergence itself, but prevents arrogant behaviour by ensuring convergence only to fixed points. Formally, the first criterion is given by ξ p , ξ 0 ≥ 0. If -α χ , ξ 0 ≥ 0 then ξ p , ξ 0 ≥ 0 automatically, so we choose p = 1 for maximal shaping. Otherwise choose p = min 1, -a ξ 0 2 -α χ , ξ 0 with any hyperparameter 0 < a < 1. This guarantees a positive inner product

ξ p , ξ 0 = p -α χ , ξ 0 + ξ 0 2 ≥ -a ξ 0 2 + ξ 0 2 = ξ 0 2 (1 -a) > 0 .

We complement this with a second criterion ensuring local convergence. The idea is to scale p by a function of ξ if ξ is small enough, which certainly holds in neighbourhoods of fixed points. Let 0 < b < 1 be a hyperparameter and take p = ξ 2 if ξ < b, otherwise p = 1. Choosing p 1 and p 2 according to these criteria, the two-part criterion is p = min{p 1 , p 2 }. SOS is obtained by combining p-LOLA with this criterion, as summarised in Algorithm 1. Crucially, all theoretical results in the next section are independent from the choice of hyperparameters a and b.

Algorithm 1: Stable Opponent Shaping 1 Initialise θ randomly and fix hyperparameters a, b ∈ (0, 1).

2 while not done do

3 Compute ξ 0 = (I -αH o )ξ and χ = diag(H o ∇L) at θ. 4 if -α χ , ξ 0 > 0 then p 1 = 1 else p 1 = min 1, -a ξ0 2 -α χ ,ξ0 5 if ξ < b then p 2 = ξ 2 else p 2 = 1 6 Let p = min{p 1 , p 2 }, compute ξ p = ξ 0 -pα χ and assign θ ← θ -αξ p . 7 end

4. THEORETICAL RESULTS:

Our central theoretical contribution is that LookAhead and SOS converge locally to SFP and avoid strict saddles in all differentiable games. Since the learning gradients involve second-order Hessian terms, our results assume thrice continuously differentiable losses (omitted hereafter). Losses which are C 2 but not C 3 are very degenerate, so this is a mild assumption. Statements made about SOS crucially hold for any hyperparameters a, b ∈ (0, 1). See Appendices D and E for detailed proofs.

4.1. LOCAL CONVERGENCE TO STABLE FIXED POINTS:

Convergence is proved using Ostrowski's Theorem. This reduces convergence of a gradient adjustment g to positive stability (eigenvalues with positive real part) of ∇g at stable fixed points. Theorem 2. Let H 0 be invertible with symmetric diagonal blocks. Then there exists > 0 such that (I -αH o )H is positive stable for all 0 < α < . This type of result would usually be proved either by analytical means showing positive definiteness and hence positive stability, or direct eigenvalue analysis. We show in Appendix D that (I -αH o )H is not necessarily positive definite, while there is no necessary relationship between eigenpairs of H and H o . This makes our theorem all the more interesting and non-trivial. We use a similarity transformation trick to circumvent the dual obstacle, allowing for analysis of positive definiteness with respect to a new inner product. We obtain positive stability by invariance under change of basis. Corollary 3. LookAhead converges locally to stable fixed points for α > 0 sufficiently small. Using the second criterion for p, we prove local convergence of SOS in all differentiable games despite the presence of a shaping term (unlike LOLA). Theorem 4. SOS converges locally to stable fixed points for α > 0 sufficiently small.

4.2. AVOIDING STRICT SADDLES:

Using the first criterion for p, we prove that SOS only converges to fixed points (unlike LOLA). Proposition 5. If SOS converges to θ and α > 0 is small then θ is a fixed point of the game. Now assume that θ is initialised randomly (or with arbitrarily small noise), as is standard in ML. Let F (θ) = θ -αξ p (θ) be the SOS iteration. Using both the second criterion and the Stable Manifold Theorem from dynamical systems, we can prove that every strict saddle θ has a neighbourhood U such that {θ ∈ U | F n (θ) → θ as n → ∞} has measure zero for α > 0 sufficiently small. Since θ is initialised randomly, we obtain the following result. Theorem 6. SOS locally avoids strict saddles almost surely, for α > 0 sufficiently small. This also holds for LookAhead, and could be strenghtened to global initialisations provided a strong boundedness assumption on H 2 . This is trickier for SOS since p(θ) is not globally continuous. Altogether, our results for LookAhead and the correct criterion for p-LOLA lead to some of the strongest theoretical guarantees in multi-agent learning. Furthermore, SOS retains all of LOLA's opponent shaping capacity while LookAhead does not, as shown experimentally in the next section.

5. EXPERIMENTS AND DISCUSSION:

We evaluate the performance of SOS in three differentiable games. We first showcase opponent shaping and superiority over LA/CO/SGA/NL in the Iterated Prisoner's Dilemma (IPD). This leaves SOS and LOLA, which have differed only in theory up to now. We bridge this gap by showing that SOS always outperforms LOLA in the tandem game, avoiding arrogant behaviour by decaying p while LOLA overshoots. Finally we test SOS on a more involved GAN learning task, with results similar to dedicated methods like Consensus Optimisation.

5.1. EXPERIMENTAL SETUP:

IPD: This game is an infinite sequence of the well-known Prisoner's Dilemma, where the payoff is discounted by a factor γ ∈ [0, 1) at each iteration. Agents are endowed with a memory of actions at the previous state. Hence there are 5 parameters for each agent i: the probability P i (C | state) of cooperating at start state s 0 = ∅ or state s t = (a 1 t-1 , a 2 t-1 ) for t > 0. One Nash equilibrium is to always defect (DD), with a normalised loss of 2. A better equilibrium with loss 1 is named tit-for-tat (TFT), where each player begins by cooperating and then mimicks the opponent's previous action.

We run 300 training episodes for SOS, LA, CO, SGA and NL. The parameters are initialised following a normal distribution around 1/2 probability of cooperation, with unit variance. We fix α = 1 and γ = 0.96, following Foerster et al. (2018) . We choose a = 0.5 and b = 0.1 for SOS. The first is a robust and arbitrary middle ground, while the latter is intentionally small to avoid poor SFP. Tandem: Though local convergence is guaranteed for SOS, it is possible that SOS diverges from poor initialisations. This turns out to be impossible in the tandem game since the Hessian is globally positive semi-definite. We show this explicitly by running 300 training episodes for SOS and LOLA. Parameters are initialised following a normal distribution around the origin. We found performance to be robust to hyperparameters a, b. Here we fix a = b = 0.5 and α = 0.1.

Gaussian mixtures::

We reproduce a setup from Balduzzi et al. (2018) . The game is to learn a Gaussian mixture distribution using GANs. Data is sampled from a highly multimodal distribution designed to probe the tendency to collapse onto a subset of modes during training -see ground truth in Appendix F. The generator and discriminator networks each have 6 ReLU layers of 384 neurons, with 2 and 1 output neurons respectively. Learning rates are chosen by grid search at iteration 8k, with a = 0.5 and b = 0.1 for SOS, following the same reasoning as the IPD. SOS and LOLA mostly succeed in playing tit-for-tat, displayed by the accumulation of points in the correct corners of (A) plots. For instance, CC and CD points are mostly in the top right and left corners so agent 2 responds to cooperation with cooperation. Agents also cooperate at the start state, represented by ∅ points all hidden in the top right corner. Tit-for-tat strategy is further indicated by the losses close to 1 in part (B). On the other hand, most points for LA/CO/SGA/NL are accumulated at the bottom left, so agents mostly defect. This results in poor losses, demonstrating the limited effectiveness of recent proposals like SGA and CO. Finally note that trained parameters and losses for SOS are almost identical to those for LOLA, displaying equal capacity in opponent shaping while also inheriting convergence guarantees and outperforming LOLA in the next experiment.

5.2. RESULTS AND:

Tandem: Results are given in Figure 3 . SOS always succeeds in decreasing p to reach the correct equilibria, with losses averaging at 0. LOLA fails to preserve fixed points, overshooting with losses averaging at 4/9. The criterion for SOS is shown in action in part (B), decaying p to avoid overshooting. This illustrates that purely theoretical guarantees descend into practical outperfor- mance. Note that SOS even gets away from the LOLA fixed points if initialised there (not shown), converging to improved losses using the alignment criterion with LookAhead.

Gaussian mixtures: The generator distribution and KL divergence are given at {2k, 4k, 6k, 8k} iterations for NL, CO and SOS in Figure 4 . Results for SGA, LOLA and LA are in Appendix F. SOS achieves convincing results by spreading mass across all Gaussians, as do CO/SGA/LOLA. LookAhead is significantly slower, while NL fails through mode collapse and hopping. Only visual inspection was used for comparison by Balduzzi et al. (2018) , while KL divergence gives stronger numerical evidence here. SOS and CO are slightly superior to others with reference to this metric. However CO is aimed specifically toward two-player zero-sum GAN optimisation, while SOS is widely applicable with strong theoretical guarantees in all differentiable games.

6. CONCLUSION:

Theoretical results in machine learning have significantly helped understand the causes of success and failure in applications, from optimisation to architecture. While gradient descent on single losses has been studied extensively, algorithms dealing with interacting goals are proliferating, with little grasp of the underlying dynamics. The analysis behind CO and SGA has been helpful in this respect, though lacking either in generality or convergence guarantees. The first contribution of this paper is to provide a unified framework and fill this theoretical gap with robust convergence results for LookAhead in all differentiable games. Capturing stable fixed points as the correct solution concept was essential for these techniques to apply.

Furthermore, we showed that opponent shaping is both a powerful approach leading to experimental success and cooperative behaviour -while at the same time preventing LOLA from preserving fixed points in general. This conundrum is solved through a robust interpolation between LookAhead and LOLA, giving birth to SOS through a robust criterion. This was partially enabled by choosing to preserve the 'middle' term in LOLA, and using it to inherit stability from LookAhead. This results in convergence guarantees stronger than all previous algorithms, but also in practical superiority over LOLA in the tandem game. Moreover, SOS fully preserves opponent shaping and outperforms SGA, CO, LA and NL in the IPD by encouraging tit-for-tat policy instead of defecting. Finally, SOS convincingly learns Gaussian mixtures on par with the dedicated CO algorithm.

APPENDIX A STABLE FIXED POINTS:

In the main text we showed that Nash equilibria are inadequate in multi-agent learning, exemplified by the simple game given by L 1 = L 2 = xy, where the origin is a global Nash equilibrium but a saddle point of the losses. It is not however obvious that SFP are a better solution concept. We begin by pointing out that for single losses, invertibility and symmetry of the Hessian imply positive definiteness at SFP. These are exactly local minima of L detected by the second partial derivative test, namely those points provably attainable by gradient descent.

To emphasise this, note that gradient descent does not converge locally to all local minima. This can be seen by considering the example L(x, y) = y 2 and the local (global) minimum (0, 0). There is no neighbourhood for which gradient descent converges to (0, 0), since initialising at (x 0 , y 0 ) will always converge to (x 0 , 0) for appropriate learning rates, with x 0 = 0 almost surely. This occurs precisely because the Hessian is singular at (0, 0). Though a degenerate example, this suggests an important difference to make between the ideal solution concept (local minima) and that for which local convergence claims are possible to attain (local minima with invertible H 0).

Accordingly, the definition of SFP is the immediate generalisation of 'fixed points with positive semi-definite Hessian', or in other words, 'second-order-tractable local minima'. It is important to impose only positive semi-definiteness to keep the class as large as possible, despite strict positive definiteness holding for single losses due to symmetry. Imposing strict positivity would for instance exclude the origin in the cyclic game L 1 = xy = -L 2 , a point certainly worthy of convergence.

Note also that imposing a weaker condition than H 0 would be incorrect. Invertibility aside, local convergence of gradient descent on single functions cannot be guaranteed if H 0, since such points are strict saddles. These are almost always avoided by gradient descent, as proven by Lee et al. (2016) and Panageas & Piliouras (2017) . It is thus necessary to impose H 0 as a minimal requirement in optimisation methods attempting to generalise gradient descent. Remark A.1. A matrix H is positive semi-definite iff the same holds for its symmetric part S = (H + H )/2, so SFP could equivalently be defined as S( θ) 0. This is the original formulation given by part of the authors (Balduzzi et al., 2018) , who also imposed the extra requirement S(θ) 0 in a neighbourhood of θ. After discussion we decided to drop this assumption, pointing out that it is 1) more restrictive, 2) superficial to all theoretical results and 3) weakens the analogy with tractable local minima. The only thing gained by imposing semi-positivity in a neighbourhood is that SFP become a subset of Nash equilibria.

Regarding unstable fixed points and strict saddles, note that H( θ) 0 implies H(θ) 0 in a neighbourhood, hence being equivalent to the definition in Balduzzi et al. (2018) . It follows also that unstable points are a subset of strict saddles: if H( θ) ≺ 0 then all eigenvalues are negative since any eigenpair (v, λ) satisfies

0 > Re(v Hv) = Re(λv v) = Re(λ) .

We introduced strict saddles in this paper as a generalisation of unstable FP, which are more difficult to handle but nonetheless tractable using dynamical systems. The name is chosen by analogy to the definition in Lee et al. (2016) for single losses.

B LOLA VECTORIAL FORM:

Proposition B.1. The LOLA gradient adjustment is LOLA = (I -αH o )ξ -α diag(H o ∇L) .

in the usual assumption of equal learning rates.

Proof. Recall the modified objective

L 1 (θ 1 , θ 2 -α∇ 2 L 2 , . . . , θ n -α∇ n L n )

for agent 1, and so on for each agent. First-order Taylor expansion yields

L 1 -α j =1 (∇ j L 1 ) ∇ j L j

and similarly for each agent. Differentiating with respect to θ i , the adjustment for player i is

LOLA i = ∇ i   L i -α j =i (∇ j L i ) ∇ j L j   = ∇ i L i -α j =i (∇ ji L i ) ∇ j L j + (∇ ji L j ) ∇ j L i = ∇ i L i -α j =i ∇ ij L i ∇ j L j -α j =i (∇ ji L j ) ∇ j L i = ξ i -α j (H o ) ij ξ j -α j (H o ) ij (∇L) ji = ξ i -α(H o ξ) i -α(H o ∇L) ii = ξ -αH o ξ -α diag(H o ∇L) i and thus LOLA = (I -αH o )ξ -α diag(H o ∇L) as required.

C TANDEM GAME:

We provide a more detailed exposition of the tandem game in this section, including computation of fixed points for NL/LOLA and corresponding losses. Recall that the game is given by

L 1 (x, y) = (x + y) 2 -2x and L 2 (x, y) = (x + y) 2 -2y .

Intuitively, agents wants to have x ≈ -y since (x + y) 2 is the leading loss, but would also prefer to have positive x and y. These are incompatible, so the agents must not be 'arrogant' and instead make concessions. The fixed points are given by ξ = 2(x + y -1) 1 1 = 0 , namely any pair (x, 1 -x). The corresponding losses are L 1 = 1 -2x = -L 2 , summing to 0 for any x. We have

H = 2 1 1 1 1 0

everywhere, so all fixed points are SFP. LOLA fails to preserve these, since

χ = diag(H o ∇L) = 4 diag 0 1 1 0 x + y -1 x + y x + y x + y -1 = 4(x + y) 1 1

which is non-zero for any SFP (x, 1 -x). Instead, LOLA can only converge to points such that

LOLA = ξ -αH o ξ -α χ = 0 .

We solve this explicitly as follows:

LOLA = 2(x + y -1) 1 1 -4α(x + y -1) 0 1 1 0 1 1 -4α(x + y) 1 1 = 2 [(1 -4α)(x + y) -(1 -2α)] 1 1 .

The fixed points for LOLA are thus pairs (x, y) such that

x + y = 1 -2α 1 -4α ,

noting that (1 -2α)/(1 -4α) > 1 for all α > 0. This leads to worse losses

L 1 = 1 -2α 1 -4α 2 -2x > 1 -2x = L 1 (x, 1 -x)

for agent 1 and similarly for agent 2. In particular, losses always sum to something greater than 0. This becomes negligible as the learning rate becomes smaller, but is always positive nonetheless Taking α arbitrarily small is not a viable solution since convergence will in turn be arbitrarily slow.

LOLA is thus not a strong algorithm candidate for all differentiable games.

D CONVERGENCE PROOFS:

We use Ostrowski's theorem as a unified framework for proving local convergence of gradient-based methods. This is a standard result on fixed-point iterations, adapted from (Ortega & Rheinboldt, 2000, 10.1.3) . We also invoke and prove a topological result of our own, Lemma D.9, at the end of this section. This is useful in deducing local convergence, though not central to intuition. Theorem D.1 (Ostrowski). Let F : Ω → R d be continuously differentiable on an open subset Ω ⊆ R d , and assume x ∈ Ω is a fixed point. If all eigenvalues of ∇F (x) are strictly in the unit circle of C, then there is an open neighbourhood U of x such that for all x 0 ∈ U , the sequence F (k) (x 0 ) converges to x. Moreover, the rate of convergence is at least linear in k. Definition D.2. A matrix M is called positive stable if all its eigenvalues have positive real part.

Recall the simultaneous gradient ξ and the Hessian H defined for differentiable games. Let X be any matrix with continuously differentiable entries. Corollary D.3. Assume x is a fixed point of a differentiable game such that XH(x) is positive stable. Then the iterative procedure

F (x) = x -αXξ(x)

converges locally to x for α > 0 sufficiently small.

Proof. By definition of fixed points, ξ(x) = 0 and so

∇[Xξ](x) = ∇X(x)ξ(x) + X(x)∇ξ(x) = XH(x)

is positive stable by assumption, namely has eigenvalues a k + ib k with a k > 0. It follows that

∇F (x) = I -α∇[Xξ](x)

has eigenvalues 1 -αa k -iαb k , which are in the unit circle for small α. More precisely,

|1 -αa k -iαb k | 2 < 1 ⇐⇒ 1 -2αa k + α 2 a 2 k + α 2 b 2 k < 1 ⇐⇒ 0 < α < 2a k a 2 k + b 2 k

which is always possible for a k > 0. Hence ∇F (x) has eigenvalues in the unit circle for 0 < α < min k 2a k /(a 2 k + b 2 k ), and we are done by Ostrowski's Theorem since x is a fixed point of F .

We apply this corollary to LookAhead, which is given by

F (θ) = θ -αXξ(θ)

where X = (I -αH o ). It is thus sufficient to prove the following result. Theorem D.4. Let H 0 invertible with symmetric diagonal blocks. Then there exists > 0 such that (I -αH o )H is positive stable for all 0 < α < .

Remark D.5. Note that (I -αH o )H may fail to be positive definite, though true in the case of 2 × 2 matrices. This no longer holds in higher dimensions, exemplified by the Hessian

H =       9 -4 -3 -3 -2 1 2 1 -3 0 1 0 -3 1 2 1      

.

By direct computation (symbolic in α), one can show that G = (I -αH o )H always has positive eigenvalues for small α > 0, whereas its symmetric part S always has a negative eigenvalue with magnitude in the order of α. This implies that S and in turn G is not positive definite. As such, an analytical proof of the theorem involving bounds on the corresponding bilinear form will fail.

This makes the result all the more interesting, but more involved. Central to the proof is a similarity transformation proving positive definiteness with respect to a different inner product, a novel technique we have not found in the multi-agent learning literature.

Proof. We cannot study the eigenvalues of G directly, since there is no necessary relationship between eigenpairs of H and H o . In the aim of using analytical tools, the trick is to find a positive definite matrix which is similar to G, thus sharing the same positive eigenvalues. First define

G 1 = (I + αH d )H and G 2 = -αH 2 ,

where H d is the sub-matrix of diagonal blocks,and rewrite

G = (I -αH o )H = (I -α(H -H d ))H = (I + αH d )H -αH 2 = G 1 + G 2 .

Note that H d is block diagonal with symmetric blocks ∇ ii L i 0, so (I + αH d ) is symmetric and positive definite for all α ≥ 0. In particular its principal square root

M = (I + αH d ) 1/2

is unique and invertible. Now note that

M -1 G 1 M = M -1 M 2 HM = M HM ,

which is positive semi-definite since u M HM u = (M u) H(M u) ≥ 0 for all non-zero u. In particular M provides a similarity transformation which eliminates H d from G 1 while simultaneously delivering positive semi-definiteness. We can now prove that

M -1 GM = M -1 G 1 M + M -1 G 2 M

is positive definite, establishing positive stability of G by similarity. Let m = d -1 where d is the vector space dimension, namely H ∈ R d×d . Recall that the m-sphere S m ⊂ R d is the space of unit vectors in R d . Take any u ∈ S m and consider the quantity

u M -1 GM u .

First note that a Taylor expansion of M in α yields

M = (I + αH d ) 1/2 = I + O(α) and M -1 = (I + αH d ) -1/2 = I + O(α) .

This implies in turn that u M -1 GM u = u Gu + O(α) .

There are two cases to distinguish. If u Hu > 0 then

u M -1 GM u = u Gu + O(α) = u G 1 u + O(α) = u Hu + O(α) > 0

for α sufficiently small. Otherwise, u Hu = 0 and consider decomposing H into symmetric and antisymmetric parts S = (H + H )/2 and A = (H -H )/2, so that H = S + A. By antisymmetry of A we have u Au = 0 and hence u Hu = 0 = u Su. Now H 0 implies S 0, so by Cholesky decomposition of S there exists a matrix T such that S = T T . In particular 0 = u Su = T u 2 implies T u = 0, and in turn Su = 0. Since H is invertible and u = 0, we have 0 = Hu = Au and so Au 2 > 0. It follows in particular that

-αu H 2 u = -αu (S -A )(S + A)u = αu A Au = α Au 2 > 0 .

Using positive semi-definiteness of

M -1 G 1 M , u M -1 GM u = u M -1 G 1 M u + u M -1 G 2 M u ≥ -αu M -1 H 2 M u = -αu H 2 u + O(α 2 ) = α Au 2 + O(α 2 ) > 0

for α > 0 small enough. We conclude that for any u ∈ S m there is u > 0 such that u M -1 GM u > 0 for all 0 < α < u , where g(α, u) = u M -1 GM u is a function g : R + × S m → R with S m compact. By Lemma D.9, this can be extended uniformly with some > 0 such that u M -1 GM u > 0 for all u ∈ S m and 0 < α < . It follows that M -1 GM is positive definite for all 0 < α < and thus G is positive stable for α in the same range, by similarity.

Corollary D.6. LookAhead converges locally to stable fixed points for α > 0 sufficiently small.

Proof. For any SFP θ we have ξ( θ) = 0 and H( θ) 0 invertible by definition, with diagonal blocks ∇ ii L i symmetric by twice continuous differentiability. We are done by the result above and Corollary D.3.

We now prove that local convergence results descend to SOS. The following lemma establishes the crucial claim that our criterion for p is C 1 in neighbourhoods of fixed points. This is necessary to invoke analytical arguments including Ostrowski's Theorem, and would be untrue globally. Lemma D.7. If θ is a fixed point and α is sufficiently small then p = ξ 2 in a neighbourhood of θ.

Proof. First note that ξ( θ) = 0, so there is a (bounded) neighbourhood V of θ such that ξ(θ) < b for all θ ∈ V , for any choice of hyperparameter b ∈ (0, 1). In particular p 2 (θ) = ξ(θ) 2 by definition of the second criterion. We want to show that p(θ) = p 2 (θ) near θ, or equivalently p

1 (θ) ≥ p 2 (θ). Since p 2 (θ) = ξ(θ) 2 < b 2 < 1 in V , it remains only to show that -a ξ 0 2 -α χ , ξ 0 ≥ ξ(θ) 2

in some neighbourhood U ⊆ V of θ, for any choice of hyperparameter a ∈ (0, 1). Now by boundedness of V and continuity of χ , there exists c > 0 such that -α χ (θ) = α 2 χ (θ) < c for all θ ∈ V and bounded α. It follows by Cauchy-Schwartz that

-a ξ 0 2 -α χ , ξ 0 ≥ a ξ 0 -α χ > a ξ 0 /c in V . Now note that ξ 0 = (I -αH o )ξ ≥ d ξ

in V , for some d > 0 and α sufficiently small, by boundedness of V and continuity of H o . Finally there is a sub-neighbourhood U ⊂ V such that ξ(θ) < ad/c for all θ ∈ U , so that ad ξ /c > ξ(θ) 2 and hence

-a ξ 0 2 -α χ , ξ 0 > ξ 2 = p 2 in U . Hence p(θ) = min{p 1 (θ), p 2 (θ)} = p 2 (θ) = ξ(θ)

2 for all θ ∈ U , as required.

Proof. Assume for contradiction that L > L , then there exists δ > 0 such that L > L + δ. By definition of limits, there exist M, N ∈ N such that

|c k -L| < δ/2 and |a k -L | < δ/2

for all k ≥ M , k ≥ N . Expanding the absolute value, this implies

L -δ/2 < c k < L + δ/2 and L -δ/2 < a k < L + δ/2 for all k ≥ max{M, N }. Now c k ≤ a k for all k, hence L -δ/2 < c k ≤ a k < L + δ/2

which implies the contradiction L < L + δ .

Proposition E.2. If SOS converges to θ and α > 0 is small then θ is a fixed point of the game.

Proof. The iterative procedure is given by

θ k+1 = F (θ k ) = θ k -αξ p (θ k ) .

If θ k → θ as k → ∞ then taking limits on both sides of the iteration yields

θ = θ -α lim k→∞ ξ p (θ k ) and so lim k ξ p (θ k ) = 0, omitting k → ∞ for convenience. It follows by continuity that ξ 0 ( θ) + lim k p(θ k ) -α χ ( θ) = 0 ,

noting that p(θ) is not a globally continuous function. Assume for contradiction that ξ 0 ( θ) = 0.

There are two cases to distinguish for clarity.

(i) First assume -α χ , ξ 0 ( θ) ≥ 0. Note that lim k p(θ k ) ≥ 0 since p(θ) ≥ 0 for all θ, and so

lim k ξ p (θ k ), ξ 0 ( θ) = lim k p(θ k ) -α χ , ξ 0 ( θ) + ξ 0 ( θ) 2 > 0 .

This is a contradiction since lim k ξ p (θ k ) = 0.

(ii) Otherwise, -α χ , ξ 0 ( θ) < 0 and hence -α χ , ξ 0 (θ) < 0 in a neighbourhood. In particular there exists N ∈ N such that -α χ , ξ 0 (θ k ) < 0 for all k ≥ N . In particular

p 1 (θ k ) = min 1, -a ξ 0 (θ k ) 2 -α χ , ξ 0 (θ k ) for all k ≥ N . Now notice that lim k p(θ k ) = lim k min 1, -a ξ 0 (θ k ) 2 -α χ , ξ 0 (θ k ) , p 2 (θ k ) , which implies lim k p(θ k ) ≤ lim k -a ξ 0 (θ k ) 2 -α χ , ξ 0 (θ k ) = -a ξ 0 ( θ) 2

-α χ , ξ 0 ( θ) by continuity and Lemma E.1. Finally we conclude

lim k ξ p , ξ 0 (θ k ) = lim k p(θ k ) -α χ , ξ 0 ( θ) + ξ 0 ( θ) 2 ≥ -a ξ 0 ( θ) 2 + ξ 0 ( θ) 2 > 0

for any a ∈ (0, 1), a contradiction.

In both cases a contradiction is obtained, hence ξ 0 ( θ) = 0 = (I -αH o )ξ( θ). Now note that (I -αH o )( θ) is singular iff H o ( θ) has an eigenvalue 1/α, which is impossible for α sufficiently small. Hence (I -αH o )ξ( θ) = 0 implies ξ( θ) = 0, as required. Now assume that θ is initialised randomly (or with arbitrarily small noise around a point), as is standard in ML. We prove that SOS locally avoids strict saddles using the Stable Manifold Theorem, inspired from Lee et al. (2017) . Theorem E.3 (Stable Manifold Theorem). Let x be a fixed point for the C 1 local diffeomorphism F : U → R d , where U is a neighbourhood of x in R d . Let E s ⊕ E u be the generalised eigenspaces of ∇F (x) corresponding to eigenvalues with |λ| ≤ 1 and |λ| > 1 respectively. Then there exists a local stable center manifold W with tangent space E s at x and a neighbourhood B of x such that

F (W ) ∩ B ⊂ W and ∩ ∞ n=0 F -n (B) ⊂ W .

In particular, if ∇F (x) has at least one eigenvalue |λ| > 1 then E u has dimension at least 1. Since W has tangent space E s at x, with codimension at least one, it follows that W has measure zero in R d . This is central in proving that the set of initial points in a neighbourhood which converge through SOS to a given strict saddle θ has measure zero. Theorem E.4. SOS locally avoids strict saddles almost surely, for α > 0 sufficiently small.

Proof. Let θ a strict saddle and recall that SOS is given by F (θ) = θ -α(I -αH o )ξ(θ) + α 2 p(θ) χ (θ) .

Recall by Lemma D.7 that p(θ) = ξ(θ) 2 for all θ in a neighbourhood U of θ. Restricting F to U , all terms involved are continuously differentiable and

∇F ( θ) = I -α(I -αH o )H( θ)

by assumption that ξ( θ) = 0. Since all terms except I are of order at least α, ∇F ( θ) is invertible for all α sufficiently small. By the inverse function theorem, there exists a neighbourhood V of θ such that F is has a continuously differentiable inverse on V . Hence F restricted to U ∩ V is a C 1 diffeomorphism with fixed point θ.

By definition of strict saddles, H( θ) has a negative eigenvalue. It follows by continuity that (I -αH o )H( θ) also has a negative eigenvalue a + ib with a < 0 for α sufficiently small. Finally, ∇F ( θ) = I -α(I -αH o )H( θ) has an eigenvalue λ = 1 -αa -iαb with |λ| = 1 -2αa + α 2 (a 2 + b 2 ) ≥ 1 -2αa > 1 .

It follows that E s has codimension at least one, implying in turn that the local stable set W has measure zero. We can now prove that

Z = {θ ∈ U ∩ V | lim n→∞ F n (θ) = θ}

has measure zero, or in other words, that local convergence to θ occurs with zero probability. Let B the neighbourhood guaranteed by the Stable Manifold Theorem, and take any θ ∈ Z. By definition of convergence there exists N ∈ N such that F N +n (θ) ∈ B for all n ≥ 0, so that F N (θ) ∈ ∩ ∞ n=0 F -n (B) ⊂ W by the Stable Manifold Theorem. This implies that θ ∈ F -N (W ), and finally θ ∈ ∪ n∈N F -n (W ). Since θ was arbitrary, we obtain the inclusion Z ⊆ ∪ n∈N F -n (W ) . Now F -1 is C 1 , hence locally Lipschitz and thus preserves sets of measure zero, so that F -n (W ) has measure zero for each n. Countable unions of measure zero sets are still measure zero, so we conclude that Z also has measure zero. In other words, SOS converges to θ with zero probability upon random initialisation of θ in U . 

F FURTHER GAUSSIAN MIXTURE EXPERIMENTS:

In the Gaussian mixture experiment, data is sampled from a highly multimodal distribution designed to probe the tendency to collapse onto a subset of modes during training, given in Figure 5 .

The generator distribution and KL divergence are given at {2k, 4k, 6k, 8k} iterations for LA, LOLA and SGA in Figure 6 . LOLA and SGA successfully spread mass across all Gaussians. LookAhead displays mode collapse and hopping in early stages, but begins to incorporate further mixtures near 8k iterations. We ran further iterations and discovered that LookAhead eventually spreads mass across all mixtures, though very slowly. Comparing with results for NL/CO/SOS in the main text, we see that CO/SOS/LOLA/SGA are equally successful in qualitative terms.

Note that SOS/CO are slightly superior with respect to KL divergence after 6-8k iterations, though LOLA is initially faster. This may be due only to random sampling. We also noticed experimentally that LOLA often moves away from the correct distribution after 8-10k iterations (not shown), while SOS stays stable in the long run. This may occur thanks to the two-part criterion encouraging convergence, while LOLA continually attempts to exploit opponent learning.

Finally we plot ξ at all iterations up to 12k for SOS, LA and NL in Figure 7 (other algorithms are omitted for visibility). This gives further evidence of SOS converging quite rapidly to the correct distribution, while NL perpetually suffers from mode hopping and LA lags behind significantly. 

:

This operator is implemented in TensorFlow as stop gradient and in PyTorch as detach.

Figures and Tables:

Figure 1: Illustration of the tandem game.

DISCUSSION IPD: Results are given in Figure 2. Parameters in part (A) are the end-run probabilities of cooperating for each memory state, encoded in different colours. Only 50 runs are shown for visibility. Losses at each step are displayed in part (B), averaged across 300 episodes with shaded deviations.

Figure 3: Results in the tandem game. (A) Average loss and (B) average p at each learning step, across 300 runs, with shaded deviations. SOS decays p to avoid arrogance and outperforms LOLA. Iteration = 2k 4k 6k 8k

Figure 5: Ground truth in the Gaussian mixture experiment.

Figure6: Generator distribution at sampled iterations for LA/LOLA/SGA. LA suffers in the early stages from mode collapse and hopping, but incorporates more mixtures later on. LOLA and SGA learn the correct mixture of Gaussians. Below each plot: KL divergence D KL (P || Q) from generator P to ground truth Q, estimated from 25600 samples. Best result at each iteration in bold.

First we formulate n-player LOLA in vectorial form. Let H d and H o be the matrices of diagonal and anti-diagonal blocks of H, so that H = H d + H o . Also define L = (L 1 , . . . , L n ) and the operator diag : R d×n → R d constructing a vector from the block matrix diagonal, namely diag(M ) i = M ii .

Table content: <table><tr><td>Proposition 1 (Appendix B). Writing χ = diag(H o ∇L), the LOLA gradient adjustment is</td></tr><tr><td>LOLA</td></tr></table>

Results in the IPD. (A) Probability that agents cooperate, given memory state, at the end of 50 training runs. SOS and LOLA mostly play tit-for-tat, while others mostly defect. (B) Average loss at each step, across 300 runs, with shaded deviations. SOS and LOLA outperform all others.

Table content: <table><tr><td>1.0</td><td/><td/><td>(A) SOS</td><td/><td/><td/><td/><td>(A) LOLA</td><td/><td/><td/><td/><td/><td>(A) LA</td><td/><td/><td/><td/><td>(A) CO</td><td/></tr><tr><td>0.8</td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/></tr><tr><td>0.4 0.6 | state) P 2 (C</td><td/><td/><td>CC CD DC DD</td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/></tr><tr><td>0.2</td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/></tr><tr><td>0.0 1.0</td><td>0.0</td><td>0.2</td><td>0.4 P 1 (C | state) 0.6 (A) SGA</td><td>0.8</td><td>1.0</td><td>0.0</td><td>0.2</td><td>0.4 P 1 (C | state) 0.6 (A) NL</td><td>0.8</td><td>1.0</td><td colspan="2">0.0 2.2</td><td>0.2</td><td>0.4 P 1 (C | state) 0.6</td><td>0.8</td><td>1.0</td><td>0.0 (B)</td><td>0.2</td><td>0.4 P 1 (C | state) 0.6</td><td>0.8</td><td>1.0</td></tr><tr><td>0.8</td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td>2.0</td><td/><td/><td/><td/><td/><td/><td/><td/></tr><tr><td>0.2 0.4 0.6 | state) P 2 (C</td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td>Average Loss</td><td>1.2 1.4 1.6 1.8</td><td/><td/><td/><td/><td>SOS LOLA LA CO SGA NL</td><td/><td/><td/></tr><tr><td>0.0</td><td>0.0</td><td>0.2</td><td>0.4 P 1 (C | state) 0.6</td><td>0.8</td><td>1.0</td><td>0.0</td><td>0.2</td><td>0.4 P 1 (C | state) 0.6</td><td>0.8</td><td>1.0</td><td/><td>1.0</td><td>0</td><td>20</td><td/><td colspan="3">40 Learning Step 60</td><td>80</td><td/><td>100</td></tr><tr><td colspan="3">Figure 2:</td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/></tr></table>

Peer review (iclr2019):
 This paper focuses on the problem of convergence in multi-objective optimisation with differentiable losses. This topic is timely and relevant, given the increasing amount of recent work on multi-objective architectures, e.g. GANs, adversarial learning, multi-agent reinforcement learning. The authors focus on stable fixed points (SFP), rather than Nash equilibria, as the solution concept in the entirety of their analysis. Casting the recently proposed LOLA gradient adjustment into a general matrix form, they diagnose an example where the shaping term in LOLA prevents convergence to SFP. They also find that discarding the shaping term leads to an earlier method (which they name **LA**) with convergence guarantees in two-player two-action games. However, this also loses the opponent shaping ability of LOLA. To address these limitation, the authors propose SOS, which interpolates between LA and LOLA, and dynamically chooses the interpolation coefficient so that their adjusted gradient preserves LOLA*s shaping ability only to the extent allowed by the constraint of moving in LA*s direction. The main goal of the paper is to show that SOS converges locally to SFP, and to fixed points only, while avoiding strict saddles. Experiments on synthetic games show that SOS preserves the benefit of LOLA while avoiding its theoretically-predicted issues, and a more complex Gaussian mixture GAN experiment shows SOS is empirically competitive with other gradient adjustment methods. The main conceptual novelty consists of the dynamic interpolation term to combine advantages of LOLA and LA while avoiding pitfalls of both. The major strength of the paper lies in the clear justification for this interpolation approach. The paper contains strong theoretical results for general differentiable games, and deserves the notice of the ICLR community if valid. However, I have major concerns with the proof of Theorem 2 (i.e. Theorem D.4 in the appendix), which affects the validity of Corollary 3 and Theorem 4. In the proof of Theorem D.4: 1. How does the expression have conformable dimensions, when while ? Was any assumption made about the matrix ? 2. In the middle of page 14, a unit vector is defined, but it is not clear what vector space is meant by . 3. In the second-to-last line of page 14, a quantity is used but not defined clearly in any preceding part of the proof. Remark D.5 refers to as the symmetric part of , and asserts that S is not positive definite. If the quantity used in the proof is the same non-PD quantity, then does not have a Cholesky factorisation. So how is Cholesky decomposition conducted at end of page 14? 4. In the first line of page 15, a quantity is used but not defined anywhere else in the entire paper. 5. From the subsequent line, it appears to be the anti-symmetric part of H. Is it correct assumption? If so, is not . If you replace it with correct form, whole quantity does not compute to be positive or becomes meaningless. As Theorem 2 is the crux for all the theoretical advancement presented in the paper, clarifications on above correctness questions is very important for clear acceptance of this work. While Definition 1 precisely defines differentiable games to have *twice* differentiable losses, why do the authors assume *thrice* differentiable losses at the start of Section 4? In Section 2.2, the authors make a broad statement that **Nash equilibria cannot be the right solution concept for multi-agent learning.** They provide one example where Nash is undesirable (L^1 = L^2 = xy). However, since this example can be viewed as a fully cooperative game with joint loss L = 2xy, it does not support the broader statement that Nash is undesirable in all games. Because this statement directly motivates the authors to focus on stable fixed points, rather than Nash, as the solution concept in their subsequent analysis, it is very important to provide better justification for the claim. Minor comments: 1. Under Proposition 1, the authors suddenly speak of **...the policy being optimal**. Since the author*s work pertains to general multi-objective settings, not solely multi-agent reinforcement learning, the word **policy** sounds strange in context. 2. The statement of Proposition B.1, and the concluding line of the derivation, left out a coefficient that is present in Proposition 1 in the main text. 3. While the authors claim and prove independence of theoretical results from choice of a and b, are there any practical implications in terms of performance or convergence?