Title: DYNAMIC MODEL PRUNING WITH FEEDBACK

Abstract: Deep neural networks often have millions of parameters. This can hinder their deployment to low-end devices, not only due to high memory requirements but also because of increased latency at inference. We propose a novel model compression method that generates a sparse trained model without additional overhead: by allowing (i) dynamic allocation of the sparsity pattern and (ii) incorporating feedback signal to reactivate prematurely pruned weights we obtain a performant sparse model in one single training pass (retraining is not needed, but can further improve the performance). We evaluate our method on CIFAR-10 and ImageNet, and show that the obtained sparse models can reach the state-of-the-art performance of dense models. Moreover, their performance surpasses that of models generated by all previously proposed pruning schemes.

1. INTRODUCTION:

Highly overparametrized deep neural networks show impressive results on machine learning tasks. However, with the increase in model size comes also the demand for memory and computer power at inference stage-two resources that are scarcely available on low-end devices. Pruning techniques have been successfully applied to remove a significant fraction of the network weights while preserving test accuracy attained by dense models. In some cases, the generalization of compressed networks has even been found to be better than with full models (Han et al., 2015; 2017; Mocanu et al., 2018) .

The sparsity of a network is the number of weights that are identically zero, and can be obtained by applying a sparsity mask on the weights. There are several different approaches to find sparse models. For instance, one-shot pruning strategies find a suitable sparsity mask by inspecting the weights of a pretrained network (Mozer & Smolensky, 1989; LeCun et al., 1990; Han et al., 2017) . While these algorithms achieve a substantial size reduction of the network with little degradation in accuracy, they are computationally expensive (training and refinement on the dense model), and they are outperformed by algorithms that explore different sparsity masks instead of a single one. In dynamic pruning methods, the sparsity mask is readjusted during training according to different criteria (Mostafa & Wang, 2019; Mocanu et al., 2018) . However, these methods require fine-tuning of many hyperparameters.

We propose a new pruning approach to obtain sparse neural networks with state-of-the-art test accuracy. Our compression scheme uses a new saliency criterion that identifies important weights in the network throughout training to propose candidate masks. As a key feature, our algorithm not only evolves the pruned sparse model alone, but jointly also a (closely related) dense model that is used in a natural way to correct for pruning errors during training. This results in better generalization properties on a wide variety of tasks, since the simplicity of the scheme allows us further to study it from a theoretical point of view, and to provide further insights and interpretation. We do not require tuning of additional hyperparameters, and no retraining of the sparse model is needed (though can further improve performance).

Contributions.:

• A novel dynamic pruning scheme, that incorporates an error feedback in a natural way Sec. 3 and finds a trained sparse model in one training pass.

Sec. 5 • We demonstrate state-of-the-art performance (in accuracy and sparsity), Sec. 5 ourperforming all previously proposed pruning schemes.

Sec. 5 • We complement our results by an ablation study that provides further insights.

Sec. 6 and convergence analysis for convex and non-convex objectives.

Sec. 4

2. RELATED WORK:

Previous works on obtaining pruned networks can (loosely) be divided into three main categories.

Pruning after training.:

Training approaches to obtain sparse networks usually include a three stage pipeline-training of a dense model, one-shot pruning and fine-tuning-e.g., (Han et al., 2015) . Their results (i.e., moderate sparsity level with minor quality loss) made them the standard method for network pruning and led to several variations (Guo et al., 2016; Carreira-Perpinán & Idelbayev, 2018) .

Pruning during training. Zhu & Gupta (2017) propose the use of magnitude-based pruning and to gradually increase the sparsity ratio while training the model from scratch. A pruning schedule determines when the new masks are computed (extending and simplifying (Narang et al., 2017) ). He et al. (2018) (SFP) prune entire filters of the model at the end of each epoch, but allow the pruned filters to be updated when training the model. Deep Rewiring (DeepR) (Bellec et al., 2018) allows for even more adaptivity by performing pruning and regrowth decisions periodically. This approach is computationally expensive and challenging to apply to large networks and datasets. Sparse evolutionary training (SET) (Mocanu et al., 2018) simplifies prune-regrowth cycles by using heuristics for random growth at the end of each training epoch and NeST (Dai et al., 2019) by inspecting gradient magnitudes.

Dynamic Sparse Reparameterization (DSR) (Mostafa & Wang, 2019) implements a pruneredistribute-regrowth cycle where target sparsity levels are redistributed among layers, based on loss gradients (in contrast to SET, which uses fixed, manually configured, sparsity levels). Sparse Momentum (SM) (Dettmers & Zettlemoyer, 2019) follows the same cycle but instead using the mean momentum magnitude of each layer during the redistribute phase. SM outperforms DSR on ImageNet for unstructured pruning by a small margin but has no performance difference on CIFAR experiments. Our approach also falls in the dynamic category but we use error compensation mechanisms instead of hand crafted redistribute-regrowth cycles.

Pruning before training. Recently-spurred by the lottery ticket hypothesis (LT) (Frankle & Carbin, 2019) -methods which try to find a sparse mask that can be trained from scratch have attracted increased interest. For instance, Lee et al. (2019) propose SNIP to find a pruning mask by inspecting connection sensitivities and identifying structurally important connections in the network for a given task. Pruning is applied at initialization, and the sparsity mask remains fixed throughout training. 2018) learn gating variables (e.g. through 0 regularization) that minimize the number of nonzero weights, recent parallel work studies filter pruning for pre-trained models (You et al., 2019) . Gal et al. (2017) ; Neklyudov et al. (2017) ; Molchanov et al. (2017) prune from Bayesian perspectives to learn dropout probabilities during training to prune and sparsify networks as dropout weight probabilities reach 1. Gale et al. (2019) extensively study recent unstructured pruning methods on large-scale learning tasks, and find that complex techniques (Molchanov et al., 2017; Louizos et al., 2018) perform inconsistently. Simple magnitude pruning approaches achieve comparable or better results (Zhu & Gupta, 2017) .

3. METHOD:

We consider the training of a non-convex loss function f : R d → R. We assume for a weight vector w ∈ R d to have access to a stochastic gradient g(w) ∈ R d such that E[g(w)] = ∇f (w). This corresponds to the standard machine learning setting with g(w) representing a (mini-batch) gradient of one (or several) components of the loss function. Stochastic Gradient Descent (SGD) computes a sequence of iterates by the update rule

EQUATION

for some learning rate γ t . To obtain a sparse model, a general approach is to prune some of the weights of w t , i.e., to set them to zero. Such pruning can be implemented by applying a mask m ∈ {0, 1} d to the weights, resulting in a sparse model w t := m w t , where denotes the entry-wise (Hadamard) product. The mask could potentially depend on the weights w t (e.g., smallest magnitude pruning), or depend on t (e.g., the sparsity is incremented over time).

Before we introduce our proposed dynamic pruning scheme, we formalize the three main existing types of pruning methodologies (summarized in Figure 1 ). These approaches differ in the way the mask is computed, and the moment when it is applied.1 

Pruning before training. A mask m 0 (depending on e.g. the initialization w 0 or the network architecture of f ) is applied and (SGD) is used for training on the resulting subnetwork f (w) := f (m 0 w) with the advantage that only pruned weights need to be stored and updated2 , and that by training with SGD a local minimum of the subnetwork f (but not of f -the original training target) can be reached.

In practice however, it remains a challenge to efficiently determine a good mask m 0 and a wrongly chosen mask at the beginning strongly impacts the performance.

Pruning after training (one-shot pruning). A dense model is trained, and pruning is applied to the trained model w T . As the pruned model w T = m T w T is very likely not at a local optimum of f , fine-tuning (retraining with the fixed mask m T ) is necessary to improve performance.

