Title: SPATIAL-WINOGRAD PRUNING ENABLING SPARSE WINOGRAD CONVOLUTION

Abstract: Deep convolutional neural networks (CNNs) are deployed in various applications but demand immense computational requirements. Pruning techniques and Winograd convolution are two typical methods to reduce the CNN computation. However, they cannot be directly combined because Winograd transformation fills in the sparsity resulting from pruning. Li et al. (2017) propose sparse Winograd convolution in which weights are directly pruned in the Winograd domain, but this technique is not very practical because Winograd-domain retraining requires low learning rates and hence significantly longer training time. Besides, Liu et al. (2018) move the ReLU function into the Winograd domain, which can help increase the weight sparsity but requires changes in the network structure.
To achieve a high Winograd-domain weight sparsity without changing network structures, we propose a new pruning method, spatial-Winograd pruning. As the first step, spatial-domain weights are pruned in a structured way, which efficiently transfers the spatial-domain sparsity into the Winograd domain and avoids Winograd-domain retraining. For the next step, we also perform pruning and retraining directly in the Winograd domain but propose to use an importance factor matrix to adjust weight importance and weight gradients. This adjustment makes it possible to effectively retrain the pruned Winograd-domain network without changing the network structure. For the three models on the datasets of CIFAR-10, CIFAR-100, and ImageNet, our proposed method can achieve the Winograddomain sparsities of 63%, 50%, and 74%, respectively.

1. INTRODUCTION:

Deep convolutional neural networks (CNNs) have been ubiquitously utilized in various application domains. However, their performance comes at the cost of a significant amount of computation which keeps growing over time. As an example, for the ImageNet challenge (Russakovsky et al., 2015) , Krizhevsky et al. (2012) proposed AlexNet which requires more than 1.1 × 10 9 multiplications. Later, in 2016, the ResNet-152 model (He et al., 2016) increased the computation cost to 11.3 × 10 9 multiplications. This high computation cost limits the deployment of larger and deeper CNN models.

There are two primary methods to reduce the required computation of CNN models: pruning techniques and Winograd/FFT convolution. Pruning removes redundant weight parameters, inducing sparsity into the network. On the other hand, Winograd convolution (Lavin & Gray, 2016) and FFT convolution (Mathieu et al., 2013) transform the computation into different domains. The convolution operations can then be replaced by element-wise multiplications. For the typical convolution kernel size of 3 × 3, Winograd convolution can achieve more than twofold speedup over highly optimized spatial convolution algorithms, and typically requires fewer flops than FFT-based approaches (Li et al., 2017) . Therefore, in this paper, we focus on the Winograd convolution.

The pruning techniques and Winograd convolution are not directly compatible with each other. Sparse weight matrices, which are generated by pruning, lose most of the sparsity after the Winograd transformation from the spatial (original) domain to the Winograd domain. The remaining sparsity is much lower than what we need for improving computation performance.

To increase the Winograd-domain sparsity, Li et al. (2017) propose to perform pruning and retraining directly on Winograd-domain weights. However, it requires using an extremely small learning rate, e.g., 200x smaller for AlexNet, in retraining and is difficult to be applied to deep networks. Besides, Winograd-ReLU pruning (Liu et al., 2018) moves ReLU function into the Winograd domain, which helps increase Winograd-domain sparsity but requires changes in the network structure.

In this paper, to further improve the sparsity of Winograd-domain weights without changing the network structure, we propose a new pruning method, spatial-Winograd pruning. It includes two parts: spatial structured pruning and Winograd direct pruning. In spatial structured pruning, we prune the spatial-domain weights in a structured way, in which the structures are designed to transfer the spatial-domain sparsity into the Winograd domain efficiently. After spatial structured pruning, weights of the pruned layers will be converted to and kept in the Winograd domain. Then, for Winograd direct pruning, we perform pruning and retraining entirely in the Winograd domain to improve the sparsity further. This paper makes the following contributions:

• We propose a new pruning method, spatial-Winograd pruning. Without changing the network structure, it can achieve higher sparsity in Winograd-domain weights compared with previous methods. • As the first part of spatial-Winograd pruning, we provide a structured pruning method to transfer the spatial-domain sparsity into the Winograd domain efficiently. It can help avoid Winograd-domain retraining in this part and accelerate the pruning process. • In the second part, to perform pruning directly in the Winograd domain, we present a new approach to measuring the importance of each Winograd-domain weight based on its impact on output activations. Also, we propose to use an importance factor matrix to adjust the gradients of Winograd-domain weights, which makes it much faster to retrain deep networks directly in the Winograd domain without changing the network structure.

2. PRELIMINARY AND RELATED WORK:

Winograd convolution (Lavin & Gray, 2016 ) is a typical algorithm to reduce the arithmetic complexity of CNNs. It transforms the computation into the Winograd domain, and the convolution operations can then be replaced by element-wise multiplications. We call the domain, in which the conventional convolution operation is executed, to be the spatial domain.

The basic block of Winograd convolution works on a 2D input tile, I, with a size of m × m and a 2D weight filter, W , with a size of n × n. In this case, the 2D output tile generated, O, will have a size of (m -n + 1) × (m -n + 1). For a typical convolutional layer, the input feature maps are first disassembled into input tiles and, after the Winograd convolution, the output tiles will be reassembled into the output feature maps.

Figure 1a shows how conventional Winograd convolution works. As the first step, the weight filter W and the input tile I are converted into the Winograd domain using the predefined matrices G and B. Element-wise multiplication is then applied to the Winograd-domain weight filter, GW G , and input tile, B I B, to generate the Winograd-domain output tile with a size of m × m. In the last step, the output tile is converted back into the spatial domain with another predefined matrix A.

With as the Hadamard product (element-wise multiplication), the entire process can be written as The transform/inverse-transform matrices A, B and G are only determined by m and n. These matrices contain many repeating elements and applying them requires only few multiplications. In this case, considering only the element-wise multiplication between GW G and B I B, the Winograd convolution can reduce the number of multiplications from (m -n + 1) 2 n 2 to m 2 .

O = A [(GW G ) (B I B)]A (1)

In addition to Winograd convolution, pruning is also a well-explored method to reduce CNN computation. Han et al. (2015b; a) propose to perform pruning and retraining iteratively, which can help reduce the computation by up to 5×. To fully utilize the sparsity incurred by pruning to accelerate CNN computation, Wen et al. (2016) and Yu et al. (2017) prune networks in a structured way: weights are clustered into groups with hardware-friendly structures and then get pruned in groups.

However, Winograd convolution is not directly compatible with conventional pruning algorithms.

The transformation GW G fills in the zeros in the sparse weight filters generated by pruning.

There have been several research attempts to solve this problem. Liu & Turakhia (2016) propose to directly mask out Winograd-domain weights and use backpropagation to train the spatial-domain weights. However, compared with spatial-domain weights, Winograd-domain weights are in a higher-dimensional space. Directly setting Winograd-domain weights to zero will cause an inconsistency between the spatial domain and the Winograd domain. This inconsistency will lead to a significant accuracy loss or a low sparsity on networks, e.g., AlexNet, for large datasets (Li et al., 2017) .

To address the inconsistency between the spatial and Winograd domain, Li et al. (2017) 

EQUATION

The Winograd-domain inputs also become sparse, which helps further reduce the required computation. Besides, a higher sparsity in Winograd-domain weights can be achieved. However, with weight filters being sparse, it is challenging to utilize both the weight and input sparsity for CNN acceleration on general-purpose processors due to more irregularity in the access pattern and control flow. Also, Winograd-ReLU pruning cannot be applied to conventional CNN models since the new computation in Equation 2does not correspond to the original convolution operation. It requires changing the network structure and retraining the network from scratch.

3. SPATIAL-WINOGRAD PRUNING:

In this paper, to achieve a high Winograd-domain weight sparsity on deep CNN models without changing network structures, we propose the spatial-Winograd pruning. As shown in Figure 2 , it consists of two parts: spatial structured pruning and Winograd direct pruning.

