B IDIRECTIONALLY S ELF -N ORMALIZING N EURAL N ETWORKS
Anonymous authors
Paper under double-blind review
A BSTRACT
The problem of exploding and vanishing gradients has been a long-standing obstacle that hinders the effective training of neural networks.
Despite various tricks and techniques that have been employed to alleviate the problem in practice, there still lacks satisfactory theories or provable solutions.
In this paper, we address the problem from the perspective of high-dimensional probability theory.
We provide a rigorous result that shows, under mild conditions, how the exploding/vanishing gradient problem disappears with high probability if the neural networks have sufﬁcient width.
Our main idea is to constrain both forward and backward signal propagation in a nonlinear neural network through a new class of activation functions, namely Gaussian-Poincar ´
e normalized functions, and orthogonal weight matrices.
Experiments on both synthetic and real-world data validate our theory and conﬁrm its effectiveness on very deep neural networks when applied in practice.
I NTRODUCTION
Neural networks have brought unprecedented performance in various artiﬁcial intelligence tasks (Graves et al., 2013; Krizhevsky et al., 2012; Silver et al., 2017).
However, despite decades of research, training neural networks is still mostly guided by empirical observations and successful training often requires various heuristics and extensive hyperparameter tuning.
It is therefore desirable to understand the cause of the difﬁculty in neural network training and to propose theoretically sound solutions.
A major difﬁculty is the gradient exploding/vanishing problem (Glorot & Bengio, 2010; Hochreiter, 1991; Pascanu et al., 2013; Philipp et al., 2018).
That is, the norm of the gradient in each layer is either growing or shrinking at an exponential rate as the gradient signal is propagated from the top layer to bottom layer.
For deep neural networks, this problem might cause numerical overﬂow and make the optimization problem intrinsically difﬁcult, as the gradient in each layer has vastly different magnitude and therefore the optimization landscape becomes pathological.
One might attempt to solve the problem by simply normalizing the gradient in each layer.
Indeed, the adaptive gradient optimization methods (Duchi et al., 2011; Kingma & Ba, 2015; Tieleman & Hinton, 2012) implement this idea and have been widely used in practice.
However, one might also wonder if there is a solution more intrinsic to deep neural networks, whose internal structure if well-exploited would lead to further advances.
To enable the trainability of deep neural networks, batch normalization (Ioffe & Szegedy, 2015) was proposed in recent years and achieved widespread empirical success.
Batch normalization is a differentiable operation which normalizes its inputs based on mini-batch statistics and inserted between the linear and nonlinear layers.
It is reported that batch normalization can accelerate neural network training signiﬁcantly (Ioffe & Szegedy, 2015).
However, batch normalization does not solve the gradient exploding/vanishing problem (Philipp et al., 2018).
Indeed it is proved that batch normalization can actually worsen the problem (Yang et al., 2019).
Besides, batch normalization requires separate training and testing phases and might be ineffective when the mini-batch size is small (Ioffe, 2017).
The shortcomings of batch normalization motivate us to search for a more principled and generic approach to solve the gradient exploding/vanishing problem.
Alternatively, self-normalizing neural networks (Klambauer et al., 2017) and dynamical isometry theory (Pennington et al., 2017) were proposed to combat gradient exploding/vanishing problem.
In self-normalizing neural networks, the output of each network unit is constrained to have zero mean and unit variance.
Based on this motivation, a new activation function, scaled exponential linear unit (SELU), was devised.
In dynamical isometry theory, all singular values of the input-output
Jacobian matrix are constrained to be close to one at initialization.
This amounts to initializing the functionality of a network to be close to an orthogonal matrix.
While the two theories dispense batch normalization, it is shown that SELU still suffers from exploding/vanishing gradient problem and dynamical isometry restricts the functionality of the network to be close to linear (pseudolinearity) (Philipp et al., 2018).
In this paper, we follow the above line of research to investigate neural network trainability.
Our contributions are three-fold: First, we introduce bidirectionally self-normalizing neural network (BSNN) that consist of orthogonal weight matrices and a new class of activation functions which we call Gaussian-Poincar´e normalized (GPN) functions.
We show many common activation functions can be easily transformed into their respective GPN versions.
Second, we rigorously prove that the gradient exploding/vanishing problem disappears with high probability in BSNNs if the width of each layer is sufﬁciently large.
Third, with experiments on synthetic and real-world data, we conﬁrm that BSNNs solve the gradient vanishing/exploding problem to large extent while maintaining nonlinear functionality.
T HEORY
In this section, we introduce bidirectionally self-normalizing neural network (BSNN) formally and analyze its properties.
All the proofs of our results are left to Appendix.
To simplify the analysis, we deﬁne a neural network in a restricted sense as follows:
Deﬁnition 1 ( Neural Network ) .
A neural network is a function from R n to R n composed of layerwise operations for l = 1 , . . .
, L as
[Equation]
where W ( l ) ∈ R n × n , φ : R → R is a differentiable function applied element-wise to a vector, x (1) is the input and x ( L +1) is the output.
Under this deﬁnition, φ is called the activation function, { W ( l ) } Ll =1 are called the parameters, n is called the width and L is called the depth.
Superscript ( l ) denotes the l -th layer of a neural network.
The above formulation is similar to (Pennington et al., 2017) but we omit the bias term in (1) for simplicity as it plays no role in our analysis.
Let E be the objective function of { W ( l ) } Ll =1 and D ( l ) be a diagonal matrix with diagonal elements
D ( l ) ii = φ (cid:48)
( h ( l ) i ) , where φ (cid:48) denotes the derivative of φ .
Now, the error signal is back propagated via
[Equation]
and the gradient of the weight matrix for layer l can be computed as
[Equation]
To solve the gradient exploding/vanishing problem, we constrain the forward signal x ( l ) and the backward signal d ( l ) in order to constrain the norm of the gradient.
This leads to the following deﬁnition and proposition.
Deﬁnition 2 ( Bidirectional Self-Normalization ) .
A neural network is bidirectionally selfnormalizing if
[Equation] [Equation]
Proposition 1.
If a neural network is bidirectionally self-normalizing, then
[Equation] [Equation]
In the rest of this section, we derive the conditions under which bidirectional self-normalization is achievable for a neural network.
2. C ONSTRAINTS ON W EIGHT M ATRICES
We constrain the weight matrices to be orthogonal since multiplication by an orthogonal matrix preserves the norm of a vector.
For linear networks, this guarantees bidirectionally selfnoramlizingnormalization and its further beneﬁts are discussed in (Saxe et al., 2014).
Even for nonlinear neural networks, orthogonal constraints are shown to improve the trainability with proper scaling (Mishkin & Matas, 2016; Pennington et al., 2017).
2. C ONSTRAINTS ON A CTIVATION F UNCTIONS
To achieve bidirectionally self-noramlizingnormalization for a nonlinear network, it is not enough only to constrain the weight matrices.
We also need to constrain the activation function in such a way that both forward and backward signals are normalized.
To this end, we propose the following constraint that captures the relationship between a function and its derivative.
Deﬁnition 3 ( Gaussian-Poincar´e Normalization ) .
Function φ : R → R is Gaussian-Poincar ´
e normalized if it is differentiable and
[Equation]
The deﬁnition is inspired by the following theorem which shows the fundamental relationship between a function and its derivative under Gaussian measure.
Theorem 1 ( Gaussian-Poincar´e Inequality (Bogachev, 1998)) .
If function φ : R → R is differentiable with bounded E x
∼N (0 , 1) [ φ ( x ) 2 ] and E x ∼N (0 , 1) [ φ (cid:48) ( x ) 2 ] , then [Equation]
[Equation]
Note that there is an implicit assumption that the input is approximately Gaussian for a GaussianPoincar´e normalized (GPN) function.
Even though this is standard in the literature (Klambauer et al., 2017; Pennington et al., 2017; Schoenholz et al., 2017), we will rigorously prove that this assumption is valid when orthogonal weight matrices are used in equation 1.
Next, we state a property of GPN functions.
Proposition 2.
Function φ : R → R is Gaussian-Poincar´e normalized and E x ∼N (0 , 1) [ φ ( x )]
=
if and only if φ ( x ) = x or φ ( x ) = − x .
This result indicates that any nonlinear function with zero mean under Gaussian distribution ( e.g .,
Tanh and SELU) is not GPN.
Now we show that a large class of activation functions can be converted into their respective GPN versions using an afﬁne transformation.
Proposition 3.
For any differentiable function φ : R → R with non-zero and bounded E x ∼N (0 , 1) [ φ ( x ) 2 ] and E x ∼N (0 , 1) [ φ (cid:48)
( x ) 2 ] , there exist two constants a and b such that aφ ( x ) + b is Gaussian-Poincar´e normalized.
To obtain a and b , one can use numerical procedure to compute the values of E x ∼N (0 , 1) [ φ (cid:48)
( x ) 2 ] ,
E x ∼N (0 , 1) [ φ ( x ) 2 ] and E x ∼N (0 , 1) [ φ ( x )] and then solve the quadratic equations
[Equation]
We computed a and b (not unique) for several common activation functions with their default hyperparameters
and the results are listed in Table 1.
Note that ReLU, LeakyReLU and SELU are not differentiable at x =
but they can be regarded as approximations of their smooth counterparts.
We ignore such point and evaluate the integrals for x ∈ ( −∞ , 0) ∪ (0 , ∞ ) .
With the orthogonal constraint on the weight matrices and the Gaussian-Poincar´e normalization on the activation function, we prove that bidirectionally self-noramlizingnormalization is achievable with high probability under mild conditions in the next subsection.
[Table]
2. N ORM -P RESERVATION T HEOREMS
The bidirectionally self-noramlizingnormalization may not be achievable precisely in general unless the neural network is a linear one.
Therefore, we investigate the properties of neural networks in a probabilistic framework.
The random matrix theory and the high-dimensional probability theory allow us to characterize the behaviors of a large class of neural networks by its mean behavior, which is signiﬁcantly simpler to analyze.
Therefore, we study neural networks of random weights whose properties may shed light on the trainability of neural networks in practice.
First, we need a probabilistic version of the vector norm constraint.
Deﬁnition 4 ( Thin-Shell Concentration ) .
Random vector x ∈ R n is thin-shell concentrated if for any (cid:15) >
[Equation]
(11) as n → ∞ .
The deﬁnition is modiﬁed from the one in (Bobkov, 2003).
Examples of thin-shell concentrated distributions include standard multivariate Gaussian and any distribution on the n -dimensional sphere of radius √ n .
Assumptions.
To prove the main results, i.e ., the norm-preservation theorems, we require the following assumptions:
Random vector x ∈ R n is thin-shell concentrated.
Random orthogonal matrix W = ( w 1 , w 2 , ..., w n ) T is uniformly distributed.
Function φ : R → R is Gaussian-Poincar´e normalized.
Function φ and its derivative are Lipschitz continuous.
The above assumptions are not restrictive.
For Assumption 1, one can always normalize the input vectors of a neural network.
For Assumption 2, orthogonal constraint or its relaxation has already been employed in neural network training (Brock et al., 2017).
Note, in Assumption 2, uniformly distributed means that W is distributed under Haar measure, which is the unique rotation invariant probability measure on orthogonal matrix group.
We refer the reader to (Meckes, 2019) for details.
Furthermore, all the activation functions or their smooth counterparts listed in Table 1 satisfy
Assumptions 3 and 4.
With the above assumptions, we can prove the following norm-preservation theorems.
Theorem 2 ( Forward Norm-Preservation ) .
Random vector
[Equation]
is thin-shell concentrated.
This result shows the transformation (orthogonal matrix followed by the GPN activation function) can preserve the norm of its input with high probability.
Since the output is thin-shell concentrated, it serves as the input for the next layer and so on.
Hence, the forward pass can preserve the norm of its input in each layer along the forward path when n is sufﬁciently large.
Theorem 3 ( Backward Norm-Preservation ) .
Let D be the diagonal matrix whose diagonal elements are D ii = φ (cid:48) ( w Ti x ) and y ∈ R n be a ﬁxed vector with bounded (cid:107) y (cid:107) ∞ .
Then for any (cid:15) >
[Equation]
(13) as n → ∞ .
This result shows that the multiplication by D preserves the norm of its input with high probability.
Since orthogonal matrix W also preserves the norm of its input, when the gradient error signal is propagated backwards as in (2), the norm is preserved in each layer along the backward path when n is sufﬁcient large.
Hence, combining Theorems 2 and 3, we proved that bidirectionally self-noramlizingnormalization is achievable with high probability if the neural network is wide enough and the conditions in the Assumptions are satisﬁed.
Then by Proposition 1, the gradient exploding/vanishing problem disappears with high probability.
Sketch of the proofs.
The proofs of Theorems 2 and 3 are mainly based on a phenomenon in high-dimensional probability theory, concentration of measure.
We refer the reader to (Vershynin, 2018) for an introduction to the subject.
Brieﬂy, it can be shown that for some high-dimensional probability distributions, most mass is concentrated around a certain range.
For example, while most mass of a low-dimensional standard multivariate Gaussian distribution is concentrated around the center, most mass of a high-dimensional standard multivariate Gaussian distribution is concentrated around a thin-shell.
Furthermore, the random variables transformed by Lipschitz functions are also concentrated around certain values.
Using this phenomenon, it can be shown that rows { w i } of a random orthogonal matrix W in high dimension are approximately independent random unit vectors and the inner product w Ti x for thin-shell concentrated vector x can be shown to be approximately
Gaussian.
Then from the assumptions that φ is GPN and φ and φ (cid:48) are Lipschitz continuous, the proofs follow.
Each of these steps is rigorously proved in Appendix.
E XPERIMENTS
We verify our theory on both synthetic and real-world data.
More experimental results can be found in Appendix.
In short, while very deep neural networks with non-GPN activations show vanishing/exploding gradients, GPN versions show stable gradients and improved trainability in both synthetic and real data.
Furthermore, compared to dynamical isometry theory, BSNNs do not exhibit pseudo-linearity and maintain nonlinear functionality.
[Equation]
We create synthetic data to test the norm-preservation properties of the neural networks.
The input x 1 is 500 data points of random standard Gaussian vectors of 500 dimension.
The gradient error ∂E/∂ x L +1 is also random standard Gaussian vector of 500 dimension.
All the neural networks have depth 200.
All the weight matrices are random orthogonal matrices uniformly generated.
No training is performed.
In Figure 1, we show the norm of inputs and gradient of the neural networks of width 500.
From the results, we can see that with GPN, the gradient exploding/vanishing problem is eliminated to large extent.
The neural network with Tanh activation function does not show gradient exploding/vanishing problem either.
However, (cid:107) x ( l ) (cid:107) is close to zero for large l and each layer is close to a linear one since Tanh ( x ) ≈ x
when x ≈ 0 (pseudo-linearity), for which dynamical isometry is achieved.
One might wonder if bidirectionally self-noramlizingnormalization has the same effect as dynamical isometry in solving the gradient exploding/vanishing problem, that is, to make the neural network close to an orthogonal matrix.
To answer this question, we show the histogram of φ (cid:48)
( h ( l ) i ) in Figure 2.
If the functionality of a neural network is close to an orthogonal matrix, since the weight matrices are orthogonal, then the values of φ (cid:48) ( h ( l ) i ) would concentrate around one (Figure 2 (a)), which is not the case for BSNNs (Figure 2 (b)).
This shows that BSNNs do not suffer from the gradient vanishing/explosion problem while exhibiting nonlinear functionality.
In Figure 7 in Appendix, we show the gradient norm of BSNNs with varying width.
There we note that as the width increases, the norm of gradient in each layer of the neural network becomes more equalized, as predicted by our theory.
[Figure] [Figure]
[Figure] [Figure]
[Table]
3. R EAL -W ORLD D ATA
We run experiments on real-world image datasets MNIST and CIFAR-10.
The neural networks have width 500 and depth 200 (plus one unconstrained layer at bottom and one at top to ﬁt the dimensionality of the input and output).
We use stochastic gradient descent of momentum 0.
with mini-batch size 64 and learning rate 0.0001.
The training is run for 50 epochs for MNIST and 100 epochs for CIFAR-10.
We do not use data augmentation.
Since it is computationally expensive to enforce the orthogonality constraint, we simply constrain each row of the weight matrix to have l 2 norm one as a relaxation of orthogonality by the following parametrization
W = ( v 1 / (cid:107) v 1 (cid:107) 2 , v 2 / (cid:107) v 2 (cid:107) 2 , ..., v n / (cid:107) v n (cid:107) 2 ) T and optimize V =
( v 1 , v 2 , ..., v n ) T as an unconstrained problem.
We summarize the results in Table 2.
We can see that, for activation functions ReLU, LeakyReLU and GELU, the neural networks are not trainable.
But once these functions are GPN, the neural network can be trained.
GPN activation functions consistently outperform their unnormalized counterparts in terms of the trainability, as the training accuracy is increased, but not necessarily generalization ability.
We show the test accuracy during training in Figure 3, from which we can see the training is accelerated when SELU is GPN.
ReLU, LeakyReLU and GELU, if not GPN, are completely untrainable due to the vanished gradient (see Appendix).
We observe that batch normalization leads to gradient explosion when combining with any of the activation functions.
This conﬁrms the claim of (Philipp et al., 2018) and (Yang et al., 2019) that batch normalization does not solve the gradient exploding/vanishing problem.
On the other hand, without batch normalization the neural network with any GPN activation function has stable gradient magnitude throughout training (see Appendix).
This indicates that BSNNs can dispense with batch normalization and therefore avoid its shortcomings.
R ELATED W ORK
We compare our theory to several most relevant theories in literature.
A key distinguishing feature of our theory is that we provide rigorous proofs of the conditions under which the exploding/vanishing problem disappears.
To the best of our knowledge, this is the ﬁrst time that the problem is provably solved for nonlinear neural networks.
4. S ELF -N ORMALIZING N EURAL N ETWORKS
Self-normalizing neural networks enforce zero mean and unit variance for the output of each unit with the SELU activation function (Klambauer et al., 2017).
However, as pointed out in (Philipp et al., 2018) and conﬁrmed in our experiments, only constraining forward signal propagation does not solve the gradient exploding/vanishing problem since the norm of the backward signal can grow or shrink.
The signal propagation in both directions needs to be constrained as in our theory.
4. D EEP S IGNAL P ROPAGATION
Our theory is developed from the deep signal propagation theory (Poole et al., 2016; Schoenholz et al., 2017).
Both theories require E x ∼N (0 , 1) [ φ (cid:48)
( x ) 2 ] = 1 .
However, ours also requires the quantity
E x ∼N (0 , 1) [ φ ( x ) 2 ] to be one while in Poole et al .
(Poole et al., 2016; Schoenholz et al., 2017) it can be an arbitrary positive number.
We emphasize that it is desirable to enforce E x ∼N (0 , 1) [ φ ( x ) 2 ]
=to avoid trivial solutions.
For example, if φ ( x ) =
Tanh ( (cid:15)x ) with (cid:15) ≈ 0 , then φ ( (cid:15)x ) ≈
(cid:15)x and the neural network becomes essentially a linear one for which depth is unnecessary (pseudo-linearity (Philipp et al., 2018)).
This is observed in Figure 1 (a).
Moreover, in (Poole et al., 2016; Schoenholz et al., 2017)
the signal propagation analysis is done based on random weights under i.i.d.
Gaussian distribution whereas we proved how one can solve gradient vanishing/exploding problem assuming the weight matrices are orthogonal and uniformly distributed under Haar measure.
4. D YNAMICAL I SOMETRY
Dynamical isometry theory (Pennington et al., 2017) enforces the Jacobian matrix of the inputoutput function of a neural network to have all singular values close to one.
Since the weight matrices are constrained to be orthogonal, it is equivalent to enforce each D ( l ) in (2) to be close to the identity matrix, which implies the functionality of neural network at initialization is close to an orthogonal matrix (pseudo-linearity).
This indeed enables trainability since linear neural networks with orthogonal weight matrices do not suffer from the gradient exploding/vanishing problem.
As neural networks need to learn a nonlinear input-output functionality to solve certain tasks, during training the weights of a neural network are unconstrained so that the neural network would move to a nonlinear region where the gradient exploding/vanishing problem might return.
In our theory, although the orthogonality of weight matrices is also required, we approach the problem from a different perspective.
We do not encourage the linearity at initialization.
The neural network can be initialized to be nonlinear and stay nonlinear during the training even when the weights are constrained.
C ONCLUSION
In this paper, we have introduced bidirectionally self-normalizing neural network (BSNN) which constrains both forward and backward signal propagation using a new class of Gaussian-Poincar´e normalized activation functions and orthogonal weight matrices.
BSNNs are not restrictive in the sense that many commonly used activation functions can be Gaussian-Poincar´e normalized.
We have rigorously proved that gradient vanishing/exploding problem disappears in BSNNs with high probability under mild conditions.
Experiments on synthetic and real-world data conﬁrm the validity of our theory and demonstrate that BSNNs have excellent trainability without batch normalization.
Currently, the theoretical analysis is limited to same width, fully-connected neural networks.
Future work includes extending our theory to more sophisticated networks such as convolutional architectures as well as investigating the generalization capabilities of BSNNs.
A PPENDIX A P ROOFS
Proposition 1.
If a neural network is bidirectionally self-noramlizingnormalizing, then
[Equation]
Proof.
For each l , we have
[Equation]
By the deﬁnition of bidirectional self-normalization, we have (cid:107) ∂E
∂ W (1) (cid:107) F
= ...
= (cid:107) ∂E
∂ W ( L ) (cid:107) F .
Proposition 2.
Function φ : R → R is Gaussian-Poincar´e normalized and E x ∼N (0 , 1) [ φ ( x )]
=
if and only if φ ( x ) = x or φ ( x ) = − x .
Proof.
Since E x ∼N (0 , 1) [ φ ( x ) 2 ] < ∞ and E x ∼N (0 , 1)
[ φ (cid:48)
( x ) 2 ] < ∞ , φ ( x ) and φ (cid:48)
( x ) can be expanded in terms of Hermite polynomials.
Let the Hermite polynomial of degree
k be
[Equation]
and due to H (cid:48) k
( x ) = √ kH k − 1 ( x ) , we have
[Equation] [Equation]
Since E x ∼N (0 , 1) [ φ ( x )]
= 0 , we have a
[Equation]
[Equation]
[Equation]
(24) = 0 .
Since
E x ∼N (0 , 1)
[Equation]
[Equation]
[Equation]
(26) and Hermite polynomials are orthonormal, we have
E x ∼N (0 , 1) [
[Equation]
( x [Equation]
[Equation]
∞ (cid:88) k
= a
[Equation]
k
= E x ∼N (0 , 1) [
[Equation]
(cid:48) ( x ) 2 ] = ∞ (cid:88) k = ka 2 k
= 1 .
Therefore, we have ∞ (cid:88) k
= ka 2 k − ∞ (cid:88) k
= a 2 k =(28) that is ∞ (cid:88)
[Equation]
Since each term in (cid:80) ∞ k =2 ( k − 1) a 2 k is nonnegative, the only solution is a k = 0 for k ≥ 2 .
And since E x ∼N (0 , 1) [ φ ( x ) 2 ] = a 2 1 = 1 , we have a 1 = ± 1 .
Hence, φ ( x ) = ± H 1 ( x ) = ± x .
Proposition 3.
For any differentiable function φ : R → R with non-zero and bounded E x ∼N (0 , 1) [ φ ( x ) 2 ] and E x ∼N (0 , 1) [ φ (cid:48)
( x ) 2 ] , there exist two constants a and b such that aφ ( x ) + b is Gaussian-Poincar´e normalized.
Proof.
Let ϕ ( x ) =
φ ( x ) + c .
Then let
[Equation]
Therefore, ψ ( c ) is a quadratic function of c .
We also have ψ ( c ) > 0 as c → ∞ and ψ ( − E x ∼N (0 , 1) [ φ ( x )]) ≤ 0 due to Gaussian-Poincar´e inequality.
Hence, there exists c for which ψ ( c ) =
such that E x ∼N (0 , 1) [( φ ( x )+ c ) 2 ] =
E x ∼N [Equation] 0 , 1) [ φ
[Equation]
( x ) 2 ] .
Let a =
( E x ∼N [Equation] [Equation] 1) [ φ
[Equation] ( x ) 2 ]) − [Equation]
and b = ac , we have E x ∼N (0 , 1) [( aφ ( x ) + b ) 2 ] =
E x ∼N (0 , 1) [( aφ (cid:48) ( x )) 2 ] = 1.
The proof is largely due to (Eldredge, 2020) with minor modiﬁcation in here.
Assumptions.
Random vector x ∈ R n is thin-shell concentrated.
Random orthogonal matrix W = ( w 1 , w 2 , ..., w n ) T is uniformly distributed.
Function φ : R → R is Gaussian-Poincar´e normalized.
Function φ and its derivative are Lipschitz continuous.
Theorem 2 ( Forward Norm-Preservation ) .
Random vector
[Equation]
[Equation]
33) is thin-shell concentrated.
Theorem 3 ( Backward Norm-Preservation ) .
Let D be the diagonal matrix whose diagonal elements are D ii = φ (cid:48) ( w Ti x ) and y ∈ R n be a ﬁxed vector with bounded (cid:107) y (cid:107) ∞ .
Then for any (cid:15) >
[Equation]
34) as n → ∞ .
Notations.
S n − 1 = { x ∈ R n : (cid:107) x (cid:107) 2 = 1 } .
O ( n ) is the orthogonal matrix group of size n .
1 {·} denotes the indicator function.
0 n denotes the vector of dimension n and all elements equal to zero.
I n denotes the identity matrix of size n × n .
Lemma 1.
If random variable x ∼ N (0 , 1) and function f : R → R is Lipschitz continuous,
then random variable f
( x ) is sub-gaussian.
Proof.
Due to the Gaussian concentration theorem (Theorem 5.2.2 in (Vershynin, 2018)), we have
[Equation]
(35) where (cid:107) · (cid:107) ψ 2 denotes sub-gaussian norm, C is a constant and K is the Lipschitz constant of f .
This implies f ( x ) − E [ f ( x )] is sub-gaussian (Proposition 2.5.2 in (Vershynin, 2018)).
Therefore f ( x ) is sub-gaussian (Lemma 2.6.8 in (Vershynin, 2018)).
Lemma 2.
Let x = ( x 1 , x 2 , ..., x n ) ∈ R n be a random vector that each coordinate x i is independent and sub-gaussian and E [ x 2 i ] = 1 .
Let y =
( y 1 , y 2 , ..., y n ) ∈ R n be a ﬁxed vector with bounded (cid:107) y (cid:107) ∞ .
Then
P
[Equation]
as n → ∞ .
Proof.
Since y i x i is sub-gaussian, then
y 2 i x
i is sub-exponential (Lemma 2.7.6 in (Vershynin, 2018)).
Since E [ y 2 i x 2 i ]
=
y 2 i E
[ x 2 i ] =
y 2 i , y 2 i x
i − y
i is sub-exponential with zero mean (Exercise 2.7.10 in (Vershynin, 2018)).
Applying Bernsteins inequality (Corollary 2.8.3 in (Vershynin, 2018)), we proved the lemma.
Lemma 3.
Let z ∼ N ( 0 n , I n ) .
Then for any 0 < δ < [Equation] See (Alberts & Khoshnevisan, 2018) (Theorem 1.2) for a proof.
Lemma 4.
Let z ∼ N ( 0 n , I n ) .
Then z / (cid:107) z (cid:107) 2 is uniformly distributed on S n − 1 . See (Dawkins, 2016) for a proof.
Lemma 5. Let z =
( z 1 , z 2 , ..., z n ) ∼ N ( 0 n , I n ) , a = ( a 1 , a 2 , ..., a n ) be a ﬁxed vector with bounded (cid:107) a (cid:107) ∞ and f : R → R be a continuous function.
Then for any (cid:15) >
[Equation]
as n → ∞ . Proof.
Since
[Equation]
if, as n → ∞ ,
[Equation]
n
[Equation]
then
[Equation]
For 0 < δ < 1 , let A = (cid:110)
[Equation]
i
[Equation]
U
[Equation]
[Equation] Then P
[Equation]
n
[Equation]
Let δ = n − 1 [-√ n (cid:12)(cid:12)(cid:12) ≤ 1 √ n (cid:15) 1 − (cid:15) (54) and (cid:12) (cid:12)(cid:12) y T v (cid:107) v (cid:107) 2 − y T v √ n (cid:12) (cid:12)(cid:12) = (cid:12) (cid:12)(cid:12) 
-]/ 4 .
From Lemma 3, we have, as n → ∞ ,
[Equation]
[Equation]
For z ∈ U δ and δ = n − 1 / 4 , we have (cid:107) z (cid:107) 2 → √ n , √ {+n/ (cid:107) z (cid:107) 
+}
z
i
→
{+z i and therefore f ( √ n/ (cid:107) z (cid:107) 
+}
z i ) →
f ( z i ) as n → ∞ .
Hence, (cid:82) U δ 1 { z ∈ A } d z → 0 , as n → ∞ .
Lemma 6.
Let random matrix W be uniformly distributed on O ( n ) random vector θ be uniformly distributed on S n − 1 and random vector x ∈ R n be thin-shell concentrated. Then Wx → √ n θ as n → ∞ .
Proof.
Let y ∈ R n be any vector with (cid:107) y (cid:107) 2 = √ n and e = ( √ n, 0 , ..., 0)
∈ R n .
Since W is uniformly distributed, Wy has the same distribution as We .
We is the ﬁrst row of √ n W , which is equivalent to random vector √ n θ .
Since x is thin-shell concentrated, x → √ n (cid:107) x (cid:107) 2 x = y and therefore
Wx → √ n θ as n → ∞ .
Proof of Theorem 2.
Let z
=
( z 1 , z 2 , ..., z n ) ∼ N ( 0 , I ) .
Due to Lemma 1, random variable φ ( z i ) is sub-gaussian.
Since φ is Gaussian-Poincar´e normalized, E z i ∼N (0 , 1) [ φ ( z i ) 2 ] = 1 .
Applying Lemma 2 with each y i = 1 , we have for (cid:15) >
[Equation]
as n → ∞ .
Due to Lemma 4 and 5 (with each a i = 1 ), for random vector θ = ( θ 1 , θ 2 , ..., θ n ) uniformly distributed on S n − 1 , we have
[Equation]
and therefore P
[Equation]
as n → ∞ .
Then from Lemma 6, we have Wx → √ n θ and therefore
P
[Equation]
as n → ∞ .
Proof of Theorem 3.
Let z = ( z 1 , z 2 , ..., z n ) be the standard multivariate Gaussian random vectors.
Due to Lemma 1, random variable φ (cid:48) ( z i ) is sub-gaussian.
Since φ is Gaussian-Poincar ´ e normalized, E z i ∼N (0 , 1) [ φ (cid:48)
( z i ) 2 ] = 1 .
Applying Lemma 2, we have
[Equation]
as n → ∞ .
Due to Lemma 4 and 5 (with each a i = y 2 i ), for random vector θ = ( θ 1 , θ 2 , ..., θ n ) uniformly distributed on S n − 1 , we have
[Equation]
52) and therefore
[Equation]
53) as n → ∞ .
Then from Lemma 6, we have Wx → √ n θ and therefore
[Equation]
54) as n → ∞ .
A PPENDIX B A DDITIONAL E XPERIMENTS
Due to the space limitation
[Figure]
we only showed the experiments with Tanh and SELU activation functions in the main text.
In this section
[Figure]
we show the experiments with ReLU, LeakyReLU, ELU and SELU.
Additionally, we also measure the gradient exploding/vanishing during training on the real-world data
[Figure]
B. S YNTHETIC D ATA
We show the ﬁgures of the experimental results in addition to the ones in the main text
[Figure]
[Figure] [Figure]
[Figure]
B. R EAL -W ORLD D ATA
In Figure 8, we show the test accuracy during training on MNIST.
In Figure 9, we show the test accuracy during training on CIFAR-10.
In Figure 10, 11, 12 and 13, we show a measure of gradient exploding/vanishing during training for different activation functions.
The measure is deﬁned as the ratio of the maximum gradient norm and the minimum gradient norm across layers.
Since we use the parametrization
[Equation]
(55) with V =
( v 1 , v 2 , ..., v n ) T , the gradient norm ratio is deﬁned on the unconstrained weights V , that is,
[Equation]
Note that for ReLU, LeakyReLU and GELU, the gradient vanishes during training in some experiments and therefore the plots are empty.
From the ﬁgures, we can see that batch normalization leads to gradient explosion especially at the early stage of training.
On the other hand, without batch normalization, the gradient is stable throughout training for GPN activation functions.
[Figure]
[Figure] [Figure]
[Figure] [Figure]
[Figure] [Figure]