Title: SUPER-CONVERGENCE: VERY FAST TRAINING OF RESIDUAL NETWORKS USING LARGE LEARNING RATES

Abstract: In this paper, we show a phenomenon, which we named "super-convergence", where residual networks can be trained using an order of magnitude fewer iterations than is used with standard training methods. The existence of super-convergence is relevant to understanding why deep networks generalize well. One of the key elements of super-convergence is training with cyclical learning rates and a large maximum learning rate. Furthermore, we present evidence that training with large learning rates improves performance by regularizing the network. In addition, we show that super-convergence provides a greater boost in performance relative to standard training when the amount of labeled training data is limited. We also derive a simplification of the Hessian Free optimization method to compute an estimate of the optimal learning rate. The architectures to replicate this work will be made available upon publication.

1. INTRODUCTION:

While deep neural networks have achieved amazing successes in a range of applications, understanding why stochastic gradient descent (SGD) works so well remains an open and active area of research. This paper provides unique empirical evidence that supports the theories in some papers but not others. Specifically, we show that, for certain datasets, residual network architectures (He et al., 2016) , and hyper-parameter values, using very large learning rates with the cyclical learning rate (CLR) method (Smith, 2015; 2017) can speed up training by an order of magnitude. Analogous to the phenomenon of super-conductivity that only happens in limited circumstances and provides theoretical insights of materials, we named this phenomenon "super-convergence." While super-convergence might be of some practical value, the primary purpose of this paper is to provide empirical support and theoretical insights to the active discussions in the literature on SGD and understanding generalization.

Figure 1a provides a comparison of test accuracies from a super-convergence example and the result of a typical (piecewise constant) training regime for Cifar-10, both using a 56 layer residual network architecture. Piecewise constant training reaches a peak accuracy of 91.2% after approximately 80,000 iterations, while the super-convergence method reaches a higher accuracy (92.4%) after only 10,000 iterations. Figure 1b shows the results for a range of CLR stepsize values, where training lasted only one cycle. This modified learning rate schedule achieves a higher final test accuracy (92.1%) than typical training (91.2%) after only 6,000 iterations. In addition, as the total number of iterations increases from 2,000 to 20,000, the final accuracy improves from 89.7% to 92.7%.

The contributions of this paper are:

1. We demonstrate a new training phenomenon and systematically investigate it. 2. We show evidence that large learning rates regularize the trained network and hypothesize that this allows SGD to find a flat local minima that generalizes well. 3. We derive a simplification of the second order, Hessian-free optimization method to estimate optimal learning rates which demonstrates that large learning rates find wide, flat minima. 4. We demonstrate that the effects of super-convergence are increasingly dramatic when less labeled training data is available. 

2. BACKGROUND:

In this paper, when we refer to a typical, standard, or a piecewise-constant training regime, it means the practice of using a global learning rate, (≈ 0.1), for many epochs, until the test accuracy plateaus, and then continuing to train with a lower learning rate decreased by a factor of 0.1. This process of reducing the learning rate and continuing to train is often repeated two or three times.

There exists extensive literature on stochastic gradient descent (SGD) (see Goodfellow et al. (2016) and Bottou (2012) ) which is relevant to this work. Also, there exists a significant amount of literature on the loss function topology of deep networks (see Chaudhari et al. (2016) for a review of the literature). This paper contains a discussion of the loss function topology which follows from the work of Goodfellow et al. (2014) on characterizing the landscape. Our use of large learning rate values is in contrast to suggestions in the literature of a maximum learning rate value Bottou et al. (2016) . Methods for adaptive learning rates have also been an active area of research. This paper uses a simplification of the second order Hessian-Free optimization (Martens, 2010) to estimate optimal values for the learning rate. In addition, we utilize some of the techniques described in Schaul et al. (2013) and Gulcehre et al. (2017) . Also, we show that adaptive learning rate methods such as Nesterov momentum (Sutskever et al., 2013; Nesterov, 1983) , AdaDelta (Duchi et al., 2011 ), AdaGrad (Zeiler, 2012) , and Adam (Kingma & Ba, 2014) do not use sufficiently large learning rates when they are effective nor do they lead to super-convergence without using CLR. A warmup learning rate strategy (He et al., 2016; Goyal et al., 2017) could be considered a discretization of CLR, which was also recently suggested in (Jastrzębski et al., 2017) . We note that Loshchilov & Hutter (2016) subsequently proposed a similar method to CLR, which they call SGDR. The SGDR method uses a sawtooth pattern with a cosine followed by a jump back up to the original value. Experiments show that it is not possible to observe the super-convergence phenomenon when using their pattern.