In spatial structured pruning, spatial-domain weights are pruned in a structured way and then retrained to regain the original accuracy. After spatial structured pruning, the weights of the pruned model will be transferred into and kept in the Winograd domain. The Winograd direct pruning then performs pruning and retraining directly onto the weights in the Winograd domain. The pruning and retraining steps in both spatial structured pruning and Winograd direct pruning will be iteratively executed until we achieve the desired sparsity or the produced model loses much accuracy.

3.1. SPATIAL STRUCTURED PRUNING:

The first part of the spatial-Winograd pruning is spatial structured pruning. Spatial-domain weights which affect the same Winograd-domain weight are clustered into the same group. Less important weight groups are removed, and the pruned network will be retrained to regain the accuracy. This structured pruning method can help transfer more spatial-domain sparsity into the Winograd domain.

Spatial Pruning For each spatial-domain filter W , we need to generate a mask M spatial to indicate the redundant weights. Assuming W has a size of n × n and the Winograd-

domain filter Q is m × m, we have Q = GW G . Each element of the Winograd-domain filter, Q i,j

, is the weighted sum of the spatial-domain weights:

Q i,j = 0 u,v n-1 (S i,j,u,v • W u,v ) 0 i, j m -1 (3)

where S is a 4D tensor containing the weight coefficients of the spatial-domain weights and is only determined by m and n. Details about the calculation of S can be found in Appendix A.1.

For each Winograd-domain weight Q i,j , we can create a set D i,j containing the spatial-domain weights which affect the value of Q i,j . D i,j is defined as

EQUATION

In this case, for each weight group D i,j , we use a function h(D i,j ) to measure its importance. In this paper, we use the maximum norm function as h

EQUATION

With a specific threshold t spatial , if h(D i,j ) < t spatial , then D i,j is considered as redundant and all weights included need to be removed. In this case, the corresponding Q i,j will be fixed to 0 and also removed. The set of redundant weights for entire W is the union of all redundant D i,j and can be calculated as

EQUATION

Here we define D i,j in a structured way based on the relation between spatial-domain weights and Winograd-domain weights. It helps transfer as much spatial-domain sparsity into Winograd-domain sparsity as possible.

The mask matrix M spatial can be generated by

M spatial u,v = 0 W u,v ∈ D 1 W u,v / ∈ D 0 u, v n -1 (7)

Spatial Retraining After spatial pruning, we can perform the spatial retraining with conventional training algorithms, e.g., stochastic gradient descent (SGD). The removed weights are fixed to 0 by applying W = W M spatial after each training iteration. is element-wise multiplication. The steps of spatial pruning and spatial retraining will be iteratively performed until the retrained model loses much accuracy. The threshold t spatial is gradually increased to incur more sparsity into the network.

In spatial structured pruning, both pruning and retraining steps are performed in the spatial domain.

It helps avoid the Winograd-domain retraining to accelerate the pruning process but, at the same time, incurs high Winograd-domain sparsity.

3.2. WINOGRAD DIRECT PRUNING:

After spatial structured pruning, as in sparse Winograd convolution, weights of the pruned model will be transferred into and kept in the Winograd domain. In Winograd direct pruning, we measure the importance of each weight based on its impact on output activations, and unimportant weights are removed. The pruned network is then retrained in the Winograd domain, and an importance factor matrix is deployed to adjust the weight gradients.

Winograd Pruning Similar to spatial pruning, in Winograd pruning, we need to generate a mask matrix M W inograd for each Winograd-domain filter Q to indicate the redundant weights. With the weight filter Q in the Winograd domain, the output tile O is calculated as

EQUATION

Each output element can be considered as the weighted sum of the products of weights and inputs

O x,y = 0 i,j,s,t m-1 (H x,y,i,j,s,t • Q i,j • I s,t ) 0 x, y m -n (9)

where H is a 6D tensor containing the weight coefficients of different products (Q i,j • I s,t ) and is only determined by m and n. Details about the calculation of H can be found in Appendix A.2.

