arXiv:1705.02894v2 [stat.ML] 9 May 2017

Geometric GAN
Jae Hyun Lim1 , Jong Chul Ye2,3
1
ETRI, South Korea
jaehyun.lim@etri.re.kr
2
Dept. of Bio and Brain engineering, KAIST, South Korea
3
Dept. of Mathematical Sciences, KAIST, South Korea
jong.ye@kaist.ac.kr

Abstract
Generative Adversarial Nets (GANs) represent an important milestone for effective generative models, which has inspired numerous variants seemingly different
from each other. One of the main contributions of this paper is to reveal a unified
geometric structure in GAN and its variants. Specifically, we show that the adversarial generative model training can be decomposed into three geometric steps:
separating hyperplane search, discriminator parameter update away from the separating hyperplane, and the generator update along the normal vector direction of
the separating hyperplane. This geometric intuition reveals the limitations of the
existing approaches and leads us to propose a new formulation called geometric
GAN using SVM separating hyperplane that maximizes the margin. Our theoretical analysis shows that the geometric GAN converges to a Nash equilibrium
between the discriminator and generator. In addition, extensive numerical results
show that the superior performance of geometric GAN.

1

Introduction

Recently, inspired by the success of the deep discriminative models, Goodfellow et al [1] proposed
a novel generative model training method called generative adversarial nets (GAN). GAN is formulated as a minimax game between a generative network (generator) that maps a random vector into
the data space and a discriminative network (discriminator) trying to distinguish the generated samples from real samples. Unlike the classical generative models such as Variational Auto-Encoders
(VAEs) [2], the minimax formation of GAN can transfer the success of deep discriminative models
to generative models, resulting in significant improvement in generative model performance [1].
Specifically, the original form of the GAN solves the following minmax game:
min max LGAN (D, G)
G

D

(1)

where
LGAN (D, G)

:= Ex∼PS [log D(x)] + Ez∼PZ [log(1 − D(G(z))] ,

(2)

where PS is the sample distribution; D(x) is the discriminator that takes x ∈ X as input and outputs
a scalar between [0, 1]; G(z) is the generator that maps a sample z drawn from a distribution PZ
to the input space X . The meaning of (1) is that the generator tries to fool out the discriminator
while the discriminator wants to maximize the differentiation power between the true and generated
samples. The authors further showed that the GAN training is indeed to approximate the minimization of the symmetric Jensen-Shannon divergence [1]. This idea has been generalized by the authors
in [3] for all f -divergences. Moreover, Maximum Mean Discrepancy objective (MMD) for GAN
training was also proposed in [4, 5].
It is well-known that the training GAN is difficult. In particular, the authors in [6] have identified the following sources of the difficulties: 1) when the discriminator becomes accurate, the

gradient for generator vanishes, 2) a popular fixation using a generator gradient updating with
Ez∼PZ [− log D(G(z))] is unstable because of the singularity at the denominator when the discriminator is accurate. The main motivation of Wasserstein GAN (W-GAN) [7] was, therefore, to
introduce the weight clipping to address the above-described limitations. In fact, Wasserstein GAN
is a special instance of minimizing the integral probability metric (IPM) [8], and Mroueh et al [9]
recently generalized the W-GAN for wider function classes and proposed the mean feature matching
and/or covariance feature matching GAN (McGAN) using the IPM minimization framework [9].
Inspired by McGAN, here we propose a novel geometric generalization called geometric GAN.
Specifically, geometric GAN is inspired by our novel observation that McGAN is composed of
three geometric operations in feature space:
• Separating hyperplane search: finding the separating hyperplane for a linear classifier
[10, 11, 12]
• Discriminator update away from the hyperplane: discriminator parameter update away
from the separating hyperplane using stochastic gradient direction (SGD).
• Generator update toward the hyperplane: generator parameter update along the normal
vector direction of the separating hyperplane using stochastic gradient direction (SGD).
This geometric interpretation proves to be very general, so it can be applied to most of the existing
GAN and its variants. Indeed, the main differences between the algorithms come from the construction of the separating hyperplanes for a linear classifier on feature space and the geometric scaling
factors for the feature vectors. Based on this observation, we provide new geometric interpretations
of GAN [1], f -GAN [3], EB-GAN [13], and W-GAN [7] in terms of separating hyperplanes and
geometric scaling factors, and discuss their limitations. Furthermore, we propose a novel geometric GAN using the support vector machine (SVM) separating hyperplane that has maximal margin
between two classes of separable data [14, 11]. Our numerical experiments clearly showed that the
proposed geometric GAN outperforms the existing GANs in all data set.

2

Related Approaches

In order to introduce the geometric interpretation of GAN and its variants, we begin with the review
of the mean feature matching GAN (McGAN) [9].
2.1

Mean feature matching GAN

Let F be a set of bounded real valued functions on the sample space X . Suppose that P and Q are
two probability distributions on X . Then, the integral probability metric (IPM) between P and Q
on the function space F is defined as follows [8, 15, 16]:
Z
Z
f dP − f dQ
dF (P, Q) = sup
f ∈F

=

sup |Ex∼P [f (x)] − Ex∼Q [f (x)]|

(3)

f ∈F

where Ex∼P [·] denotes the expectation with respect to the probability distribution P (x). We can
easily show that dF is non-negative, symmetric and satisfies the triangle inequality. So dF can be
used as a distance measure in the probability space. For example, when F is defined as a collection
of functions with a finite Lipschitz constant, the IPM is the Wasserstein distance or the earth mover’s
distance [15, 16] that forms the basis of the Wasserstein GAN [7].
In McGAN, the generator network gθ : Z → X that maps a random input z ∈ Z to a target x ∈ X
is shown as a block in Fig. 1(a), and the function space F under study is defined as follows [9]:
Fw,ζ = {f (x) = hw, Φζ (x)i | w ∈ Ξ, kwk2 ≤ 1}
where Φζ : X → Ξ is a bounded map from X to a (often higher-dimensional) feature space Ξ. Note
that Fw,ζ is a symmetric function spaces because if f ∈ Fw,ζ , then −f ∈ Fw,ζ . Then, the IPM
between Px and Pg is given by
dF (Px , Pg )