Pruning during training (incremental and dynamic pruning). Dynamic schemes change the mask m t every (few) iterations based on observations during training (i.e. by observing the weights and stochastic gradients). Incremental schemes monotonically increase the sparsity pattern, fully dynamic schemes can also reactivate previously pruned weights. In contrast to previous dynamic schemes that relied on elaborated heuristics to adapt the mask m t , we propose a simpler approach:

Dynamic pruning with feedback (DPF, Algorithm 1). Our scheme evaluates a stochastic gradient at the pruned model w t = m t w t and applies it to the (simultaneously maintained) dense model w t :

w t+1 := w t -γ t g(m t w t ) = w t -γ t g( w t ) . (DPF)

Applying the gradient to the full model allows to recover from "errors", i.e. prematurely masking out important weights: when the accumulated gradient updates from the following steps drastically change a specific weight, it can become activated again (in contrast to incremental pruning approaches that have to stick to sub-optimal decisions). For illustration, observe that (DPF) can equivalently be written as

w t+1 = w t -γ t g(w t + e t ),

where e t := w t -w t is the error produced by the compression. This provides a different intuition of the behavior of (DPF), and connects it with the concept of error-feedback (Stich et al., 2018; Karimireddy et al., 2019) . 3 We illustrate this principle in Figure 2 and give detailed pseudocode and further implementation details in Appendix A.1. The DPF scheme can also be seen as an instance of a more general class of schemes that apply (arbitrary) perturbed gradient updates to the dense model. For instance straight-through gradient estimators (Bengio et al., 2013) that are used to empirically simplify the backpropagation can be seen as such perturbations. Our stronger assumptions on the structure of the perturbation allow to derive non-asymptotic convergence rates in the next section, though our analysis could also be extended to the setting in (Yin et al., 2019) if the perturbations can be bounded.

4. CONVERGENCE ANALYSIS:

We now present convergence guarantees for (DPF). For the purposes of deriving theoretical guarantees, we assume that the training objective is smooth, that is ∇f (w) -∇f (v) ≤ L wv , ∀w, v ∈ R d , for a constant L > 0, and that the stochastic gradients are bounded E g( w t ) 2 ≤ G 2 for every pruned model w t = m t (w t ) w t . The quality of this pruning is defined as the parameter

δ t ∈ [0, 1] such that δ t := w t -w t 2 w t 2 .

(1)

Pruning without information loss corresponds to w t = w t , i.e., δ t = 0, and in general δ t ≤ 1.

Convergence on Convex functions. We first consider the case when f is in addition µ-strongly

convex, that is ∇f (v), w -v ≤ f (w) -f (v) -µ 2 w -v 2 , ∀w, v ∈ R d .

While it is clear that this assumption does not apply to neural networks, it eases the presentation as strongly convex functions have a unique (global) minimizer w := arg min w∈R d f (w).

Theorem 4.1. Let f be µ-strongly convex and learning rates given as γ t = 4 µ(t+2) . Then for a randomly chosen pruned model u of the iterates { w 0 , . . . , w T } of DPF, concretely u = w t with probability p t = 2(t+1) (T +1)(T +2) , it holds that-in expectation over the stochasticity and the selection of u:

EQUATION

)

The rightmost term in (2) measures the average quality of the pruning. However, unless δ t → 0 or w t → 0 for t → ∞, the error term never completely vanishes, meaning that the method converges only to a neighborhood of the optimal solution (this not only holds for the pruned model, but also for the jointly maintained dense model, as we will show in the appendix). This behavior is expected, as the global optimal model w might be dense and cannot be approximated well by a sparse model.

For one-shot methods that only prune the final (SGD) iterate w T at the end, we have instead:

Ef ( w t ) -f (w ) ≤ 2E (f (w T ) -f (w )) + Lδ T E w T 2 = O LG 2 µ 2 T + LE δ T w T 2 ,

as we show in the appendix. First, we see from this expression that the estimate is very sensitive to δ T and w T , i.e. the quality of the pruning the final model. This could be better or worse than the average of the pruning quality of all iterates. Moreover, one looses also a factor of the condition number L µ in the asymptotically decreasing term, compared to (2). This is due to the fact that standard convergence analysis only achieves optimal rates for an average of the iterates (but not the last one). This shows a slight theoretical advantage of DPF over rounding at the end.

Convergence on Non-Convex Functions to Stationary Points. Secondly, we consider the case when f is a non-convex function and show convergence (to a neighborhood) of a stationary point.

Theorem 4.2. Let learning rate be given as

γ = c √ T , for c = f (w0)-f (w )

LG 2

. Then for pruned model u chosen uniformly at random from the iterates { w 0 , . . . , w T } of DPF, concretely u := w t with probability p t = 1 T +1 , it holds-in expectation over the stochasticity and the selection of u:

EQUATION

)

Extension to Other Compression Schemes. So far we put our focus on simple mask pruning schemes to achieve high model sparsity. However, the pruning scheme in Algorithm 1 could be replaced by an arbitrary compressor C : R d → R d , i.e., w t = C(w t ). Our analysis extends to compressors as e.g. defined in (Karimireddy et al., 2019; Stich & Karimireddy, 2019) , whose quality is also measured in terms of the δ t parameters as in (1). For example, if our objective was not to obtain a sparse model, but to produce a quantized neural network where inference could be computed faster on low-precision numbers, then we could define C as a quantized compressor. One variant of this approach is implemented in the Binary Connect algorithm (BC) (Courbariaux et al., 2015) with prominent results, see also (Li et al., 2017) for further insights and discussion.

5. EXPERIMENTS:

We evaluated DPF together with its competitors on a wide range of neural architectures and sparsity levels. DPF exhibits consistent and noticeable performance benefits over its competitors.

5.1. EXPERIMENTAL SETUP:

Datasets. We evaluated DPF on two image classification benchmarks: (1) CIFAR-10 (Krizhevsky & Hinton, 2009) (50K/10K training/test samples with 10 classes), and (2) ImageNet (Russakovsky et al., 2015) (1.28M/50K training/validation samples with 1000 classes). We adopted the standard data augmentation and preprocessing scheme from He et al. (2016a) ; Huang et al. (2016) ; for further details refer to Appendix A.2.

Models.:

Following the common experimental setting in related work on network pruning (Liu et al., 2019; Gale et al., 2019; Dettmers & Zettlemoyer, 2019; Mostafa & Wang, 2019) , our main experiments focus on ResNet (He et al., 2016a) and WideResNet (Zagoruyko & Komodakis, 2016) . However, DPF can be effectively extended to other neural architectures, e.g., VGG (Simonyan & Zisserman, 2015) , DenseNet (Huang et al., 2017) . We followed the common definition in He et al. (2016a) ; Zagoruyko & Komodakis (2016) and used ResNet-a, WideResNet-a-b to represent neural network with a layers and width factor b.

Baselines. We considered the state-of-the-art model compression methods presented in the table below as our strongest competitors. We omit the comparison to other dynamic reparameterization methods, as DSR can outperform DeepR (Bellec et al., 2018) and SET (Mocanu et al., 2018 ) by a noticeable margin (Mostafa & Wang, 2019) . overhead. We perform pruning across all neural network layers (no layer-wise pruning) using magnitude-based unstructured weight pruning (inherited from Han et al. (2015) ). We found the best preformance when updating the mask every 16 iterations (see also Table 11 ) and we keep this value fixed for all experiments (independent of the architecture or task).

Implementation

