Title: ON LARGE-BATCH TRAINING FOR DEEP LEARNING: GENERALIZATION GAP AND SHARP MINIMA

Abstract: The stochastic gradient descent (SGD) method and its variants are algorithms of choice for many Deep Learning tasks. These methods operate in a small-batch regime wherein a fraction of the training data, say 32-512 data points, is sampled to compute an approximation to the gradient. It has been observed in practice that when using a larger batch there is a degradation in the quality of the model, as measured by its ability to generalize. We investigate the cause for this generalization drop in the large-batch regime and present numerical evidence that supports the view that large-batch methods tend to converge to sharp minimizers of the training and testing functions-and as is well known, sharp minima lead to poorer generalization. In contrast, small-batch methods consistently converge to flat minimizers, and our experiments support a commonly held view that this is due to the inherent noise in the gradient estimation. We discuss several strategies to attempt to help large-batch methods eliminate this generalization gap.

1. INTRODUCTION:

Deep Learning has emerged as one of the cornerstones of large-scale machine learning. Deep Learning models are used for achieving state-of-the-art results on a wide variety of tasks including computer vision, natural language processing and reinforcement learning; see (Bengio et al., 2016) and the references therein. The problem of training these networks is one of non-convex optimization. Mathematically, this can be represented as:

EQUATION

where f i is a loss function for data point i ∈ {1, 2, • • • , M } which captures the deviation of the model prediction from the data, and x is the vector of weights being optimized. The process of optimizing this function is also called training of the network. Stochastic Gradient Descent (SGD) (Bottou, 1998; Sutskever et al., 2013) and its variants are often used for training deep networks.

These methods minimize the objective function f by iteratively taking steps of the form:

EQUATION

where B k ⊂ {1, 2, • • • , M } is the batch sampled from the data set and α k is the step size at iteration k. These methods can be interpreted as gradient descent using noisy gradients, which and are often referred to as mini-batch gradients with batch size |B k |. SGD and its variants are employed in a small-batch regime, where |B k | M and typically |B k | ∈ {32, 64, • • • , 512}. These configurations have been successfully used in practice for a large number of applications; see e.g. (Simonyan & Zisserman, 2014; Graves et al., 2013; Mnih et al., 2013) . Many theoretical properties of these methods are known. These include guarantees of: (a) convergence to minimizers of strongly-convex functions and to stationary points for non-convex functions (Bottou et al., 2016) , (b) saddle-point avoidance (Ge et al., 2015; Lee et al., 2016) , and (c) robustness to input data (Hardt et al., 2015) .

Stochastic gradient methods have, however, a major drawback: owing to the sequential nature of the iteration and small batch sizes, there is limited avenue for parallelization. While some efforts have been made to parallelize SGD for Deep Learning (Dean et al., 2012; Das et al., 2016; Zhang et al., 2015) , the speed-ups and scalability obtained are often limited by the small batch sizes. One natural avenue for improving parallelism is to increase the batch size |B k |. This increases the amount of computation per iteration, which can be effectively distributed. However, practitioners have observed that this leads to a loss in generalization performance; see e.g. (LeCun et al., 2012) . In other words, the performance of the model on testing data sets is often worse when trained with largebatch methods as compared to small-batch methods. In our experiments, we have found the drop in generalization (also called generalization gap) to be as high as 5% even for smaller networks.

In this paper, we present numerical results that shed light into this drawback of large-batch methods. We observe that the generalization gap is correlated with a marked sharpness of the minimizers obtained by large-batch methods. This motivates efforts at remedying the generalization problem, as a training algorithm that employs large batches without sacrificing generalization performance would have the ability to scale to a much larger number of nodes than is possible today. This could potentially reduce the training time by orders-of-magnitude; we present an idealized performance model in the Appendix C to support this claim.

The paper is organized as follows. In the remainder of this section, we define the notation used in this paper, and in Section 2 we present our main findings and their supporting numerical evidence. In Section 3 we explore the performance of small-batch methods, and in Section 4 we briefly discuss the relationship between our results and recent theoretical work. We conclude with open questions concerning the generalization gap, sharp minima, and possible modifications to make large-batch training viable. In Appendix E, we present some attempts to overcome the problems of large-batch training.

1.1. NOTATION:

We use the notation f i to denote the composition of loss function and a prediction function corresponding to the i th data point. The vector of weights is denoted by x and is subscripted by k to denote an iteration. We use the term small-batch (SB) method to denote SGD, or one of its variants like ADAM (Kingma & Ba, 2015) and ADAGRAD (Duchi et al., 2011) , with the proviso that the gradient approximation is based on a small mini-batch. In our setup, the batch B k is randomly sampled and its size is kept fixed for every iteration. We use the term large-batch (LB) method to denote any training algorithm that uses a large mini-batch. In our experiments, ADAM is used to explore the behavior of both a small or a large batch method.

2. DRAWBACKS OF LARGE-BATCH METHODS:



2.1. OUR MAIN OBSERVATION:

As mentioned in Section 1, practitioners have observed a generalization gap when using large-batch methods for training deep learning models. Interestingly, this is despite the fact that large-batch methods usually yield a similar value of the training function as small-batch methods. One may put forth the following as possible causes for this phenomenon: (i) LB methods over-fit the model; (ii) LB methods are attracted to saddle points; (iii) LB methods lack the explorative properties of SB methods and tend to zoom-in on the minimizer closest to the initial point; (iv) SB and LB methods converge to qualitatively different minimizers with differing generalization properties. The data presented in this paper supports the last two conjectures.

The main observation of this paper is as follows:

The lack of generalization ability is due to the fact that large-batch methods tend to converge to sharp minimizers of the training function. These minimizers are characterized by a significant number of large positive eigenvalues in ∇ 2 f (x), and tend to generalize less well. In contrast, small-batch methods converge to flat minimizers characterized by having numerous small eigenvalues of ∇ 2 f (x). We have observed that the loss function landscape of deep neural networks is such that large-batch methods are attracted to regions with sharp minimizers and that, unlike small-batch methods, are unable to escape basins of attraction of these minimizers.

The concept of sharp and flat minimizers have been discussed in the statistics and machine learning literature. (Hochreiter & Schmidhuber, 1997) (informally) define a flat minimizer x as one for which the function varies slowly in a relatively large neighborhood of x. In contrast, a sharp minimizer x is such that the function increases rapidly in a small neighborhood of x. A flat minimum can be described with low precision, whereas a sharp minimum requires high precision. The large sensitivity of the training function at a sharp minimizer negatively impacts the ability of the trained model to generalize on new data; see Figure 1 for a hypothetical illustration. This can be explained through the lens of the minimum description length (MDL) theory, which states that statistical models that require fewer bits to describe (i.e., are of low complexity) generalize better (Rissanen, 1983) . Since flat minimizers can be specified with lower precision than to sharp minimizers, they tend to have better generalization performance. Alternative explanations are proffered through the Bayesian view of learning (MacKay, 1992) , and through the lens of free Gibbs energy; see e.g. Chaudhari et al. (2016) .

Flat Minimum Sharp Minimum:

Training Function Testing Function 

f (x)

2.2. NUMERICAL EXPERIMENTS:

In this section, we present numerical results to support the observations made above. To this end, we make use of the visualization technique employed by (Goodfellow et al., 2014b ) and a proposed heuristic metric of sharpness (Equation ( 4)). We consider 6 multi-class classification network configurations for our experiments; they are described in Table 1 . The details about the data sets and network configurations are presented in Appendices A and B respectively. As is common for such problems, we use the mean cross entropy loss as the objective function f .

The networks were chosen to exemplify popular configurations used in practice like AlexNet (Krizhevsky et al., 2012) and VGGNet (Simonyan & Zisserman, 2014) . Results on other networks and using other initialization strategies, activation functions, and data sets showed similar behavior.

Since the goal of our work is not to achieve state-of-the-art accuracy or time-to-solution on these tasks but rather to characterize the nature of the minima for LB and SB methods, we only describe the final testing accuracy in the main paper and ignore convergence trends.

For all experiments, we used 10% of the training data as batch size for the large-batch experiments and 256 data points for small-batch experiments. We used the ADAM optimizer for both regimes.