=

max hw, Ex∼Px Φζ (x) − Egθ ∼Pg Φζ (gθ )i

kwk2 ≤1

2

Figure 1: Structure of the mean feature matching GAN and its extension to geometric GAN.
under the constraints that Φζ is bounded and kwk2 ≤ 1.
Now, given a finite sequence of mini-batch training data set S = {(z1 , x1 ), · · · , (zn , xn )}, an empirical estimate of the minmax game that minimizes the dF (Px , Pg ) is given by [9]:
min max

L̂(w, ζ, θ)

kwk2 ≤1,ζ

θ

(4)
(5)

where

*
L̂(w, ζ, θ) :=

+
n
n
1X
1X
Φζ (xi ) −
Φζ (gθ (zi ))
w,
n i=1
n i=1

Then, the discriminator update can be done using a stochastic gradient descent (SGD) [9]:


(w, ζ) ← (w, ζ) + η ∇w L̂(w, ζ, θ), ∇ζ L̂(w, ζ, θ)

(6)

where η is a learning rate. The authors in [9] also used the projection onto the unit lp ball and
weight clipping for w and ζ updates, respectively, to meet the constraints. Using the updated w, the
generator update is then given by [9]:
θ

← θ+η

n
X

hw, ∇θ Φ(gθ (zi ))i /n .

(7)

i=1

2.2

Geometric interpretation of the mean feature matching GAN

The primal form of the w update in (6) is geometrically less informative, so here we use the CauchySwartz inequality to obtain a closed-form update for w:
w∗ = c

n
X

(Φζ (xi ) − Φζ (gθ (zi ))) /n

(8)

i=1

Pn
1
where the constant c is given by c = k i=1 (Φζ (xi ) − Φζ (gθ (zi ))) /nk 2 . Given w∗ , the corresponding discriminator and generator updates are represented by
ζ
θ

← ζ +η
← θ+η

n
X
i=1
n
X

hw∗ , ∇ζ Φζ (xi ) − ∇ζ Φζ (gθ (zi ))i/n
hw∗ , ∇θ Φζ (gθ (zi ))i /n

(9)
(10)

i=1

Note that the update equations (8) to (10) are equivalent to the dual form of the McGAN [9], where
the authors derived the following minmax problem by using (8):
n
n
X
1 X
Φζ (xi ) −
Φζ (gθ (zi ))
min max
θ
ζ
2 i=1
i=1

2

.

(11)

However, we notice that the explicit representation by Eqs. (8) to (10) gives clearer geometric
intuition that plays the key role in designing a geometric GAN, ass will become clear soon.
3

More specifically, in designing a linear classifier for two class classification problems, (8) is known
as the normal vector for the separating hyperplane for the mean difference (MD) classifier [10, 11,
12]. As shown in Fig. 2, once the separating hyperplane is defined, the SGD udpate (9) is to update
the discriminator parameters such that the true and fake samples are maximally separately away
from the separating hyperplane parallel to the normal vector. On the other hand, the SGD udpate
using (10) is to update the generator parameters to make the fake samples approach the separating
hyperplane along the normal vector direction (see Fig. 2).

Figure 2: Geometry of the mean feature matching GAN.
Recall that a linear classifier is defined via the normal vector to the separating hyperplane and an
offset [14]. Thus, comparing the direction between two classifiers means comparing their normal
vector directions. As shown in Appendix, aside from the geometric scaling factors, the existing GAN
and its variants mainly differ in their definition of the normal vector for the separating hyperplane.
Based on this observation, in the next section, we discuss an optimal separating hyperplane for
generative model training that has the maximal margin.

3
3.1

Geometric GAN
Linear classifiers in high-dimension low-sample size

In adversarial training in feature space, a discriminator is interested in discriminating true samples
{Φζ (xi )}ni=1 and the fake samples {Φζ (gθ (zi ))}ni=1 . In practice, the minibatch size n is much
smaller than the dimension of the feature space d, and this type of classification problem is often
called the high-dimension low-sample size (HDLSS) problem [10, 11, 12].
In fact, the mean difference (MD) classifier is one of the popular methods for HDLSS. Specifically,
the MD classifier selects the hyperplane that lies half way between the two class means. In particular
the normal vector w for the seperating hyperplane is given by the difference of the class means:
n

wM D

=

n

1X
1X
Φζ (xi ) −
Φζ (gθ (zi ))
n i=1
n i=1

(12)

Note that if the variables are first mean centered then scaled by the standard deviation, then the
mean difference is equivalent to the naive Bayes classifier [10, 11, 12]. Furthermore, in HDLSS,
there always exists a maximal data piling direction (MDP) [12], where mulitple points in each class
have the identical projection on the line spanned by the normal vector.
On the other hand, the Support Vector Machine (SVM) [14] and its many variants is one of the most
widely used and well studied classification algorithms and its robustness has been also proven for
HDLSS setup [10, 11, 12]. Although the aforementioned classification algorithms are motivated
4