Unlike our competitors that may ignore some layers (e.g. the first convolution and downsampling layers in DSR), we applied DPF (as well as the One-shot P+FT and Incremental baselines) to all convolutional layers while keeping the last fully-connected layer4 , biases and batch normalization layers dense. Lastly, our algorithm gradually increases the sparsity s t of the mask from 0 to the desired sparsity using the same scheduling as in (Zhu & Gupta, 2017) ; see Appendix A.2.

Training schedules. For all competitors, we adapted their open-sourced code and applied a consistent (and standard) training scheme over different methods to ensure a fair comparison. Following the standard training setup for CIFAR-10, we trained ResNet-a for 300 epochs and decayed the learning rate by 10 when accessing 50% and 75% of the total training samples (He et al., 2016a; Huang et al., 2017) ; and we trained WideResNet-a-b as Zagoruyko & Komodakis (2016) for 200 epochs and decayed the learning rate by 5 when accessing 30%, 60% and 80% of the total training samples. For ImageNet training, we used the training scheme in (Goyal et al., 2017) for 90 epochs and decayed learning rate by 10 at 30, 60, 80 epochs. For all datasets and models, we used mini-batch SGD with Nesterov momentum (factor 0.9) with fine-tuned learning rate for DPF. We reused the tuned (or recommended) hyperparameters for our baselines (DSR and SM), and fine-tuned the optimizer and learning rate for One-shot P+FT, Incremental and SNIP. The mini-batch size is fixed to 128 for CIFAR-10 and 1024 for ImageNet regardless of datasets, models and methods.

5.2. EXPERIMENT RESULTS:

CIFAR-10. Figure 3 shows a comparison of different methods for WideResNet-28-2. For low sparsity level (e.g. 50%), DPF outperforms even the dense baseline, which is in line with regularization properties of network pruning. Furthermore, DPF can prune the model up to a very high level (e.g. 99%), and still exhibit viable performance. This observation is also present in Table 1 , where the results of training different state-of-the-art DNN architectures with higher sparsities are depicted. DPF shows reasonable performance even with extremely high sparsity level on large models (e.g. WideResNet-28-8 with sparsity ratio 99.9%), while other methods either suffer from significant quality loss or even fail to converge.

Because simple model pruning techniques sometimes show better performance than complex techniques (Gale et al., 2019) , we further consider these simple models in Table 2 . While DPF outperforms them in almost all settings, it faces difficulties pruning smaller models to extremely high sparsity ratios (e.g. ResNet-20 with sparsity ratio 95%).This seems however to be an artifact of fine-tuning, as DPF with extra fine-tuning convincingly outperforms all other methods regardless of the network size. This comes to no surprise as schemes like One-shot P+FT and Incremental do not benefit from extra fine-tuning, since it is already incorporated in their training procedure and they might become stuck in local minima. On the other hand, dynamic pruning methods, and in particular DPF, work on a different paradigm, and can still heavily benefit from fine-tuning. 5Figure 13 (in Appendix A.3.4) depicts another interesting property of DPF. When we search for a subnetwork with a (small) predefined number of parameters for a fixed task, it is much better to run DPF on a large model (e.g. WideResNet-28-8) than on a smaller one (e.g. WideResNet-28-2). That is, DPF performs structural exploration more efficiently in larger parametric spaces. ImageNet. We compared DPF to other dynamic reparameterization methods as well as the strong Incremental baseline in Table 3 . For both sparsity levels (80% and 90%), DPF shows a significant improvement of top-1 test accuracy with fewer or equal number of parameters. We also directly (and fairly) compare with the results of Incremental (Zhu & Gupta, 2017) reproduced and fine-tuned by Gale et al. (2019) , where they consider layer-wise sparsity ratio and fine-tune both the sparsity warmup schedule and label-smoothing for better performance. 

6. DISCUSSION:

Besides the theoretical guarantees, a straightforward benefit of DPF over one-shot pruning in practice is its fine-tuning free training process. Figure 12 in the appendix (Section A.3.3) demonstrates the trivial computational overhead (considering the dynamic reparameterization cost) of involving DPF to train the model from scratch. Small number of hyperparameters compared to other dynamic reparameterization methods (e.g. DSR and SM) is another advantage of DPF and Figure 11 further studies how different setups of DPF impact the final performance. Notice also that for DPF, inference is done only at sparse models, an aspect that could be leveraged for more efficient computations.

Empirical difference between one-shot pruning and DPF. From the Figure 2 one can see that DPF tends to oscillate among several local minima, whereas one-shot pruning, even with fine tuning, converges to a single solution, which is not necessarily close to the optimum. We believe that the wider exploration of DPF helps to find a better local minima (which can be even further improved by fine-tuning, as shown in Table 2 ). We empirically analyzed how drastically the mask changes between each reparameterization step, and how likely it is for some pruned weight to become non-zero in the later stages of training. Figure 4 shows at what stage of the training each element of the final mask becomes fixed. For each epoch, we report how many mask elements were flipped starting from this epoch. As an example, we see that for sparsity ratio 95%, after epoch 157 (i.e. for 43 epochs left), only 5% of the mask elements were changing. This suggests that, except for a small percentage of weights that keep oscillating, the mask has converged early in the training. In the final epochs, the algorithm keeps improving accuracy, but the masks are only being fine-tuned. A similar mask convergence behavior can be found in Appendix (Figure 7 ) for training ResNet-20 on CIFAR-10. DPF does not find a lottery ticket. The LT hypothesis (Frankle & Carbin, 2019) conjectures that for every desired sparsity level there exists a sparse submodel that can be trained to the same or better performance as the dense model. In Figure 5 we show that the mask found by DPF is not a LT, i.e., training the obtained sparse model from scratch does not recover the same performance. The (expensive) procedure proposed in Frankle & Carbin (2019) ; Frankle et al. (2019) finds different masks and achieves the same performance as DPF for mild sparsity levels, but DPF is much better for extremely sparse models (99% sparsity). Extension to structured pruning. The current state-of-the-art dynamic reparameterization methods only consider unstructured weight pruning. Structured filter pruning6 is either ignored (Bellec et al., 2018; Mocanu et al., 2018; Dettmers & Zettlemoyer, 2019) or shown to be challenging (Mostafa & Wang, 2019) even for the CIFAR dataset. In Figure 6 below we presented some preliminary results on CIFAR-10 to show that our DPF can also be applied to structured filter pruning schemes. DPF outperforms the current filter-norm based state-of-the-art method for structured pruning (e.g. SFP (He et al., 2018) ) by a noticeable margin. Figure 16 in Appendix A.4.3 displays the transition procedure of the sparsity pattern (of different layers) for WideResNet-28-2 with different sparsity levels. DPF can be seen as a particular neural architecture search method, as it gradually learns to prune entire layers under the guidance of the feedback signal.

We followed the common experimental setup as mentioned in Section 5 with 2 norm based filter selection criteria for structured pruning extension on CIFAR-10. We do believe a better filter selection scheme (Ye et al., 2018; He et al., 2019; Lym et al., 2019) could further improve the results but we leave this exploration for the future work. We trigger the mask update every p = 16 iterations (see also Figure 11 ) and we keep this parameter fixed throughout all experiments, independent of architecture or task.

We perform pruning across all neural network layers (no layer-wise pruning) using magnitude-based unstructured weight pruning (inherited from Han et al. (2015) ). Pruning is applied to all convolutional layers while keeping the last fully-connected layer, biases and batch normalization layers dense.

We gradually increases the sparsity s t of the mask from 0 to the desired sparsity using the same scheduling as in Zhu & Gupta (2017) ; see Appendix A.2 below.

A.2 IMPLEMENTATION DETAILS:

We implemented our DPF in PyTorch (Paszke et al., 2017) . All experiments were run on NVIDIA Tesla V100 GPUs. Sparse tensors in our implementation are respresented as the dense tenors multiplied by the corresponding binary masks.

Datasets We evaluate all methods on the following standard image classification tasks:

• Image classification for CIFAR-10 (Krizhevsky & Hinton, 2009) . Dataset consists of a training set of 50K and a test set of 10K color images of 32 × 32 pixels, as well as 10 target classes. We adopt the standard data augmentation and preprocessing scheme (He et al., 2016a; Huang et al., 2016) .

• Image classification for ImageNet (Russakovsky et al., 2015) . The ILSVRC 2012 classification dataset consists of 1.28 million images for training, and 50K for validation, with 1K target classes. We use ImageNet-1k (Deng et al., 2009) and adopt the same data preprocessing and augmentation scheme as in He et al. (2016a; b) ; Simonyan & Zisserman (2015) .

Gradual Pruning Scheme For Incremental baseline, we tuned their automated gradual pruning scheme

s t = s f + (s i -s f ) 1 -t-t0

n∆t 3 to gradually adjust the pruning sparsity ratio s t for t ∈ {t 0 , . . . , t 0 + n∆t}. That is, in our setup, we increased from an initial sparsity ratio s i = 0 to the desired target model sparsity ratio s f over the epoch (n) when performing the second learning rate decay, from the training epoch t 0 = 0 and with pruning frequency ∆t = 1 epoch. In our experiments, we used this gradual pruning scheme over different methods, except One-shot P+FT, SNIP, and the methods (DSR, SM) that have their own fine-tuned gradual pruning scheme.

Hyper-parameters tuning procedure We grid-searched the optimal learning rate, starting from the range of {0.05, 0.10, 0.15, 0.20}. More precisely, we will evaluate a linear-spaced grid of learning rates. If the best performance was ever at one of the extremes of the grid, we would try new grid points so that the best performance was contained in the middle of the parameters.

We trained most of the methods by using mini-batch SGD with Nesterov momentum. For baselines involving fine-tuning procedure (e.g. Table 2 ), we grid-searched the optimal results by tuning the optimizers (i.e. mini-batch SGD with Nesterov momentum, or Adam) and the learning rates.

The optimal hyper-parameters for DPF The mini-batch size is fixed to 128 for CIFAR-10 and 1024 for ImageNet regardless of datasets and models.

For CIFAR-10, we trained ResNet-a and VGG for 300 epochs and decayed the learning rate by 10 when accessing 50% and 75% of the total training samples (He et al., 2016a; Huang et al., 2017) ; and we trained WideResNet-a-b as Zagoruyko & Komodakis (2016) for 200 epochs and decayed the learning rate by 5 when accessing 30%, 60% and 80% of the total training samples. The optimal learning rate for ResNet-a, WideResNet-a-b and VGG are 0.2, 0.1 and 0.2 respectively; the corresponding weight decays are 1e-4, 5e-4 and 1e-4 respectively.

For ImageNet training, we used the training scheme in Goyal et al. (2017) for 90 epochs, where we gradually warmup the learning rate from 0.1 to 0.4 and decayed learning rate by 10 at 30; 60; 80 epochs. The used weight decay is 1e-4. and Incremental (Zhu & Gupta, 2017) , demonstrating the fact that our scheme enables a drastical reparameterization over the dense parameter space for better generalization performance. sizes, as well as whether layer-wise pruning or not. We can witness that the optimal reparameterization (i.e p = 16) is quite consistent over different sparsity ratios and different mini-batch sizes, and we used it in all our experiments. The global-wise unstructured weight pruning (instead of layer-wise weight pruning) allows our DPF more flexible to perform dynamic parameter reallocation, and thus can provide better results especially for more aggressive pruning sparsity ratios. However, we also need to note that, for the same number of compressed parameters (layerwise or globalwise unstructured weight pruning), using global-wise pruning leads to a slight increase in the amount of MACs, as illustrated Table 5 . Figure 12 demonstrates the trivial computational overhead of involving DPF to gradually train a compressed model (ResNet-50) from scratch (on ImageNet). Note that we evaluated the introduced reparameterization cost for dynamic pruning, which is independent of (potential) significant system speedup brought by the extreme high model sparsity. Even though our work did not estimate the practical speedup, we do believe we can have a similar training efficiency as the values reported in Dettmers & Zettlemoyer (2019) . Figure 13 below demonstrates that for the same pruned model size (i.e. any point in the x-axis), we can always perform "architecture search" to get a better (in terms of generalization) pruned model, from a larger network (e.g. WideResNet-28-8) rather than the one searched from a relatively small network (e.g. WideResNet-28-4). (e) Structured pruning with 80% sparsity.

A.3 ADDITIONAL RESULTS FOR UNSTRUCTURED PRUNING:



A.3.1 COMPLETE RESULTS OF UNSTRUCTURED PRUNING ON CIFAR-10:

Figure 16 : The element-wise sparsity of each layer for WideResNet28-2 (trained on CIFAR-10 via DPF), under different structured pruning sparsity ratios. DPF for model compression (with structured pruning) performs implicit as a neural architecture search.

B MISSING PROOFS:

In this section we present the proofs for the claims in Section 4.

First, we give the proof for the stongly convex case. Here we follow Lacoste-Julien et al. ( 2012) for the general structure, combined with estimates from the error-feedback framework (Stich et al., 2018; Stich & Karimireddy, 2019) to control the pruning errors.

Proof of Theorem 4.1. By definition of (DPF), w t+1 = w t -γ t g( w t ), hence, et al. (2012) . Combining these two estimates (with w t = w T ) shows the claim. Furthermore, we also claimed that also the dense model converges to a neighborhood of optimal solution. This follows by L-smoothness and (4): For any fixed model w t we have the estimate (4), hence for a randomly chosen (dense) model u (from the same distribution as the sparse model in Theorem 4.1) we have Ef (u) -f (w ) Summing these inequalities from t = 0 to t = T gives

E

1 T + 1 T t=0 E ∇f ( w t ) 2 ≤ 2 γ(T + 1) T t=0 (E [f (w t )] -E [f (w t+1 )]) + γLG 2 + L 2 T + 1 T t=0 E w t -w t 2 ≤ 2 (f (w 0 ) -f (w )) γ(T + 1) + LγG 2 + L 2 T + 1 T t=0 E e t 2 .

Finally, using w t -w t 2 = δ t w t 2 by (1), and plugging in the stepsize γ that minimizes the right hand side shows the claim.

:

The method introduced in Section

typically follow one of these broad themes loosely, with slight variations in detail. For the sake of clarity we omit a too technical and detailed discussion here.2 When training on f (w), it suffices to access stochastic gradients of f (w), denoted by g(w), which can potentially be cheaper be computed than by naively applying the mask to g(w) (note g(w) = m0 g(w)).

Our variable wt corresponds to xt in the notation ofKarimireddy et al. (2019). Their error-fixed SGD algorithm evaluates gradients at perturbed iterates xt := xt + et, which correspond precisely to wt = wt + et in our notation. This shows the connection of these two methods.

The last fully-connected layer normally makes up only a very small faction of the total MACs, e.g. 0.05% for ResNet-50 on ImageNet and 0.0006% for WideResNet-28-2 on CIFAR-10.

Besides that a large fraction of the mask elements already converge during training (see e.g. Figure4below), not all mask elements converge. Thus DPF can still benefit from fine-tuning on the fixed sparsity mask.

Lym et al. (2019) consider structured filter pruning and reconfigure the large (but sparse) model to small (but dense) model during the training for better training efficiency. Note that they perform model update on a gradually reduced model space, and it is completely different from the dynamic reparameterization methods (e.g. DSR, SM and our scheme) that perform reparameterization under original (full) model space.

