Title: A CLOSER LOOK AT THE OPTIMIZATION LANDSCAPES OF GENERATIVE ADVERSARIAL NETWORKS

Abstract: Generative adversarial networks have been very successful in generative modeling, however they remain relatively challenging to train compared to standard deep neural networks. In this paper, we propose new visualization techniques for the optimization landscapes of GANs that enable us to study the game vector field resulting from the concatenation of the gradient of both players. Using these visualization techniques we try to bridge the gap between theory and practice by showing empirically that the training of GANs exhibits significant rotations around Local Stable Stationary Points (LSSP), similar to the one predicted by theory on toy examples. Moreover, we provide empirical evidence that GAN training converge to a stable stationary point which is a saddle point for the generator loss, not a minimum, while still achieving excellent performance.

1. INTRODUCTION:

Deep neural networks have exhibited remarkable success in many applications (Krizhevsky et al., 2012) . This success has motivated many studies of their non-convex loss landscape (Choromanska et al., 2015; Kawaguchi, 2016; Li et al., 2018b) , which, in turn, has led to many improvements, such as better initialization and optimization methods (Glorot and Bengio, 2010; Kingma and Ba, 2015) .

While most of the work on studying non-convex loss landscapes has focused on single objective minimization, some recent class of models require the joint minimization of several objectives, making their optimization landscape intrinsically different. Among these models is the generative adversarial network (GAN) (Goodfellow et al., 2014) which is based on a two-player game formulation and has achieved state-of-the-art performance on some generative modeling tasks such as image generation (Brock et al., 2019) .

On the theoretical side, many papers studying multi-player games have argued that one main optimization issue that arises in this case is the rotation due to the adversarial component of the game (Mescheder et al., 2018; Balduzzi et al., 2018; Gidel et al., 2019b) . This has been extensively studied on toy examples, in particular on the so-called bilinear example (Goodfellow, 2016) (a.k.a Dirac GAN (Mescheder et al., 2018) ). However, those toy examples are very far from the standard realistic setting of image generation involving deep networks and challenging datasets. To our knowledge it remains an open question if this rotation phenomenon actually occurs when training GANs in more practical settings.