by fitting a statistical distribution to the data, SVM is motivated by a geometric heuristic that leads
directly to an optimization problem: maximize the margin between two classes of separable data.
In addition, soft-margin SVM balances two competing objectives to maximize the margin while
penalizing points on the wrong side of the margin.
Recently, Carmichael et al [11] investigated the Karush-Kuhn-Tucker conditions to provide rigorous
mathematical proof for new insights into the behaviour of soft-margin SVM in the large and small
tuning parameter regimes in HDLSS. They revealed that for small tuning parameter, if the number
of data in two classes are the same (which is the case in our problem), then the SVM direction
becomes exactly the MD direction. In addition, for sufficiently large tuning parameter, the authors
showed that soft margin SVM is equivalent to hard margin SVM if the data are separable, and the
hard-margin SVM has data piling. Due to this generality of the soft-margin SVM, the proposed
geometric GAN is designed based on soft-margin SVM linear classifier.
3.2

Geometric GAN with SVM hyperplane

Note that soft-margin SVM is designed by adding a tunning parameter C and slack variables ξi
which allows points to be on the wrong side of the margin [14]. In our problem classifying the true
samples versus fake samples, the primal form of soft-margin SVM can be formulated by
P
1
2
0
min
i (ξi + ξi )
2 kwk + C
w,b

subject to

hw, Φζ (xi )i + b ≥ 1 − ξi ,
i = 1, · · · , n
0
hw, Φζ (gθ (zi ))i + b ≤ ξi − 1, i = 1, · · · , n
ξi , ξi0 ≥ 0, i = 1, · · · , n

Equivalently, the primal form of the soft-margin SVM can be represented using loss + penalty form
[17, 14]:
min Rθ (w, b; ζ)

(13)

w,b

where
n

Rθ (w, b; ζ)

=

1X
1
kwk2 +
max (0, 1 − hw, Φζ (xi )i − b)
2Cn
n i=1
n

1X
+
max (0, 1 + hw, Φζ (gθ (zi ))i + b)
n i=1

(14)

The goal of the SVM optimization (13) is to maximize the margin between the two classes. This
implies that the discriminator update can be also easily incorporated with SVM update, because
the goal of the discriminator update can be also regarded to maximize the margin between the two
classes. More specifically, our optimization problem is given by
min Rθ (w, b; ζ)

(15)

w,b,ζ

for a given generator parameter θ.
Specifically, in SVM, the normal vector for the optimal separating hyperplane wSV M from (13) is
given by [11, 14]:
wSV M :=

n
X

αi Φζ (xi ) −

i=1

n
X

βi Φζ (gθ (zi ))

(16)

i=1

where (αi , βi ) will be nonzero only for the support vectors, where the set of support vectors now
includes all data points on the margin boundary as well as those on the wrong side of the margin
boundary (see Fig. 3). More specifically, we define the region M between the margin boudaries as
shown in Fig. 3(a):
M

= {φ ∈ Ξ | |hwSV M , φi + b| ≤ 1} .
5

(17)

Then, for given wSV M and b, the cost function (14) then becomes
1 X SV M
1 X SV M
hw
hw
Rθ (w, b; ζ) :=
, Φζ (gθ (zi ))i −
, Φζ (xi )i + constant
n
n
=

1
n

i∈IS
n
X

i∈IT

hwSV M , si Φζ (gθ (zi )) − ti Φζ (xi )i + constant

(18)

i=1

where (ti , si ) are geometric scaling factors defined by


1, Φζ (xi ) ∈ M
1, Φζ (gθ (zi )) ∈ M
ti =
, si =
0, otherwise
0, otherwise

(19)

This is because the SVM cost function value is not dependent on the feature vectors outside of the
margin boundaries and is now fully determined by the supporting vectors in M.
Accordingly, the discriminator update is given by following SGD updates:
ζ

←

ζ +η

n
X
hwSV M , ti ∇ζ Φζ (xi ) − si ∇ζ Φζ (gθ (zi ))/ni

(20)

i=1

In another word, to update the discriminator parameters, we only need to push out the supporting
vectors toward the margin boundaries.
On the other hand, generator update requires more geometric intuition. As shown in Fig. 3, the
generator update tries to move the fake feature vectors toward the normal vector direction of the
separating hyperplane so that they can be classified as the true feature vectors. This means that the
generator update should be given by the following minimization problem:
min Lw,b,ζ (θ)
θ

where
n

Lw,b,ζ (θ)

:= −

1X
Dw,b,ζ (gθ (zi ))
n i=1

(21)

with the linear classifier
Dw,b,ζ (x) := hw, Φζ (x)i + b .

(22)

This is because −Dw,b,ζ (x) gets smaller as Φζ (gθ (z)) moves toward the upper-left side in Fig. 3
along the normal vector direction wSV M . This results in the following SGD updates:
θ

← θ+η

n
X
hwSV M , ∇θ Φζ (gθ (zi ))i/n .

(23)

i=1

Note that the SGD updates (20) and (23) are strikingly similar to (9) and (10) of McGAN. Aside from
the different choice of separating hyperplane by (16), the discriminator update (20) has additional
geometric scaling factors (ti , si ). As will be shown in Appendix A, the appearance of the geometric
scaling factors is a recurrent theme in geometric interpretation of GAN and its variants, which we
believe is fundamental to account for the geometry of the classifiers.
3.3

Convergence of Geometric GAN

In order to show the convergence of the geometric GAN to a Nash equilibrium, we investigate the
behaviour at large sample limit. Specifically, as n → ∞ for a fixed C, the soft margin SVM cost in
(14) becomes
R(D, g)

=
=
=

Ex∼Px [max (0, 1 − Dw,b,ζ (x))] + Ez∼Pz [max (0, 1 + Dw,b,ζ (gθ (z)))]
Z
Z




1 − Dw,b,ζ (x) + dPx +
1 + Dw,b,ζ (gθ (z)) + dPz
Z
Z




px (x) 1 − Dw,b,ζ (x) + dx + pgθ (x) 1 + Dw,b,ζ (x) + dx
6