Our work is intertwined with several active lines of research in the deep learning research community, including a lively discussion on stochastic gradient descent (SGD) and understanding why solutions generalize so well, research on SGD and the importance of noise for generalization, and the generalization gap between small and large mini-batches. We defer our discussion of these lines of research to Section 7 where we can compare to our empirical results and theoretical insights.

3. SUPER-CONVERGENCE:

In this work, we use cyclical learning rates (CLR) and the learning rate range test (LR range test) which were first suggested by Smith (2015) and later updated in Smith (2017) . To use CLR, one specifies minimum and maximum learning rate boundaries and a stepsize. The stepsize is the number of iterations used for each step and a cycle consists of two such steps -one in which the learning rate increases and the other in which it decreases. Smith (2015) tested numerous ways to vary the learning rate between the two boundary values, found them to be equivalent and therefore recommended the simplest, which is letting the learning rate change linearly (similarly Jastrzębski et al. (2017) suggest discrete jumps). Please note that cyclical noise (via learning rates, batch size, etc.) is obtained by combining curriculum learning (Bengio et al., 2009) and simulated annealing (Aarts & Korst, 1988) , both of which have a long history of use in deep learning. The LR range test can be used to determine if super-convergence is possible for an architecture. In the LR range test, training starts with a zero or very small learning rate which is slowly increased linearly throughout a pre-training run. This provides information on how well the network can be trained over a range of learning rates. Figure 2a shows a typical curve from a LR range test, where the test accuracy has a distinct peak.1 When starting with a small learning rate, the network begins to converge and, as the learning rate increases, it eventually becomes too large and causes the training/test accuracy to decrease. The learning rate at this peak is the largest value to use as the maximum learning rate bound when using CLR. The minimum learning rate can be chosen by dividing the maximum by a factor of 3 or 4. The optimal initial learning rate for a typical (piecewise constant) training regime usually falls between these minimum and maximum values.

If one runs the LR range test for Cifar-10 on a 56 layer residual networks, one obtains the curves shown in Figure 2b . Please note that learning rate values up to 3.0 were tested, which is at least an order of magnitude lager than typical values of the learning rate. The test accuracy remains consistently high over this unusually long range of large learning rates. This unusual behavior motivated our experimentation with much higher learning rates, and we believe that such behavior during a LR range test is indicative of potential for super-convergence. The three curves in this figure are for runs with a maximum number of iterations of 5,000, 20,000, and 100,000, where the test accuracy remains consistently high for a long range of large learning rate values. In early iterations, the learning rate must be small in order for the training to make progress in an appropriate direction. The Figure also

shows that significant progress is made in those early iterations. However, as the slope decreases so does the amount of progress per iteration and little improvement occurs over the bulk of the iterations. Figure 3b shows a close up of the final parts of the training where the solution maneuvers through a valley to the local minimum within a trough.

Cyclical learning rates are well suited for training when the loss topology takes this form. The learning rate initially starts small to allow convergence to begin. As the network traverses the flat valley, the learning rate is large, allowing for faster progress through the valley. In the final stages of the training, when the training needs to settle into the local minimum (as seen in Figure 3b ), the learning rate is once again reduced to its original small value. 

4. LARGE LEARNING RATE REGULARIZATION:

In our experiments, we only observed super-convergence with large learning rates. The LR range test reveals evidence of regularization through results shown in Figure 4a . There is additional evidence that large learning rates are regularizing the training. As shown in Figure 1a , the final test accuracy results from a super-convergence training is demonstrably better than the accuracy results from a typical training method. In the literature, this type of improvement in the final training accuracy is often taken as evidence of regularization. Furthermore, others show that large learning rates leads to larger gradient noise, which leads to better generalization (i.e., Jastrzębski et al. (2017) ; Smith et al. (2017) ).

5. ESTIMATING OPTIMAL LEARNING RATES:

Gradient or steepest descent is an optimization method that uses the slope as computed by the derivative to move in the direction of greatest negative gradient to iteratively update a variable. That is, given an initial point x 0 , gradient descent proposes the next point to be:

x = x 0 - x f (x) (1)

where is the step size or learning rate . If we denote the parameters in a neural network (i.e., weights) as θ ∈ R N and f (θ) is the loss function, we can apply gradient descent to learn the weights of a network; i.e., with input x, a solution y, and non-linearity σ:

EQUATION

where W l ∈ θ are the weights for layer l and b l ∈ θ are biases for layer l.

The Hessian-free optimization method (Martens, 2010) suggests a second order solution that utilizes the slope information contained in the second derivative (i.e., the derivative of the gradient θ f (θ)).

From Martens (2010) , the main idea of the second order Newton's method is that the loss function can be locally approximated by the quadratic as:

f (θ) ≈ f (θ 0 ) + (θ -θ 0 ) T θ f (θ 0 ) + 1 2 (θ -θ 0 ) T H(θ -θ 0 ) (3)

where H is the Hessian, or the second derivative matrix of f (θ 0 ). Writing Equation 1 to update the parameters at iteration i as:

EQUATION

) allows Equation 3 to be re-written as:

f (θ i - θ f (θ i )) ≈ f (θ i ) + (θ i+1 -θ i ) T θ f (θ i ) + 1 2 (θ i+1 -θ i ) T H(θ i+1 -θ i ) (5)

In general it is not feasible to compute the Hessian matrix, which has Ω(N 2 ) elements, where N is the number of parameters in the network, but it is unnecessary to compute the full Hessian. The Hessian expresses the curvature in all directions in a high dimensional space, but the only relevant curvature direction is in the direction of steepest descent that SGD will traverse. This concept is contained within Hessian-free optimization, as Martens (2010) suggests a finite difference approach for obtaining an estimate of the Hessian from two gradients:

H(θ) = lim δ→0 f (θ + δ) -f (θ) δ (6)

where δ should be in the direction of the steepest descent. The AdaSecant method (Gulcehre et al., 2014; 2017) builds an adaptive learning rate method based on this finite difference approximation as:

* ≈ θ i+1 -θ i f (θ i+1 ) -f (θ i ) (7)

where * represents the optimal learning rate for each of the neurons. Utilizing Equation 4, we rewrite Equation 7 in terms of the differences between the weights from three sequential iterations as:

* = θ i+1 -θ i 2θ i+1 -θ i -θ i+2 (8)

where on the right hand side is the learning rate value actually used in the calculations to update the weights. Equation 8 is an expression for an adaptive learning rate for each weight update. We borrow the method in Schaul et al. (2013) to obtain an estimate of the global learning rate from the weight specific rates by summing over the numerator and denominator, with one minor difference. In Schaul et al. (2013) their expression is squared, leading to positive values -therefore we sum the absolute values of each quantity to maintain positivity (using the square root of the sum of squares of the numerator and denominator of Equation 8 leads to similar results).

For illustrative purposes, we computed the optimal learning rates from the weights of every iteration using Equation 8for two runs: first when the learning rate was a constant value of 0.1 and second with CLR in the range of 0.1 -3 with a stepsize of 5,000 iterations. Since the computed learning rate exhibited rapid variations, we computed a moving average of the estimated learning rate as LR = α * + (1 -α)LR with α = 0.1 and the results are shown in Figure 5a for the first 300 iterations. This curve qualitatively shows that the optimal learning rates should be in the range of 2 to 4 for this architecture. In Figure 5b , we used the weights as computed every 10 iterations and ran the BN MAF is the value use for the moving_average_f raction parameter with batch normalization. PC-LR is a standard piecewise constant learning rate policy described in Section 2 with an initial learning rate of 0.35. learning rate estimation for 10,000 iterations. An interesting divergence happens here: when keeping the learning rate constant, the learning rate estimate initially spikes to a value of about 3 but then drops down near 0.2. On the other hand, the learning rate estimate using CLR remains high until the end where it settles down to a value of about 0.5. The large learning rates indicated by these Figures is caused by small values of our Hessian approximation and small values of the Hessian implies that SGD is finding flat and wide local minima.

In this paper we do not perform a full evaluation of the effectiveness of this technique as it is tangential to the theme of this work. We only use this method here to demonstrate that training with large learning rates are indicated by this approximation. We leave a full assessment and tests of this method to estimate optimal adaptive learning rates as future work.

(a) Estimated learning rates computed from the weights at every iteration for 300 iterations.

(b) Estimated learning rates computed from the weights at every 10 iterations for 10,000 iterations. 

6. EXPERIMENTS AND ANALYSIS:

This section highlights a few of our more significant experiments. Additional experiments and details of our architecture are illustrated in the Supplemental Section. Figure 6a provides a comparison of super-convergence with a reduced number of training samples.

When the amount of training data is limited, the gap in performance between the result of standard training and super-convergence increases. With a piecewise constant learning rate schedule the training encounters difficulties and diverges along the way. On the other hand, a network trained with specific CLR parameters exhibits super-convergence and trains without difficulties. The highest accuracies attained using standard learning rate schedules are listed in Table 1 and super-convergence test accuracy is 1.2%, 5.2%, and 9.2% better for 50,000, 20,000, and 10,000 training cases, respectively. Hence, super-convergence becomes more beneficial when training data is more limited.