Figures and Tables:

Figure 1: Schematic view of different pruning methodologies and their properties.

Figure2: Left: One-shot pruning (red) computes a stochastic gradient at w and takes a step towards the best dense model. In contrast, DPF (blue) computes a stochastic gradient at the pruned model w (here obtained by smallest magnitude pruning), and takes a step that best suits the compressed model. Right: One-shot pruning commits to a single sparsity mask and might obtain sparse models that generalize poorly (without retraining). DPF explores different available sparsity patterns and finds better sparse models.

Figure4: Convergence of the pruning mask mt of DPF for different target sparsity levels (see legend). The y-axis represent the percentage of mask elements that still change after a certain epoch (x-axis). The illustrated example are from WideResNet-28-2 on CIFAR-10. We decayed the learning rate at 60,120,160 epochs.

Figure6: MAC v.s. top-1 test accuracy, for training WideResNet-28 (with different width) on CIFAR-10. The reported results are averaged over three runs. The WideResNet-28-2 has 216M MACs, WideResNet-28-4 has 848M MACs and WideResNet-28-8 has 3366M MACs. Other detailed information refers to the Appendix A.4.1, e.g., the # of params v.s. top-1 test accuracy in Figure14, and the numerical test accuracy score in Table6.

Figure 7 and Figure 8 complements Figure 4, and details the training dynamics (e.g. the converge of δ and masks) of DPF from the other aspect. Figure 9 compares the training dynamics between DPF

Figure 7: Convergence of the pruning mask mt of DPF for different target sparsity levels (see legend). The y-axis represent the percentage of mask elements that still change after a certain epoch (x-axis). The illustrated example are from ResNet-20 on CIFAR-10. We decayed the learning rate at 150 and 225 epochs.

Figure 8: Training dynamics of DPF (WideResNet-28-2 on CIFAR-10) for unstructured pruning with different sparsity ratios. IoU stands for Intersection over Union for the non-masked elements of two consecutive masks; the smaller value the more fraction of the masks will flip.

Figure9: of the pruning mask mt of DPF V.S. Incremental(Zhu & Gupta, 2017) for different target sparsity levels (see legend). The y-axis represent the percentage of mask elements that still change after a certain epoch (x-axis). The illustrated example are from WideResNet-28-2 on CIFAR-10. We decayed the learning rate at 60; 120; 160 epochs. These two schemes use the same gradual warmup schedule (and hyper-parameters) for the pruning ratio during the training.

Figure 10 in addition to the Figure 5 (in the main text) further studies the lottery ticket hypothesis under different training budgets (same epochs or same total flops). The results of DPF also demonstrate the importance of training-time structural exploration as well as the corresponding implicit regularization effects. Note that we do not want to question the importance of the weight initialization or the existence of the lottery ticket. Instead, our DPF can provide an alternative training scheme to compress the model to an extremely high compression ratio without sacrificing the test accuracy, where most of the existing methods still meet severe quality loss (including Frankle & Carbin (2019); Liu et al. (2019); Frankle et al. (2019)).A.3.3 COMPUTATIONAL OVERHEAD AND THE IMPACT OF HYPER-PARAMETERSIn Figure11, we evaluated the top-1 test accuracy of a compressed model trained by DPF under different setups, e.g., different reparameterization period p, different sparsity ratios, different mini-batch

Figure 10: Investigate the effect of lottery ticket for model compression (unstructured weight pruning for WideResNet28-2 on CIFAR-10). It complements the observations in Figure 5 by retraining the model for the same amount of computation budget (i.e. flops).

Reparameterization period vs. sparsity ratio. The heatmap value is the test accuracy minus the best accuracy of the corresponding sparsity. We use mini-batch size 128 w/o layerwise reparameterization.

Figure 11: Investigate how the reparameterization period/scheme and mini-batch size impact the generalization performance (test top-1 accuracy), for dynamically training (and reparameterizing) a compressed model from scratch (ResNet-20 with CIFAR-10).

Figure12: The learning curve of our DPF (with unstructured magnitude pruning) and the standard mini-batch SGD for training ResNet50 on ImageNet. Our proposed DPF has trivial computational overhead. We trained ResNet50 on 4 NVIDIA V100 GPUs with 1024 mini-batch size. The target sparsity ratio is 80%.

Figure 13: Test top-1 accuracy vs. the compressed model size, for training WideResNet-28 (with different widths) on CIFAR-10. The compressed model is searched from WideResNet-28 (fixed depth) with different width (number of filters per layer).

Figure 14: # of params v.s. top-1 test accuracy, for training WideResNet28 (with different width) on CIFAR-10. Structured filter-wise pruning is used here. The reported results are averaged over three runs.

2E [f ( u) -f (w )] + LE δ t w t 2 (Thm 4.1)

WideResNet-28-2 on CIFAR-10 for unstructured weight pruning. The original model has 1.47M parameters with 216M MACs (Multiplier-ACcumulator). We varied the sparsity ratio from 50% to 99%. The complete numerical test accuracy values refer to Table4in Appendix A.3.1. The lower # of params and MACs the model has, the higher sparsity ratio it uses. All results are averaged over three runs. Note that different methods might consider different types of layers and thus the same pruning sparsity ratio might result in the slight difference for both of # of params and MACs. The DSR cannot converge when using the extreme high sparsity ratio (99%). Top-1 test accuracy of SOTA DNNs on CIFAR-10 for unstructured weight pruning. We considered unstructured pruning and the indicates the method cannot converge. The results are averaged for three runs. The results we presented for each model consider some reasonable pruning ratios (we prune more aggressively for deeper and wider neural networks), and readers can refer to Table4(in Appendix A.3.1) for a complete overview.