Figure 3: Geometric GAN using SVM hyperplane. Discriminator and generator update directions
are shown.
where [x]+ = max{0, x} and Dw,b,ζ (x) is a linear discriminator in (22) parameterized by (w, b, ζ).
Here, px (x) and pgθ (x) denote the probability density functions (pdf) for the distribution Px and
Pz (gθ (z)), respectively. Similarly, the generator cost function in (21) becomes
L(D, g)

−Ez∼Pz [Dw,b,ζ (gθ (z))]
Z
− pgθ (x)Dw,b,ζ (x)dx

=
=

Then, the adversarial training between discriminator and generator can be achieved by the following
alternating minimization:
min R(D, g)

:=

min L(D, g)

:=

D
g

min R(Dw,b,ζ , g)

(24)

min L(D, gθ )

(25)

w,b,ζ
θ

Suppose that the optimal solution of the aforementioned adversarial training is a pair (D∗ , g ∗ ). Then,
we can prove the following key convergence result.
Theorem 3.1. Suppose that (D∗ , g ∗ ) is a minimizer of the alternating minimization of (24) and
(25). Then, pg∗ (x) = px (x) almost everywhere, and R(D∗ , G∗ ) = 2.
Proof. See Appendix B.
In the following example, we provide a specific example where the discriminator and generator cost
function has close form expressions, which also has intuitive meaning of the minimum value 2 in
Theorem 3.1. In particular, we consider an example of learning parallel lines as in the original
Wasserstein GAN paper [7].
Example 3.1 (Learning parallel lines). Let u ∼ U [0, 1] denote the uniform distribution on the unit
interval. Let Px be the distribution of x = (0, u) ∈ R2 , uniform on a straight vertical line passing
through the origin. Suppose that the generator sample is given by gθ (z) = (θ, z) with θ a single
real parameter. We can easily see that the SVM separating hyperplane is given by
hw, xi + b = 0
7

where w = (−1, 0), b = θ/2 if θ ≥ 0 and w = (1, 0), b = −θ/2 if θ < 0. Thus, the generator cost
function becomes
L(D, g)

=
=

−E [hw, gθ (z)i + b]
|θ|/2

which achieves its minimum at θ∗ = 0. Then, the corresponding discriminator cost value at n → ∞
is given by

R(D∗ , g ∗ ) = lim E [1 − hw, xi − b]+ + E [1 + hw, gθ (z)i + b]+
θ→0

=
=

lim 2 [1 − |θ|/2]+

θ→0

2

which coincides the results by Theorem 3.1.
In this example, R(D∗ , g ∗ ) = 2 because all the true and fake samples lies on the separating hyperplane. This informs that at the Nash equilibrium of this problem, all the true samples and the fake
samples are not separable, which is the desired property of GAN training. However, Theorem 3.1
is only a necessary condition to make the true and fakes sample non-separable. The proof for the
sufficiency condition would be very interesting, which is beyond the scope of current paper.

4
4.1

Experimental Results
Mixture of Gaussians

In order to evaluate the proposed geometric GAN, we perform comparative studies with three representative types of GANs; 1) Jenson-Shannon (GAN) [1], 2) mean difference in l∞ (Wasserstein
GAN) [7], and 3) mean difference in l2 [9]. Here, the behavior of the maximum margin separating hyperplane of the geometric GAN is empirically analyzed against those of the aforementioned
approaches.
In addition, to evaluate the dependency of each variants on Lipschitz continuity constraints, the
Lipschitz constraints suggested in [7, 18] was also applied to each adversarial training approach.
More specifically, the parameters (w, b) of the final linear layer in discriminator is determined to
represent the aforementioned hyperplane properties, whereas the Lipschitz constraints are applied
for other network parameters, such as ζ in Φζ (x) and θ in gθ (z). In this paper, we only consider
Lipschitz density constraints in [18], so we follow to use weight decay on generators gθ (x) and
feature space mapping φζ (x) in discriminators.
We test the four hyperplane searching approaches for discriminators, as well as their complementary
generator losses, on two dimensional synthetic data. The synthetic data consists of 100K data points
generated from a mixture of 25 Gaussians, akin to the data that have been used for describing mode
collapsing behaviors of GANs [19, 20, 21]. Specifically, the means of the Gaussians are evenly
spaced as a 5 by 5 grid along x and y axis from -21 to 21. The standard deviation of each normal
distribution is 0.316 (so that the variance would be 0.1). The sampled data from the true distribution
can be seen in Figure 4 and 5.
For discriminator and generator, a multi-layered fully-connected neural network architecture is used,
as described below. RMSprop [22] is used to train these networks, except vanilla GAN (without any
Lipschitz constraints). For vanilla GAN, Adam [23] with momentum β1 = 0.5 is used. Base
learning rate is set to 0.001. When weight clipping is applied, parameters in feature mapping φζ (x)
is clipped within the range of [−0.01, 0.01]. When weight projection on unit l2 norm is applied, the
following rule, p = min{1, 1/kpk2 } × p described in [9] is used to update any parameter p for every
iteration. For weight decay, weight decaying parameter is set to 0.001. Batch size is set to 500 for
all experiment. For the number of discriminator update Kd and the one of generator update Kg , we
set them as 1, i.e. (Kd = 1, Kg = 1).
• Discriminator: FC(2, 128)-ReLU-FC(128, 128)-ReLU-FC(128, 128)-ReLU-FC(128, 1)
• Generator:
FC(4, 128)-BN-ReLU-FC(128, 128)-BN-ReLU-FC(128, 128)-BN-ReLUFC(128, 2)
8

(a) True data

(b) GAN [1]

(c) Proposed

(d) WGAN [7]

(e) meanGAN +
wproj on φζ [9]

Figure 4: Generated samples of GAN variants for the mixture of 25 Gaussians, and the ones from
the true data distribution.

(a) True data

(b) GAN + wdecay (c) Proposed + (d) WGAN + wde- (e) meanGAN +
on Φζ and gθ
wdecay on Φζ and cay on Φζ and gθ
wdecay on Φζ and
gθ
gθ