We also ran experiments with Resnets with a number of layers in the range of 20 to 110 layers; that is, we ran experiments on residual networks with l layers, where l = 20 + 9n; for n = 0, 1, ...10. Figure 6b illustrates the results for Resnet-20 and Resnet-110, for both a typical (piecewise constant) training regime with a standard initial learning rate of 0.35 and for CLR with a stepsize of 10,000 iterations. For this entire range of architecture depths, super-convergence was possible. The accuracy increase due to super-convergence appears to be greater for the shallower architectures (Resnet-20: CLR 90.4% versus piecewise constant LR schedule 88.6%) than for the deeper architectures (Resnet-110: CLR 92.1% versus piecewise constant LR schedule 91.0%).

Recently, there are discussions in the deep learning literature on the effects of larger batch size and the generalization gap (Keskar et al., 2016; Jastrzębski et al., 2017; Chaudhari & Soatto, 2017; Hoffer et al., 2017) . Hence, we investigated the effects total mini-batch size3 used in super-convergence training and found a small improvement in performance with larger batch sizes, as can be seen in Figure 7a . In addition, Figure 7b shows that the generalization gap (the difference between the training and test accuracies) are approximately equivalent for small and large mini-batch sizes. This result differs than results reported elsewhere (Keskar et al., 2016) and illustrates that a benefit of training with large learning rates is the ability to use large batch sizes (Goyal et al., 2017) . We ran a series of experiments with a variety of ranges for CLR. Table 1 shows the results for maximum learning rate bounds from 1.0 to 3.5. These experiments show that a maximum learning rate of approximately 3 performed well. Although our experiments used a minimum learning rate of 0.1, training with a single cycle of CLR and minimum learning rate of 0 performs quite well.

We tested the effect of batch normalization on super-convergence. Initially, we found that having use_global_stats : true in the test phase prevents super-convergence. However, we realized this was due to using the default value of moving_average_f raction = 0.999 that is only appropriate for the typical, long training times. However, when a network is trained very quickly, as with super-convergence, a value of 0.999 does not update the accumulated global statistics quickly enough. Hence, we found that the smaller values listed in Table 1 (column BN MAF) were more appropriate.

We ran a variety of other experiments in which super-convergence continued to occur. We ran experiments with Cifar-100 that demonstrated the super-convergence phenomenon. We found that adding dropout to the architecture still permitted super-convergence and improved the results a small amount. Our experiments also investigated whether adaptive learning rate methods in a piecewise constant training regime would learn to adaptively use a large learning rate to improve performance. We tested Nesterov momentum (Sutskever et al., 2013; Nesterov, 1983) , AdaDelta (Duchi et al., 2011 ), AdaGrad (Zeiler, 2012) , and Adam (Kingma & Ba, 2014) on Cifar-10 with the Resnet-56 architecture but none of these methods speed up the training process in a similar fashion to super-convergence. This is important because it is indicative of the failing in the theory behind these approaches. We also tried these methods with CLR with appropriate bounds and found super-convergence happens for AdaDelta, AdaGrad and Nesterov momentum but we were unable to achieve super-convergence with Adam. In addition, we investigated the effect of momentum. We ran experiments with momentum in the range of 0.80 -0.95 and found momentums between 0.8 -0.9 yield a higher final test accuracy.

We illustrate the results of these experiments in the supplemental materials.

7. DISCUSSION:

There is substantial discussion in the literature on stochastic gradient descent (SGD) and understanding why solutions generalize so well (i.e, Chaudhari et al. ( 2016 2017)). Super-convergence provides empirical evidence that supports some theories, contradicts some others, and points to the need for further theoretical understanding. We hope the response to super-convergence is similar to the reaction to the initial report of network memorization (Zhang et al., 2016) , which sparked an active discussion within the deep learning research community on better understanding of the factors in SGD leading to solutions that generalize well (i.e., (Arpit et al., 2017) ).

Our work impacts the line of research on SGD and the importance of noise for generalization. In this paper we focused on the use of CLR with very large learning rates, which adds noise in the middle part of training. Recently, Jastrzębski et al. (2017) stated that higher levels of noise lead SGD to solutions with better generalization. Specifically, they showed that the ratio of the learning rate to the batch size, along with the variance of the gradients, controlled the width of the local minima found by SGD. Independently, Chaudhari & Soatto (2017) show that SGD performs regularization by causing SGD to be out of equilibrium, which is crucial to obtain good generalization performance, and derive that the ratio of the learning rate to batch size alone controls the entropic regularization term. They also state that data augmentation increases the diversity of SGD's gradients, leading to better generalization. These two papers provide theoretical support for the super-convergence phenomenon. In addition there are several other papers in the literature which state that wide, flat local minima produce solutions that generalize better than sharp minima (Hochreiter & Schmidhuber, 1997; Keskar et al., 2016; Wu et al., 2017) . Our super-convergence results align with these results.