Experiments with other optimizers for the large-batch experiments, including ADAGRAD (Duchi et al., 2011) , SGD (Sutskever et al., 2013) and adaQN (Keskar & Berahas, 2016) , led to similar results. All experiments were conducted 5 times from different (uniformly distributed random) starting points and we report both mean and standard-deviation of measured quantities. The baseline performance for our setup is presented Table 2 . From this, we can observe that on all networks, both approaches led to high training accuracy but there is a significant difference in the generalization performance. The networks were trained, without any budget or limits, until the loss function ceased to improve. We emphasize that the generalization gap is not due to over-fitting or over-training as commonly observed in statistics. This phenomenon manifest themselves in the form of a testing accuracy curve that, at a certain iterate peaks, and then decays due to the model learning idiosyncrasies of the training data. This is not what we observe in our experiments; see Figure 2 for the training-testing curve of the F 2 and C 1 networks, which are representative of the rest. As such, early-stopping heuristics aimed at preventing models from over-fitting would not help reduce the generalization gap. The difference between the training and testing accuracies for the networks is due to the specific choice of the network (e.g. AlexNet, VGGNet etc.) and is not the focus of this study. Rather, our goal is to study the source of the testing performance disparity of the two regimes, SB and LB, on a given network model.

2.2.1. PARAMETRIC PLOTS:

We first present parametric 1-D plots of the function as described in (Goodfellow et al., 2014b ). Let x s and x indicate the solutions obtained by running ADAM using small and large batch sizes respectively. We plot the loss function, on both training and testing data sets, along a line-segment containing the two points. Specifically, for α ∈ [-1, 2], we plot the function f (αx + (1 -α)x s ) and also superimpose the classification accuracy at the intermediate points; see Figure 3 1 . For this experiment, we randomly chose a pair of SB and LB minimizers from the 5 trials used to generate the data in Table 2 . The plots show that the LB minima are strikingly sharper than the SB minima in this one-dimensional manifold. The plots in Figure 3 only explore a linear slice of the function, but in Figure 7 in Appendix D, we plot f (sin( απ 2 )x + cos( απ 2 )x s ) to monitor the function along a curved path between the two minimizers . There too, the relative sharpness of the minima is evident.

2.2.2. SHARPNESS OF MINIMA:

So far, we have used the term sharp minimizer loosely, but we noted that this concept has received attention in the literature (Hochreiter & Schmidhuber, 1997) . Sharpness of a minimizer can be characterized by the magnitude of the eigenvalues of ∇ 2 f (x), but given the prohibitive cost of this computation in deep learning applications, we employ a sensitivity measure that, although imperfect, is computationally feasible, even for large networks. It is based on exploring a small neighborhood of a solution and computing the largest value that the function f can attain in that neighborhood. We use that value to measure the sensitivity of the training function at the given local minimizer. Now, since the maximization process is not accurate, and to avoid being mislead by the case when a large value of f is attained only in a tiny subspace of R n , we perform the maximization both in the entire space R n as well as in random manifolds. For that purpose, we introduce an n × p matrix A, whose columns are randomly generated. Here p determines the dimension of the manifold, which in our experiments is chosen as p = 100. Specifically, let C denote a box around the solution over which the maximization of f is performed, and let A ∈ R n×p be the matrix defined above. In order to ensure invariance of sharpness to problem dimension and sparsity, we define the constraint set C as:

EQUATION

where A + denotes the pseudo-inverse of A. Thus controls the size of the box. We can now define our measure of sharpness (or sensitivity). Metric 2.1. Given x ∈ R n , > 0 and A ∈ R n×p , we define the (C , A)-sharpness of f at x as:

EQUATION

Unless specified otherwise, we use this metric for sharpness for the rest of the paper; if A is not specified, it is assumed to be the identity matrix, I n . (We note in passing that, in the convex optimization literature, the term sharp minimum has a different definition (Ferris, 1988) , but that concept is not useful for our purposes.)

In Tables 3 and 4 , we present the values of the sharpness metric (4) for the minimizers of the various problems. Table 3 explores the full-space (i.e., A = I n ) whereas Table 4 uses a randomly sampled n × 100 dimensional matrix A. We report results with two values of , (10 -3 , 5 • 10 -4 ). In all experiments, we solve the maximization problem in Equation ( 4) inexactly by applying 10 iterations of L-BFGS-B (Byrd et al., 1995) . This limit on the number of iterations was necessitated by the large cost of evaluating the true objective f . Both tables show a 1-2 order-of-magnitude difference between the values of our metric for the SB and LB regimes. These results reinforce the view that the solutions obtained by a large-batch method defines points of larger sensitivity of the training function. In Appedix E, we describe approaches to attempt to remedy this generalization problem of LB methods. These approaches include data augmentation, conservative training and adversarial training. Our preliminary findings show that these approaches help reduce the generalization gap but still lead to relatively sharp minimizers and as such, do not completely remedy the problem.