Figure 5: Generated samples of GAN variants under trained with Lipschitz density constraints suggested in [18] for the mixture of 25 Gaussians, and the sample from the true data distribution. During
training, weight decay was applied on ζ in Φζ (x) and θ in gθ (z) in addition to their own hyperplane
constraints.

The results of the experiment with the mixture of 25 Gaussians are illustrated in Figure 4 and 5.
Amongst all GAN variants in this experiment, geometric GAN demonstrated the least mode collapsing behavior independently with Lipschitz continuity regularization constraints.
As shown in Fig. 5, under the same Lipschitz density constraints, linear hyperplane approaches
demonstrated less mode collapsing behaviors by virtue of consistent gradients unlike nonlinear separating hyperplane of original GAN. However, mean difference-driven hyperplanes in Wasserstein
GAN or McGAN led generators to the mean of arbitrary number of modes in true distributions since
the characteristics of mean difference. One the other hand, geometric GAN generally showed robust
and consistent convergence behavior towards true distributions.
4.2

Image Datasets

In order to analyze the proposed method on large-scale dataset, the geometric GAN is empirically
analyzed on well-studied datasets in the context of adversarial training; MNIST, CelebA, and LSUN
datasets. Since consistent quantitative measures are still under debate, we only perform qualitative
comparisons of generated samples from the learned generators of the propsed method against the
results of previous literatures. In favor of fair comparisons with other adversarial training methods, we adopt the settings from the previous literatures except the hyperparameters of stochastic
optimizations and the tuning parameter of the proposed method.
The DCGAN neural network architectures [24] was used, including batch normalization for generator. Note that the currently known adversarial training methods that demonstrated stable learning
without batch normalization [7, 9, 18] have resorted to Lipschitz constraints; therefore, it can also
be applied to other adversarial training criterions, including geometric GAN, in order to train batch
normalization-free generators.
Each pixel value in input image was rescaled to [−1, 1] for all dataset, including MNIST dataset.
During all training, mini-batch size was set to 64. For stochastic gradient update during training,
RMSprop was used [22]. The number of generator’s updates per each discriminator’s update is set
to 10 (Kd = 1, Kg = 10). Learning rate is set to 0.0002 for both discriminators and generators, and
a tuning parameter C for discriminator is set to 1.
Specifically for MNIST dataset, input images were resized to 64 by 64 pixels in order to use the
same DCGAN network architecture, and the number of epochs for training was set to 20. For
CelebA dataset, input images were resized to 96 by 96 pixels and center-cropped with 64 by 64
pixels, and the number of epochs for training is set to 50. For LSUN dataset, only bedroom dataset
9

is used, and an input image is resized to 64 by 64 pixels. The number of epochs for training is set to
2 for LSUN dataset.
The results in Figure 6, 7, and 8 clearly show that the geometric GAN generates very realistic
images without mode collapsing or divergent behaviours.

Figure 6: Generated samples of Geometric GAN trained for MNIST dataset.

5

Conclusion

This paper proposed a novel geometric GAN using SVM separating hyperplane, based on geometric
intuitions revealed from previous adversarial training approaches. The geometric GAN was based
on SVM separating hyperplanes that has the maximal margins between the two classes. Compared
to the most of the existing approaches that are based on statistical design criterion, the geometric
GAN is derived based on geometric intuition similar to the derivation of SVM. Extensive numerical experiments showed that the proposed method has demonstrated less mode collapsing and more
stable training behavior. Moreover, our theoretical results showed that the proposed algorithm converges to the Nash equilibrium between the discriminator and generator, which has also geometric
meaning.

Acknowledgement
This work is supported by Korea Science and Engineering Foundation, Grant number
NRF2016R1A2B3008104. The first author would like to thank Yunhun Jang for helpful discussions.

References
[1] I. Goodfellow, J. Pouget-Abadie, M. Mirza, B. Xu, D. Warde-Farley, S. Ozair, A. Courville,
and Y. Bengio, “Generative adversarial nets,” in Advances in Neural Information Processing
Systems, 2014, pp. 2672–2680.
10

Figure 7: Generated samples of Geometric GAN trained for CelebA dataset.

Figure 8: Generated samples of Geometric GAN trained for LSUN dataset.

11