By removing one weight Q i,j , the change on each output O x,y is

EQUATION

In Winograd pruning, we need to remove a certain amount of weights while minimizing the change of the output activations ||∆O|| 2 . Removing an important weight will lead to a larger change in output activations. Therefore, we propose to measure the importance of each weight Q i,j by the expected value of ||∆O| Qi,j || 2 2 . In this case, we have

EQUATION

For simplicity, we can assume input values are independent and identically distributed (i.i.d.), and have expected values of 0. With this assumption, we have

EQUATION

Since the importance of weights are relative numbers, we can assume E(I 2 s,t ) = 1. In this case,

EQUATION

Based on Equation. 13, we can generate an importance factor matrix F , where

F i,j = E(||∆O| Qi,j || 2 2 ) Q 2 i,j = 0 x,y m-n,0 s,t m-1 H 2 x,y,i,j,s,t 0 i, j m -1 (14)

Therefore, F is only determined by m and n, and keeps the same for all 2D Winograd-domain filters Q in a specific layer. Then Equation. 13 can be simplified to

EQUATION

In this case, with a specific threshold t W inograd , we can generate the mask matrix M W inograd as

EQUATION

For a specific weight Q i,j , conventional pruning algorithms (Han et al., 2015b; Guo et al., 2016 ) use its absolute value |Q i,j | as the weight importance, which is equivalent to using Q 2 i,j . Therefore, in Equation 16, the employed weight importance, Q 2 i,j •F 2 i,j , can be considered as using the importance factor matrix F to adjust the conventional weight importance Q 2 i,j . Winograd Retraining As the same with the spatial retraining, we fix the removed Winograddomain weights to 0 by applying Q = Q M W inograd after each training iteration.

However, using conventional SGD to retrain the Winograd-domain parameters will lead to divergence. This is because, as shown in Equation. 14, different locations of Winograd-domain weights have different importance and, therefore, require different learning speeds. Using an extremely small learning rate can avoid the divergence but makes the retraining much slower.

To address this problem, in Winograd retraining, we propose to adjust the gradients of Winograddomain weights with the importance factor matrix F . Assume loss to be the loss value. At the training step k, after the backward computation, the gradients of Q, ∂ loss ∂Q | k , will be adjusted by

( ∂ loss ∂Q k ) adjusted = ∂ loss ∂Q k F •α (17)

where and •α are the Hadamard division (element-wise division) and Hadamard power (elementwise power of α) function, respectively. In this paper, based on empirical results, α is fixed to 1.5. In this case, with the learning rate of η, the SGD update for the Winograd-domain weights Q at the training step k becomes

EQUATION

4. EXPERIMENTS:

To evaluate the spatial-Winograd pruning, we perform the experiments on three datasets: CIFAR-10, CIFAR-100 (Krizhevsky, 2009) and ImageNet (ILSVRC-2012) (Russakovsky et al., 2015) . Py-Torch (Paszke et al., 2017) is used to implement the pruning framework.

We use the Winograd-ReLU pruning (Liu et al., 2018) as the baseline pruning technique. To show the effectiveness of our proposed method, we test the same models as in Winograd-ReLU pruning: VGG-nagadomi (Nagadomi, 2014), ConvPool-CNN-C (Springenberg et al., 2014) and ResNet-18 (He et al., 2016) on the three datasets tested, respectively. Those models are chosen since the majority of the included convolutional layers use 3 × 3 kernels.

For 3×3 kernels, we set the input tile size m to 6 instead of 4. A larger input tile size can help achieve higher computation speedup. With our proposed method, we expect that lower input tile sizes can lead to a similar or higher sparsity. This is because, with lower input tile sizes, the spatial-domain weights have less correlation between each other and the spatial structured pruning can achieve a higher sparsity.

4.1. CIFAR-10: VGG-NAGADOMI:

For the CIFAR-10 dataset, we test the VGG-nagadomi model (Nagadomi, 2014) . It contains 8 convolutional layers with 3 × 3 kernels. We use batch normalization instead of dropout to regularize the convolutional layers. The original model has a prediction accuracy of 93.96%. We prune the first convolutional layer with a fixed Winograd-domain sparsity of 20%. For the remaining convolutional layers, we incur a uniform Winograd-domain sparsity, increasing from 20% to 80%, for simplicity.

Figure 3a shows the pruning results. The baseline result reported in (Liu et al., 2018) prediction accuracy of 69.95%. We prune the first convolutional layer with a fixed Winograd-domain sparsity of 20%. Similar to the VGG-nagadomi model, the remaining 6 convolutional layers with 3 × 3 kernels are iteratively pruned and retrained with uniform Winograd-domain sparsities.

Figure . 3b shows the result of the relative accuracy against the Winograd-domain sparsity. The baseline result reported in (Liu et al., 2018) is shown as the dashed line. With <0.1% accuracy loss, it achieves a sparsity of 40%. Winograd direct pruning is applied to the model pruned by spatial structured pruning with 30% sparsity. With no accuracy loss, spatial-Winograd pruning can reach a sparsity of 50%, which is 10% higher than Winograd-ReLU pruning.

4.3. IMAGENET: RESNET-18:

We test the ResNet-18 model on the ImageNet (ILSVRC-2012) dataset. As the same with Winograd-ReLU pruning, we replace each 2×2-stride 3×3 convolutional layer with a 2×2-stride max-pooling layer followed by a 1 × 1-stride 3 × 3 convolutional layer. This change makes it easier to apply Winograd convolution on most of the convolutional layers.

The original model has a top-1/top-5 prediction accuracy of 69.82%/89.55%. However, for Winograd-ReLU pruning, Liu et al. (2018) use the model with the original top-1/top-5 accuracy of only 66.67%/87.42%. Despite this, we still use the relative accuracies reported in (Liu et al., 2018) as the baseline.

We prune the 16 convolutional layers in the residual blocks with the same Winograd-domain sparsity.

The first convolutional layer and the downsample layers are kept intact. Figure . 4 shows the results of the relative accuracy against the Winograd-domain sparsity. As the dashed line show, the Winograd-ReLU pruning achieves a sparsity of 70%/65% with <0.1% top-1/top-5 accuracy loss.

We apply Winograd direct pruning to the models pruned by spatial structured pruning with 65% and 70% Winograd-domain sparsity, annotated as Winograd direct pruning -0.65 and 0.70, respectively. As shown in the figure, with <0.1% top-1/top-5 accuracy loss, applying Winograd direct pruning to the model with 70% Winograd-domain sparsity can achieve a higher sparsity of 74%/72%. This is because, with the sparsity increasing, Winograd direct pruning makes the prediction accuracy drop much faster than spatial structured pruning. Although we can use the importance factor matrix F to adjust the weight gradients to accelerate the Winograd-domain retraining, the learning rate still needs to be much lower than in the spatial retraining. In this case, the accuracy loss recovered through Winograd retraining is limited, which makes the accuracy drop much faster when applying Winograd direct pruning.

4.4. EFFECTIVENESS OF IMPORTANCE FACTOR MATRIX:

In Winograd direct pruning, we use the importance factor matrix F to adjust the weight importance and gradients for different locations of Winograd-domain weights. Here we test the effectiveness of employing the importance factor matrix in both the Winograd pruning and retraining.

We first test how the importance factor matrix F helps in Winograd pruning. Winograd pruning without retraining is applied to the model pruned by spatial structured pruning with 70% sparsity. Figure 5a shows the relative accuracy against the sparsity when pruning with weight importance unadjusted or adjusted with F . As shown in the figure, adjusting the weight importance with the importance factor matrix can dramatically reduce the accuracy loss when performing Winograd pruning. When pruning the model to 76% sparsity, using the absolute value as the weight importance will cause a 22% accuracy loss. In comparison, using the importance factor matrix to adjust the weight importance can help reduce the accuracy loss to 10%.