In addition, there are several recent papers on the generalization gap between small and large minibatches and the relationship between gradient noise, learning rate, and batch size. Our results here supplements this other work by illustrating the possibility of time varying high noise levels during training. As mentioned above, Jastrzębski et al. (2017) showed that SGD noise is proportional to the learning rate, the variance of the loss gradients, divided by the batch size. Similarly Smith and Le (Smith & Le, 2017) derived the noise scale as g ≈ N/B(1 -m), where g is the gradient noise, the learning rate, N the number of training samples, and m is the momentum coefficient. Furthermore, Smith et al. (2017) showed an equivalence of increasing batch sizes instead of a decreasing learning rate schedule. Importantly, these authors demonstrated that the noise scale g is relevant to training and not the learning rate or batch size. Keskar et al. (2016) study the generalization gap between small and large mini-batches, stating that small mini-batch sizes lead to wide, flat minima and large batch sizes lead to sharp minima. They also suggest a batch size warm start for the first few epochs, then using a large batch size, which amounts to training with large gradient noise for a few epochs and then removing it. Our results contradicts this suggestion as we found it preferable to start training with little or no noise and let it increase (i.e., curriculum training), reach a noise peak, and reduce the noise level in the last part of training (i.e., simulated annealing). Goyal et al. (2017) use a very large mini-batch size of up to 8,192 and adjust the learning rate linearly with the batch size. They also suggest a gradual warmup of the learning rate, which is a discretized version of CLR and matches our experience with an increasing learning rate. They make a point relevant to adjusting the batch size; if the network uses batch normalization, different mini-batch sizes leads to different statistics, which must be handled. Hoffer et al. (2017) made a similar point about batch norm and suggested using ghost statistics. Also, Hoffer et al. (2017) show that longer training lengths is a form of regularization that improves generalization. On the other hand, our results show a different form of regularization that comes from training with very large learning rates, which permits much shorter training lengths.

Furthermore, this paper points the way towards new research directions, such as the following three:

1. Characterizing "good noise" that improves the trained network's ability to generalize versus "bad noise" that interferes with finding a good solution (i.e., Zhang et al. (2016) ). We find that there is a lack of a unified framework for treating SGD noise/diversity, such as architectural noise (e.g., dropout (Srivastava et al., 2014) , dropconnect (Wan et al., 2013) ), noise from hyper-parameter settings (e.g., learning rate, mini-batch size), adding gradient noise, adding noise to weights (Fortunato et al., 2017) , and input diversity (e.g., data augmentation, noise).

Gradient diversity has been shown to lead to flatter local minimum and better generalization.

A unified framework should resolve conflicting claims in the literature on the value of each of these, such as for architectural noise (Srivastava et al. (2014) versus Hoffer et al. (2017) ). Furthermore, many papers study each of these factors independently and by focusing on the trees, one might miss the forest. 2. Time dependent application of good noise during training: As described above, combining curriculum learning with simulated annealing leads to cyclical application. To the best of our knowledge, this has only been applied sporadically in a few methods such as CLR (Smith, 2017) or cyclical batch sizes (Jastrzębski et al., 2017) but these all fall under a single umbrella of time dependent gradient diversity (i.e., noise). Also, one might learn an optimal noise level while training the network. 3. Discovering new ways to stabilize optimization (i.e., SGD) with large noise levels: Our evidence indicates that normalization (and batch normalization in particular) is the catalyst enabling super-convergence in the face of destabilizing noise from the large learning rates. from a piecewise constant training regime. Furthermore, the final accuracy for the super-convergence curve is 68.6%, while the accuracy for the piecewise constant method is 59.8%, which is an 8.8% improvement.

As discussed in the main text, we tested various adaptive learning rate methods with Resnet-56 training on Cifar-10 to determine if they are capable of recognizing the need for using very large learning rates. Figure 9a shows the results of this training for Nesterov momentum (Sutskever et al., 2013; Nesterov, 1983) , AdaDelta (Duchi et al., 2011) , AdaGrad (Zeiler, 2012), and Adam (Kingma & Ba, 2014). We found no sign that any of these methods discovered the utility of large learning rates nor any indication of super-convergence-like behavior. We also ran CLR with these adaptive learning methods and found that Nesterov, AdaDelta, and AdaGrad allowed super-convergence to occur, but we were unable to create this phenomenon with Adam. For example, Figure 9b shows a comparison of super-convergence to a piecewise constant training regime with the Nesterov momentum method. Here super-convergence yields a final test accuracy after 10,000 iterations of 92.1% while the piecewise constant training regime at iteration 80,000 has an accuracy of 90.9%.