[2] D. P. Kingma and M. Welling, “Auto-encoding variational bayes,” arXiv preprint
arXiv:1312.6114, 2013.
[3] S. Nowozin, B. Cseke, and R. Tomioka, “f-GAN: Training generative neural samplers using
variational divergence minimization,” in Advances in Neural Information Processing Systems,
2016, pp. 271–279.
[4] Y. Li, K. Swersky, and R. Zemel, “Generative moment matching networks,” in Proceedings of
The 32nd International Conference on Machine Learning, 2015, pp. 1718–1727.
[5] G. Dziugaite, D. Roy, and Z. Ghahramani, “Training generative neural networks via maximum
mean discrepancy optimization,” in Uncertainty in Artificial Intelligence-Proceedings of the
31st Conference, UAI 2015, 2015, pp. 258–267.
[6] M. Arjovsky and L. Bottou, “Towards principled methods for training generative adversarial
networks,” in NIPS 2016 Workshop on Adversarial Training. In review for ICLR, vol. 2016,
2017.
[7] M. Arjovsky, S. Chintala, and L. Bottou, “Wasserstein GAN,” arXiv preprint
arXiv:1701.07875, 2017.
[8] A. Müller, “Integral probability metrics and their generating classes of functions,” Advances in
Applied Probability, vol. 29, no. 02, pp. 429–443, 1997.
[9] Y. Mroueh, T. Sercu, and V. Goel, “McGan: Mean and covariance feature matching GAN,”
arXiv preprint arXiv:1702.08398, 2017.
[10] J. S. Marron, M. J. Todd, and J. Ahn, “Distance-weighted discrimination,” Journal of the
American Statistical Association, vol. 102, no. 480, pp. 1267–1271, 2007.
[11] I. Carmichael and J. Marron, “Geometric insights into support vector machine behavior using
the KKT conditions,” arXiv preprint arXiv:1704.00767, 2017.
[12] J. Ahn and J. Marron, “The maximal data piling direction for discrimination,” Biometrika,
vol. 97, no. 1, pp. 254–259, 2010.
[13] J. Zhao, M. Mathieu, and Y. LeCun, “Energy-based generative adversarial network,” arXiv
preprint arXiv:1609.03126, 2017.
[14] B. Schölkopf and A. J. Smola, Learning with kernels: support vector machines, regularization,
optimization, and beyond. MIT press, 2002.
[15] B. K. Sriperumbudur, A. Gretton, K. Fukumizu, B. Schölkopf, and G. R. Lanckriet, “Hilbert
space embeddings and metrics on probability measures,” Journal of Machine Learning Research, vol. 11, no. Apr, pp. 1517–1561, 2010.
[16] B. K. Sriperumbudur, K. Fukumizu, A. Gretton, B. Schölkopf, and G. R. Lanckriet,
“On integral probability metrics,φ-divergences and binary classification,” arXiv preprint
arXiv:0901.2698, 2009.
[17] T. Hastie, S. Rosset, R. Tibshirani, and J. Zhu, “The entire regularization path for the support
vector machine,” Journal of Machine Learning Research, vol. 5, no. Oct, pp. 1391–1415, 2004.
[18] G.-J. Qi, “Loss-sensitive generative adversarial networks on lipschitz densities,” arXiv preprint
arXiv:1701.06264, 2017.
[19] V. Dumoulin, I. Belghazi, B. Poole, A. Lamb, M. Arjovsky, O. Mastropietro, and A. C.
Courville, “Adversarially learned inference,” arXiv preprint arXiv:1606.00704, 2016.
[20] L. Metz, B. Poole, D. Pfau, and J. Sohl-Dickstein, “Unrolled generative adversarial networks,”
arXiv preprint arXiv:1611.02163, 2016.
[21] T. Che, Y. Li, A. P. Jacob, Y. Bengio, and W. Li, “Mode regularized generative adversarial
networks,” arXiv preprint arXiv:1612.02136, 2016.
[22] T.
Tieleman
and
G.
Hinton,
“RMSprop
Gradient
Optimization,”
http://www.cs.toronto.edu/t̃ijmen/csc321/slides/lecture-slides-lec6.pdf.
[23] D. P. Kingma and J. Ba, “Adam: A method for stochastic optimization,” arXiv preprint
arXiv:1412.6980, 2014.
[24] A. Radford, L. Metz, and S. Chintala, “Unsupervised representation learning with deep convolutional generative adversarial networks,” arXiv preprint arXiv:1511.06434, 2015.

12

A

Geometric interpretation of GAN and its variants

This appendix provides geometric interpretation of GAN and its variants. In particular, we consider
a specific form of the discriminator Dw,ζ (xi ) given by
Dw,ζ (xi ) = Sf (Vw,ζ (xi )),

where Vw,ζ (xi ) := hw, Φζ (xi )i

(26)

where Sf is an output activation function, Vw,ζ (xi ) is the output layer composed of linear layer
w and the convolutional neural network below corresponding to Φζ (x). Under this choice of the
discriminator, we will show that the differences between existing approaches come from the choice
of the separating hyperplanes and geometric scaling factors.
A.1

GAN

Recall that the empirical estimate of the GAN cost in (2) is given by :
n

L̂GAN (w, ζ, θ) :=

n