Table content: <table><tr><td>Top-1 test accuracy</td><td>84 86 88 90 92 94</td><td colspan="2">Dense model SM DPF (ours)</td><td>DSR SNIP</td><td/><td/><td>Top-1 test accuracy</td><td>84 86 88 90 92 94</td><td/><td>Dense model SM DPF (ours)</td><td>DSR SNIP</td><td>Top-1 test accuracy</td><td>84 86 88 90 92 94</td><td>Dense model SM DPF (ours)</td><td>DSR SNIP</td></tr><tr><td/><td>0.5</td><td>0.6</td><td colspan="2">0.7 Target sparsity ratio 0.8</td><td>0.9</td><td>1.0</td><td/><td>0.0</td><td>0.2</td><td>0.4 # of params (M) 0.6</td><td>0.8</td><td>0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 MAC</td><td>1e8</td></tr><tr><td colspan="7">(a) Sparsity ratio v.s. top-1 test acc.</td><td colspan="5">(b) # of params v.s. top-1 test acc.</td><td>(c) MAC v.s. top-1 test acc.</td></tr><tr><td colspan="12">Figure 3: Top-1 test accuracy of Methods</td></tr><tr><td colspan="2">Model</td><td/><td/><td colspan="2">Baseline on dense model</td><td colspan="3">SNIP (L + , 2019)</td><td/><td>SM (DZ, 2019)</td><td colspan="2">DSR (MW, 2019)</td><td>DPF</td><td>Target Pr. ratio</td></tr><tr><td colspan="3">VGG16-D</td><td/><td colspan="2">93.74 ± 0.13</td><td colspan="3">93.04 ± 0.26</td><td/><td>93.59 ± 0.17</td><td>-</td><td>93.87 ± 0.15</td><td>95%</td></tr><tr><td/><td/><td/><td/><td/><td/><td colspan="3">91.10 ± 0.22</td><td/><td>91.98 ± 0.01</td><td colspan="2">92.00 ± 0.19</td><td>92.42 ± 0.14</td><td>70%</td></tr><tr><td colspan="3">ResNet-20</td><td/><td colspan="2">92.48 ± 0.20</td><td colspan="3">90.53 ± 0.27 88.50 ± 0.13</td><td/><td>91.54 ± 0.16 89.76 ± 0.40</td><td colspan="2">91.78 ± 0.28 87.88 ± 0.04</td><td>92.17 ± 0.21 90.88 ± 0.07</td><td>80% 90%</td></tr><tr><td/><td/><td/><td/><td/><td/><td colspan="3">84.91 ± 0.25</td><td/><td>83.03 ± 0.74</td><td/><td>88.01 ± 0.30</td><td>95%</td></tr><tr><td colspan="3">ResNet-32</td><td/><td colspan="2">93.83 ± 0.12</td><td colspan="3">90.40 ± 0.26 87.23 ± 0.29</td><td/><td>91.54 ± 0.18 88.68 ± 0.22</td><td colspan="2">91.41 ± 0.23 84.12 ± 0.32</td><td>92.42 ± 0.18 90.94 ± 0.35</td><td>90% 95%</td></tr><tr><td colspan="3">ResNet-56</td><td/><td colspan="2">94.51 ± 0.20</td><td colspan="3">91.43 ± 0.34</td><td/><td>92.73 ± 0.21 90.96 ± 0.40</td><td colspan="2">93.78 ± 0.20 92.57 ± 0.09</td><td>93.95 ± 0.11 92.74 ± 0.08</td><td>90% 95%</td></tr><tr><td/><td/><td/><td/><td/><td/><td colspan="3">92.58 ± 0.22</td><td/><td>93.41 ± 0.22</td><td colspan="2">93.88 ± 0.08</td><td>94.36 ± 0.24</td><td>90%</td></tr><tr><td colspan="6">WideResNet-28-2 95.01 ± 0.04</td><td colspan="3">90.80 ± 0.04</td><td/><td>92.24 ± 0.14</td><td colspan="2">92.74 ± 0.17</td><td>93.62 ± 0.05</td><td>95%</td></tr><tr><td/><td/><td/><td/><td/><td/><td colspan="3">83.45 ± 0.38</td><td/><td>85.36 ± 0.80</td><td/><td>88.92 ± 0.29</td><td>99%</td></tr><tr><td/><td/><td/><td/><td/><td/><td colspan="3">93.62 ± 0.17</td><td/><td>94.45 ± 0.14</td><td colspan="2">94.63 ± 0.08</td><td>95.38 ± 0.04</td><td>95%</td></tr><tr><td colspan="6">WideResNet-28-4 95.69 ± 0.10</td><td colspan="3">92.06 ± 0.38</td><td/><td>93.80 ± 0.24</td><td colspan="2">93.92 ± 0.16</td><td>94.98 ± 0.08</td><td>97.5%</td></tr><tr><td/><td/><td/><td/><td/><td/><td colspan="3">89.49 ± 0.20</td><td/><td>92.18 ± 0.04</td><td colspan="2">92.50 ± 0.07</td><td>93.86 ± 0.20</td><td>99%</td></tr><tr><td/><td/><td/><td/><td/><td/><td colspan="3">95.49 ± 0.21</td><td/><td>95.67 ± 0.14</td><td colspan="2">95.81 ± 0.10</td><td>96.08 ± 0.15</td><td>90%</td></tr><tr><td/><td/><td/><td/><td/><td/><td colspan="3">94.92 ± 0.13</td><td/><td>95.64 ± 0.07</td><td colspan="2">95.55 ± 0.12</td><td>95.98 ± 0.10</td><td>95%</td></tr><tr><td colspan="6">WideResNet-28-8 96.06 ± 0.06</td><td colspan="3">94.11 ± 0.19</td><td/><td>95.31 ± 0.20</td><td colspan="2">95.11 ± 0.07</td><td>95.84 ± 0.04</td><td>97.5%</td></tr><tr><td/><td/><td/><td/><td/><td/><td colspan="3">92.04 ± 0.11</td><td/><td>94.38 ± 0.12</td><td colspan="2">94.10 ± 0.12</td><td>95.63 ± 0.16</td><td>99%</td></tr><tr><td/><td/><td/><td/><td/><td/><td colspan="3">74.50 ± 2.23</td><td/><td/><td colspan="2">88.65 ± 0.36</td><td>91.76 ± 0.18</td><td>99.9%</td></tr></table>

Top-1 test accuracy of SOTA DNNs on CIFAR-10 for unstructured weight pruning via some simple pruning techniques. This table complements Table1and evaluates the performance of model compression under One-shot P+FT and Incremental, as well as how extra fine-tuning (FT) impact the performance of Incremental and our DPF. Note that One-shot P+FT prunes the dense model and uses extra fine-tuning itself. The Dense, Incremental and DPF train with the same number of epochs from scratch. The extra fine-tuning procedure considers the model checkpoint at the end of the normal training, uses the same number of training epochs (60 epochs in our case) with tuned optimizer and learning rate. Detailed hyperparameters tuning procedure refers to Appendix A.2.

Table content: <table><tr><td>Methods</td></tr></table>

Top-1 test accuracy of ResNet-50 on ImageNet for unstructured weight pruning. The # of parameters for the full model is 25.56 M. We used the results of DSR fromMostafa & Wang (2019) as we use the same (standard) training/data augmentation scheme for the same neural architecture. Note that different methods prune different types of layers and result in the different # of parameters for the same target pruning ratio.

Table content: <table/>

details the numerical results for training SOTA DNNs on CIFAR-10. Some results of it reconstruct the Table1 and Figure 3.

Table content: <table/>

Top-1 test accuracy for training (compressed) SOTA DNNs on CIFAR-10 from scratch. We considered unstructured pruning and the indicates the method cannot converge. The results are averaged for three runs.

Table content: <table><tr><td>Methods</td></tr></table>

Investigate how the reparameterization scheme (layer-wise or not) impact the MACs (for the same number of compressed parameters), for using DPF on ResNet-20 with CIFAR-10.

Table content: <table><tr><td>Target sparsity</td><td>50%</td><td>60%</td><td>70%</td><td>80%</td><td>90%</td><td>95%</td></tr><tr><td colspan="7">w/o layerwise reparameterization 22.80M 19.13M 15.18M 11.12M 6.54M 4.02M</td></tr><tr><td colspan="7">w/ layerwise reparameterization 20.99M 16.91M 12.83M 8.75M 4.67M 2.63M</td></tr></table>

Performance evaluation of DPF (and other baseline methods) for training (Wide)ResNet variants on CIFAR-10. We use the norm-based criteria for filter selection (as in SFP(He et al., 2018)) to estimate the output channel-wise pruning threshold. We follow the gradual pruning warmup scheme (as inZhu & Gupta (2017)) from 0 epoch to the epoch when performing the second learning rate decay. Note that SFP prunes filters within the layer by a given ratio while our DPF prunes filters across layers. Due to the difference between filters for different layers, the # of parameters pruned by DPF might slight different from the one pruned by SFP. The pruning ratio refers to either prune filters within the layer or across the layers.