Note that Metric 2.1 is closely related to the spectrum of ∇ 2 f (x). Assuming to be small enough, when A = I n , the value (4) relates to the largest eigenvalue of ∇ 2 f (x) and when A is randomly sampled it approximates the Ritz value of ∇ 2 f (x) projected onto the column-space of A. We conclude this section by noting that the sharp minimizers identified in our experiments do not resemble a cone, i.e., the function does not increase rapidly along all (or even most) directions. By sampling the loss function in a neighborhood of LB solutions, we observe that it rises steeply only along a small dimensional subspace (e.g. 5% of the whole space); on most other directions, the function is relatively flat.

3. SUCCESS OF SMALL-BATCH METHODS:

It is often reported that when increasing the batch size for a problem, there exists a threshold after which there is a deterioration in the quality of the model. This behavior can be observed for the F 2 and C 1 networks in Figure 4 . In both of these experiments, there is a batch size (≈ 15000 for F 2 and ≈ 500 for C 1 ) after which there is a large drop in testing accuracy. Notice also that the upward drift in value of the sharpness is considerably reduced around this threshold. Similar thresholds exist for the other networks in Table 1 .

Let us now consider the behavior of SB methods, which use noisy gradients in the step computation.

From the results reported in the previous section, it appears that noise in the gradient pushes the iterates out of the basin of attraction of sharp minimizers and encourages movement towards a flatter minimizer where noise will not cause exit from that basin. When the batch size is greater than the threshold mentioned above, the noise in the stochastic gradient is not sufficient to cause ejection from the initial basin leading to convergence to sharper a minimizer.

To explore that in more detail, consider the following experiment. We train the network for 100 epochs using ADAM with a batch size of 256, and retain the iterate after each epoch in memory. Using these 100 iterates as starting points we train the network using a LB method for 100 epochs and receive a 100 piggybacked (or warm-started) large-batch solutions. We plot in Figure 5 the testing accuracy and sharpness of these large-batch solutions, along with the testing accuracy of the small-batch iterates. Note that when warm-started with only a few initial epochs, the LB method does not yield a generalization improvement. The concomitant sharpness of the iterates also stays high. On the other hand, after certain number of epochs of warm-starting, the accuracy improves and sharpness of the large-batch iterates drop. This happens, apparently, when the SB method has ended its exploration phase and discovered a flat minimizer; the LB method is then able to converge towards it, leading to good testing accuracy.

It has been speculated that LB methods tend to be attracted to minimizers close to the starting point x 0 , whereas SB methods move away and locate minimizers that are farther away. Our numerical experiments support this view: we observed that the ratio of x s -x 0 2 and x -x 0 2 was in the range of 3-10.

In order to further illustrate the qualitative difference between the solutions obtained by SB and LB methods, we plot in Figure 6 our sharpness measure (4) against the loss function (cross entropy) for one random trial of the F 2 and C 1 networks. For larger values of the loss function, i.e., near the initial point, SB and LB method yield similar values of sharpness. As the loss function reduces, the sharpness of the iterates corresponding to the LB method rapidly increases, whereas for the SB method the sharpness stays relatively constant initially and then reduces, suggesting an exploration phase followed by convergence to a flat minimizer.

4. DISCUSSION AND CONCLUSION:

In this paper, we present numerical experiments that support the view that convergence to sharp minimizers gives rise to the poor generalization of large-batch methods for deep learning. To this end, we provide one-dimensional parametric plots and perturbation (sharpness) measures for a variety of deep learning architectures. In Appendix E, we describe our attempts to remedy the problem, including data augmentation, conservative training and robust optimization. Our preliminary investigation suggests that these strategies do not correct the problem; they improve the generalization of large-batch methods but still lead to relatively sharp minima. Another prospective remedy includes the use of dynamic sampling where the batch size is increased gradually as the iteration progresses (Byrd et al., 2012; Friedlander & Schmidt, 2012) . The potential viability of this approach is suggested by our warm-starting experiments (see Figure 5 ) wherein high testing accuracy is achieved using a large-batch method that is warm-start with a small-batch method.