1X
1X
log Dw,ζ (xi ) +
log(1 − Dw,ζ (gθ (zi )),
n i=1
n i=1

We now define geometric scaling factors for true and synthetic (or fake) feature vectors:
ti :=

S 0 (hw, Φζ (xi )i)
,
D(xi )

si :=

S 0 (hw, Φζ (gθ (zi ))i)
1 − D(gθ (zi ))

In particular, if the activation function is the sigmoid, i.e. S(u) = 1/(1 + e−u ), then we can easily
see that
ti := 1 − D(xi ),
si := D(gθ (zi )).
Then, the separating hyperplane update is given by:
wGAN

← wGAN + η

n
X

(ti Φζ (xi ) − si Φζ (gθ (xi ))) /n

(27)

i=1

Using another application of chain rules,
X
ζ ← ζ +η
hwGAN , ti ∇ζ Φζ (xi ) − si ∇ζ Φζ (gθ (zi ))i/n

(28)

i∈I

θ

← θ+η

X
hwGAN , si ∇θ Φζ (gθ (zi ))i/n

(29)

i∈I

Aside from different choice of separating hyperplane, the only difference is that the features vectors
needs to be scaled appropriated using geometric scaling parameters. In fact, the scale parameter
is directly related to the geometry of the underlying curved feature spaces due to the log(·) and
nonlinear activations.
From (29), we can easily see that as discriminator becomes accurate, we have si = D (gθ (zi )) ' 0,
so the update of the generator becomes more difficult. This is the main technical limitation of the
GAN training.
A.2

f -GAN

The f -GAN formulation is given by the minmax game of the following empirical cost:
n

F (w, ζ, θ)

=

n

1X ∗
1X
Sf (Vw,ζ (xi )) −
f (Sf (Vw,ζ (gθ (zi )))
n i=1
n i=1

where f ∗ is the convex conjugate of the divergence function f . We again define a geometric scale
factors for true and fake feature vectors:
ti := Sf0 (hw, Φζ (xi )i) ,

si := (f ∗ )0 (Sf (hw, Φζ (xi )i)) Sf0 (hw, Φζ (xi )i) .

The explicit forms of the geometric scaling factors for different f -divergences are‘ summarized in
Table 1,
13

Name
Sf (v)
f ∗ (v)
ti (u)
si (u)
1
1
1
Total variation
tanh(v)
v
coth(u)
coth(u)
2
2
2
Kullback-Leiber (KL)
v
exp(v − 1)
1
exp(u − 1)
Reverse KL
− exp(v)
−1 − log(−v)
− exp(u)
1
Pearson χ2
v
v 2 /4 + v
1
u/2 + 1
e−u
1
Jensen-Shannon
log(2) − log(1 + exp(−v)) − log(2 − exp(v))
1+e−u
1+e−u
e−u
1
GAN
− log(1 + exp(−v))
− log(1 − exp(v))
1+e−u
1+e−u
Table 1: Recommended final layer activation functions for f -GAN [3] and their geometric scaling
factors.

Then, the separating hyperplane update is given by:
wf GAN

← wf GAN + η

n
X

(ti Φζ (xi ) − si Φζ (gθ (xi ))) /n

(30)

From the chain rules, we have discriminator and generator update rules:
X
ζ ← ζ +η
hwf GAN , ti ∇ζ Φζ (xi ) − si ∇ζ Φζ (gθ (zi ))i/n

(31)

i=1

i∈I

← θ+η

θ

X

hwf GAN , si ∇θ Φζ (gθ (zi ))i/n

(32)

i∈I

Note that f -GAN is only different from each other in their construction of the weight coefficient
(ti , si ) (see Table 1) that reflects the underlying geometry of the curved feature space. Other than the
total variation-based divergence, the scaling factors are asymmetric. Thus, controlling the balance
between discriminator and generator updates are one of the important technical issues of f -GAN
training.
A.3

Wasserstein GAN

Wasserstein GAN [7] minimizes the following IPM:
n

n

1X
1X
f (xi ) −
f (gθ (zi ))
n i=1
kf kL ≤1 n i=1

dF (P, Q) = sup

where kf kL := sup{|f (x) − f (y)|/ρ(x, y) : x 6= y ∈ M } is called the Lipschitz seminorm of a
real-valued function f on M . Using the discriminator model (26), the Wasserstein GAN update can
be written by:
Pn
Pn
w, n1 i=1 Φζ (xi ) − n1 i=1 Φζ (gθ (zi ))
min max
(33)
θ

kwk∞ ≤1,ζ

Therefore, other than the mean difference on l∞ ball for the hyperplane normal vector w update, the
W-GAN update is same as the mean matching GAN update with geometric scaling factor ti = si =
1, ∀i.
A.4

Energy-based GAN

For a given a positive margin m, the energy-based GAN (EBGAN) is given by the alternating minimization of the discriminator and generator cost functions [13]:
n

LD (w, ζ)

=


1 X
Dw,ζ (xi ) + [m − Dw,ζ (gθ (zi ))]+
n i=1

=

1X
Dw,ζ (gθ (zi ))
n i=1

n

LG (θ)

14

where [x]+ = max{0, x}. For a function ψ(y) = ay + b[m − y]+ with y, a, b ≥ 0, its subgradient
is given by:

y ∈ [0, m]
a − b,
ψ 0 (y) = a,
y ∈ (m, ∞)

[a − b, a], y ∈ m
Due to the margin, geometric scale factors for true and fake feature vectors should be defined accordingly. More specifically, we have

s , Dw,ζ (gθ (zi )) ∈ [0, m]
0
G
0
ti := Sf (hw, Φζ (xi )i) ,
si := Sf (hw, Φζ (xi )i) , si := i
0, otherwise
Then, the separating hyperplane update is given by:
wf GAN

← wf GAN + η

n
X

(ti Φζ (xi ) − si Φζ (gθ (xi ))) /n

(34)

hwf GAN , ti ∇ζ Φζ (xi ) − si ∇ζ Φζ (gθ (zi ))i/n

(35)

hwf GAN , sG
i ∇θ Φζ (gθ (zi ))i/n

(36)

i=1

Similarly, we have
←

ζ

ζ +η

X
i∈I

← θ+η

θ

X
i∈I

It is worthy to note that the introduction of margin appears similar to our geometric GAN with SVM
hyperplane. In particular, when a linear activation function is used, we have ti = sG
i = 1, the update
equations (35) and (36) appears very similar to (20) and (23), respectively. However, there exists
fundamental differences. First, in EB-GAN, only the fake samples outside the margins are excluded
for the hyperplane and discriminator updates. On the other hand, in geometric GAN, both the true
and fake samples outside the margins are excluded for the hyperplane and discriminator updates.
The symmetric exclusion in geometric GAN is observed to make the algorithm more robust to
outliers. Second, in EBGAN, the margin is defined for the discriminator values. On the other hand,
in geometric GAN, the margin is determined by the geometric distance between the feature vectors.
Therefore, it is much easier to rely on geometric intuition in designing the geometric GAN.
A.5

Empirical risk minimization

The empirical risk minimization (ERM) with l2 cost is one of the standard method for regression
problems. Although the empirical risk minimization (ERM) is rarely used for generator model, our
analysis also provides the geometric intuition of ERM update.
Specifically, for a given mini-batch training data set S = {(z1 , x1 ), · · · , (zn , xn )}, recall that the
empirical risk minimization (ERM) in the feature space [14] is given by
n

min
θ

1X
2
kΦζ (xi ) − Φζ (gθ (zi ))k
2 i=1

Then, the stochastic gradient for (37) can be represented in the identical form to (7):
X
ζ ← ζ +η
hwiERM , ∇ζ Φζ (xi ) − ∇ζ Φζ (gθ (zi ))i/n
θ

← θ+η

i∈I
n
X

(37)

(38)

wiERM , ∇θ Φζ (gθ (zi )) /n

(39)

Φζ (xi ) − Φζ (gθ (zi ))

(40)

i=1

where wiERM is now defined as
wiERM

=

which is dependent on the sample index i. Therefore, aside from the geometric scaling factors, the
main difference comes from the separating hyperplane for linear classifiers. More specifically, the
hyperplane for geometric GAN is obtained for samples within each mini-batch, while the classifier
for regression is optimally designed for each pair of samples.
15

B

Proof for Theorem 3.1

The proof technique is inspired from that of EB-GAN [13]. We first need the following two lemmas
as the extensions of Lemma 1 in [13].
Lemma B.1. Let ϕ(y) = (m − y) + [m + y]+ . The minimum of ϕ(y) is 2m and is reached at all
y ≥ −m.
Proof. If y ≥ −m, then ϕ(y) = 2m. If y ≤ −m, the ϕ(y) = −2y, whose minimum 2m is achieved
at y = −m.
Lemma B.2. For given α, β ≥ 0, The minimum of ϕ(y) = α [m − y]+ + β [m + y]+ exist if
y ∈ [−m, m]. More specifically, the minimum of ϕ(y) is 2βm at y = m if α > β, or 2αm at
y = −m if α ≤ β.
Proof. If y ≥ m, then ϕ(y) = β [m + y]+ = β(m + y) ≥ 2βm. Thus, inf y∈[m,∞) ϕ(y) = 2βm
at y = m. Similarly, if y ≤ −m, then inf y∈(−∞,−m] ϕ(y) = 2αm at y = −m. For y ∈ [−m, m],
ϕ(y) = α(m−y)+β(m+y) = (α+β)m+(β −α)y. If α > β, inf y∈[−m,m] ϕ(y) = 2βm at y = m
since ϕ(y) is a decreasing function on y ∈ [−m, m]. Similarly, if α ≤ β, inf y∈[−m,m] ϕ(y) = 2αm
at y = −m since ϕ(y) is increasing at y ∈ [−m, m].
Now we are ready for the proof.
Proof of Theorem 3.1. Since R(D, g) and L(D, g) are lower semi-continuous functions, R(D∗ , g ∗ )
has a finite value for the optimal solution (D∗ , g ∗ ). Moreover, due to the alternating minimization,
the pair satisfies:
R(D∗ , g ∗ ) ≤ R(D, g ∗ ), ∀D
(41)
L(D∗ , g ∗ ) ≤ L(D∗ , g), ∀g
(42)
First, define the set
A = {x|px (x) ≤ pg∗ (x)}
and observe that
Z




R(D, g ∗ ) =
px (x) 1 − Dw,b,ζ (x) + + pgθ∗ (x) 1 + Dw,b,ζ (x) + dx
Z


=
1A (x)px (x) 1 − Dw,b,ζ (x) + dx
Z


+ 1A (x)pgθ∗ (x) 1 + Dw,b,ζ (x) + dx
Z


+ 1Ac (x)px (x) 1 − Dw,b,ζ (x) + dx
Z


+ 1Ac (x)pgθ∗ (x) 1 + Dw,b,ζ (x) + dx
where Ac denotes the complementary set and 1A (x) is an indicator function

1, x ∈ A
1A (x) =
.
0, x ∈
/A
From Lemma B.2, we know that 1) when px (x) < pgθ∗ (x), the term within the integral achieves its
minimum value of 2px (x) at D∗ (x) = −1, or 2) when px (x) > pgθ∗ (x), the term within the integral
achieves its minimum value of 2pgθ∗ (x) at D∗ (x) = 1. Therefore,
Z
Z
∗ ∗
R(D , g ) = 2 1A (x)px (x)dx + 2 1Ac (x)pgθ∗ (x)dx
Z