In this paper, we aim at closing this gap between theory and practice. Following Mescheder et al. (2017) and Balduzzi et al. (2018) , we argue that instead of studying the loss surface, we should study the game vector field (i.e., the concatenation of each player's gradient), which can provide better insights to the problem. To this end, we propose a new visualization technique that we call Path-angle which helps us observe the nature of the game vector field close to a stationary point for high dimensional models, and carry on an empirical investigation of the properties of the optimization landscape of GANs. The core questions we want to address may be summarized as the following:

Is rotation a phenomenon that occurs when training GANs on real world datasets, and do existing training methods find local Nash equilibria?

To answer this question we conducted extensive experiments by training different GAN formulations (NSGAN and WGAN-GP) with different optimizers (Adam and ExtraAdam) on three datasets (MoG, MNIST and CIFAR10) . Based on our experiments and using our visualization techniques we observe that the landscape of GANs is fundamentally different from the standard loss surfaces of deep networks. Furthermore, we provide evidence that existing GAN training methods do not converge to a local Nash equilibrium.

Contributions More precisely, our contributions are the following: (i) We propose studying empirically the game vector field (as opposed to studying the loss surfaces of each player) to understand training dynamics in GANs using a novel visualization tool, which we call Path-angle and that captures the rotational and attractive behaviors near local stationary points (ref. 

2. RELATED WORK:

Improving the training of GANs has been an active research area in the past few years. Most efforts in stabilizing GAN training have focused on formulating new objectives (Arjovsky et al., 2017) , or adding regularization terms (Gulrajani et al., 2017; Mescheder et al., 2017; 2018) . In this work, we try to characterize the difference in the landscapes induced by different GAN formulations and how it relates to improving the training of GANs.

Recently, Nagarajan and Kolter (2017) ; Mescheder et al. (2018) show that a local analysis of the eigenvalues of the Jacobian of the game can provide guarantees on local stability properties. However, their theoretical analysis is based on some unrealistic assumptions such as the generator's ability to fully capture the real distribution. In this work, we assess experimentally to what extent these theoretical stability results apply in practice.

Rotations in differentiable games has been mentioned and interpreted by (Mescheder et al., 2018; Balduzzi et al., 2018) and Gidel et al. (2019b) . While these papers address rotations in games from a theoretical perspective, it was never shown that GANs, which are games with highly non-convex losses, suffered from these rotations in practice. To our knowledge, trying to quantify that GANs actually suffer from this rotational component in practice for real world dataset is novel.

The stable points of the gradient dynamics in general games have been studied independently by Mazumdar and Ratliff (2018) and Adolphs et al. (2018) . They notice that the locally stable stationary point of some games are not local Nash equilibria. In order to reach a local Nash equilibrium, Adolphs et al. (2018) ; Mazumdar et al. (2019) develop techniques based on second order information. In this work, we argue that reaching local Nash equilibria may not be as important as one may expect and that we do achieve good performance at a locally stable stationary point.

Several works have studied the loss landscape of deep neural networks. Goodfellow et al. (2015) proposed to look at the linear path between two points in parameter space and show that neural networks behave similarly to a convex loss function along this path. Draxler et al. (2018) proposed an extension where they look at nonlinear paths between two points and show that local minima are connected in deep neural networks. Another extension was proposed by (Li et al., 2018a) where they use contour plots to look at the 2D loss surface defined by two directions chosen appropriately. In this paper, we use a similar approach of following the linear path between two points to gain insight about GAN optimization landscapes. However, in this context, looking at the loss of both players along that path may be uninformative. We propose instead to look, along a linear path from initialization to best solution, at the game vector field, particularly at its angle w.r.t. the linear path, the Path-angle.

Another way to gain insight into the landscape of deep neural networks is by looking at the Hessian of the loss; this was done in the context of single objective minimization by (Dauphin et al., 2014; Sagun et al., 2016; 2017; Alain et al., 2019) . Compared to linear path visualizations which can give global information (but only along one direction), the Hessian provides information about the loss landscape in several directions but only locally. The full Hessian is expensive to compute and one often has to resort to approximations such has computing only the top-k eigenvalues. While, the Hessian is symmetric and thus has real eigenvalues, the Jacobian of a game vector field is significantly different since it is in general not symmetric, which means that the eigenvalues belong to the complex plane. In the context of GANs, Mescheder et al. (2017) introduced a gradient penalty and use the eigenvalues of the Jacobian of the game vector field to show its benefits in terms of stability. In our work, we compute these eigenvalues to assess that, on different GAN formulations and datasets, existing training procedures find a locally stable stationary point that is a saddle point for the loss function of the generator.

3. FORMULATIONS FOR GAN OPTIMIZATION AND THEIR PRACTICAL IMPLICATIONS:



3.1. THE STANDARD GAME THEORY FORMULATION:

From a game theory point of view, GAN training may be seen as a game between two players: the discriminator D ϕ and the generator G θ , each of which is trying to minimize its loss L D and L G , respectively. Using the same formulation as Mescheder et al. (2017) , the GAN objective takes the following form (for simplicity of presentation, we focus on the unconstrained formulation):

EQUATION

The solution (θ * , ϕ * ) is called a Nash equilibrium (NE). In practice, the considered objectives are non-convex and we typically cannot expect better than a local Nash equilibrium (LNE), i.e. a point at which (1) is only locally true (see e.g. (Adolphs et al., 2018) for a formal definition). Ratliff et al. (2016) derived some derivative-based necessary and sufficient conditions for being a LNE. They show that, for being a local NE it is sufficient to be a differential Nash equilibrium:

Definition 1 (Differential NE). A point (θ * , ϕ * ) is a differential Nash equilibrium (DNE) iff ∇ θ L G (θ * , ϕ * ) = ∇ ϕ L D (θ * , ϕ * ) = 0 , ∇ 2 θ L G (θ * , ϕ * ) 0 and ∇ 2 ϕ L D (θ * , ϕ * ) 0 (2)

where S 0 if and only if S is positive definite.

Being a DNE is not necessary for being a LNE because a local Nash equilibrium may have Hessians that are only semi-definite. NE are commonly used in GANs to describe the goal of the learning procedure (Goodfellow et al., 2014) : in this definition, θ * (resp. ϕ * ) is seen as a local minimizer of

L G (•, ϕ * ) (resp. L D (θ * , •)).

Under this view, however, the interaction between the two networks is not taken into account. This is an important aspect of the game stability that is missed in the definition of DNE (and Nash equilibrum in general). We illustrate this point in the following section, where we develop an example of a game for which gradient methods converge to a point which is a saddle point for the generator's loss and thus not a DNE for the game.

3.2. AN ALTERNATIVE FORMULATION BASED ON THE GAME VECTOR FIELD:

In practice, GANs are trained using first order methods that compute the gradients of the losses of each player. Following Gidel et al. (2019a) , an alternative point of view on optimizing GANs is to jointly consider the players' parameters θ and ϕ as a joint state ω := (θ, ϕ), and to study the vector field associated with these gradients,2 which we call the game vector field

v(ω) := ∇ θ L G (ω) ∇ ϕ L D (ω)

where ω := (θ, ϕ) .

(3) With this perspective, the notion of DNE is replaced by the notion of locally stable stationary point (LSSP). Verhulst (1989, Theorem 7 .1) defines a LSSP ω * using the eigenvalues of the Jacobian of the game vector field ∇v(ω * ) at that point.

Definition 2 (LSSP). A point ω * is a locally stable stationary point (LSSP) iff

EQUATION

)

where denote the real part of the eigenvalue λ belonging to the spectrum of ∇v(ω * ).

This definition is not easy to interpret but one can intuitively understand a LSSP as a stationary point (a point ω * where v(ω * ) = 0) to which all neighbouring points are attracted. We will formalize this intuition of attraction in Proposition 1. In our two-player game setting, the Jacobian of the game vector field around the LSSP has the following block-matrices form:

EQUATION

)