Figure 10 shows a comparison of runs of the super-convergence phenomenon, both with and without dropout. The LR range test with and without dropout is shown in Figure 10a . Figure 10b shows the results of training for 10,000 iterations. In both cases, the dropout ratio was set to 0.2 and the Figure shows a small improvement with dropout. We also ran with other values for the dropout ratio and consistently saw similar improvements.

The effect of mini-batch size is discussed in the main paper but here we present a table containing the final accuracies of super-convergence training with various mini-batch sizes. One can see in Table 5 the final test accuracy results and this table shows that the larger the mini-batch size, the better the final accuracy, which differs from the results shown in the literature4 . Most of results reported in this paper are with a total mini-batch size of 1,000.

In addition, we ran experiments on Resnet-56 on Cifar-10 with modified values for momentum and weight decay to determine if they might hinder the super-convergence phenomenon. Figure 11a shows the results for momentum set to 0.8, 0.85, 0.9, and 0.95 and the final test accuracies are listed in Table 6 . These results indicate only a small change in the results, with a setting of 0.9 being a bit better than the other values. In Figure 11b are the results for weight decay values of 10 -3 , 10 -4 , 10 -5 , and 10 -6 . In this case, a weight decay value of 10 -3 prevents super-convergence, while the smaller values do not. This Figure also shows that a weight decay value of 10 -4 performs well. We did not observe super-convergence in many experiments, however this does not necessarily mean that the super-convergence phenomenon is impossible in those cases but only that we did not find a way to induce it. Figure 12a shows an example of the LR range test with a bottleneck version of Resnet-56. There appears to be a peak from a learning rate of approximately 0.15 but even more pronounced are the large swings in the accuracy from around 0.1 and onward. These oscillations and the decay of accuracy at learning rates greater than 0.1 indicate that one cannot train this network with large learning rates and our experiments confirmed this to be true. We also tried Cifar-10 with the ResNeXt architecture and were unable to produce the super-convergence phenomenon. Figure 12b shows an example of the LR range test with ResNeXt. This result is similar to Figure 12a . There appears to be a peak and large swings in the accuracy from a learning rate of approximately 0.1. These oscillations and the decay of accuracy at learning rates greater than 0.1 are indications that one cannot train this network with large learning rates.

Our experiments with Densenets and all our experiments with the Imagenet dataset for a wide variety of architectures failed to produce the super-convergence phenomenon. Here we will list some of the experiments we tried that failed to produce the super-convergence behavior. Super-convergence did 

:

Figure reproduced from Smith (2017) with permission.

Figure reproduced from Goodfellow et al. (2014) with permission.

Most of our reported results are with a total mini-batch size of 1,000 and we primarily used 8 GPUs and split the total mini-batch size 8 ways over the GPUs.

GPU memory limitations prevented our testing a total batch size greater than 1,530

Figures and Tables:

(a) Comparison of test accuracies of superconvergence example to a typical (piecewise constant) training regime. (b) Comparison of test accuracies of superconvergence for a range of stepsizes.

Figure 1: Examples of super-convergence with Resnet-56 on Cifar-10.

(a) Typical learning rate range test result where there is a peak to indicate max_lr. (b) Learning rate range test result with the Resnet-56 architecture on Cifar-10.

Figure 2: Comparison of learning rate range test results.

Figure 3a provides an example of transversing the loss function topology. 2 This figure helps give an intuitive understanding of how super-convergence happens. The blue line in the Figure represents the trajectory of the training while converging and the x's indicate the location of the solution at each iteration and indicates the progress made during the training.In early iterations, the learning rate must be small in order for the training to make progress in an appropriate direction. The Figurealsoshows that significant progress is made in those early iterations. However, as the slope decreases so does the amount of progress per iteration and little improvement occurs over the bulk of the iterations. Figure3bshows a close up of the final parts of the training where the solution maneuvers through a valley to the local minimum within a trough.

(a) Visualization of how training transverses a loss function topology. (b) A close up of the end of the training for the example in Figure 3a.

Figure 3: The 3-D visualizations from Goodfellow et al. (2014). The z axis represents the loss potential.

Figure 4: Evidence of regularization with large learning rates: decreasing generalization error as the learning rate increases from 0.3 to 1.5.

Figure 4a shows an increasing training loss and decreasing test loss while the learning rate increases from approximately 0.2 to 2.0 when training with the Cifar-10 dataset and a Resnet-56 architecture, which implies that regularization is occurring while training with these large learning rates. Similarly, Figure 4b presents the training and test loss curves for Resnet-20 and Resnet-110, where one can see the same decreasing generalization error. In addition, we ran the LR range test on residual networks with l layers, where l = 20 + 9n; for n = 0, 1, ...10 and obtained similar results.

