G RADIENT D ESCENT P ROVABLY O PTIMIZES O VER - PARAMETERIZED N EURAL N ETWORKS
Anonymous authors
Paper under double-blind review
A BSTRACT
One of the mysteries in the success of neural networks is randomly initialized ﬁrst order methods like gradient descent can achieve zero training loss even though the objective function is non-convex and non-smooth.
This paper demystiﬁes this surprising phenomenon for two-layer fully connected ReLU activated neural networks.
For an m hidden node shallow neural network with ReLU activation and n training data, we show as long as m is large enough and no two inputs are parallel, randomly initialized gradient descent converges to a globally optimal solution at a linear convergence rate for the quadratic loss function.
Our analysis relies on the following observation: over-parameterization and random initialization jointly restrict every weight vector to be close to its initialization for all iterations, which allows us to exploit a strong convexity-like property to show that gradient descent converges at a global linear rate to the global optimum.
We believe these insights are also useful in analyzing deep models and other ﬁrst order methods.
I NTRODUCTION
Neural networks trained by ﬁrst order methods have achieved a remarkable impact on many applications, but their theoretical properties are still mysteries.
One of the empirical observation is even though the optimization objective function is non-convex and non-smooth, randomly initialized ﬁrst order methods like stochastic gradient descent can still ﬁnd a global minimum.
Surprisingly, this property is not correlated with labels.
In Zhang et al.
(2016), authors replaced the true labels with randomly generated labels, but still found randomly initialized ﬁrst order methods can always achieve zero training loss.
A widely believed explanation on why a neural network can ﬁt all training labels is that the neural network is over-parameterized.
For example, Wide ResNet (Zagoruyko and Komodakis) uses 100x parameters than the number of training data. Thus there must exist one such neural network of this architecture that can ﬁt all training data.
However, the existence does not imply why the network found by a randomly initialized ﬁrst order method can ﬁt all the data.
The objective function is neither smooth nor convex, which makes traditional analysis technique from convex optimization not useful in this setting.
To our knowledge, only the convergence to a stationary point is known (Davis et al., 2018).
In this paper we demystify this surprising phenomenon on two-layer neural networks with rectiﬁed linear unit (ReLU) activation.
Formally, we consider a neural network of the following form.
[Equation]
where x ∈ R
d is the input, w r ∈ R d is the weight vector of the ﬁrst layer, a r ∈ R is the output weight and σ ( · ) is the ReLU activation function: σ ( z )
=
z if z ≥ 0 and σ ( z ) =
if z < 0 .
We focus on the empirical risk minimization problem with a quadratic loss.
Given a training data set { ( x i , y i ) } ni =1 , we want to minimize
[Equation]
Our main focus of this paper is to analyze the following procedure. We ﬁx the second layer and apply gradient descent (GD) to optimize the ﬁrst layer
[Equation]
where η > 0 is the step size.
Here the gradient formula for each weight vector is
[Equation]
Though this is only a shallow fully connected neural network, the objective function is still nonsmooth and non-convex due to the use of ReLU activation function.
Even for this simple function, why randomly initialized ﬁrst order method can achieve zero training error is not known.
Many previous works have tried to answer this question or similar ones.
Attempts include landscape analysis (Soudry and Carmon, 2016), partial differential equations (Mei et al.)
, analysis of the dynamics of the algorithm (Li and Yuan, 2017), optimal transport theory (Chizat and Bach, 2018), to name a few.
These results often make strong assumptions on the labels and input distributions or do not imply why randomly initialized ﬁrst order method can achieve zero training loss.
See Section 2 for detailed comparisons between our result and previous ones.
In this paper, we rigorously prove that as long as no two inputs are parallel and m is large enough, with randomly initialized a and W (0) , gradient descent achieves zero training loss at a linear convergence rate, i.e., it ﬁnds a solution W ( K ) with L ( W ( K )) ≤ (cid:15) in K
= O (log (1 /(cid:15) )) iterations.
Thus, our theoretical result not only shows the global convergence but also gives a quantitative convergence rate in terms of the desired accuracy.
Analysis Technique Overview Our proof relies on the following insights.
First we directly analyze the dynamics of each individual prediction f
( W , a , x i ) for i = 1 , . . . , n .
This is different from many previous work (Du et al., 2017b; Li and Yuan, 2017) which tried to analyze the dynamics of the parameter ( W ) we are optimizing.
Note because the objective function is non-smooth and non-convex, analysis of the parameter space dynamics is very difﬁcult.
In contrast, we ﬁnd the dynamics of prediction space is governed by the spectral property of a Gram matrix (which can vary in each iteration, c.f. Equation (6)) and as long as this Gram matrix’s least eigenvalue is lower bounded, gradient descent enjoys a linear rate. It is easy to show as long as no two inputs are parallel, in the initialization phase, this Gram matrix has a lower bounded least eigenvalue.
(c.f. Theorem 3.1).
Thus the problem reduces to showing the Gram matrix at later iterations is close to that in the initialization phase.
Our second observation is this Gram matrix is only related to the activation patterns ( I (cid:8) w (cid:62) r x i ≥ 0 (cid:9) ) and we can use matrix perturbation analysis to show if most of the patterns do not change, then this Gram matrix is close to its initialization.
Our third observation is we ﬁnd over-parameterization, random initialization, and the linear convergence jointly restrict every weight vector w r to be close to its initialization.
Then we can use this property to show most of the patterns do not change.
Combining these insights we prove the ﬁrst global quantitative convergence result of gradient descent on ReLU activated neural networks for the empirical risk minimization problem.
Notably, our proof only uses linear algebra and standard probability bounds so we believe it can be easily generalized to analyze deep neural networks.
Notations We let [ n ] = { 1 , 2 , . . . , n } .
Given a set S , we use unif { S } to denote the uniform distribution over S .
Given an event E , we use I { A } to be the indicator on whether this event happens.
We use N ( 0 , I ) to denote the standard Gaussian distribution.
For a matrix A , we use A ij to denote its ( i, j ) -th entry.
We use (cid:107)·(cid:107) 2 to denote the Euclidean norm of a vector, and use (cid:107)·(cid:107) F to 1 In Section 3.2, we also extend our technique to analyze the setting where we train both layers jointly
We remark that if one ﬁxes the ﬁrst layer and only optimizes the output layer, then the problem becomes a convex and smooth one.
If m is large enough, one can show the global minimum has zero training loss (Nguyen and Hein, 2018).
Though for both cases (ﬁxing the ﬁrst layer and ﬁxing the output layer), gradient descent achieves zero training loss, the learned prediction functions are different.
denote the Frobenius norm of a matrix. If a matrix A is positive semi-deﬁnite, we use λ min ( A ) to denote its smallest eigenvalue.
We use (cid:104)· , ·(cid:105) to denote the standard Euclidean inner product between two vectors.
C OMPARISON WITH P REVIOUS R ESULTS
In this section, we survey an incomplete list of previous attempts in analyzing why ﬁrst order methods can ﬁnd a global minimum.
Landscape Analysis A popular way to analyze non-convex optimization problems is to identify whether the optimization landscape has some good geometric properties.
Recently, researchers found if the objective function is smooth and satisﬁes (1) all local minima are global and (2) for every saddle point, there exists a negative curvature, then the noise-injected (stochastic) gradient descent (Jin et al., 2017;
Ge et al., 2015; Du et al., 2017a) can ﬁnd a global minimum in polynomial time.
This algorithmic ﬁnding encouraged researchers to study whether the deep neural networks also admit these properties.
For the objective function deﬁned in Equation (2), some partial results were obtained.
Soudry and Carmon (2016) showed if md ≥ n , then at every differentiable local minimum, the training error is zero.
However, since the objective is non-smooth, it is hard to show gradient descent convergences to a differentiable local minimum.
Xie et al.
(2017) studied the same problem and related the loss to the gradient norm through the least singular value of the “extended feature matrix” D at the stationary points.
However, they did not prove the convergence rate of the gradient norm.
Interestingly, our analysis relies on the Gram matrix which is DD (cid:62) .
Landscape analyses of ReLU activated neural networks for other settings have also been studied in many previous works (Ge et al., 2017; Safran and Shamir, 2016; Zhou and Liang, 2017; Freeman and Bruna, 2016; Hardt and Ma, 2016; Nguyen and Hein, 2018).
These works establish favorable landscape properties but none of them implies that gradient descent converges to a global minimizer of the empirical risk. More recently, some negative results have also been discovered (Safran and
Shamir, 2018; Yun et al., 2018a) and new procedures have been proposed to test local optimality and escape strict saddle points at non-differentiable points (Yun et al., 2018b).
However, the new procedures cannot ﬁnd global minima as well.
For other activation functions, some previous works showed the landscape does have the desired geometric properties (Du and Lee, 2018; Soltanolkotabi et al., 2018; Nguyen and Hein, 2017; Kawaguchi, 2016; Haeffele and Vidal, 2015; Andoni et al., 2014; Venturi et al., 2018; Yun et al., 2018a).
However, it is unclear how to extend their analyses to our setting.
Analysis
of Algorithm Dynamics
Another way to prove convergence result is to analyze the dynamics of ﬁrst order methods directly.
Our paper also belongs to this category.
Many previous works assumed (1) the input distribution is Gaussian and (2) the label is generated according to a planted neural network.
Based on these two (unrealistic) conditions, it can be shown that randomly initialized (stochastic) gradient descent can learn a ReLU (Tian, 2017; Soltanolkotabi, 2017), a single convolutional ﬁlter (Brutzkus and Globerson, 2017), a convolutional neural network with one ﬁlter and one output layer (Du et al., 2018b) and residual network with small spectral norm weight matrix (Li and Yuan, 2017).
Beyond Gaussian input distribution, Du et al. (2017b) showed for learning a convolutional ﬁlter, the Gaussian input distribution assumption can be relaxed but they still required the label is generated from an underlying true ﬁlter.
Comparing with these work, our paper does not try to recover the underlying true neural network.
Instead, we focus on providing theoretical justiﬁcation on why randomly initialized gradient descent can achieve zero training loss, which is what we can observe and verify in practice.
The most related paper is by Li and Liang (2018) who observed that when training a two-layer full connected neural network, most of the patterns ( I (cid:8) w (cid:62) r x i ≥ 0 (cid:9) ) do not change over iterations, which we also use to show the stability of the Gram matrix.
They used this observation to obtain the convergence rate of GD on a two-layer over-parameterized neural network for the cross-entropy loss.
They need the number of hidden nodes m scales with poly(1 /(cid:15) ) where (cid:15) is the desired accuracy.
Thus unless the number of hidden nodes m → ∞ , their result does not imply GD can achieve zero training loss.
We improve by allowing the amount of over-parameterization to be independent of the desired accuracy and show GD can achieve zero training loss.
Furthermore, our proof is much simpler and more transparent so we believe it can be easily generalized to analyze other neural network architectures.
Other Analysis Approaches Chizat and Bach (2018) used optimal transport theory to analyze continuous time gradient descent on over-parameterized models.
They required the second layer to be inﬁnitely wide and their results on ReLU activated neural network is only at the formal level.
Mei et al. analyzed SGD for optimizing the population loss and showed the dynamics can be captured by a partial differential equation in the suitable scaling limit.
They listed some speciﬁc examples on input distributions including mixture of Gaussians.
However, it is still unclear whether this framework can explain why ﬁrst order methods can minimize the empirical risk.
Daniely (2017) built connection between neural networks with kernel methods and showed stochastic gradient descent can learn a function that is competitive with the best function in the conjugate kernel space of the network.
Again this work does not imply why ﬁrst order methods can achieve zero training loss.
C ONTINUOUS T IME A NALYSIS
In this section, we present our result for gradient ﬂow, i.e., gradient descent with inﬁnitesimal step size.
The analysis of gradient ﬂow is a stepping stone towards understanding discrete algorithms and this is the main topic of recent work (Arora et al., 2018; Du et al., 2018a).
In the next section, we will modify the proof and give a quantitative bound for gradient descent with positive step size.
Formally, we consider the ordinary differential equation 6 deﬁned by:
[Equation]
for r ∈ [ m ] .
We denote u
i ( t ) =
f ( W ( t ) , a , x i ) the prediction on input x i at time t
and we let u ( t ) =
( u 1 ( t ) , . . . ,
u n ( t ))
∈ R n be the prediction vector at time t .
We state our main assumption. Assumption 3.1. Deﬁne matrix H ∞ ∈
R n × n with
H ∞ ij =
E w ∼ N ( 0 , I ) (cid:2) x (cid:62) i x j I (cid:8) w (cid:62) x i ≥ 0 , w (cid:62) x j ≥ 0 (cid:9)(cid:3) .
We assume λ 0 (cid:44) λ min ( H ∞ ) > 0 .
H ∞ is the Gram matrix induced by the ReLU activation function and the random initialization.
Later we will show that during the training, though the Gram matrix may change (c.f. Equation (6)), it is still close to H ∞ .
Furthermore, as will be apparent in the proof (c.f. Equation (7)), H ∞ is the fundamental quantity that determines the convergence rate. Interestingly, various properties of this H ∞ matrix has been studied in previous works (Xie et al., 2017; Tsuchida et al., 2017). Now to justify this assumption, the following theorem shows if no two inputs are parallel the least eigenvalue is strictly positive. Theorem 3.1.
If for any i (cid:54) =
j , x i (cid:54)(cid:107) x j , then λ 0 > 0 .
Note for most real world datasets, no two inputs are parallel, so our assumption holds in general. Now we are ready to state our main theorem in this section.
Theorem 3.2 (Convergence Rate of Gradient Flow) . Suppose Assumption 3.1 holds and for all i ∈ [ n ] , (cid:107) x i (cid:107) 2 = 1 and | y i | ≤ C for some constant C .
Then if we set the number of hidden nodes m = Ω (cid:16) n 6 λ 40 δ 3 (cid:17) and we i.i.d. initialize w r ∼ N ( 0 , I ) , a r ∼ unif [ {− 1 , 1 } ] for r ∈ [ m ] , then with probability at least 1 − δ over the initialization, we have
[Equation]
This theorem establishes that if m is large enough, the training error converges to 0 at a linear rate.
Here we assume (cid:107) x
i (cid:107) 2 = 1 only for simplicity and it is not hard to relax this condition.
The the cleanest our technique.
bounded label condition also holds for most real world data set.
The number of hidden nodes m required is Ω (cid:16) n 6 λ 40 δ 3 (cid:17) , which depends on the number of samples n , λ 0 , and the failure probability δ .
Over-parameterization, i.e., the fact m = poly( n, 1 /λ 0 , 1 /δ ) , plays a crucial role in guaranteeing gradient descent to ﬁnd the global minimum. In this paper, we only use the simplest concentration inequalities (Hoeffding’s and Markov’s) in order to have the cleanest proof.
We believe using a more advanced concentration analysis we can further improve the dependency.
Lastly, we note the speciﬁc convergence rate depends on λ 0 but independent of the number of hidden nodes m .
3. P ROOF OF T HEOREM 3.
Our ﬁrst step is to calculate the dynamics of each prediction.
[Equation]
where H ( t ) is an n × n matrix with ( i, j ) -th entry
[Equation] [Equation]
With this H ( t ) matrix, we can write the dynamics of predictions in a compact way:
[Equation]
Remark 3.1.
Note Equation (7) completely describes the dynamics of the predictions.
In the rest of this section, we will show (1) at initialization (cid:107) H (0) − H ∞ (cid:107) 2 is O ( (cid:112) 1 /m ) and (2) for all t > 0 , (cid:107) H ( t ) − H (0) (cid:107) 2 is O ( (cid:112) 1 /m ) .
Therefore, according to Equation (7) , as m → ∞ , the dynamics of the predictions are characterized by H ∞ .
This is the main reason we believe H ∞ is the fundamental quantity that describes this optimization process.
H ( t ) is a time-dependent symmetric matrix.
We ﬁrst analyze its property when t = 0 .
The following lemma shows if m is large then H (0) has a lower bounded least eigenvalue with high probability.
The proof is by the standard concentration bound so we defer it to the appendix.
Lemma 3.1.
If m = Ω (cid:16) n 2 λ 20 log (cid:0) nδ (cid:1)(cid:17) , we have with probability at least 1 − δ , (cid:107) H (0) − H ∞ (cid:107) 2 ≤ λ 0and λ min ( H (0)) ≥ 34 λ 0 .
Our second step is to show H ( t ) is stable in terms of W ( t ) .
Formally, the following lemma shows for any W close to W (0) , the induced Gram matrix H is close to H (0) and has a lower bounded least eigenvalue.
Lemma 3.2.
If w 1 , . . . , w m are i.i.d. generated from N ( 0 , I ) , then with probability at least 1 − δ , the following holds. For any set of weight vectors w 1 , . . . , w m ∈ R d that satisfy for any r ∈ [ m ] , (cid:107) w r (0) − w r (cid:107) 2 ≤ cδλ 0 n 2 (cid:44) R for some small positive constant c , then the matrix H ∈ R n × n deﬁned by
[Equation] satisﬁes (cid:107) H − H (0) (cid:107) 2 < λ 0 4 and λ min ( H ) > λ 0 2 .
This lemma plays a crucial role in our analysis so we give the proof below.
Proof of Lemma 3.2 We deﬁne the event
[Equation]
Note this event happens if and only if (cid:12)(cid:12) w r (0) (cid:62) x i (cid:12)(cid:12) < R .
Recall w r (0) ∼ N ( 0 , I ) .
By anticoncentration inequality of Gaussian, we have P ( A ir ) =
P z ∼ N (0 , 1) ( | z | < R ) ≤ 2 R √ 2 π .
Therefore, for any set of weight vectors w 1 , . . . , w m that satisfy the assumption in the lemma, we can bound the entry-wise deviation on their induced matrix H : for any ( i, j ) ∈ [ n ] × [ n ] [Equation] [Table]
[Equation]
where the expectation is taken over the random initialization of w 1 (0) , . . , w m (0) . Summing over ( i, j ) , we have E (cid:104)(cid:80) ( n,n ) ( i,j )=(1 , 1) |
H ij − H ij (0) | (cid:105) ≤ 4 n 2 R √ 2 π . Thus by Markov’s inequality, with probability 1 − δ , we have (cid:80) ( [Equation] [Equation]
Next, we use matrix perturbation theory to bound the deviation from the initialization
[Equation]
Lastly, we lower bound the smallest eigenvalue by plugging in R
[Equation]
The next lemma shows two facts if the least eigenvalue of H ( t ) is lower bounded.
First, the loss converges to 0 at a linear convergence rate.
Second, w r ( t ) is close to the initialization for every r ∈ [ m ] .
This lemma clearly demonstrates the power of over-parameterization.
Lemma 3.3.
Suppose for 0 ≤ s ≤ t , λ min ( H ( s )) ≥ λ 0 2 .
Then we have (cid:107) y − u ( t ) (cid:107) 22 ≤ exp( − λ 0 t ) (cid:107) y − u (0) (cid:107) 22 and for any r ∈ [ m ] , (cid:107) w r ( t ) − w r (0) (cid:107) [Equation] y [Equation]
[Equation]
(cid:107) 2 √ mλ [Equation] R (cid:48) . Proof of Lemma 3.3 Recall we can write the dynamics of predictions as ddt u ( t ) =
H ( y − u ( t )) .
We can calculate the loss function dynamics
[Equation]
Thus we have ddt (cid:16) exp( λ 0 t ) (cid:107) y − u ( t ) (cid:107) 2 2 (cid:17) ≤ 0 and exp( λ 0 t ) (cid:107) y − u ( t ) (cid:107) 2 2 is a decreasing function with respect to t .
Using this fact we can bound the loss
[Equation]
Therefore, u ( t ) →
y exponentially fast.
Now we bound the gradient norm.
Recall for 0 ≤ s ≤ t , [Equation]
d
[Equation]
Integrating the gradient, we can bound the distance from the initialization
[Equation]
The next lemma shows if R (cid:48) <
R , the conditions in Lemma 3.2 and 3.3 hold for all t ≥ 0 .
The proof is by contradiction and we defer it to appendix.
Lemma 3.4.
If R (cid:48) < R , we have for all t ≥ 0 , λ min ( H ( t )) ≥ 12 λ 0 , for all r ∈ [ m ] , (cid:107) w r ( t ) − w r (0) (cid:107) 2 ≤ R (cid:48) and (cid:107) y − u ( t ) (cid:107) 22 ≤ exp( − λ 0 t ) (cid:107) y − u (0) (cid:107) 22 .
Thus it is sufﬁcient to show R (cid:48) <
R which is equivalent to m = Ω (cid:16) n 5 (cid:107) y − u (0) (cid:107) 22 λ 40 δ(cid:17) .
We bound
E (cid:104) (cid:107) y
[Equation]
u (0) (cid:107) 22 (cid:105) = n (cid:88)
i = ( y
[Equation] [Equation] [Equation] [Equation]
(cid:2)
[Equation]
[Equation]
(cid:3) )
=
n (cid:88) i =
( y 2 i + 1) = O ( n ) . Thus by Markov’s inequality, we have with probability at least 1 − δ , (cid:107) y − u (0) (cid:107)
= O ( nδ ) .
Plugging in this bound we prove the theorem. 3. J OINTLY T RAINING B OTH L AYERS
In this subsection, we showcase our proof technique can be applied to analyze the convergence of gradient ﬂow for jointly training both layers.
Formally, we consider the ordinary differential equation deﬁned by:
[Equation] d w r ( t ) [Equation] for r = 1 , . . .
, m .
The following theorem shows using gradient ﬂow to jointly train both layers, we can still enjoy linear convergence rate towards zero loss.
Theorem 3.3 (Convergence Rate of Gradient Flow for Training Both Layers) .
Under the same assumptions as in Theorem 3.2, if we set the number of hidden nodes m = Ω (cid:16) n 6 log( m/δ ) λ 40 δ(cid:17) and we i.i.d.
initialize w r ∼ N ( 0 , I ) , a r ∼ unif [ {− 1 , 1 } ] for r ∈ [ m ] , with probability at least 1 − δ over the initialization we have
[Equation]
Theorem 3.3 shows under the same assumptions as in Theorem 3.2, we can achieve the same convergence rate as that of only training the ﬁrst layer.
The proof of Theorem 3.3 relies on the same arguments as the proof of Theorem 3.2.
Again we consider the dynamics of the predictions and this dynamics is characterized by a Gram matrix.
We can show for all t > 0 , this Gram matrix is close to the Gram matrix at the initialization phase.
We refer readers to appendix for the full proof.  D ISCRETE T IME A NALYSIS In this section, we show randomly initialized gradient descent with a constant positive step size converges to the global minimum at a linear rate.
We ﬁrst present our main theorem.
Theorem 4.1 (Convergence Rate of Gradient Descent) .
Under the {+same assumptions as in Theorem 3.2, if we set the number of hidden nodes m = Ω (cid:16) n 6 λ 
+}
δ 3 (cid:17) , we i.i.d.
initialize w r ∼ N ( 0 , I ) , a r ∼ unif [ {− 1 , 1 } ] for r ∈ [ m ] , and we set the step size η = O (cid:0) λ 0 n 2 (cid:1) then with probability at least 1 − δ over the random initialization we have for k = 0 , 1 , 2 , . . .
[Equation]
(cid:107) [Table] Theorem 4.1 shows even though the objective function is non-smooth and non-convex, gradient descent with a constant step size still enjoys a linear convergence rate.
Our assumptions on the least eigenvalue and the number of hidden nodes are exactly the same as the theorem for gradient ﬂow.
4. P ROOF OF T HEOREM 4.
We prove Theorem 4.1 by induction.
Our induction hypothesis is just the following convergence rate of the empirical loss.
Condition 4.1.
At the k -th iteration, we have (cid:107) y − u ( k ) (cid:107)
[Equation] [Equation]
u
[Equation]
[Equation]
A directly corollary of this condition is the following bound of deviation from the initialization.
The proof is similar to that of Lemma 3.3 so we defer it to appendix.
Corollary 4.1.
If Condition 4.1 holds for k (cid:48)
= 0 , . .
, k , then we have for every r ∈ [ m ]
[Equation]
Now we show Condition 4.1 holds for every k = 0 , 1 , . . . .
For the base case k = 0 , by deﬁnition Condition 4.1 holds.
Suppose for k (cid:48)
= 0 , . .
, k , Condition 4.1 holds and we want to show
Condition 4.1 holds for k (cid:48)
=
k + 1 .
Our strategy is similar to the proof of Theorem 3.2.
We deﬁne the event
[Equation]
where R = cλ 0 n 2 for some small positive constant c .
Different from gradient ﬂow, for gradient descent we need a more reﬁned analysis.
We let S i = { r ∈
m ] : I { A ir } = 0 } and S ⊥ i =
[ m ] \ S i .
The following lemma bounds the sum of sizes of S ⊥ i .
The proof is similar to the analysis used in
Lemma 3.2.
See Section A for the whole proof.
Lemma 4.1.
With probability at least 1 − δ over the initialization, we have (cid:80) ni
=1 (cid:12)(cid:12)
S ⊥ i (cid:12)(cid:12) ≤ CmnRδ for some positive constant C > 0 .
Next, we calculate the difference of predictions between two consecutive iterations, analogue to du i ( t ) dt term in Section 3.
[Equation]
Here we divide the right hand side into two parts.
I i 1 accounts for terms that the pattern does not change
and I i 2 accounts for terms that pattern may change.
[Equation]
We view I i 2 as a perturbation and bound its magnitude.
Because ReLU is a 1 -Lipschitz function and | a r | = 1 , we have
[Equation]
To analyze I i 1 , by Corollary 4.1, we know (cid:107) w r ( k ) − w r (0) (cid:107) ≤ R (cid:48) and (cid:107) w r ( k )
− w r (0) (cid:107) ≤ R (cid:48) for all r ∈ [ m ] .
Furthermore, because R (cid:48) <
R , we know I (cid:8) w r ( k + 1) (cid:62) x i ≥ 0 (cid:9) =
I (cid:8) w r ( k ) (cid:62) x i ≥ 0 (cid:9) for r ∈ S i .
Thus we can ﬁnd a more convenient expression of I i 1 for analysis
[Equation]
where H ij ( k ) = 1 m (cid:80) mr =1 x (cid:62) i x j I (cid:8) w r ( k ) (cid:62) x i ≥ 0 , w r ( k ) (cid:62) x j ≥ 0 (cid:9) is just the ( i, j ) -th entry of a discrete version of Gram matrix deﬁned in Section 3 and H ⊥ ij ( k )
[Equation]
1 m (cid:80) r ∈ S [Equation]
i x (cid:62) i x j I (cid:8) w r ( k ) (cid:62) x i ≥ 0 , w r ( k ) (cid:62) x j ≥ 0 (cid:9) is a perturbation matrix.
Let H ⊥ ( k ) be the n × n matrix with ( i, j ) -th entry being H ⊥ ij ( k ) .
Using Lemma 4.1, we obtain an upper bound of the operator norm
[Equation]
Similar to the classical analysis of gradient descent, we also need bound the quadratic term.
[Equation]
With these estimates at hand, we are ready to prove the induction hypothesis.
[Equation]
The third equality we used the decomposition of u ( k + 1) − u ( k ) .
The ﬁrst inequality we used the
Lemma 3.2, the bound on the step size, the bound on I 2 , the bound on (cid:13)(cid:13) H ( k ) ⊥
(cid:13)(cid:13) 2 and the bound on (cid:107) u ( k + 1)
− u ( k ) (cid:107) 22 .
The last inequality we used the bound of the step size and the bound of R .
Therefore Condition 4.1 holds for k (cid:48)
=
k + 1 .
Now by induction, we prove Theorem 4.1.  E XPERIMENTS
In this section, we use synthetic data to corroborate our theoretical ﬁndings.
We use the initialization and training procedure described in Section 1. For all experiments, we run 100 epochs of gradient descent and use a ﬁxed step size. We uniformly generate n = 1000 data points from a d =dimensional unit sphere and generate labels from a one-dimensional standard Gaussian distribution. We test three metrics with different widths ( m ).
First, we test how the amount of overparameterization affects the convergence rates.
Second, we test the relation between the amount of over-parameterization and the number of pattern changes.
Formally, at a given iteration k , we check
[Equation]
(there are mn patterns).
This aims to verify Lemma 3.2.
Last, we test the relation between the amount of over-parameterization and the maximum of the distances between weight vectors and their initializations.
Formally, at a given iteration k , we check max r
[Equation]
[ m ]
[Equation]
r
[Equation]
k
[Equation]
[Equation]
This aims to verify Lemma 3.3 and Corollary 4.1.
Figure 1a shows as m becomes larger, we have better convergence rate.
We believe the reason is as m becomes larger, H ( t ) matrix becomes more stable, and thus has larger least eigenvalue.
Figure 1b and Figure 1c show as m becomes larger, the percentiles of pattern changes and the maximum distance from the initialization become smaller.
These empirical ﬁndings are consistent with our theoretical results.  C ONCLUSION AND D ISCUSSION In this paper we show with over-parameterization, gradient descent provable converges to the global minimum of the empirical loss at a linear convergence rate.
The key proof idea is to show the [Figure] over-parameterization makes Gram matrix remain positive deﬁnite for all iterations, which in turn guarantees the linear convergence.
Here we list some future directions.
First, we believe our approach can be generalized to deep neural networks. We elaborate the main idea here for gradient ﬂow. Consider a deep neural network of the form [Equation]
where x ∈ R
d is the input, W (1)
∈ R m ×
d is the ﬁrst layer,
W ( h ) ∈
R m × m for h = 2 , . .
, H are the middle layers and a ∈ R m is the output layer.
Recall u
i is the i -th prediction.
If we use the quadratic loss, we can compute
[Equation]
Similar to Equation (5), we can calculate
[Equation]
G
[Equation]
( t ) where G ( h ) ( t ) ∈ R n × n with G ( h ) ij ( t ) = (cid:104)
∂u
i ( t ) W ( h ) ( t ) , ∂u j ( t ) W ( h ) ( t ) (cid:105) .
Therefore, similar to Equation (7), we can write
[Equation]
Note for every h ∈ [ H ] , G ( h ) is a Gram matrix and thus it is positive semideﬁnite. If (cid:80) Hh =1 G ( h ) ( t ) has a lower bounded least eigenvalue for all t , then similar to Section 3, gradient ﬂow converges to zero training loss at a linear convergence rate. Based on our observations in Remark 3.1, we conjecture that if m is large enough, (cid:80) Hh =1 G ( h ) (0) is close to a ﬁxed matrix (cid:80)
Hh =1 G ( h ) ∞ and (cid:80) Hh =1 G ( h ) ( t ) is close its initialization (cid:80)
Hh =1 G ( h ) (0) for all t > 0 .
Therefore, using the same arguments as we used in Section 3, as long as (cid:80) Hh =1 G ( h ) ∞ has a lower bounded least eigenvalue, gradient ﬂow converges to zero training loss at a linear convergence rate. Second, we believe the number of hidden nodes m required can be reduced.
For example, previous work (Soudry and Carmon, 2016) showed m ≥ nd is enough to make all differentiable local minima global. In our setting, using advanced tools from probability and matrix perturbation theory to analyze H ( t ) , we may be able to tighten the bound.
Lastly, in our paper, we used the empirical loss as a potential function to measure the progress.
If we use another potential function, we may be able to prove the convergence rates of accelerated methods.
This technique has been exploited in Wilson et al.
(2016) for analyzing convex optimization.
It would be interesting to bring their idea to analyze other ﬁrst order methods for optimizing neural networks.
A T ECHNICAL P ROOFS FOR S ECTION
Proof of Theorem 3.1.
The proof of this lemma just relies on standard real and functional analysis.
Let H be the Hilbert space of integrable d -dimensional vector ﬁelds on R d : f ∈
H if E w ∼ N
[Equation]
I
[Equation]
< ∞ .
The inner product of this space is then (cid:104) f, g (cid:105) H =
E w ∼ N ( 0 , I ) (cid:2) f ( w ) (cid:62) g ( w ) (cid:3) .
ReLU activation induces an inﬁnite-dimensional feature map φ which is deﬁned as for any x ∈ R d , ( φ ( x ))( w ) =
x I (cid:8) w (cid:62) x ≥ 0 (cid:9) where w can be viewed as the index.
Now to prove H ∞ is strictly positive deﬁnite, it is equivalent to show φ ( x 1 ) , . .
, φ ( x n ) ∈ H are linearly independent.
Suppose that there are α 1 , · · · , α n ∈ R such that
[Equation]
This means that
[Equation]
Now we prove α i = 0 for all i .
We deﬁne D i = (cid:8) w ∈ R d : w (cid:62) x i = 0 (cid:9) .
This is set of discontinuities of φ ( x i ) .
The following lemma characterizes the basic property of these discontinuity sets.
Lemma A.1.
If for any i (cid:54) =
j , x i (cid:54)(cid:107) x j , then for any i ∈ [ m ] , D i (cid:54)⊂ (cid:83) j
(cid:54) =
i D j .
Now for a ﬁxed i ∈ [ n ] , since D i (cid:54)⊂ (cid:83) j (cid:54) =
i D j , we can choose z ∈ D i \ (cid:83)
j (cid:54) = i D j .
Note D j , j (cid:54) =
i are closed sets.
We can pick r 0 > 0 small enough such that B ( z , r ) ∩ D j = ∅ , ∀ j (cid:54) =
i, r ≤ r 0 .
Let
B ( z , r ) =
B + r (cid:116)
B − r where
[Equation]
For j (cid:54) =
{+i , φ ( x j )( w ) is continuous in a neighborhood of z , then for any (cid:15) > 
+}
there is a small enough r > 0 such that
[Equation]
Let µ be the Lebesgue measure on R d .
We have
[Equation]
and similarly
[Equation]
[Equation]
Thus, we have lim
[Equation]
Therefore, as r → 0+ , by continuity, we have
[Equation]
i
[Equation]
[Equation]
Next recall that ( φ ( x ))( w ) =
x I (cid:8) x (cid:62) w > 0 (cid:9) , so for w ∈ B + r and x i , ( φ ( x i ))
( w ) =
x i .
Then, we have
[Equation]
For w ∈ B − r and x i , we know ( φ ( x i ))
( w ) = 0 .
Then we have
[Equation]
Now recall (cid:80) i
α i φ ( x i ) ≡ 0 .
Using Equation (9), (10) and (11), we have
[Equation]
Since x i (cid:54) = 0 , we must have α i = 0 .
We complete the proof.
Proof of Lemma A.1. Let µ be the canonical Lebesgue measure on D i . We have (cid:80) j (cid:54) = i µ ( D i ∩ D j ) {+= 
+}
because D i ∩ D j is a hyperplane in D i . Now we bound [Equation] This implies our desired result.
Proof of Lemma 3.1. For every ﬁxed ( i, j ) pair, H ij (0) is an average of independent random variables. Therefore, by Hoeffding inequality, we have with probability 1 − δ (cid:48) ,
[Equation]
Setting δ (cid:48)
= n 2 δ and applying union bound over ( i, j ) pairs, we have for every ( i, j ) pair with probability at least 1 − δ
[Equation]
Thus we have
[Equation]
Thus if m = Ω (cid:16) n 2 log( n/δ ) λ 2(cid:17) we have the desired result.
Proof of Lemma 3.4.
Suppose the conclusion does not hold at time t .
If there exists r ∈ [ m ] , (cid:107) w r ( t ) − w r (0) (cid:107) ≥ R (cid:48) or (cid:107) y − u ( t ) (cid:107) 22 > exp( − λ 0 t ) (cid:107) y − u (0) (cid:107) 22 , then by Lemma 3.3 we know there exists s ≤ t such that λ min ( H ( s )) < 12 λ 0 . By Lemma 3.2 we know there exists t 0 = inf
[Equation]
(cid:27) . Thus at t 0 , there exists r ∈ [ m ] , (cid:107) w r ( t 0 ) − w r (0) (cid:107)
= R .
Now by Lemma 3.2, we know H ( t 0 ) ≥ 12 λ 0 for t (cid:48) ≤ t 0 .
However, by Lemma 3.3, we know (cid:107) w r
t 0 ) − w r (0) (cid:107) 2 < R (cid:48) <
R .
Contradiction.
For the other case, at time t , λ min ( H ( t )) < 12 λ
we know there exists
[Equation]
The rest of the proof is the same as the previous case.
A. P ROOF OF T HEOREM 3.
In this section we show using gradient ﬂow to jointly train both the ﬁrst layer and the output layer we can still achieve 0 training loss.
We follow the same approach we used in Section 3.
Recall the gradient for a .
[Equation]
We compute the dynamics of an individual prediction. [Equation] Recall we have found a convenient expression for the ﬁrst term. [Equation]
where
[Equation]
For the second term
[Equation]
it easy to derive
[Equation]
where
[Equation]
Therefore we have
[Equation]
First, the following lemma shows as long as H ( t ) has lower bounded least eigenvalue, gradient ﬂow enjoys a linear convergence rate.
The proof is analogue to the ﬁrst part of the proof of Lemma 3.3.
Lemma A.2.
If for≤ s ≤ t , λ min ( H ( s )) ≥ λ 0 2 , we have (cid:107) y − u ( t ) (cid:107)≤ exp( − λ 0 t ) (cid:107) y − u (0) (cid:107) 2 2 .
Proof of Lemma A.2.
We can calculate the loss function dynamics
[Equation]
where in the ﬁrst inequality we use the fact that G ( t ) is Gram matrix thus it is positive.
Thus we have ddt (cid:16) exp( λ 0 t ) (cid:107) y − u ( t ) (cid:107) 22 (cid:17) ≤ 0 and exp( λ 0 t ) (cid:107) y − u ( t ) (cid:107) 22 is a decreasing function with respect to t .
Using this fact we can bound the loss
[Equation]
22 ≤
[Equation]
y
[Equation]
We continue to follow the analysis in Section 3.
For convenience, we deﬁne
[Equation]
The ﬁrst lemma characterizes how the perturbation in a and W affect the Gram matrix.
Lemma A.3.
With probability at least 1 − δ over initialization, if a set of weight vectors { w r } mr = and the output weight a satisfy for all r ∈ [ m ] , (cid:107) w r − w r
0) (cid:107) 2 ≤ R w and | a r − a r (0) | ≤ R a , then the matrix H ∈ R n × n deﬁned by
[Equation]
{+satisﬁes (cid:107) H − H (0) (cid:107) 2 ≤ λ 0 
+}
and λ min ( H ) > λ 0 2 .
Proof of Lemma A.3.
Deﬁne a surrogate Gram matrix,
[Equation]
Using the same analysis for Lemma 3.3, we know (cid:107) H (cid:48)
− H (0) (cid:107) 2 ≤ 4 n 2 R w √ 2 πδ .
Now we bound H − H (cid:48) .
For ﬁxed ( i, j ) ∈ [ n ] × [ n ] , we have
[Equation]
a .
Therefore, we have (cid:107) H − H (cid:48) (cid:107) 2 ≤ (cid:80) ( i,j ) ∈ [ n ] × [ n ] (cid:12)(cid:12) H ij − H (cid:48) ij (cid:12)(cid:12) ≤ 2 n 2 R a .
Combing these two inequalities we have (cid:107) H − H (0) (cid:107) 2 ≤ 4 n 2 R w √ 2 πδ + 2 n 2 R a ≤ λ 0 4 .
Lemma A.4.
Suppose for 0 ≤ s ≤ t , λ min ( H ( s )) ≥ λ 0 2 and | a r ( s ) − a r (0) | ≤ R a .
Then we have (cid:107) w r ( t ) − w r (0) (cid:107) 2 ≤ R (cid:48) w .
Proof of Lemma A.4.
We bound the gradient.
Recall for 0 ≤ s ≤ t ,
[Equation]
Integrating the gradient and using Lemma A.2, we can bound the distance from the initialization
[Equation]
Lemma A.5.
With probability at least 1 − δ over initialization, the following holds.
Suppose for 0 ≤ s ≤ t , λ min ( H ( s )) ≥ λ 0 2 and (cid:107) w r ( s ) − w r (0) (cid:107) 2 ≤ R w .
Then we have | a ( t ) − a r (0)
| ≤ R (cid:48) a for all r ∈ [ m ] .
Proof of Lemma A.5.
Note for any i ∈ [ n ] and r ∈ [ m ] , w r (0) (cid:62) x i ∼ N (0 , 1) .
Therefore applying
Gaussian tail bound and union bound we have with probability at least 1 − δ , for all i ∈ [ n ] and
[Table]
[Table]
[Equation]
[Table]
[Table]
[Equation]
√
[Equation]
Lemma A.6.
If R (cid:48) w < R w and R (cid:48) a < R a , we have for all t ≥ 0 , λ min ( H ( t )) ≥ 12 λ 0 , for all r ∈ [ m ] , (cid:107) w r ( t ) − w r (0) (cid:107) 2 ≤ R (cid:48) w , | a r ( t ) − a r (0) | ≤ R (cid:48) a and (cid:107) y − u ( t ) (cid:107) 22 ≤ exp( − λ 0 t ) (cid:107) y − u (0) (cid:107) 22 .
Proof of Lemma A.6.
We prove by contradiction.
Let t > 0 be the smallest time that the conclusion does not hold.
Then either λ min ( H ( t )) <
λ 0 2 or there exists r ∈
[ m ] , (cid:107) w r ( t )
− w r (0) (cid:107) 2 ≤ R (cid:48) w or there exists r ∈ [ m ] , | a r − a r (0) | ≤ R (cid:48) a .
If λ min ( H ( t )) < λ 0 2 , by Lemma A.3, we know there exists s < t such that either there exists r ∈ [ m ] , (cid:107) w r ( s ) − w r (0) (cid:107) 2 ≤ R w
or there exists r ∈
[ m ] , | a r ( s ) − a r (0) | ≤ R a .
However, since R (cid:48) w < R w and R (cid:48) a < R a .
This contradicts with the minimality of t .
If there exists r ∈ [ m ] , (cid:107) w r − w r
0) (cid:107) 2 ≤ R (cid:48) w , then by Lemma
A.4, we know there exists s < t such that r ∈
[ m ] , | a r ( s ) − a r (0) | ≤ R a or λ min ( H ( s )) <
λ 0 2 .
However, since
R (cid:48)
w < R w and R (cid:48) a < R a .
This contradicts with the minimality of t .
The last case is similar for which we can simply apply Lemma A.5.
Based on Lemma A.2, we only need to ensure R (cid:48)
w < R w and R a < R (cid:48) a .
By the proof in Section 3, we know with probability at least δ , (cid:107) y − u (0) (cid:107) 2 ≤ C √ n δ for some large constant C .
Note our section on m in Theorem 3.3 sufﬁces to ensure R (cid:48)
w < R w and R a < R (cid:48) a .
We now complete the proof.
B T ECHNICAL P ROOFS FOR S ECTION
Proof of Corollary 4.1. We use the norm of gradient to bound this distance.
[Equation]
Proof of Lemma 4.1.
For a ﬁxed i ∈ [ n ] and r ∈ [ m ] , by anti-concentration inequality, we know
P ( A ir ) ≤ 2 R √ 2 π .
Thus we can bound the size of S ⊥ i in expectation.
[Equation]
Summing over i = 1 , . .
. , n , we have
[Equation]
Thus by Markov’s inequality, we have with probability at least 1 − δ
[Equation]
for some large positive constant C > 0 .