We also test the effectiveness of the importance factor matrix in Winograd retraining. For the model pruned with spatial structured pruning (70% sparsity), Winograd pruning is applied to increase the sparsity to 74%. We then perform Winograd retraining for 10 epochs. Figure 5b shows the relative accuracy against the retraining epochs with unadjusted and adjusted gradients. For unadjusted gradients, we try three learning rates of 1e-7, 1e-8 and 1e-9. Higher learning rates, e.g., 1e-6, will lead to an accuracy drop through retraining. As shown in the figure, adjusting the gradients with the importance factor matrix can substantially accelerate the convergence. With retraining of only 10 epochs, it reduces the accuracy loss to 0.2% while retraining without gradient adjustment only reduces the accuracy loss to 0.7%.

5. CONCLUSION:

In this paper, we present a new pruning method, spatial-Winograd pruning, to improve the Winograd-domain weight sparsity without changing network structures. It includes two steps: spatial structured pruning and Winograd direct pruning. In spatial structured pruning, we prune the spatial-domain weights based on the internal structure in the Winograd transformation. It can help efficiently transfer the spatial-domain sparsity into the Winograd domain. For Winograd direct pruning, we perform both pruning and retraining in the Winograd domain. An importance factor matrix is proposed to adjust the weight gradients in Winograd retraining, which makes it possible to effectively retrain the Winograd-domain network to regain the original accuracy without changing the network structure. We evaluate spatial-Winograd pruning on three datasets, CIFAR-10, CIFAR-100, ImageNet, and it can achieve the Winograd-domain sparsities of 63%, 50%, and 74%, respectively.

A COEFFICIENT TENSORS S AND H:

The coefficient tensors S and H contain the weight coefficients in Equation 3 and Equation 9. They are only determined by the input tile size m and the weight filter size n.

A.1 COEFFICIENT TENSOR S:

To calculate the tensor S, we first introduce an equivalent transformation: for two vectors a and b with a size of m, and a matrix C with a size of m × m, we have

a Cb = 1 [(ab ) C] 1 (19)

where 1 is a vector of size m and all its entries are 1.

For matrix S, with the weight transform matrix G, we have Q = GW G . Each element in Winograd-domain weight filter Q is calculated as

EQUATION

where 0 i, j m -1. In this case, compared with Equation 3, each element in the coefficient tensor S can be calculated as

EQUATION

where 0 i, j m -1, 0 u, v n -1.

A.2 COEFFICIENT TENSOR H:

With the Winograd-domain weight filter Q, the output tile O is calculated as

EQUATION

Each element O x,y is calculated as

EQUATION

where 0 x, y m -n. Based on Equation 19, we have

EQUATION

where 0 i, j m -1. The element (B I B) i,j can be calculated as

EQUATION

Based on Equation 24, 25 and 26, we have Therefore, compared with Equation 9, each element in the coefficient tensor H is calculated as

EQUATION

EQUATION

where 0 x, y m -n, 0 i, j, s, t m -1.

B RESNET-18 PRUNING WITH VARIED SPARSITIES ACROSS LAYERS:

In addition to pruning ResNet-18 with the same sparsity across all targeting layers, we experiment incurring different sparsities into different layers with spatial-Winograd pruning.

For spatial structured pruning, we first test the pruning sensitivity of each convolutional layer to decide which layers need to be pruned and the corresponding thresholds. To choose the targeting layers, we measure the accuracy loss when 60% of Winograd-domain weights are pruned for each layer. Only one layer is pruned at one time, and other layers are kept intact. Figure . 6 shows the results. In ResNet-18, the i-th residual block contains two convolutional layers, i-a and i-b. As shown in the figure, the first layer in each residual block is much more sensitive to pruning than the second layer. Therefore, we will only prune the second convolutional layer, i-b, in each residual block.