When B = -A , being a DNE is a sufficient condition for being of LSSP (Mazumdar and Ratliff, 2018) . However, some LSSP may not be DNE (Adolphs et al., 2018) , meaning that the optimal generator θ * could be a saddle point of L G (•, ϕ * ), while the optimal joint state (θ * , ϕ * ) may be a LSSP of the game. We summarize these properties in Table 1 . In order to illustrate the intuition behind this counter-intuitive fact, we study a simple example where the generator is 2D and the discriminator is 1D.

Example 1. Let us consider L G as a hyperbolic paraboloid (a.k.a., saddle point function) centered in (1, 1) where (1, ϕ) is the principal descent direction and (-ϕ, 1) is the principal ascent direction, while L D is a simple bilinear objective. 1b . Note that the discriminator ϕ controls the principal descent direction of L G .

L G (θ 1 , θ 2 , ϕ) = (θ 2 -ϕθ 1 -1) 2 -1 2 (θ 1 + ϕθ 2 -1) 2 , L D (θ 1 , θ 2 , ϕ) = ϕ(5θ 1 + 4θ 2 -9) We plot L G in Fig.

We show (see:

§ A.2) that (θ * 1 , θ * 2 , ϕ * ) = (1, 1, 0

) is a locally stable stationary point but is not a DNE: the generator loss at the optimum

(θ 1 , θ 2 ) → L G (θ 1 , θ 2 , ϕ * ) = θ 2 2 -1 2 θ 2

1 is not at a DNE because it has a clear descent direction, (1, 0). However, if the generator follows this descent direction, the dynamics will remain stable because the discriminator will update its parameter, rotating the saddle and making (1, 0) an ascent direction. We call this phenomenon dynamic stability: the loss L G (•, ϕ * ) is unstable for a fixed ϕ * but becomes stable when ϕ dynamically interacts with the generator around ϕ * .

A mechanical analogy for this dynamic stability phenomenon is a ball in a rotating saddle-even though the gravity pushes the ball to escape the saddle, a quick enough rotation of the saddle would trap the ball at the center (see (Thompson et al., 2002) for more details). This analogy has been used to explain Paul's trap (Paul, 1990 ): a counter-intuitive way to trap ions using a dynamic electric field. In Example 1, the parameter ϕ explicitly controls the rotation of the saddle. This example illustrates the fact that the DNE corresponds to a notion of static stability: it is the stability of one player's loss given the other player is fixed. Conversely, LSSP captures a notion of dynamic stability that considers both players jointly. By looking at the game vector field we capture these interactions. Fig. 1b only captures a snapshot of the generator's loss surface for a fixed ϕ and indicates static instability (the generator is at a saddle point of its loss). In Fig. 1a , however, one can see that, starting from any point, we will rotate around the stationary point (ϕ * , θ * 1 ) = (0, 1) and eventually converge to it. The visualization of the game vector field reveals an interesting behavior that does not occur in single objective minimization: close to a LSSP, the parameters rotate around it. Understanding this phenomenon is key to grasp the optimization difficulties arising in games. In the next section, we formally characterize the notion of rotation around a LSSP and in §4 we develop tools to visualize it in high dimensions. Note that gradient methods may converge to saddle points in single objective minimization, but these are not stable stationary points, unlike in our game example.