Table content: <table><tr><td/><td/><td/><td colspan="2">Methods</td></tr><tr><td colspan="2">Model</td><td>Baseline on dense model</td><td>SFP (H + , 2018)</td><td colspan="2">DPF</td><td>Target Pr. ratio</td></tr><tr><td colspan="2">ResNet-20</td><td>92.48 ± 0.20</td><td>92.18 ± 0.31</td><td colspan="2">92.54 ± 0.07</td><td>10%</td></tr><tr><td colspan="2">ResNet-20</td><td>92.48 ± 0.20</td><td>91.12 ± 0.20</td><td colspan="2">91.90 ± 0.06</td><td>20%</td></tr><tr><td colspan="2">ResNet-20</td><td>92.48 ± 0.20</td><td>90.32 ± 0.25</td><td colspan="2">91.07 ± 0.40</td><td>30%</td></tr><tr><td colspan="2">ResNet-20</td><td>92.48 ± 0.20</td><td>89.60 ± 0.46</td><td colspan="2">90.28 ± 0.26</td><td>40%</td></tr><tr><td colspan="2">ResNet-32</td><td>93.52 ± 0.13</td><td>92.07 ± 0.22</td><td colspan="2">92.18 ± 0.16</td><td>30%</td></tr><tr><td colspan="2">ResNet-32</td><td>93.52 ± 0.13</td><td>91.14 ± 0.45</td><td colspan="2">91.50 ± 0.21</td><td>40%</td></tr><tr><td colspan="2">ResNet-56</td><td>94.51 ± 0.20</td><td>93.99 ± 0.27</td><td colspan="2">94.53 ± 0.13</td><td>30%</td></tr><tr><td colspan="2">ResNet-56</td><td>94.51 ± 0.20</td><td>93.57 ± 0.16</td><td colspan="2">94.03 ± 0.38</td><td>40%</td></tr><tr><td colspan="3">WideResNet-28-2 95.01 ± 0.04</td><td>94.02 ± 0.24</td><td colspan="2">94.52 ± 0.08</td><td>40%</td></tr><tr><td colspan="3">WideResNet-28-2 95.01 ± 0.04</td><td>93.34 ± 0.14</td><td colspan="2">94.11 ± 0.12</td><td>50%</td></tr><tr><td colspan="3">WideResNet-28-2 95.01 ± 0.04</td><td>92.07 ± 0.09</td><td colspan="2">93.74 ± 0.25</td><td>60%</td></tr><tr><td colspan="3">WideResNet-28-2 95.01 ± 0.04</td><td>90.66 ± 0.62</td><td colspan="2">92.89 ± 0.16</td><td>70%</td></tr><tr><td colspan="3">WideResNet-28-2 95.01 ± 0.04</td><td>86.00 ± 1.09</td><td colspan="2">90.53 ± 0.17</td><td>80%</td></tr><tr><td colspan="3">WideResNet-28-4 95.69 ± 0.10</td><td>95.15 ± 0.11</td><td colspan="2">95.50 ± 0.05</td><td>40%</td></tr><tr><td colspan="3">WideResNet-28-4 95.69 ± 0.10</td><td>94.86 ± 0.10</td><td colspan="2">95.43 ± 0.16</td><td>50%</td></tr><tr><td colspan="3">WideResNet-28-4 95.69 ± 0.10</td><td>94.47 ± 0.10</td><td colspan="2">95.20 ± 0.05</td><td>60%</td></tr><tr><td colspan="3">WideResNet-28-4 95.69 ± 0.10</td><td>93.37 ± 0.18</td><td colspan="2">94.67 ± 0.08</td><td>70%</td></tr><tr><td colspan="3">WideResNet-28-4 95.69 ± 0.10</td><td>91.88 ± 0.59</td><td colspan="2">93.79 ± 0.09</td><td>80%</td></tr><tr><td colspan="3">WideResNet-28-8 96.06 ± 0.06</td><td>95.62 ± 0.04</td><td colspan="2">96.06 ± 0.12</td><td>40%</td></tr><tr><td colspan="3">WideResNet-28-8 96.06 ± 0.06</td><td>95.59 ± 0.09</td><td colspan="2">96.03 ± 0.02</td><td>50%</td></tr><tr><td colspan="3">WideResNet-28-8 96.06 ± 0.06</td><td>95.40 ± 0.14</td><td colspan="2">95.88 ± 0.16</td><td>60%</td></tr><tr><td colspan="3">WideResNet-28-8 96.06 ± 0.06</td><td>94.99 ± 0.22</td><td colspan="2">95.71 ± 0.16</td><td>70%</td></tr><tr><td colspan="3">WideResNet-28-8 96.06 ± 0.06</td><td>94.22 ± 0.21</td><td colspan="2">95.15 ± 0.03</td><td>80%</td></tr><tr><td colspan="5">A.4.2 UNDERSTANDING THE LOTTERY TICKET EFFECT</td></tr><tr><td colspan="6">Similar to the observations in Section 6 (for unstructured pruning), Figure 15 instead considers</td></tr><tr><td colspan="6">structured pruning and again we found DPF does not find a lottery ticket. The superior generalization</td></tr><tr><td colspan="6">performance of DPF cannot be explained by the found mask or the weight initialization scheme.</td></tr><tr><td>Top-1 test accuracy</td><td>91 92 93 94</td><td colspan="4">DPF re-train w/ original init (mask from DPF) re-train w/o original init (mask from DPF)</td></tr><tr><td/><td>0.40</td><td>0.50</td><td colspan="2">0.60 Sparsity ratio</td><td>0.70</td><td>0.80</td></tr><tr><td colspan="6">Figure 15: Investigate the effect of lottery ticket for model compression (WideResNet28-2 with CIFAR-10) for</td></tr><tr><td colspan="6">structured pruning. We retrained the model with the mask from the model trained by DPF, by using the same</td></tr><tr><td>epoch budget.</td><td/><td/><td/><td/></tr><tr><td colspan="4">A.4.3 MODEL SPARSITY VISUALIZATION</td><td/></tr></table>

w t+1 -w 2 | w t = w t -w 2 -2γ t w t -w , Eg( w t ) + γ 2 t E g( w t ) 2 ≤ w t -w 2 -2γ t w t -w , ∇f ( w t ) + γ 2 t G 2 = w t -w 2 -2γ t w t -w , ∇f ( w t ) + γ 2 t G 2 + 2γ t w t -w t , ∇f ( w t ) . By strong convexity, -2 w t -w , ∇f ( w t ) ≤ -µ w t -w 2 -2 (f ( w t ) -f (w )) , + γ t λ t G 2 + 3λ t LE w t -w t 2 .By plugging in the learning rate, γ t = 4 µ(t+2) and setting λ t = (t + 1) we obtainλ t E (f ( w t ) -f (w )) ≤ µ 4 t(t + 1)E w t -w 2 -(t + 1)(t + 2)E w t+1 -wBy summing from t = 0 to t = T these λ t -weighted inequalities, we obtain a telescoping sum:Finally, using w t -w t 2 = δ t w t 2 by (1), shows the theorem.Before giving the proof of Theorem 4.2, we first give a justification for the remark just below Theorem 4.1 on the one-shot pruning of the final iterate.We have by L-smoothness and a, b ≤ 1 2α a 2 + α 2 b 2 for a, b ∈ R d and α > 0 for any iterate w t : f ( w t ) -f (w ) ≤ f (w t ) -f (w ) + ∇f (w t ), w t -w t +