Figure 5: Estimated learning rate from the simplified Hessian-free optimization (see text for additional information).

(a) Comparison of test accuracies for Cifar-10 with limited training samples. (b) Comparison of test accuracies for Resnet-20 and Resnet-110.

Figure 6: Comparisons of super-convergence to typical training outcome with piecewise constant learning rate schedule.

(a) Comparison of test accuracies for Cifar-10, Resnet-56 while varying the total batch sizes. (b) Comparison of the generalization gap (trainingtest accuracy) for a small and a large batch size.

Figure 7: Comparisons of super-convergence to over a range of batch sizes. These results show that a large batch size is more effective than a small batch size for super-convergence training.We ran a series of experiments with a variety of ranges for CLR. Table1shows the results for maximum learning rate bounds from 1.0 to 3.5. These experiments show that a maximum learning

); Chaudhari & Soatto (2017); Im et al. (2016); Jastrzębski et al. (2017); Smith & Le (2017); Kawaguchi et al. (

(a) Learning rate range test. (b) Comparison of test accuracies for superconvergence to a piecewise constant training regime.

Figure 8: Comparisons for Cifar-100, Resnet-56 of super-convergence to typical (piecewise constant) training regime.

(a) Comparison of test accuracies with various adaptive learning rate methods. (b) Comparison of test accuracies With Nesterov method.

Figure 9: Comparisons for Cifar-10, Resnet-56 of super-convergence to piecewise constant training regime.

(a) Comparison of LR range test for superconvergence with and without dropout. (b) Comparison of test accuracies for superconvergence with and without dropout.

Figure 10: Comparisons of super-convergence with and without dropout (dropout ratio=0.2).

Figure 11: Comparisons for Cifar-10, Resnet-56 of super-convergence to typical training regime.

(a) Learning rate range test result with the bottleneck Resnet-56 architecture. (b) Learning rate range test result with the ResNeXt-56 architecture.

Figure 12: Comparison of learning rate range test results on Cifar-10 with alternate architectures. Mini-batch size Accuracy (%) 1536 92.1 1000 92.4 512 91.7 256 89.5



Table content: <table><tr><td colspan="5"># training samples Policy (Range) BN MAF Total Iterations Accuracy (%)</td></tr><tr><td>40,000</td><td>PC-LR=0.35</td><td>0.999</td><td>80,000</td><td>89.1</td></tr><tr><td>40,000</td><td>CLR (0.1-3)</td><td>0.95</td><td>10,000</td><td>91.1</td></tr><tr><td>30,000</td><td>PC-LR=0.35</td><td>0.999</td><td>80,000</td><td>85.7</td></tr><tr><td>30,000</td><td>CLR (0.1-3)</td><td>0.95</td><td>10,000</td><td>89.6</td></tr><tr><td>20,000</td><td>PC-LR=0.35</td><td>0.999</td><td>80,000</td><td>82.7</td></tr><tr><td>20,000</td><td>CLR (0.1-3)</td><td>0.95</td><td>10,000</td><td>87.9</td></tr><tr><td>10,000</td><td>PC-LR=0.35</td><td>0.999</td><td>80,000</td><td>71.4</td></tr><tr><td>10,000</td><td>CLR (0.1-3)</td><td>0.95</td><td>10,000</td><td>80.6</td></tr><tr><td>50,000</td><td>CLR (0.1-3.5)</td><td>0.95</td><td>10,000</td><td>92.1</td></tr><tr><td>50,000</td><td>CLR (0.1-3)</td><td>0.95</td><td>10,000</td><td>92.4</td></tr><tr><td>50,000</td><td>CLR (0.1-2.5)</td><td>0.95</td><td>10,000</td><td>92.3</td></tr><tr><td>50,000</td><td>CLR (0.1-2)</td><td>0.95</td><td>10,000</td><td>91.7</td></tr><tr><td>50,000</td><td>CLR (0.1-1.5)</td><td>0.95</td><td>10,000</td><td>90.9</td></tr><tr><td>50,000</td><td>CLR (0.1-1)</td><td>0.95</td><td>10,000</td><td>91.3</td></tr><tr><td>50,000</td><td>CLR (0.1-3)</td><td>0.97</td><td>20,000</td><td>92.7</td></tr><tr><td>50,000</td><td>CLR (0.1-3)</td><td>0.95</td><td>10,000</td><td>92.4</td></tr><tr><td>50,000</td><td>CLR (0.1-3)</td><td>0.93</td><td>8,000</td><td>91.7</td></tr><tr><td>50,000</td><td>CLR (0.1-3)</td><td>0.90</td><td>6,000</td><td>92.1</td></tr><tr><td>50,000</td><td>CLR (0.1-3)</td><td>0.85</td><td>4,000</td><td>91.1</td></tr><tr><td>50,000</td><td>CLR (0.1-3)</td><td>0.80</td><td>2,000</td><td>89.7</td></tr></table>

Normalization methods (batch norm, layer normalization Ba et al. (2016), streaming normalizationLiao et al. (2016)) and new techniques (i.e., cyclical gradient clipping) to stabilize training need further investigation to discover better ways to keep SGD stable in the presence of enormous good noise.Classical physics was insufficient for explaining super-conductivity when it was discovered and pointed to the need for new theories, such as quantum mechanics. Similarly, super-convergence indicates a need for new theories of SGD and generalization. Overall architecture for ResNet-56.

Table content: <table><tr><td>Layer/Block Type</td><td>Parameters</td></tr><tr><td>Conv Layer</td><td>padding = 1</td></tr><tr><td/><td>kernel = 3x3</td></tr><tr><td/><td>stride = 2</td></tr><tr><td/><td>channels = 16</td></tr><tr><td>Batch Norm</td><td>moving_average_fraction=0.95</td></tr><tr><td>ReLU</td><td>-</td></tr><tr><td>ResNet Standard Block x9</td><td>numChannels = 16</td></tr><tr><td>ResNet Downsample Block</td><td>numChannels = 16</td></tr><tr><td>ResNet Standard Block x8</td><td>numChannels = 32</td></tr><tr><td>ResNet Downsample Block</td><td>numChannels = 32</td></tr><tr><td>ResNet Standard Block x8</td><td>numChannels = 64</td></tr><tr><td>Average Pooling</td><td>padding = 0</td></tr><tr><td/><td>kernel = 8x8</td></tr><tr><td/><td>stride = 1</td></tr><tr><td>Fully Connected Layer</td><td>-</td></tr></table>

Comparison of accuracy results for various momentum values with Resnet-56 on Cifar-10 using CLR=0.1-3 and stepsize=5,000.

Table content: <table><tr><td colspan="2">Momentum Accuracy (%)</td></tr><tr><td>0.80</td><td>92.1</td></tr><tr><td>0.85</td><td>91.9</td></tr><tr><td>0.90</td><td>92.4</td></tr><tr><td>0.95</td><td>90.7</td></tr></table>

Peer review (iclr2018):
 In this paper, the authors analyze training of residual networks using large cyclic learning rates (CLR). The authors demonstrate (a) fast convergence with cyclic learning rates and (b) evidence of large learning rates acting as regularization which improves performance on test sets – this is called “super-convergence”. However, both these effects are only shown on a specific dataset, architecture, learning algorithm and hyper parameter setting. Some specific comments by sections: 2. Related Work: This section loosely mentions other related works on SGD, topology of loss function and adaptive learning rates. The authors mention Loshchilov & Hutter in next section but do not compare it to their work. The authors do not discuss a somewhat contradictory claim from NIPS 2017 (as pointed out in the public comment): http://papers.nips.cc/paper/6770-train-longer-generalize-better-closing-the-generalization-gap-in-large-batch-training-of-neural-networks.pdf 3. Super-convergence: This is a well explained section where the authors describe the LR range test and how it can be used to understand potential for super-convergence for any architecture. The authors also provide sufficient intuition for super-convergence. Since CLRs were already proposed by Smith (2015), the originality of this work would be specifically tied to their application to residual units. It would be interesting to see a qualitative analysis on how the residual error is impacting super-convergence. 4. Regularization: While Fig 4 demonstrates the regularization property, the reference to Fig 1a with better test error compared to typical training methods could simply be a result of slower convergence of typical training methods. 5. Optimal LRs: Fig.5b shows results for 1000 iterations whereas the text says 10000 (seems like a typo in scaling the plot). Figs 1 and 5 illustrate only one cycle (one increase and one decrease) of CLR. It would be interesting to see cases where more than one cycle is required and to see what happens when the LR increases the second time. 6. Experiments: This is a strong section where the authors show extensive reproducible experimentation to identify settings under which super-convergence works or does not work. However, the fact that the results only applies to CIFAR-10 dataset and could not be observed for ImageNet or other architectures is disappointing and heavily takes away from the significance of this work. Overall, the work is presented as a positive result in very specific conditions but it seems more like a negative result. It would be more appealing if the paper is presented as a negative result and strengthened by additional experimentation and theoretical backing.