3.3. ROTATION AND ATTRACTION AROUND LOCALLY STABLE STATIONARY POINTS IN GAMES:

In this section, we formalize the notions of rotation and attraction around LSSP in games, which we believe may explain some difficulties in GAN training. The local stability of a LSSP is characterized by the eigenvalues of the Jacobian ∇v(ω * ) because we can linearize v(ω) around ω * :

EQUATION

If we assume that ( 6) is an equality, we have the following theorem.

Proposition 1. Let us assume that (6) is an equality and that ∇v(ω * ) is diagonalizable, then there exists a basis P such that the coordinates ωj (t) := [P (ω(t)ω * )] j where ω(t) is a solution of (6) have the following behavior: for λ j ∈ Sp ∇v(ω * ) we have, 1. If λ j ∈ R, we observe pure attraction: ωj (t) = e -λj t ωj (0) .

2. If (λ j ) = 0, we observe pure rotation: ωj (t) ωj+1 (t) = cos |λ j t| sin |λ j t| sin |λ j t| cos |λ j t| ωj (0) ωj+1 (0) .

3. Otherwise, we observe both: ωj (t) ωj+1 (t) = e -Re(λj )t cos Im(λ j t) sin Im(λ j t) sin Im(λ j t) cos Im(λ j t) ωj (0) ωj+1 ( 0) .

Note that we re-ordered the eigenvalues such that the complex conjugate eigenvalues form pairs: if

λ j / ∈ R then λ j+1 = λj .

Matrices in 2. and 3. are rotations matrices. They induce a rotational behavior illustrated in Fig 1a .

This proposition shows that the dynamics of ω(t) can be decomposed in a particular basis into attractions and rotations over components that do not interact between each other. Rotation does not appear in single objective minimization around a local minimum, because the eigenvalues of the Hessian of the objective are always real. Mescheder et al. (2017) discussed that difficulties in training GANs may be a result of the imaginary part of the eigenvalues of the Jacobian of the game vector field and Gidel et al. (2019b) mentioned that games have a natural oscillatory behavior. This cyclic behavior has been explained in (Balduzzi et al., 2018 ) by a non-zero Hamiltonian component in the Helmholtz decomposition of the Jacobian of the game vector field. All these explanations are related to the spectral properties of this Jacobian. The goal of Proposition 1 is to provide a formal definition to the notions of rotation and attraction we are dealing with in this paper.

In the following section, we introduce a new tool in order to assess the magnitude of the rotation around a LSSP compared to the attraction to this point.

4. VISUALIZATION FOR THE VECTOR FIELD LANDSCAPE:

Neural networks are parametrized by a large number of variables and visualizations are only possible using low dimensional plots (1D or 2D). We first present a standard visualization tool for deep neural network loss surfaces that we will exploit in §4.2.

4.1. STANDARD VISUALIZATIONS FOR THE LOSS SURFACE:

One way to visualize a neural network's loss landscape is to follow a parametrized path ω(α) that connects two parameters ω, ω (often one is chosen early in learning and another one is chosen late in learning, close to a solution). A path is a continuous function ω(•) such that ω(0) = ω and ω(1) = ω . Goodfellow et al. (2015) considered a linear path ω(α) = αω + (1α)ω . More complex paths can be considered to assess whether different minima are connected (Draxler et al., 2018) .

4.2. PROPOSED VISUALIZATION: PATH-ANGLE:

We propose to study the linear path between parameters early in learning and parameters late in learning. We illustrate the extreme cases for the game vector field along this path in simple examples in Figure 2(a-c ): pure attraction occurs when the vector field perfectly points to the optimum (Fig. 2a ) and pure rotation when the vector field is orthogonal to the direction to the optimum (Fig. 2b ). In practice, we expect the vector field to be in between these two extreme cases (Fig. 2c ). In order to determine in which case we are, around a LSSP, in practice, we propose the following tools.