For each targeting layer i-b, we determine the corresponding pruning threshold t spatial i based on its pruning sensitivity. We gradually increase the threshold until the validation accuracy drops by 2% and the threshold is recorded as t spatial, 2% loss i . Then in spatial structured pruning, we can calculate the threshold used for layer i-b as t spatial i = β • t spatial, 2% loss i where β is a multiplier shared across all targeting layers. With a larger β, the threshold and, therefore, the sparsity will be higher. Also, in Winograd direct pruning, we use the same strategy to choose the thresholds used for different layers. 

C SPARSITY DISTRIBUTION:

For the pruned ResNet-18 model, we analyze more detailed sparsity distribution across and inside 2D weight filters. Here each Winograd-domain weight matrix Q is considered as a 2D filter. We use the last convolutional layer (7-b) as an example. The model with a uniform sparsity of 74% across all pruned layers, which corresponds to point P in Figure 4a , is tested. Figure 7 shows the sparsity distribution across the filters. As shown in the figure, more than half (62%) of the filters have all weights removed. An interesting observation is that a large portion of the filters have exact 20 weights removed.

To explain why many filters have exact 20 weights removed, we visualize the sparsity distribution inside the filters. Figure 8a shows the sparsity of different locations for the filters with 20 weights removed. Darker locations have higher sparsities where more weights are removed. The border part of the 6 × 6 filter, which includes 20 weights, has much higher sparsity than the central part.

It means the border part of the Winograd-domain weights is much less important than weights in the central part. A potential reason is that the weights in the central part are correlated to more spatial-domain weights and, therefore, removing them will lead to a larger difference in the output activations. In Figure 8b , we also visualize the sparsity distribution inside the filters with at least one weight remaining, and it shows a similar pattern.

:



Figures and Tables:

Figure 1: Conventional Winograd convolution and sparse Winograd convolution (m = 4, n = 3).

Figure 2: Overview of the spatial-Winograd pruning.

ConvPool-CNN-C on CIFAR-100.

Figure 3: Pruning of (a) VGG-nagadomi on CIFAR-10 (b) ConvPool-CNN-C on CIFAR-100 with uniform sparsity across layers.

Figure 4: Pruning of ResNet-18 on ImageNet with uniform sparsity across the pruned layers.

retraining: unadjusted or adjusted gradients for different locations..

Figure 5: Effectiveness of employing importance factor matrix F in (a) Winograd pruning and (b) Winograd retraining.

Figure 7: Distribution of 2D filters with different numbers of weights pruned.

The sparsity for the pruned convolutional layers when pruning the second convolutional layer in each residual block of ResNet-18.

Table content: <table><tr><td colspan="4">Figure 6: Accuracy loss of ResNet-18 when incurring 60% Winograd-domain sparsity into different</td></tr><tr><td colspan="4">layers. Spatial structured pruning is applied with no retraining.</td></tr><tr><td/><td colspan="2">Spatial Structured</td><td>Winograd Direct</td></tr><tr><td>Layer</td><td colspan="2">Pruning Winograd Sparsity Corr. Spatial Sparsity</td><td>Pruning Winograd Sparsity</td></tr><tr><td>0-b</td><td>78.0 %</td><td>88.8 %</td><td>85.0 %</td></tr><tr><td>1-b</td><td>77.1 %</td><td>87.9 %</td><td>84.0 %</td></tr><tr><td>2-b</td><td>76.4 %</td><td>88.5 %</td><td>84.4 %</td></tr><tr><td>3-b</td><td>85.2 %</td><td>93.1 %</td><td>90.4 %</td></tr><tr><td>4-b</td><td>76.9 %</td><td>88.9 %</td><td>84.7 %</td></tr><tr><td>5-b</td><td>88.6 %</td><td>95.1 %</td><td>93.0 %</td></tr><tr><td>6-b</td><td>74.4 %</td><td>88.3 %</td><td>82.4 %</td></tr><tr><td>7-b</td><td>82.6 %</td><td>87.8 %</td><td>92.3 %</td></tr><tr><td>Average</td><td>79.4 %</td><td>88.9 %</td><td>87.6 %</td></tr><tr><td>Top-1 Acc.</td><td colspan="2">69.92 %</td><td>69.94 %</td></tr><tr><td>Top-5 Acc.</td><td colspan="2">89.34 %</td><td>89.51 %</td></tr></table>