Table content: <table><tr><td colspan="3">Hence, for Λ T :=</td><td/><td colspan="5">T t=0 λ t = (T +1)(T +2) 2</td><td>,</td></tr><tr><td/><td>1 Λ T</td><td colspan="2">T t=0</td><td colspan="6">λ t E (f ( w t ) -f (w )) ≤</td><td>4(T + 1) µΛ T</td><td>G 2 +</td><td>3L Λ T</td><td>T t=0</td><td>λ t E w t -w t</td><td>2</td></tr><tr><td/><td/><td/><td/><td/><td/><td/><td/><td/><td>= O</td><td>G 2 µT</td><td>+</td><td>L Λ T</td><td>T t=0</td><td>λ t E w t -w t</td><td>2</td><td>.</td></tr><tr><td/><td/><td/><td/><td/><td/><td/><td/><td/><td>L 2</td><td>w t -w t</td><td>2</td></tr><tr><td colspan="3">and with a + b</td><td colspan="7">2 ≤ 2 a -w t -w 2 + 2 b ≤ f (w t ) -f (w ) + 2 further 2 ≤ -2 w t -w 1 1 2L ∇f (w t ) 2 + w t -w t</td><td>2</td></tr><tr><td colspan="5">and with a, b ≤ 1 2α a</td><td>2 + α 2 b</td><td/><td/><td/></tr><tr><td/><td/><td/><td/><td/><td/><td/><td/><td/><td>2</td></tr><tr><td/><td/><td/><td/><td/><td colspan="3">≤ 1 -</td><td colspan="2">µγ t 2</td><td>w t -w</td></tr><tr><td/><td/><td/><td/><td/><td/><td colspan="2">γ t</td><td/><td>2)</td><td>E w t -w</td><td>2 -</td><td>λ t γ t</td><td>E w t+1 -w</td></tr><tr><td/><td/><td/><td/><td/><td/><td/><td/><td/><td>2</td></tr><tr><td/><td/><td/><td/><td/><td>+</td><td colspan="4">4(t + 1) µ(t + 2)</td><td>G 2 + 3(t + 1)LE w t -w t</td><td>2 .</td></tr><tr><td>T t=0</td><td colspan="5">λ t E (f ( w t ) -f (w )) ≤</td><td>µ 4</td><td colspan="3">0 -(T + 1)(T + 2)E w t+1 -w</td><td>2 +</td><td>4(T + 1) µ</td><td>G 2</td></tr><tr><td/><td/><td/><td/><td/><td/><td/><td/><td/><td>T</td></tr><tr><td/><td/><td/><td/><td/><td/><td/><td colspan="3">+ 3L</td><td>λ t E w t -w t</td><td>2</td></tr><tr><td/><td/><td/><td/><td/><td/><td/><td/><td/><td>t=0</td></tr><tr><td/><td/><td/><td/><td/><td>≤</td><td colspan="4">4(T + 1) µ</td><td>G 2 + 3L</td></tr></table>

Lastly, we give the proof of Theorem 4.2, followingKarimireddy et al. (2019).E [f (w t+1 ) | w t ] ≤ f (w t ) -γ ∇f (w t ), Eg( w t ) + γ 2 L 2 E g( w t ) (w t ) -γ ∇f (w t ), ∇f ( w t ) + γ 2 LG 2 2 = f (w t ) -γ ∇f ( w t ), ∇f ( w t ) + γ 2 LG 2 2 + γ ∇f ( w t ) -∇f (w t ), ∇f ( w t ) ≤ f (w t ) -γ ∇f ( w t )

Table content: <table><tr><td/><td/><td/><td/><td/><td/><td/><td>= O</td><td>G 2 µT</td><td>+ LE δ t w t</td><td>2</td></tr><tr><td colspan="8">Proof of Theorem 4.2. By smoothness, and a, b ≤ 1 2 a</td><td>2 + 1 2 b</td></tr><tr><td/><td/><td/><td>+</td><td>γ 2</td><td colspan="3">2 + γ 2 LG 2 2 ∇f ( w t ) -∇f (w t ) 2 +</td><td>γ 2</td><td>∇f ( w t )</td></tr><tr><td/><td/><td/><td/><td/><td>γ 2</td><td>∇f ( w t )</td><td>2 + γ 2 LG 2 2</td><td>+</td><td>γL 2 2</td><td>w t -w t</td><td>2 ,</td></tr><tr><td>and by rearranging</td><td/><td/><td/><td/><td/><td/></tr><tr><td>E ∇f ( w t )</td><td>2 ≤</td><td>2 γ</td><td colspan="2">[Ef (w</td><td/><td/></tr></table>

Peer review (iclr2020):
This work proposes a simple pruning method that dynamically sparsifies the network during training. This is achieved by performing at fixed intervals magnitude based pruning for either individual weights or entire neurons. While similar methods have been explored before, this work proposes a slight twist; instead of updating the weights of the model by following the gradient of the parameters of the dense model, they update the parameters of the dense model according to the gradients of the sparse model. Essentially, this corresponds to a variant of the straight-through estimator [1], where in the forward pass we evaluate the compressed model, but in the backward pass we update the model as if the compression didn’t take place. The authors argue that this process allows for ``feedback” in the pruning mechanism, as the pruned weights still receive gradient updates hence they can be ``re-activated” at later stages of training. They then provide a convergence analysis about the optimization procedure with such a gradient, and show that for strongly convex functions the method converges in the vicinity of the global optimum, whereas for non-convex functions it converges to the neighbourhood of a stationary point. Finally, the authors perform extensive experimental evaluation and show that their method is better than the baselines that they considered. This work is in general well written and conveys the main idea in an effective manner. It is also a timely contribution as sparse models / compression are important topics for the deep learning community. The overall method seems simple to implement, doesn’t introduce too many hyper-parameters and seem to work very well. For this reason I tend towards recommending for acceptance, provided that the authors address /comment on a couple of issues I found in the draft. More specifically: - The connection to the error feedback is kind of loose and not well explained. After skimming Karimireddy et al. I noticed that 1. it evaluates the gradient at a point (i.e. the current estimate of the parameters), 2. compresses said gradient, 3. updates the parameters while maintaining the difference of the original w.r.t. the compressed gradient. In this sense, it seems a bit different that DPF, as your notation at the first equation of page 4 implies that you take the gradient of a different point, i.e. w_t + e_t instead of w_t. I believe that expanding a bit more about the connection would help in making the manuscript more clear. - There seems to be a typo / error on your definition of an m-strongly convex function at the “convergence of Convex functions” paragraph. I believe it should be <
abla f(v), w-v> <= f(w) - f(v) - 0.5 m ||w - v||^2, instead of <
abla f(w), w-v> <= f(w) - f(v) - 0.5 m ||w - v||^2. - The proposed gradient estimator seems to be an instance of the STE [1] estimator, that, as the authors mention, has been using at the Binary Connect algorithm. It would be interesting to see some more discussion about this similarity perhaps also expanding upon recent work that discusses the STE gradient as a form of coarse gradient [2]. - At section 5.2 the authors mention that “dynamic pruning methods, and in particular DPF, work on a different paradigm, and can still heavily benefit from fine-tuning”. This claim seems to contradict the results at Figure 4; there it seems that the masks have “converged” in the later stages of training, hence one could argue that the fine-tuning already happens thus it wouldn’t benefit DPF. I believe it would be interesting if the authors provide a similar plot as the one in Figure 4 but rather for the ResNet-20 network on CIFAR 10 (which seems to benefit heavily from FT). Do the masks still settle at the end of training (as it was the case for WideResNet-28-2) and if they do, why is fine-tuning still increasing the accuracy? - Minor: Try to use consistent coloring at Figure 6 as while (a), (b) share the same color-coding, (c) is using a different one hence could be confusing. [1] Estimating or Propagating Gradients Through Stochastic Neurons for Conditional Computation, Yoshua Bengio, Nicholas Léonard, Aaron Courville, 2013 [2] Understanding Straight-Through Estimator in Training Activation Quantized Neural Nets, Penghang Yin, Jiancheng Lyu, Shuai Zhang, Stanley Osher, Yingyong Qi, Jack Xin, 2019