Path-norm. We first ensure that we are in a neighborhood of a stationary point by computing the norm of the vector field. Note that considering independently the norm of each player may be misleading: even though the gradient of one player may be close to zero, it does not mean that we are at a stationary point since the other player might still be updating its parameters.

Path-angle. Once we are close to a final point ω , i.e., in a neighborhood of a LSSP, we propose to look at the angle between the vector field (3) and the linear path from ω to ω . Specifically, we monitor the cosine of this angle, a quantity we call Path-angle:

EQUATION

Usually [a, b] = [0, 1], but since we are interested in the landscape around a LSSP, it might be more informative to also consider further extrapolated points around ω with b > 1.

Eigenvalues of the Jacobian. Another important tool to gain insights on the behavior close to a LSSP, as discussed in §3.2, is to look at the eigenvalues of ∇v(ω * ). We propose to compute the top-k eigenvalues of this Jacobian. When all the eigenvalues have positive real parts, we conclude that we have reached a LSSP, and if some eigenvalues have large imaginary parts, then the game has a strong rotational behavior (Thm. 1). Similarly, we can also compute the top-k eigenvalues of the diagonal blocks of the Jacobian, which correspond to the Hessian of each player. These eigenvalues can inform us on whether we have converged to a LSSP that is not a LNE.

An important advantage of the Path-angle relative to the computation of the eigenvalues of ∇v(ω * ) is that it only requires computing gradients (and not second order derivatives, which may be prohibitively computationally expensive for deep networks). Also, it provides information along a whole path between two points and thus, more global information than the Jacobian computed at a single point.

In the following section, we use the Path-angle to study the archetypal behaviors presented in Thm 1.

4.3. ARCHETYPAL BEHAVIORS OF THE PATH-ANGLE AROUND A LSSP:

Around a LSSP, we have seen in ( 6) that the behavior of the vector field is mainly dictated by the Jacobian matrix ∇v(ω * ). This motivates the study of the behavior of the Path-angle c(α) where the Jacobian is a constant matrix: Notice the "bump" in path-angle (close to α = 1), characteristic of rotational dynamics.

EQUATION

Depending on the choice of S 1 , S 2 , A and B, we cover the following cases:

• S 1 , S 2 0, A = B = 0: eigenvalues are real. Thm. 1 ensures that we only have attraction.

Far from ω * , the gradient points to ω * (See Fig. 2a ) and thus c(α) = 1 for α 1 and c(α) = -1 for α 1. Since ω is not exactly ω * , we observe a quick sign switch of the Path-angle around α = 1. We plotted the average Path-angle over different approximate optima in Fig. 2a (see appendix for details).

• S 1 , S 2 = 0, A = -B : eigenvalues are pure imaginary. Thm. 1 ensures that we only have rotations. Far from the optimum the gradient is orthogonal to the direction that points to ω (See Fig. 2b ). Thus, c(α) vanishes for α 1 and α 1. Because ω is not exactly ω * , around α = 1, the gradient is tangent to the circles induced by the rotational dynamics and thus c(α) = ±1. That is why in Fig. 2b we observe a bump in c(α) when α is close to 1.

• General high dimensional LSSP (4). The dynamics display both attraction and rotation.

We observe a combination of the sign switch due to the attraction and the bump due to the rotation. The higher the bump, the closer we are to pure rotations. Since we are performing a low dimensional visualization, we actually project the gradient onto our direction of interest.

That is why the Path-angle is significantly smaller than 1 in Fig. 2c .

5. NUMERICAL RESULTS ON GANS:

Losses. We focus on two common GAN loss formulations: we consider both the original nonsaturating GAN (NSGAN) formulation proposed in Goodfellow et al. (2014) and the WGAN-GP objective described in Gulrajani et al. (2017) .

Datasets. We first propose to train a GAN on a toy task composed of a 1D mixture of 2 Gaussians (MoG) with 10,000 samples. For this task both the generator and discriminator are neural networks with 1 hidden layer and ReLU activations. We also train a GAN on MNIST, where we use the DCGAN architecture (Radford et al., 2016) with spectral normalization(see §C.2 for details). Finally we also look at the optimization landscape of a state of the art ResNet on CIFAR10 (Krizhevsky and Hinton, 2009) .