= 2
1A (x)px (x) + 1 − 1A (x) pgθ∗ (x) dx
Z

= 2 + 2 1A (x) px (x) − pgθ∗ (x) dx
≤ 2

(43)
16

where the last inequality comes from px (x) − pgθ (x) ≤ 0 for all x ∈ A.
Second, we will show that R(D∗ , g ∗ ) ≥ 2. Because (42) holds for arbitrary pdf p(x), we have
Z
Z
pgθ∗ (x) (−D∗ (x)) dx ≤
px (x) (−D∗ (x)) dx
By adding 1 to both sides, we have
Z
pgθ∗ (x) (1 − D∗ (x)) dx

Z

px (x) (1 − D∗ (x)) dx

Z

px (x) [1 − D∗ (x)]+ dx

≤
≤

where
the
comes from x ≤ [x]+ = max{0, x}.
Now, by adding
 last inequality

R
pgθ∗ (x) 1 + D∗ (x) + dx on both sides, we have:



R
R
pgθ∗ (x) 1 − D∗ (x) dx + pgθ∗ (x) 1 + D∗ (x) + dx




R
R
≤ px (x) 1 − D∗ (x) + dx + pgθ∗ (x) 1 + D∗ (x) + dx = R(D∗ , g ∗ )
(44)
 

From Lemma B.1, we know that 1 − D∗ (x) + 1 + D∗ (x) + ≥ 2. Thus, we have
Z
Z



R(D∗ , g ∗ ) ≥
pgθ∗ (x) 1 − D∗ (x) dx + pgθ∗ (x) 1 + D∗ (x) + dx
Z
≥
pgθ∗ (x)2dx ≥ 2
Thus,
2 ≤ R(D∗ , g ∗ ) ≤ 2

(45)

Finally, the equality in (43) holds if and only if
Z

1A (x) px (x) − pgθ∗ (x) dx = 0
The above equalities hold if and only if pg∗ (x) = pdata (x) almost everywhere [13]. This concludes
the proof.

17