lists the pruning results. After spatial structured pruning, we can reach an average Winograddomain sparsity of 79.4% for the pruned layers. The corresponding spatial-domain sparsity is 88.9% which is 9.5% higher. Winograd direct pruning can further improve the Winograd-domain sparsity to 87.6% and layer 5-b has the highest sparsity of 93.0%.

Table content: <table/>

Peer review (iclr2019):
 Review of Spatial-Winograd Pruning Enabling Sparse Winograd Convolution Summary: In this submission, the authors propose a new method for pruning weights in the presence in CNNs in which the convolution is expressed as a CNN. The goal of the project is to demonstrate that they can achieve a network which contains fewer weights (and runs faster) then the original network with minimal sacrifice in network performance and without altering the network architecture. Major Comments: My largest comments and concerns revolve around the degree to which the proposed pruning methods results in a model that is applicable for real world devices. 1. At the minimum, the authors should provide a table with the number of parameters in the (a) original networks, (b) network with baseline pruning method and (c) network with their pruning methods. (Are the savings entirely in convolutional filters or are there savings from fully connected layers?) 2. If the authors are indeed arguing that a goal of network pruning is to perform faster inference, then results must be shown to justify this -- as it is not obvious that speed-ups could be achieved by just pruning weights. In the case of other methods, speed-ups may be achieved by selectively pruning channel filters (as opposed to spatial positions in the convolutional filter) or pruning fully connected layers (in fact, for VGG and AlexNet, a majority of the reduction in parameters due to pruning were found in these layers). 3. Are these levels of sparsity useable in a real-world system? Often the degree of sparsity in "vanilla" CNNs are at levels (e.g. in AlexNet, 30-50% sparsity in higher layers) not high enough to be harnessed for a fast implementation. Selectively zero-ing out individual spatial components of a filter might reduce the parameter count but not achieve any speed up gains. That is, the sparsity must be structured in such a way as to permit a faster implementation [e.g. 1] 4. Considering that one of the baselines changes the network architecture itself [2], I would be curious to understand how effective this method is versus other, simple baselines that change the network architecture such as (1) decreasing the number of filter, (2) decreasing the kernel sizes, (3) swapping a convolutional filter for a separable convolution [3, 4]. All of these baselines are simple to train and experiment with and a practitioner would probably consider in many situations where speed or # parameters are a constraint. Minor Comments: - It is unclear from the presentation whether both proposed pruning methods may be trained in tandem or in series. Please clarify in the manuscript. - Why does Figure 3a, 3b focus on maintaining 20% sparsity on the 1st layer of the network systematically all other layers in sparsity? What is special about the first layer? - Why do the authors explore a different range of sparsities in Figure 3a and Figure 3b? - The authors should discuss the source of the variability (and non-monotonicity) in the plots in Figure 3 and 4 for their proposed method. How are we to interpret this? Naively, it would be appear that the method is somewhat unpredictable in performance across a range of sparsity. - Why do the blue and purple curves in Figure 4 not space the entire range of sparsities? - Figure 5b. What is the relative accuracy measured with respect to? The baseline model at that particular epoch or at final asymptotic performance? [1] Outrageously Large Neural Networks: The Sparsely-Gated Mixture-of-Experts Layer https://openreview.net/forum?id=B1ckMDqlg [2] Efficient sparse-winograd convolutional neural networks. Xingyu Liu, Jeff Pool, Song Han, and William J Dally. [3] MobileNets: Efficient Convolutional Neural Networks for Mobile Vision Applications Andrew G. Howard, Menglong Zhu, Bo Chen, Dmitry Kalenichenko, Weijun Wang, Tobias Weyand, Marco Andreetto, Hartwig Adam https://arxiv.org/abs/1704.04861 [4] Xception: Deep Learning with Depthwise Separable Convolutions François Chollet https://arxiv.org/abs/1610.02357