Optimization methods. For the mixture of Gaussian (MoG) dataset, we used the full-batch extragradient method (Korpelevich, 1976; Gidel et al., 2019a) . We also tried to use standard batch gradient descent, but this led to unstable results indicating that gradient descent might indeed be unable to converge to stable stationary points due to the rotations (see §C.4). On MNIST and CIFAR10, we tested both Adam (Kingma and Ba, 2015) and ExtraAdam (Gidel et al., 2019a) . The observations made on models trained with both methods are very similar. ExtraAdam gives slightly better performance in terms of inception score (Salimans et al., 2016) , and Adam sometimes converge to unstable points, thus we decided to only include the observations on ExtraAdam, for more details on the observations on Adam (see §C.5). As recommended by Heusel et al. (2017) , we chose different learning rates for the discriminator and the generator. All the hyper-parameters and precise details about the experiments can be found in §C.1.

5.1. EVIDENCE OF ROTATION AROUND LOCALLY STABLE STATIONARY POINTS IN GANS:

We first look, for all the different models and datasets, at the path-angles between a random initialization (initial point) and the set of parameters during training achieving the best performance (end point) (Fig. 3 ), and at the eigenvalues of the Jacobian of the game vector field for the same end point (Fig. 4 ). We're mostly interested in looking at the optimization landscape around LSSPs, so we first check if we are actually close to one. To do so we look at the gradient norm around the end point, this is shown by the orange curves in Fig. 3 , we can see that the norm of the gradient is quite small for all the models meaning that we are close to a stationary point. We also need to check that the point is stable, to do so we look at the eigenvalues of the Game in Fig. 4 , if all the eigenvalues have positive real parts then the point is also stable. We observe that most of the time, the model has reached a LSSP. However we can see that this is not always the case, for example in Fig. 4d some of the eigenvalues have a negative real part. We still include those results since although the point is unstable it gives similar performance to a LSSP.

Our first observation is that all the GAN objectives on both datasets have a non zero rotational component. This can be seen by looking at the Path-angle in Fig. 3 , where we always observe a bump, and this is also confirmed by the large imaginary part in the eigenvalues of the Jacobian in Fig. 4 . The rotational component is clearly visible in Fig. 3d , where we see no sign switch and a clear bump similar to Fig. 2b . On MNIST and CIFAR10, with NSGAN and WGAN-GP (see Fig. 3 ), we observe a combination of a bump and a sign switch similar to Fig. 2c . Also Fig. 4 clearly shows the existence of imaginary eigenvalues with large magnitude. Fig. 4c and 4e . We can see that while almost all models exhibit rotations, the distribution of the eigenvalues are very different. In particular the complex eigenvalues for NSGAN seems to be much more concentrated on the imaginary axis while WGAN-GP tends to spread the eigenvalues towards the right of the imaginary axis Fig. 4e . This shows that different GAN objectives can lead to very different landscapes, and has implications in terms of optimization, in particular that might explain why WGAN-GP performs slightly better than NSGAN. LNE. Results are shown in Fig. 5 and Fig. 6 , from which we make several observations. First, we see that the generator never reaches a local minimum but instead finds a saddle point. This means that the algorithm converges to a LSSP which is not a LNE, while achieving good results with respect to our evaluation metrics. This raises the question whether convergence to a LNE is actually needed or if converging to a LSSP is sufficient to reach a good solution. We also observe a large difference in the eigenvalues of the discriminator when using the WGAN-GP v.s. the NSGAN objective. In particular, we find that the discriminator in NSGAN converges to a solution with very large positive eigenvalues compared to WGAN-GP. This shows that the discriminator in NSGAN converges to a much sharper minimum. This is consistent with the fact that the gradient penalty acts as a regularizer on the discriminator and prevents it from becoming too sharp.

6. DISCUSSION:

Across different GAN formulations, standard optimization methods and datasets, we consistently observed that GANs do not converge to local Nash equilibria. Instead the generator often ends up being at a saddle point of the generator loss function. However, in practice, these LSSP achieve really good generator performance metrics, which leads us to question whether we need a Nash equilibrium to get a generator with good performance in GANs and whether such DNE with good performance does actually exist. Moreover, we have provided evidence that the optimization landscapes of GANs typically have rotational components specific to games. We argue that these rotational components are part of the reason why GANs are challenging to train, in particular that the instabilities observed during training may come from such rotations close to LSSP. It shows that simple low dimensional examples, such as for instance Dirac GAN, does capture some of the arising challenges for training large scale GANs, thus, motivating the practical use of method able to handle strong rotational components, such as extragradient (Gidel et al., 2019a) , averaging (Yazıcı et al., 2019) , optimism (Daskalakis et al., 2018) or gradient penalty based methods (Mescheder et al., 2017; Gulrajani et al., 2017) .