Recently, a number of researchers have described interesting theoretical properties of the loss surface of deep neural networks; see e.g. (Choromanska et al., 2015; Soudry & Carmon, 2016; Lee et al., 2016) . Their work shows that, under certain regularity assumptions, the loss function of deep learning models is fraught with many local minimizers and that many of these minimizers correspond to a similar loss function value. Our results are in alignment these observations since, in our experiments, both sharp and flat minimizers have very similar loss function values. We do not know, however, if the theoretical models mentioned above provide information about the existence and density of sharp minimizers of the loss surface.

Our results suggest some questions: (a) can one prove that large-batch (LB) methods typically converge to sharp minimizers of deep learning training functions? (In this paper, we only provided some numerical evidence. A DETAILS ABOUT DATA SETS

We summarize the data sets used in our experiments in Table 5 . TIMIT is a speech recognition data set which is pre-processed using Kaldi (Povey et al., 2011) and trained using a fully-connected network. The rest of the data sets are used without any pre-processing. For this network, we use a 784-dimensional input layer followed by 5 batch-normalized (Ioffe & Szegedy, 2015) layers of 512 neurons each with ReLU activations. The output layer consists of 10 neurons with the softmax activation.

B.2 NETWORK F 2:

The network architecture for F 2 is similar to F 1 . We use a 360-dimensional input layer followed by 7 batch-normalized layers of 512 neurons with ReLU activation. The output layer consists of 1973 neurons with the softmax activation.

B.3 NETWORKS C 1 AND C 3:

The C 1 network is a modified version of the popular AlexNet configuration (Krizhevsky et al., 2012) . For simplicity, denote a stack of n convolution layers of a filters and a Kernel size of b × c with stride length of d as n× [a, b, c, d] . The C 1 configuration uses 2 sets of [64, 5, 5, 2]-MaxPool(3) followed by 2 dense layers of sizes (384, 192) and finally, an output layer of size 10. We use batchnormalization for all layers and ReLU activations. We also use Dropout (Srivastava et al., 2014) of 0.5 retention probability for the two dense layers. The configuration C 3 is identical to C 1 except it uses 100 softmax outputs instead of 10.

B.4 NETWORKS C 2 AND C 4:

The C 2 network is a modified version of the popular VGG configuration (Simonyan & Zisserman, 2014) . The C 3 network uses the configuration: 2×[64, 3, 3, 1], 2× [128, 3, 3, 1] , 3× [256, 3, 3, 1] , 3× [512, 3, 3, 1] , 3 × [512, 3, 3, 1] which a MaxPool(2) after each stack. This stack is followed by a 512dimensional dense layer and finally, a 10-dimensional output layer. The activation and properties of each layer is as in B.3. As is the case with C 3 and C 1 , the configuration C 4 is identical to C 2 except that it uses 100 softmax outputs instead of 10.

C PERFORMANCE MODEL:

As mentioned in Section 1, a training algorithm that operates in the large-batch regime without suffering from a generalization gap would have the ability to scale to much larger number of nodes than is currently possible. Such and algorithm might also improve training time through faster convergence. We present an idealized performance model that demonstrates our goal.

For LB method to be competitive with SB method, the LB method must (i) converge to minimizers that generalize well, and (ii) do it in a reasonably number of iterations, which we analyze here. Let I s and I be number of iterations required by SB and LB methods to reach the point of comparable test accuracy, respectively. Let B s and B be corresponding batch sizes and P be number of processors being used for training. Assume that P < B , and let f s (P ) be the parallel efficiency of the SB method. For simplicity, we assume that f (P ), the parallel efficiency of the LB method, is 1.0. In other words, we assume that the LB method is perfectly scalable due to use of a large batch size.

For LB to be faster than SB, we must have

I B P < I s B s P f s (P )

.

In other words, the ratio of iterations of LB to the iterations of SB should be

I I s < B s f s (P )B .

For example, if f s (P ) = 0.2 and B s /B = 0.1, the LB method must converge in at most half as many iterations as the SB method to see performance benefits. We refer the reader to (Das et al., 2016) for a more detailed model and a commentary on the effect of batch-size on the performance.

D CURVILINEAR PARAMETRIC PLOTS:

The parametric plots for the curvilinear path from x s to x , i.e., f (sin( απ 2 )x + cos( απ 2 )x s ) can be found in Figure 7 .

E ATTEMPTS TO IMPROVE LB METHODS:

In this section, we discuss a few strategies that aim to remedy the problem of poor generalization for large-batch methods. As in Section 2, we use 10% as the percentage batch-size for large-batch experiments and 256 for small-batch methods. For all experiments, we use ADAM as the optimizer irrespective of batch-size.

E.1 DATA AUGMENTATION:

Given that large-batch methods appear to be attracted to sharp minimizers, one can ask whether it is possible to modify the geometry of the loss function so that it is more benign to large-batch methods. The loss function depends both on the geometry of the objective function and to the size and properties of the training set. One approach we consider is data augmentation; see e.g. (Krizhevsky et al., 2012; Simonyan & Zisserman, 2014) . The application of this technique is domain specific but generally involves augmenting the data set through controlled modifications on the training data. For instance, in the case of image recognition, the training set can be augmented through translations, rotations, shearing and flipping of the training data. This technique leads to regularization of the network and has been employed for improving testing accuracy on several data sets.

In our experiments, we train the 4 image-based (convolutional) networks using aggressive data augmentation and present the results in Table 6 . For the augmentation, we use horizontal reflections, random rotations up to 10 • and random translation of up to 0.2 times the size of the image. It is evident from the table that, while the LB method achieves accuracy comparable to the SB method (also with training data augmented), the sharpness of the minima still exists, suggesting sensitivity to images contained in neither training or testing set. In this section, we exclude parametric plots and sharpness values for the SB method owing to space constraints and the similarity to those presented in Section 2.2. In (Li et al., 2014) , the authors argue that the convergence rate of SGD for the large-batch setting can be improved by obtaining iterates through the following proximal sub-problem.

x k+1 = arg min

x 1 |B k | i∈B k f i (x) + λ 2 x -x k 2 2

(5)

The motivation for this strategy is, in the context of large-batch methods, to better utilize a batch before moving onto the next one. The minimization problem is solved inexactly using 3-5 iterations of gradient descent, co-ordinate descent or L-BFGS. (Li et al., 2014) report that this not only improves the convergence rate of SGD but also leads to improved empirical performance on convex machine learning problems. The underlying idea of utilizing a batch is not specific to convex problems and we can apply the same framework for deep learning, however, without theoretical guarantees. Indeed, similar algorithms were proposed in (Zhang et al., 2015) and (Mobahi, 2016) for Deep Learning. The former placed emphasis on parallelization of small-batch SGD and asynchrony while the latter on a diffusion-continuation mechanism for training. The results using the conservative training approach are presented in Figure 7 . In all experiments, we solve the problem (5) using 3 iterations of ADAM and set the regularization parameter λ to be 10 -3 . Again, there is a statistically significant improvement in the testing accuracy of the large-batch method but it does not solve the problem of sensitivity.

E.3 ROBUST TRAINING:

A natural way of avoiding sharp minima is through robust optimization techniques. These methods attempt to optimize a worst-case cost as opposed to the nominal (or true) cost. Mathematically, given an > 0, these techniques solve the problem

EQUATION

Geometrically, classical (nominal) optimization attempts to locate the lowest point of a valley, while robust optimization attempts to lower an -disc down the loss surface. We refer an interested reader to (Bertsimas et al., 2010) , and the references therein, for a review of non-convex robust optimization. A direct application of this technique is, however, not feasible in our context since each iteration is prohibitively expensive because it involves solving a large-scale second-order conic program (SOCP).

:

The code to reproduce the parametric plot on exemplary networks can be found in our GitHub repository: https://github.com/keskarnitish/large-batch-training.

Figures and Tables:

Figure 1: A Conceptual Sketch of Flat and Sharp Minima. The Y-axis indicates value of the loss function and the X-axis the variables (parameters)

Figure 2: Training and testing accuracy for SB and LB methods as a function of epochs.

Figure 3: Parametric Plots -Linear (Left vertical axis corresponds to cross-entropy loss, f , and right vertical axis corresponds to classification accuracy; solid line indicates training data set and dashed line indicated testing data set); α = 0 corresponds to the SB minimizer and α = 1 to the LB minimizer.

Figure4: Testing Accuracy and Sharpness v/s Batch Size. The X-axis corresponds to the batch size used for training the network for 100 epochs, left Y-axis corresponds to the testing accuracy at the final iterate and right Y-axis corresponds to the sharpness of that iterate. We report sharpness for two values of : 10 -3 and 5 • 10 -4 .

Figure5: Warm-starting experiments. The upper figures report the testing accuracy of the SB method (blue line) and the testing accuracy of the warm started (piggybacked) LB method (red line), as a function of the number of epochs of the SB method. The lower figures plot the sharpness measure (4) for the solutions obtained by the piggybacked LB method v/s the number of warm-starting epochs of the SB method.

Figure 6: Sharpness v/s Cross Entropy Loss for SB and LB methods.

Figure 7: Parametric Plots -Curvilinear (Left vertical axis corresponds to cross-entropy loss, f , and right vertical axis corresponds to classification accuracy; solid line indicates training data set and dashed line indicated testing data set); α = 0 corresponds to the SB minimizer while α = 1 corresponds to the LB minimizer

Network Configurations

Table content: <table><tr><td colspan="2">Name Network Type</td><td colspan="2">Architecture Data set</td></tr><tr><td>F 1</td><td>Fully Connected</td><td>Section B.1</td><td>MNIST (LeCun et al., 1998a)</td></tr><tr><td>F 2</td><td>Fully Connected</td><td>Section B.2</td><td>TIMIT (Garofolo et al., 1993)</td></tr><tr><td>C 1</td><td colspan="2">(Shallow) Convolutional Section B.3</td><td>CIFAR-10 (Krizhevsky &amp; Hinton, 2009)</td></tr><tr><td>C 2</td><td>(Deep) Convolutional</td><td>Section B.4</td><td>CIFAR-10</td></tr><tr><td>C 3</td><td colspan="2">(Shallow) Convolutional Section B.3</td><td>CIFAR-100 (Krizhevsky &amp; Hinton, 2009)</td></tr><tr><td>C 4</td><td>(Deep) Convolutional</td><td>Section B.4</td><td>CIFAR-100</td></tr></table>



Table content: <table><tr><td/><td/><td>Training Accuracy</td><td colspan="2">Testing Accuracy</td></tr><tr><td colspan="2">Name SB</td><td>LB</td><td>SB</td><td>LB</td></tr><tr><td>F 1</td><td colspan="4">99.66% ± 0.05% 99.92% ± 0.01% 98.03% ± 0.07% 97.81% ± 0.07%</td></tr><tr><td>F 2</td><td colspan="3">99.99% ± 0.03% 98.35% ± 2.08% 64.02% ± 0.2%</td><td>59.45% ± 1.05%</td></tr><tr><td>C 1</td><td colspan="2">99.89% ± 0.02% 99.66% ± 0.2%</td><td colspan="2">80.04% ± 0.12% 77.26% ± 0.42%</td></tr><tr><td>C</td><td/><td/><td/></tr></table>

Sharpness of Minima in Full Space; is defined in (3). 28.58 ± 3.13 707.23 ± 43.04 7.08 ± 0.88 227.31 ± 23.23 C 2 8.68 ± 1.32 925.32 ± 38.29 2.07 ± 0.86 175.31 ± 18.28 C 3 29.85 ± 5.98 258.75 ± 8.96 8.56 ± 0.99 105.11 ± 13.22 C 4 12.83 ± 3.84 421.84 ± 36.97 4.07 ± 0.87 109.35 ± 16.57

Table content: <table><tr><td/><td>= 10 -3</td><td/><td colspan="2">= 5 • 10 -4</td></tr><tr><td/><td>SB</td><td>LB</td><td>SB</td><td>LB</td></tr><tr><td>F 1</td><td colspan="4">1.23 ± 0.83 205.14 ± 69.52 0.61 ± 0.27 42.90 ± 17.14</td></tr><tr><td>F 2</td><td colspan="3">1.39 ± 0.02 310.64 ± 38.46 0.90 ± 0.05</td><td>93.15 ± 6.81</td></tr><tr><td>C 1</td><td/><td/><td/></tr></table>

Sharpness of Minima in Random Subspaces of Dimension 100

Table content: <table><tr><td/><td/><td>= 10 -3</td><td colspan="2">= 5 • 10 -4</td></tr><tr><td/><td>SB</td><td>LB</td><td>SB</td><td>LB</td></tr><tr><td>F 1</td><td>0.11 ± 0.00</td><td>9.22 ± 0.56</td><td>0.05 ± 0.00</td><td>9.17 ± 0.14</td></tr><tr><td>F 2</td><td>0.29 ± 0.02</td><td>23.63 ± 0.54</td><td>0.05 ± 0.00</td><td>6.28 ± 0.19</td></tr><tr><td colspan="5">C 1 2.18 ± 0.23 137.25 ± 21.60 0.71 ± 0.15 29.50 ± 7.48</td></tr><tr><td colspan="2">C 2 0.95 ± 0.34</td><td>25.09 ± 2.61</td><td>0.31 ± 0.08</td><td>5.82 ± 0.52</td></tr><tr><td colspan="5">C 3 17.02 ± 2.20 236.03 ± 31.26 4.03 ± 1.45 86.96 ± 27.39</td></tr><tr><td colspan="2">C 4 6.05 ± 1.13</td><td colspan="3">72.99 ± 10.96 1.89 ± 0.33 19.85 ± 4.12</td></tr></table>

Sixin Zhang, Anna E Choromanska, and Yann LeCun. Deep learning with elastic averaging sgd. In Advances in Neural Information Processing Systems, pp. 685-693, 2015. Stephan Zheng, Yang Song, Thomas Leung, and Ian Goodfellow. Improving the robustness of deep neural networks via stability training. arXiv preprint arXiv:1604.04326, 2016.

Table content: <table/>



Table content: <table><tr><td/><td/><td/><td colspan="2">: Data Sets</td><td/></tr><tr><td/><td colspan="2"># Data Points</td><td/><td/><td/></tr><tr><td>Data Set</td><td>Train</td><td>Test</td><td colspan="3"># Features # Classes Reference</td></tr><tr><td>MNIST</td><td>60000</td><td>10000</td><td>28 × 28</td><td>10</td><td>(LeCun et al., 1998a;b)</td></tr><tr><td>TIMIT</td><td colspan="2">721329 310621</td><td>360</td><td>1973</td><td>(Garofolo et al., 1993)</td></tr><tr><td>CIFAR-10</td><td>50000</td><td>10000</td><td>32 × 32</td><td>10</td><td>(Krizhevsky &amp; Hinton, 2009)</td></tr><tr><td colspan="2">CIFAR-100 50000</td><td>10000</td><td>32 × 32</td><td>100</td><td>(Krizhevsky &amp; Hinton, 2009)</td></tr><tr><td colspan="4">B ARCHITECTURE OF NETWORKS</td><td/><td/></tr></table>

Effect of Data Augmentation C 1 83.63% ± 0.14% 82.50% ± 0.67% 231.77 ± 30.50 45.89 ± 3.83 C 2 89.82% ± 0.12% 90.26% ± 1.15% 468.65 ± 47.86 105.22 ± 19.57 C 3 54.55% ± 0.44% 53.03% ± 0.33% 103.68 ± 11.93 37.67 ± 3.46 C 4 63.05% ± 0.5% 65.88 ± 0.13% 271.06 ± 29.69 45.31 ± 5.93

Table content: <table><tr><td colspan="2">Testing Accuracy</td><td colspan="2">Sharpness (LB method)</td></tr><tr><td>Baseline (SB)</td><td>Augmented LB</td><td>= 10 -3</td><td>= 5 • 10 -4</td></tr></table>

Effect of Conservative Training • 10 -4 F 1 98.03% ± 0.07% 98.12% ± 0.01% 232.25 ± 63.81 46.02 ± 12.58 F 2 64.02% ± 0.2% 61.94% ± 1.10% 928.40 ± 51.63 190.77 ± 25.33 C 1 80.04% ± 0.12% 78.41% ± 0.22% 520.34 ± 34.91 171.19 ± 15.13 C 2 89.24% ± 0.05% 88.495% ± 0.63% 632.01 ± 208.01 108.88 ± 47.36 C 3 49.58% ± 0.39% 45.98% ± 0.

Table content: <table><tr><td colspan="2">Testing Accuracy</td><td colspan="2">Sharpness (LB method)</td></tr><tr><td>Baseline (SB)</td><td colspan="3">Conservative LB = 5 54% = 10 -3 337.92 ± 33.09 110.69 ± 3.88</td></tr><tr><td>C 4 63.08% ± 0.10%</td><td>62.51 ± 0.67</td><td>354.94 ± 20.23</td><td>68.76 ± 16.29</td></tr><tr><td>E.2 CONSERVATIVE TRAINING</td><td/><td/><td/></tr></table>

Peer review (iclr2017):
I think that the paper is quite interesting and useful. It might benefit from additional investigations, e.g., by adding some rescaled Gaussian noise to gradients during the LB regime one can get advantages of the SB regime.