ACKNOWLEDGMENTS.

A PROOF OF THEOREMS AND PROPOSITIONS:

A.1 PROOF OF THEOREM 1

Let us recall the theorem of interest:

Proposition' 1. Let us assume that (6) is an equality and that ∇v(ω * ) is diagonalizable, then there exists a basis P such that the coordinates ω(t) := P (ω(t)ω * ) have the following behavior, 1. For λ j ∈ Sp ∇v(ω * ), λ j ∈ R, we observe pure attraction: ωj (t) = e -λj t [ ωj (0) .

2. For λ j ∈ Sp ∇v(ω * ), (λ j ) = 0, we observe pure rotation: ωj (t) ωj+1 (t) = R |λj |t ωj (0) ωj+1 (0) .

3. Otherwise, we observe both: ωj (t) ωj+1 (t) = e -Re(λj )t R Im(λj )t ωj (0) ωj+1 ( 0) .

The matrix R ϕ corresponds to a rotation of angle ϕ. Note that, we re-ordered the eigenvalues such that the complex conjugate eigenvalues form pairs: if λ j / ∈ R then λ j+1 = λj .

Proof. The ODE we consider is,

EQUATION

The solution of this ODE is ω(t) = e -(t-t0)∇v(ω * ) (ω(t 0 )ω * ) + ω * (10)

Let us now consider λ an eigenvalue of Sp(∇v(ω * )) such that Re(λ) > 0 and Im(λ) = 0. Since ∇v(ω * ) is a real matrix and Im(λ) = 0 we know that the complex conjugate λ of λ belongs to Sp(∇v(ω * )). Let u 0 be a complex eigenvector of λ, then we have that, ∇v(ω * )u 0 = λu 0 ⇒ ∇v(ω * )ū 0 = λū 0 (11) and thus ū0 is a eigenvector of λ. Now if we set u 1 := u 0 + ū0 and iu 2 := u 0 -ū0 , we have that e -t∇v(ω * ) u 1 = e -tλ u 0 + e -t λ ū0 = Re(e -tλ )u 1 + Im(e -tλ )u 2 (12)

e -t∇v(ω * ) iu 2 = e -tλ u 0e -t λ ū0 = i(Re(e -tλ )u 2 -Im(e -tλ )u 1 )

Thus if we consider the basis that diagonalizes ∇v(ω * ) and modify the complex conjugate eigenvalues in the way we described right after 11 we get the expected diagonal form in a real basis. Thus there exists P such that ∇v(ω * ) = P DP -1 (14) where D is the block diagonal matrix with the block described in Theorem 1.

A.2 BEING A DNE IS NEITHER NECESSARY OR SUFFICIENT FOR BEING A LSSP:

Let us first recall Example 1.

Example' 1. Let us consider L G as a hyperbolic paraboloid (a.k.a., saddle point function) centered in (1, 1) where (1, ϕ) is the principal descent direction and (-ϕ, 1) is the principal ascent direction, while L D is a simple bilinear objective.

L G (θ 1 , θ 2 , ϕ) = (θ 2ϕθ 1 -1) 2 -1 2 (θ 1 + ϕθ 2 -1) 2 , L D (θ 1 , θ 2 , ϕ) = ϕ(5θ 1 + 4θ 2 -9)

We want to show that (1, 1, 0) is a locally stable stationary point.

Proof. The game vector field has the following form, 15) 

EQUATION

:

Note that, in practice, the joint vector field (3) is not a gradient vector field, i.e., it cannot be rewritten as the gradient of a single function.

Figures and Tables:

2D projection of the vector field. (b) Landscape of the generator loss.

Figure 1: Visualizations of Example 1. Left: projection of the game vector field on the plane θ2 = 1. Right:Generator loss. The descent direction is (1, ϕ) (in grey). As the generator follows this descent direction, the discriminator changes the value of ϕ, making the saddle rotate, as indicated by the circular black arrow.

Figure2: Above: game vector field (in grey) for different archetypal behaviors. The equilibrium of the game is at (0, 0). Black arrows correspond to the directions of the vector field at different linear interpolations between two points: • and . Below: path-angle c(α) for different archetypal behaviors (right y-axis, in blue). The left y-axis in orange correspond to the norm of the gradients. Notice the "bump" in path-angle (close to α = 1), characteristic of rotational dynamics.

Figure 3: Path-angle for NSGAN (top row) and WGAN-GP (bottom row) trained on the different datasets, see Appendix C.3 for details on how the path-angle is computed. For MoG the ending point is a generator which has learned the distribution. For MNIST and CIFAR10 we indicate the Inception score (IS) at the ending point of the interpolation. Notice the "bump" in path-angle (close to α = 1.0), characteristic of games rotational dynamics, and absent in the minimization problem (d). Details on error bars in §C.3.

Figure 5: NSGAN. Top k-Eigenvalues of the Hessian of each player (in terms of magnitude) in descending order. Top Eigenvalues indicate that the Generator does not reach a local minimum but a saddle point (for CIFAR10 actually both the generator and discriminator are at saddle points). Thus the training algorithms converge to LSSPs which are not Nash equilibria.

Summary of the implications between Differentiable Nash Equilibrium (DNE) and a locally stable stationnary point (LSSP): in general, being a DNE is neither necessary or sufficient for being a LSSP.

Table content: <table/>

Peer review (iclr2020):
Summary: This paper proposes visualization techniques for the optimization landscape in GANs. The primary tool presented in this paper is a quantity called path-angle, which looks at the angle between the game vector field and the linear path between a point away from a stationary point and a point near a stationary point. The paper present examples of the visualization for dynamics with pure attraction, pure rotation, and a mix of attraction and rotation. Along with this, the authors propose to look at the eigenvalues of the game Jacobian and the individual player Hessian’s to evaluate convergence in GANs. The paper presents application of the tools on GANs trained with NSGAN and WGAN-GP objectives on a mixture of Gaussians, MNIST, and CIFAR10. The primary observation is that the generator performance is good, but the algorithms converge to non-Nash stable attractors. Moreover, it is shown using the path-angle plots that GANs exhibit rotational behavior around stable points. Review: There has been a lot of work in the past few years (and ongoing) on principled training approaches for GANs. The objective of the algorithms is typically to converge a differential Nash equilibrium and/or to reach a stable point of the dynamics quickly. In my view, this work fills some of the gap on the empirical side of things with respect to each goal. Notably, a main idea to speed up convergence in GANs is to change the gradient play dynamics so rotational components are neutralized. The path angle visualization provides a novel tool to evaluate the empirical ability of any dynamics proposed for GANs to cancel out rotational components. Since it is generally known that gradient play dynamics are susceptible to cycling, I would have been interesting in seeing the path angle plots for some recently proposed algorithms such as consensus, symplectic gradient adjustment, stable opponent shaping, local symplectic surgery, etc to see how they compare. This would have made the experiments using the path angle visualization stronger in my view. Nonetheless, the path angle tool is useful and I can foresee it being commonly used in the future. Aside from neutralizing rotational components, dynamics have been proposed with the goal of avoiding non-Nash stable attractors and converging only to differential Nash equilibria. However, to my knowledge, there has not been much, if any, evaluation in GANs to see if the methods are in fact converging to Nash equilibria as theory may predict. While simple, I found it interesting to evaluate the eigenvalues of the relevant quantities at convergence. I am curious why the authors evaluate the top-k eigenvalues in terms of magnitude? The scipy package referenced in the appendix can compute the largest and smallest real eigenvalues, which is what it seems like you would want to evaluate the definiteness of the game Jacobian and the individual player Hessians. The most interesting empirical result in the paper to me was that it is common to converge to non-Nash stable attractors using standard training techniques and at such stable points the generator performance is strong. This is an important observation and may cause some consideration of what points should be sought in GANs. I am not fully convinced this is always what the dynamics would always converge to depending on the network, learning rates, optimization methods, etc, but showing that it can be the case is useful. Overall, I think this paper introduces some useful tools to interpret the performance in GANs and to help understand the behavior of training dynamics. The main tool introduced was the path angle visualization and the primary empirical result was that standard GAN methods may reach non-Nash stable attractors and perform well. The paper probably be condensed in the first 4 pages, so that more experimental results could be presented and this would make the paper stronger. Post Response: Thanks for the response. I believe this paper should be accepted.