Title: PER-TENSOR FIXED-POINT QUANTIZATION OF THE BACK-PROPAGATION ALGORITHM

Abstract: The high computational and parameter complexity of neural networks makes their training very slow and difficult to deploy on energy and storage-constrained computing systems. Many network complexity reduction techniques have been proposed including fixed-point implementation. However, a systematic approach for designing full fixed-point training and inference of deep neural networks remains elusive. We describe a precision assignment methodology for neural network training in which all network parameters, i.e., activations and weights in the feedforward path, gradients and weight accumulators in the feedback path, are assigned close to minimal precision. The precision assignment is derived analytically and enables tracking the convergence behavior of the full precision training, known to converge a priori. Thus, our work leads to a systematic methodology of determining suitable precision for fixed-point training. The near optimality (minimality) of the resulting precision assignment is validated empirically for four networks on the CIFAR-10, CIFAR-100, and SVHN datasets. The complexity reduction arising from our approach is compared with other fixed-point neural network designs.
We consider a L-layer DNN deployed on a M -class classification task using the setup in Figure 1 . We denote the precision configuration as the L ×
whose l th row consists of the precision (in bits) of weight W l (B W l ), activation A l (B A l ), weight

1. INTRODUCTION:

Though deep neural networks (DNNs) have established themselves as powerful predictive models achieving human-level accuracy on many machine learning tasks (He et al., 2016) , their excellent performance has been achieved at the expense of a very high computational and parameter complexity. For instance, AlexNet (Krizhevsky et al., 2012) requires over 800 × 10 6 multiply-accumulates (MACs) per image and has 60 million parameters, while Deepface (Taigman et al., 2014) requires over 500 × 10 6 MACs/image and involves more than 120 million parameters. DNNs' enormous computational and parameter complexity leads to high energy consumption (Chen et al., 2017) , makes their training via the stochastic gradient descent (SGD) algorithm very slow often requiring hours and days (Goyal et al., 2017) , and inhibits their deployment on energy and resource-constrained platforms such as mobile devices and autonomous agents.

A fundamental problem contributing to the high computational and parameter complexity of DNNs is their realization using 32-b floating-point (FL) arithmetic in GPUs and CPUs. Reduced-precision representations such as quantized FL (QFL) and fixed-point (FX) have been employed in various combinations to both training and inference. Many employ FX during inference but train in FL, e.g., fully binarized neural networks (Hubara et al., 2016) use 1-b FX in the forward inference path but the network is trained in 32-b FL. Similarly, Gupta et al. (2015) employs 16-b FX for all tensors except for the internal accumulators which use 32-b FL, and 3-level QFL gradients were employed (Wen et al., 2017; Alistarh et al., 2017) to accelerate training in a distributed setting. Note that while QFL reduces storage and communication costs, it does not reduce the computational complexity as the arithmetic remains in 32-b FL.

Thus, none of the previous works address the fundamental problem of realizing true fixed-point DNN training, i.e., an SGD algorithm in which all parameters/variables and all computations are implemented in FX with minimum precision required to guarantee the network's inference/prediction accuracy and training convergence. The reasons for this gap are numerous including: 1) quantization Step 1: Forward Propagation

Step 2: Back Propagation

Step 3: Update ).

errors propagate to the network output thereby directly affecting its accuracy (Lin et al., 2016) ; 2) precision requirements of different variables in a network are interdependent and involve hard-toquantify trade-offs (Sakr et al., 2017) ; 3) proper quantization requires the knowledge of the dynamic range which may not be available (Pascanu et al., 2013) ; and 4) quantization errors may accumulate during training and can lead to stability issues (Gupta et al., 2015) .

Our work makes a major advance in closing this gap by proposing a systematic methodology to obtain close-to-minimum per-layer precision requirements of an FX network that guarantees statistical similarity with full precision training. In particular, we jointly address the challenges of quantization noise, inter-layer and intra-layer precision trade-offs, dynamic range, and stability. As in (Sakr et al., 2017) , we do assume that a fully-trained baseline FL network exists and one can observe its learning behavior. While, in principle, such assumption requires extra FL computation prior to FX training, it is to be noted that much of training is done in FL anyway. For instance, FL training is used in order to establish benchmarking baselines such as AlexNet (Krizhevsky et al., 2012) , VGG-Net (Simonyan and Zisserman, 2014) , and ResNet (He et al., 2016) , to name a few. Even if that is not the case, in practice, this assumption can be accounted for via a warm-up FL training on a small held-out portion of the dataset (Dwork et al., 2015) .

Applying our methodology to three benchmarks reveals several lessons. First and foremost, our work shows that it is possible to FX quantize all variables including back-propagated gradients even though their dynamic range is unknown (Köster et al., 2017) . Second, we find that the per-layer weight precision requirements decrease from the input to the output while those of the activation gradients and weight accumulators increase. Furthermore, the precision requirements for residual networks are found to be uniform across layers. Finally, hyper-precision reduction techniques such as weight and activation binarization (Hubara et al., 2016) or gradient ternarization (Wen et al., 2017) are not as efficient as our methodology since these do not address the fundamental problem of realizing true fixed-point DNN training.

We demonstrate FX training on three deep learning benchmarks (CIFAR-10, CIFAR-100, SVHN) achieving high fidelity to our FL baseline in that we observe no loss of accuracy higher then 0.56% in all of our experiments. Our precision assignment is further shown to be within 1-b per-tensor of the minimum. We show that our precision assignment methodology reduces representational, computational, and communication costs of training by up to 6×, 8×, and 4×, respectively, compared to the FL baseline and related works.

2. PROBLEM SETUP, NOTATION, AND METRICS:

gradient

G (W ) l (B G (W ) l ), activation gradient G (A) l+1 (B G (A) l+1

), and internal weight accumulator

W (acc) l (B W (acc) l

) tensors at layer l. This DNN quantization setup is summarized in Appendix A.

2.1. FIXED-POINT CONSTRAINTS & DEFINITIONS:

We present definitions/constraints related to fixed-point arithmetic based on the design of fixed-point adaptive filters and signal processing systems (Parhi, 2007) :

• A signed fixed-point scalar a with precision B A and binary representation

R A = (a 0 , a 1 , . . . , a B A -1 ) ∈ {0, 1} B A is equal to: a = r A -a 0 + B A -1 i=1

2 -i a i , where r A is the predetermined dynamic range (PDR) of a. The PDR is constrained to be a constant power of 2 to minimize hardware overhead.

• An unsigned fixed-point scalar a with precision B A and binary representation

R A = (a 0 , a 1 , . . . , a B A -1 ) ∈ {0, 1} B A is equal to: a = r A B A -1 i=0 2 -i a i . • A fixed-point scalar a is called normalized if r A = 1. • The precision B A is determined as: B A = log 2 r A ∆ A + 1

, where ∆ A is the quantization step size which is the value of the least significant bit (LSB).

• An additive model for quantization is assumed: a = ã + q a , where a is the fixed-point number obtained by quantizing the floating-point scalar ã, q a is a random variable uniformly distributed on the interval -∆ A 2 , ∆ A 2 , and the quantization noise variance is V ar(q a ) =

∆ 2 A 12

. The notion of quantization noise is most useful when there is limited knowledge of the distribution of ã.

• The relative quantization bias η A is the offset:

η A = |∆ A -µ A | µ A

, where the first unbiased quantization level

µ A = E ã ã ∈ I 1 and I 1 = ∆ A 2 , 3∆ A

2:

. The notion of quantization bias is useful when there is some knowledge of the distribution of ã.

• The reflected quantization noise variance from a tensor T to a scalar α = f (T ), for an arbitrary function f (), is :

V T →α = E T →α ∆ 2 T 12

, where ∆ T is the quantization step of T and E T →α is the quantization noise gain from T to α.

• The clipping rate β T of a tensor T is the probability: β T = Pr ({|t| ≥ r T : t ∈ T }), where r T is the PDR of T .

2.2. COMPLEXITY METRICS:

We use a set of metrics inspired by those introduced by Sakr et al. (2017) which have also been used by Wu et al. (2018a) . These metrics are algorithmic in nature which makes them easily reproducible.

• Representational Cost for weights (C W ) and activations (C A ):

C W = L l=1 |W l | B W l + B G (W ) l + B W (acc) l & C A = L l=1 |A l | B A l + B G (A) l+1

, which equals the total number of bits needed to represent the weights, weight gradients, and internal weight accumulators (C W ), and those for activations and activation gradients (

C A ). 1 • Computational Cost of training: C M = L l=1 |A l+1 | D l B W l B A l + B W l B G (A) l+1 + B A l B G (A) l+1

, where D l is the dimensionality of the dot product needed to compute one output activation at layer l. This cost is a measure of the number of 1-b full adders (FAs) utilized for all multiplications in one back-prop iteration.

2 • Communication Cost: C C = L l=1 |W l | B G (W ) l

, which represents cost of communicating weight gradients in a distributed setting (Wen et al., 2017; Alistarh et al., 2017) .

3. PRECISION ASSIGNMENT METHODOLOGY AND ANALYSIS:

We aim to obtain a minimal or close-to-minimal precision configuration C o of a FX network such that the mismatch probability p m = Pr{ Ŷfl = Ŷfx } between its predicted label ( Ŷfx ) and that of an associated FL network ( Ŷfl ) is bounded, and the convergence behavior of the two networks is similar.

Hence, we require that: (1) all quantization noise sources in the forward path contribute identically to the mismatch budget p m (Sakr et al., 2017) , (2) the gradients be properly clipped in order to limit the dynamic range (Pascanu et al., 2013) , (3) the accumulation of quantization noise bias in the weight updates be limited (Gupta et al., 2015) , (4) the quantization noise in activation gradients be limited as these are back-propagated to calculate the weight gradients, and (5) the precision of weight accumulators should be set so as to avoid premature stoppage of convergence (Goel and Shanbhag, 1998) . The above insights can be formally described via the following five quantization criteria. Criterion 1. Equalizing Feedforward Quantization Noise (EFQN) Criterion. The reflected quantization noise variances onto the mismatch probability p m from all feedforward weights ({V W l →pm } L l=1 ) and activations ({V A l →pm } L l=1 ) should be equal:

V W1→pm = . . . = V W L →pm = V A1→pm = . . . = V A L →pm

Criterion 2. Gradient Clipping (GC) Criterion. The clipping rates of weight ({β

G (W ) l } L l=1 )

and activation ({β:

G (A) l+1 } L l=1

) gradients should be less than a maximum value β 0 :

β G (W ) l < β 0 & β G (A) l+1 < β 0 for l = 1 . . . L.

Criterion 3. Relative Quantization Bias (RQB) Criterion. The relative quantization bias of weight gradients ({η

G (W ) l } L l=1

) should be less than a maximum value η 0 : reflected from quantizing

η G (W ) l < η 0 for l = 1 . . . L.

G (A) l+1 , should be less than V G (W ) l →Σ l : V G (A) l+1 →Σ l ≤ V G (W ) l →Σ l for l = 1 . . . L.

where Σ l is the total sum of element-wise variances of G , its total sum of element-wise variances.

Further explanations and motivations behind the above criteria are presented in Appendix B. The following claim ensures the satisfiability of the above criteria. This leads to closed form expressions for the precision requirements we are seeking and completes our methodology. The validity of the claim is proved in Appendix C. Claim 1. Satisfiability of Quantization Criteria. The five quantization criteria (EFQN, GC, RQB, BQN, AS) are satisfied if:

• The precisions B W l and B A l are set as follows:

EQUATION

for l = 1 . . . L, where rnd() denotes the rounding operation, E W l →pm and E A l →pm are the weight and activation quantization noise gains at layer l, respectively, B (min) is a reference minimum precision, and

E (min) = min {E W l →pm } L l=1 , {E A l →pm } L l=1 .

• The weight and activation gradients PDRs are lower bounded as follows:

EQUATION

where σ

(max)

G (W ) l and σ (max) G (A) l+1

are the largest recorded estimates of the weight and activation gradients

standard deviations σ G (W ) l and σ G (A) l+1

, respectively. • The weight and activation gradients quantization step sizes are upper bounded as follows:

∆ G (W ) l < σ (min) G (W ) l 4 & ∆ G (A) l+1 < ∆ G (W ) l λ (max) G (A) l+1 →G (W ) l   G (W ) l G (A) l+1   1/4 for l = 1 . . . L (3)

where σ

(min) G (W ) l is the smallest recorded estimate of σ G (W ) l and λ (max) G (A) l+1 →G (W ) l

is the largest singular value of the square-Jacobian (Jacobian matrix with squared entries) of

G (W ) l with respect to G (A)

l+1 .

• The accumulator PDR and step size satisfy:

EQUATION

where γ (min) is the smallest value of the learning rate used during training.

Practical considerations: Note that one of the 2L feedforward precisions will equal B (min) . The formulas to compute the quantization noise gains are given in Appendix C and require only one forward-backward pass on an estimation set. We would like the EFQN criterion to hold upon convergence; hence, (1) is computed using the converged model from the FL baseline. For backward signals, setting the values of PDR and LSB is sufficient to determine the precision using the identity , of the gradient tensors, are required. These are obtained via tensor spatial averaging, so that one estimate per tensor is required, and updated in a moving window fashion, as is done for normalization parameters in BatchNorm (Ioffe and Szegedy, 2015) . Furthermore, it might seem that computing the Jacobian in (3) is a difficult task; however, the values of its elements are already computed by the back-prop algorithm, requiring no additional computations (see Appendix C). Thus, the Jacobians (at different layers) are also estimated during training. Due to the typical very large size of modern neural networks, we average the Jacobians spatially, i.e., the activations are aggregated across channels and mini-batches while weights are aggregated across filters. This is again inspired by the work on Batch Normalization (Ioffe and Szegedy, 2015) and makes the probed Jacobians much smaller.

B A = log 2 r A ∆ A + 1,

4. NUMERICAL RESULTS:

We conduct numerical simulations in order to illustrate the validity of the predicted precision configuration C o and investigate its minimality and benefits. We employ three deep learning benchmarking datasets: CIFAR-10, CIFAR-100 (Krizhevsky and Hinton, 2009) , and SVHN (Netzer et al., 2011) . All experiments were done using a Pascal P100 NVIDIA GPU. We train the following networks:

• CIFAR-10 ConvNet: a 9-layer convolutional neural network trained on the CIFAR-10 dataset described as 2 × (64C3)

-M P 2 -2 × (128C3) -M P 2 -2 × (256C3) -2 × (512F C) -10

where C3 denotes 3 × 3 convolutions, M P 2 denotes 2 × 2 max pooling operation, and F C denotes fully connected layers. • SVHN ConvNet: the same network as the CIFAR-10 ConvNet, but trained on the SVHN dataset.

• CIFAR-10 ResNet: a wide deep residual network (Zagoruyko and Komodakis, 2016) with ResNet-20 architecture but having 8 times as many channels per layer compared to (He et al., 2016) . • CIFAR-100 ResNet: same network as CIFAR-10 ResNet save for the last layer to match the number of classes (100) in CIFAR-100.

A step by step description of the application of our method to the above four networks is provided in Appendix E. We hope the inclusion of these steps would: (1) clarify any ambiguity the reader may have from the previous section and (2) facilitate the reproduction of our results. For each network, the 5-tuple Ĉo represents the average number of bits per tensor type. For the ResNets, layer depths 21 and 22 correspond to the strided convolutions in the shortcut connections of residual blocks 4 and 7, respectively. Activation gradients go from layer 2 to L + 1 and are "shifted to the left" in order to be aligned with the other tensors.

4.1. PRECISION CONFIGURATION C o & CONVERGENCE:

The precision configuration C o , with target p m ≤ 1%, β 0 ≤ 5%, and η 0 ≤ 1%, via our proposed method is depicted in Figure 2 for each of the four networks considered. We observe that C o is dependent on the network type. Indeed, the precisions of the two ConvNets follow similar trends as do those the two ResNets. Furthermore, the following observations are made for the ConvNets:

• weight precision B W l decreases as depth increases. This is consistent with the observation that weight perturbations in the earlier layers are the most destructive (Raghu et al., 2017) .

• the precisions of activation gradients (B G (A) l:

) and internal weight accumulators (B W (acc) l ) increases as depth increases which we interpret as follows: (1) the back-propagation of gradients is the dual of the forward-propagation of activations, and (2) accumulators store the most information as their precision is the highest.

• the precisions of the weight gradients (B G (W ) l ) and activations (B A l ) are relatively constant across layers.

Interestingly, for ResNets, the precision is mostly uniform across the layers. Furthermore, the gap between B W (acc) l and the other precisions is not as pronounced as in the case of ConvNets. This suggests that information is spread equally among all signals which we speculate is due to the shortcut connections preventing the shattering of information (Balduzzi et al., 2017) . 

4.2. NEAR MINIMALITY OF C o:

To determine that C o is a close-to-minimal precision assignment, we compare it with: (a)

C +1 = C o + 1 L×5 , and (b) C -1 = C o -1 L×5

, where 1 L×5 is an L × 5 matrix with each entry equal to 13 , i.e., we perturb C o by 1-b in either direction. Figure 3 propagated via the straight through estimator (Bengio et al., 2013) as was done in (Hubara et al., 2016) , • fixed-point training with stochastic quantization (SQ). As was done in (Gupta et al., 2015) , we quantize feedforward weights and activations as well as all gradients, but accumulators are kept in floating-point. The precision configuration (excluding accumulators) is inherited from C o (hence we determine exactly how much stochastic quantization helps),

• training with ternarized gradients (TG) as was done in TernGrad (Wen et al., 2017) . All computations are done in floating-point but weight gradients are ternarized according to the instantaneous tensor spatial standard deviations {-2.5σ, 0, 2.5σ} as was suggested by Wen et al. (2017) . To compute costs, we assume all weight gradients use two bits although they are not really fixed-point and do require computation of 32-b floating-point scalars for every tensor.

The comparison is presented in Table 1 . The first observation is a massive complexity reduction compared to FL. For instance, for the CIFAR-10 ConvNet, the complexity reduction is 2.6× (= 148/56.5), 5.5× (= 9.3/1.7), 7.9× (= 94.4/11.9), and 3.5× (= 49/14) for C W , C A , C M , and C C , respectively. Similar trends are observed for the other four networks. Such complexity reduction comes at the expense of no more than 0.56% increase in test error. For the CIFAR-100 network, the accuracy when training in fixed-point is even better than that of the baseline.

The representational and communication costs of BN is significantly greater than that of FX because the gradients and accumulators are kept in full precision, which masks the benefits of binarizing feedforward tensors. However, benefits are noticeable when considering the computational cost which is lowest as binarization eliminates multiplications. Furthermore, binarization causes a severe accuracy drop for the ConvNets but surprisingly not for the ResNets. We speculate that this is due to the high dimensional geometry of ResNets (Anderson and Berg, 2017) .

As for SQ, since C o was inherited, all costs are identical to FX, save for C W which is larger due to full precision accumulators. Furthermore, SQ has a positive effect only on the CIFAR-10 ConvNet where it clearly acted as a regularizer.

TG does not provide complexity reductions in terms of representational and computational costs which is expected as it only compresses weight gradients. Additionally, the resulting accuracy is slightly worse than that of all other considered schemes, including FX. Naturally, it has the lowest communication cost as weight gradients are quantized to just 2-b.

5. DISCUSSION:



5.1. RELATED WORKS:

Many works have addressed the general problem of reduced precision/complexity deep learning.

Reducing the complexity of inference (forward path): several research efforts have addressed the problem of realizing a DNN's inference path in FX. For instance, the works in (Lin et al., 2016; Sakr et al., 2017) address the problem of precision assignment. While Lin et al. (2016) proposed a non-uniform precision assignment using the signal-to-quantization-noise ratio (SQNR) metric, Sakr et al. (2017) analytically quantified the trade-off between activation and weight precisions while providing minimal precision requirements of the inference path computations that bounds the probability p m of a mismatch between predicted labels of the FX and its FL counterpart. An orthogonal approach which can be applied on top of quantization is pruning (Han et al., 2015) . While significant inference efficiency can be achieved, this approach incurs a substantial training overhead.

A subset of the FX training problem was addressed in binary weighted neural networks (Courbariaux et al., 2015; Rastegari et al., 2016) and fully binarized neural networks (Hubara et al., 2016) , where direct training of neural networks with pre-determined precisions in the inference path was explored with the feedback path computations being done in 32-b FL.

Reducing the complexity of training (backward path): finite-precision training was explored in (Gupta et al., 2015) which employed stochastic quantization in order to counter quantization bias accumulation in the weight updates. This was done by quantizing all tensors to 16-b FX, except for the internal accumulators which were stored in a 32-b floating-point format. An important distinction our work makes is the circumvention of the overhead of implementing stochastic quantization (Hubara et al., 2016) . Similarly, DoReFa-Net (Zhou et al., 2016) stores internal weight representations in 32-b FL, but quantizes the remaining tensors more aggressively. Thus arises the need to re-scale and re-compute in floating-point format, which our work avoids. Finally, Köster et al. (2017) suggests a new number format -Flexpoint -and were able to train neural networks using slightly 16-b per tensor element, with 5 shared exponent bits and a per-tensor dynamic range tracking algorithm. Such tracking causes a hardware overhead bypassed by our work since the arithmetic is purely FX. Augmenting Flexpoint with stochastic quantization effectively results in WAGE (Wu et al., 2018b) , and enables integer quantization of each tensor.

As seen above, none of the prior works address the problem of predicting precision requirements of all training signals. Furthermore, the choice of precision is made in an ad-hoc manner. In contrast, we propose a systematic methodology to determine close-to-minimal precision requirements for FX-only training of deep neural networks.

5.2. CONCLUSION:

In this paper, we have presented a study of precision requirements in a typical back-propagation based training procedure of neural networks. Using a set of quantization criteria, we have presented a precision assignment methodology for which FX training is made statistically similar to the FL baseline, known to converge a priori. We realized FX training of four networks on the CIFAR-10, CIFAR-100, and SVHN datasets and quantified the associated complexity reduction gains in terms costs of training. We also showed that our precision assignment is nearly minimal.

The presented work relies on the statistics of all tensors being quantized during training. This necessitates an initial baseline run in floating-point which can be costly. An open problem is to predict a suitable precision configuration by only observing the data statistics and the network architecture. Future work can leverage the analysis presented in this paper to enhance the effectiveness of other network complexity reduction approaches. For instance, weight pruning can be viewed as a coarse quantization process (quantize to zero) and thus can potentially be done in a targeted manner by leveraging the information provided by noise gains. Furthermore, parameter sharing and clustering can be viewed as a form of vector quantization which presents yet another opportunity to leverage our method for complexity reduction.

B FURTHER EXPLANATIONS AND MOTIVATIONS BEHIND QUANTIZATION CRITERIA:

Criterion 1 (EFQN) is used to ensure that all feedforward quantization noise sources contribute equally to the p m budget. Indeed, if one of the 2L reflected quantization noise variances from the feedforward tensors onto p m , say V Wi→pm for i ∈ {1, . . . , L}, largely dominates all others, it would imply that all tensors but W i are overly quantized. It would therefore be necessary to either increase the precision of W i or decrease the precisions of all other tensors. The application of Criterion 1 (EFQN) through the closed form expression (1) in Claim 1 solves this issue avoiding the need for a trial-and-error approach.

Because FX numbers require a constant PDR, clipping of gradients is needed since their dynamic range is arbitrary. Ideally, a very small PDR would be preferred in order to obtain quantization steps of small magnitude, and hence less quantization noise. We can draw parallels from signal processing theory, where it is known that for a given quantizer, the signal-to-quantization-noise ratio (SQNR) is equal to SQN R(dB) = 6B + 4.78 -P AR where P AR is the peak-to-average ratio, proportional to the PDR. Thus, we would like to reduce the PDR as much as possible in order to increase the SQNR for a given precision. However, this comes at the risk of overflows (due to clipping). Criterion 2 (GC) addresses this trade-off between quantization noise and overflow errors.

Since the back-propagation training procedure is an iterative one, it is important to ensure that any form of bias does not corrupt the weight update accumulation in a positive feedback manner. FX quantization, being a uniform one, is likely to induce such bias when quantized quantities, most notable gradients, are not uniformly distributed. Criterion 3 (RQB) addresses this issue by using η as proxy to this bias accumulation a function of quantization step size and ensuring that its worst case value is small in magnitude.

Criterion 4 (BQN) is in fact an extension of Criterion 1 (EFQN), but for the back-propagation phase. Indeed, once the precision (and hence quantization noise) of weight gradients is set as per Criterion 3 (RQB), it is needed to ensure that the quantization noise source at the activation gradients would not contribute more noise to the updates. This criterion sets the quantization step of the activation gradients.

Criterion 5 (AS) ties together feedforward and gradient precisions through the weight accumulators.

It is required to increment/decrement the feedforward weights whenever the accumulated updates cross-over the weight quantization threshold. This is used to set the PDR of the weight accumulators. Furthermore, since the precision of weight gradients has already been designed to account for quantization noise (through Criteria 2-4), the criterion requires that the accumulators do not cause additional noise.

C PROOF OF CLAIM 1:

The validity of Claim 1 is derived from the following five lemmas. Note that each lemma addresses the satisfiability of one of the five quantization criteria presented in the main text and corresponds to part of Claim 1. Lemma 1. The EFQN criterion holds if the precisions B W l and B A l are set as follows: min) for l = 1 . . . L, where rnd() denotes the rounding operation, B (min) is a reference minimum precision, and E (min) is given by:

B W l = rnd log 2 E W l →pm E (min) + B (min) & B A l = rnd log 2 E A l →pm E (min) + B (

EQUATION

Proof. By definition of the reflected quantization noise variance, the EFQN, by definition, is satisfied if:

∆ 2 W1 12 E W1→pm = . . . = ∆ 2 W L 12 E W L →pm = ∆ 2 A1 12 E A1→pm = . . . = ∆ 2 A L 12 E A L →pm ,

where the quantization noise gains are given by:

EQUATION

for l = 1 . . . L, where {Z i } M i=1 are the soft outputs and Z Ŷfl is the soft output corresponding to Ŷfl . The expressions for these quantization gains are obtained by linearly expanding (across layers) those used in (Sakr et al., 2017) . Note that a second order upper bound is used as a surrogate expression for p m .

From the definition of quantization step size, the above is equivalent to:

2 -2B W 1 E W1→pm = . . . = 2 -2B W L E W L →pm = 2 -2B A 1 E A1→pm = . . . = 2 -2B A L E A L →pm .

Let E (min) be as defined in ( 5):

E (min) = min {E W l →pm } L l=1 , {E A l →pm } L l=1

. We can divide each term by E (min) :

2 -2B W 1 E W1→pm E (min) = . . . = 2 -2B W L E W L →pm E (min) = 2 -2B A 1 E A1→pm E (min) = . . . = 2 -2B A L E A L →pm E (min)

where each term is positive, so that we can take square roots and logarithms such that:

B W1 -log 2 E W1→pm E (min) = . . . = B W L -log 2 E W L →pm E (min) =B A1 -log 2 E A1→pm E (min) = . . . = B A L -log 2 E A L →pm E (min)

Thus we equate all of the above to a reference precision B (min) yielding: min) for l = 1 . . . L. Note that because E (min) is the least quantization noise gain, it is equal to one of the above quantization noise gains so that the corresponding precision actually equates B (min) . As precisions must be integer valued, each of B (min) , {B W l } L l=1 , and {B A l } L l=1 have to be integers, and thus a rounding operation is to be applied on all logarithm terms. Doing so results in (1) from Lemma 1 which completes this proof.

B W l = log 2 E W l →pm E (min) + B (min) & B A l = log 2 E A l →pm E (min) + B (

We close this appendix by discussing the approximation made by invoking the Central Limit Theorem (CLT) in the proofs of Lemmas 2 & 3. This approximation was made because, typically, a backpropagation iteration computes gradients of a loss function being averaged over a mini-batch of samples. By linearity of derivatives, the gradients themselves are averages, which warrants the invocation of the CLT. However, the CLT is an asymptotic result which might be imprecise for a finite number of samples. In typical training of neural networks, the number of samples, or mini-batch size, is in the range of hundreds or thousands (Goyal et al., 2017) . It is therefore important to quantify the preciseness, or lack thereof, of the CLT approximation. On way to do so is via the Berry-Essen Theorem which considers the average of n independent, identically distributed random variables with finite absolute third moment ρ and standard deviation σ. The worst case deviation of the cumulative distribution of the true average from the of the approximated Gaussian random variable (via the CLT), also known as the Kolmogorov-Smirnov distance, KS, is upper bounded as follows: KS < Cρ √ nσ 3 , where C < 0.4785 (Tyurin, 2010) . Observe that the quantity ρ σ 3 is data dependent. To estimate this quantity, we performed a forward-backward pass for all training samples at the start of each epoch for our four networks considered. The statistics ρ and σ were estimated by spatial (over tensors) and sample (over training samples) averages. The maximum value of the ratio ρ σ 3 for all gradient tensors was found to be 2.8097. The mini-batch size we used in all our experiments was 256. Hence, we claim that the CLT approximation in Lemmas 2 & 3 is valid in our context up to a worst case Kolmogorov-Smirnov distance of KS < 0.4785×2.8097 √ 256 = 0.084. The minimality experiments in the main paper only consider a full 1-b perturbation to the full precision configuration matrix. We further investigate the minimality of C o and its sensitivity to precision perturbation per tenor type. The results of this investigation are presented in Fig. 4 .

First, we consider random fractional precision perturbations, meaning perturbations to the precision configuration matrix where only a random fraction p of the 5L precision assignments is incremented or decremented. A fractional precision perturbation of 1 (-1) corresponds to C +1 (C -1 ). A fractional precision perturbation of 0.5 (-0.5) means that a randomly chosen half of the precision assignments is incremented (decremented). Figure 4 (a) shows the relative test error deviation compared to the test error associated with C o for various fractional precision perturbations. The error deviation is taken in a relative fashion to account for the variability of the different networks' accuracies. For instance, an absolute 1% difference in accuracy on a network trained on SVHN is significantly more severe than one on a network trained on CIFAR-100. It is observed that for negative precision perturbations the variation in test error is more important than for the case of positive perturbations. This is further encouraging evidence that C o is nearly minimal, in that a negative perturbation causes significant accuracy degradation while a positive one offers diminishing returns.

It is also interesting to study which of the 5L tensor types is most sensitive to precision reduction. To do so, we perform a similar experiment whereby we selectively decrement the precision of all tensors belonging to the same type (weights, activations, weight gradients, activation gradients, weight accumulators). The results of this experiment are found in Fig. 4 (b) . It is found that the most sensitive tensor types are weights and activation gradients while the least sensitive ones are activations and weight gradients. This is an interesting finding raising further evidence that there exists some form of duality between the forward propagation of activations and back propagation of derivatives as far as numerical precision is concerned.

E ILLUSTRATION OF METHODOLOGY USAGE:

We illustrate a step by step description of the application of our precision assignment methodology to the four networks we report results on.

E.1 CIFAR-10 CONVNET:

Feedforward Precisions: The first step in our methodology consists of setting the feedforward precisions B W l and B A l . As per Claim 1, this requires using (1). To do so, it is first needed to compute the quantization noise gains using (6). Using the converged weights from the baseline run we obtain: 

Layer Index l 1 2 3 4 5 B W l 7+B (min) 7+B (min) 8+B (min) 8+B (min) 7+B (min) B A l 4+B (min) 1+B (min) 1+B (min) 1+B (min) 2+B (min) Layer Index l 6 7 8 9 B W l 6+B (min) 5+B (min) 4+B (min) 3+B (min) B A l 1+B (min) 1+B (min) 1+B (min) 0+B (min)

The value of B (min) is swept and p m i evaluated on the validation set. It is found that the smallest value of B (min) resulting in p m < 1% is equal to 4 bits. Hence the feedforward precisions are set as follows and as illustrated in gradients. As per Claim 1, an important statistic is the spatial variance of the gradient tensors. We estimate these variances via moving window averages, where at each iteration, the running variance estimate σ2 is updated using the instantaneous variance σ2 as follows:

σ2 ← (1 -θ)σ 2 + θσ 2

where θ is the running average factor, chosen to be 0.1. The running variance estimate of each gradient tensor is dumped every epoch. Using the maximum recorded estimate and (2) we compute the PDRs of the gradient tensors (as a reminder, the PDR is forced to be a power of 2): Note that for weights, layer depths 21 and 22 correspond to the strided convolutions in the shortcut connections of residual blocks 4 and 7, respectively. The value of B (min) is again swept, and it is found that the p m < 1% for B (min) = 3. The feedforward precisions are therefore set as follows and as illustrated in Figure 2 : The value of B (min) is again swept, and it is found that the p m < 1% for B (min) = 3. The feedforward precisions are therefore set as follows and as illustrated in Figure 2 : Gradient Precisions: The spatial variance of the gradient tensors is used to determine the PDRs and the quantization step sizes of weight gradients. The singular values of the square-Jacobians are needed to determine the quantization step sizes of activation gradients. They were computed as follows:

Layer Index l 1 2 3 4 5 r G (W ) l 5.00E-01 1.25E-01 1.25E-01 1.25E-01 6.25E-02 r G (A)

Layer

Layer

:

We use the notation |T | to denote the number of elements in tensor T . Unquantized tensors are assumed to have a 32-b FL representation, which is the single-precision in a GPU.

When considering

32-b FL multiplications, we ignore the cost of exponent addition thereby favoring the FL (conventional) implementation. Boundary effects (in convolutions) are neglected.

PDRs are unchanged across configurations, except for r W (acc) l as per (4).

Figures and Tables:

Figure 1: Problem setup: FX training at layer l of a DNN showing the quantized tensors and the associated precision configuration C l = (BW l , BA l , B G (W ) l , B G (A) l+1 , B W (acc) l

Back-propagated Quantization Noise (BQN) Criterion. The reflected quantization noise variance V G (A) l+1 →Σ l , i.e., the total sum of element-wise variances of G (W ) l

Criterion 5. Accumulator Stopping (AS) Criterion. The quantization noise of the internal accumulator should be zero, equivalently:V W (acc)

as explained in Section 2.1. As per Claim 1, estimates of the second order statistics, e.g., σ G (W ) l and σ G (A) l+1

Figure 2: The predicted precision configurations Co for the CIFAR-10 ConvNet (a), SVHN ConvNet (b), CIFAR-10 ResNet (c), and CIFAR-100 ResNet (d).For each network, the 5-tuple Ĉo represents the average number of bits per tensor type. For the ResNets, layer depths 21 and 22 correspond to the strided convolutions in the shortcut connections of residual blocks 4 and 7, respectively. Activation gradients go from layer 2 to L + 1 and are "shifted to the left" in order to be aligned with the other tensors.

FX training curves in Figure3indicate that C o leads to convergence and consistently track FL curves with close fidelity. This validates our analysis and justifies the choice of C o .

Figure 3: Convergence curves for the CIFAR-10 ConvNet (a), SVHN ConvNet (b), CIFAR-10 ResNet (c), and CIFAR-100 ResNet (d) including FL training as well as FX training with precision configurations Co, C1, and C-1.

Figure 4: Additional experiments on minimality and sensitivity of C o : relative test error deviation with respect to C o as a function of (a) random fractional precision perturbations, and (b) 1-b precision reduction per tensor type.

The second step of the methodology is to determine the precisions of weight B G (W ) l and activation B G (A) l+1

-02 3.13E-02 1.56E-02 1.25E-01 r G (A) l+1 1.56E-02 7.81E-03 7.81E-03 3.13E-02 Furthermore, using the minimum recorded estimates of the weight gradient spatial variances and (3) we compute the values of the quantization step sizes of the weight tensors:

also contains the convergence curves for the two new configurations. As shown, C -1 always results in a noticeable gap compared to C o for both the loss function (except for the CIFAR-10 ResNet) and the test error. Furthermore, C +1 offers no observable improvements over C o (except for the test error of CIFAR-10 ConvNet). These results support our contention that C o is close-to-minimal in that increasing the precision above C o leads to diminishing returns while reducing precision below C o leads to a noticeable degradation in accuracy. Additional experimental results provided in Appendix D support our contention regarding the near minimality of C o . Furthermore, by studying the impact of quantizing specific tensors we determine that that the accuracy is most sensitive to the precision assigned to weights and activation gradients.4.3 COMPLEXITY VS. ACCURACYWe would like to quantify the reduction in training cost and expense in terms of accuracy resulting from our proposed method and compare them with those of other methods. Importantly, for a fair comparison, the same network architecture and training procedure are used. We report C W , C A , C M , C C , and test error, for each of the four networks considered for the following training methods: Complexity (CW , CA, CM , and CC ) and accuracy (test error) for the floating-point (FL), fixed-point (FX) with precision configuration Co, binarized network (BN), stochastic quantization (SQ), and TernGrad (TG) training schemes.

Table content: <table><tr><td/><td>CW</td><td>CA</td><td>CM</td><td>CC</td><td>Test</td><td>CW</td><td>CA</td><td>CM</td><td>CC</td><td>Test</td></tr><tr><td/><td>(10 6 b)</td><td>(10 6 b)</td><td>(10 9 FA)</td><td>(10 6 b)</td><td>Error</td><td>(10 6 b)</td><td>(10 6 b)</td><td>(10 9 FA)</td><td>(10 6 b)</td><td>Error</td></tr><tr><td/><td/><td colspan="3">CIFAR-10 ConvNet</td><td/><td/><td colspan="3">SVHN ConvNet</td><td/></tr><tr><td>FL</td><td>148</td><td>9.3</td><td>94.4</td><td>49</td><td>12.02%</td><td>148</td><td>9.3</td><td>94.4</td><td>49</td><td>2.43%</td></tr><tr><td colspan="2">FX (Co) 56.5</td><td>1.7</td><td>11.9</td><td>14</td><td>12.58%</td><td>54.3</td><td>1.9</td><td>10.5</td><td>14</td><td>2.58%</td></tr><tr><td>BN</td><td>100</td><td>4.7</td><td>2.8</td><td>49</td><td>18.50%</td><td>100</td><td>4.7</td><td>2.8</td><td>49</td><td>3.60%</td></tr><tr><td>SQ</td><td>78.8</td><td>1.7</td><td>11.9</td><td>14</td><td>11.32%</td><td>76.3</td><td>1.9</td><td>10.5</td><td>14</td><td>2.73%</td></tr><tr><td>TG</td><td>102</td><td>9.3</td><td>94.4</td><td>3.1</td><td>12.49%</td><td>102</td><td>9.3</td><td>94.4</td><td>3.1</td><td>3.65%</td></tr><tr><td/><td/><td colspan="3">CIFAR-10 ResNet</td><td/><td/><td colspan="3">CIFAR-100 ResNet</td><td/></tr><tr><td>FL</td><td>1784</td><td>96</td><td>4319</td><td>596</td><td>7.42%</td><td>1789</td><td>97</td><td>4319</td><td colspan="2">597 28.06%</td></tr><tr><td colspan="2">FX (Co) 726</td><td>25</td><td>785</td><td>216</td><td>7.51%</td><td>750</td><td>25</td><td>776</td><td colspan="2">216 27.43%</td></tr><tr><td>BN</td><td>1208</td><td>50</td><td>128</td><td>596</td><td>7.24%</td><td>1211</td><td>50</td><td>128</td><td colspan="2">597 29.35%</td></tr><tr><td>SQ</td><td>1062</td><td>25</td><td>785</td><td>216</td><td>7.42%</td><td>1081</td><td>25</td><td>776</td><td colspan="2">216 28.03%</td></tr><tr><td>TG</td><td>1227</td><td>96</td><td>4319</td><td>37.3</td><td>7.94%</td><td>1230</td><td>97</td><td>4319</td><td colspan="2">37.3 30.62%</td></tr></table>

11+B (min) 10+B (min) 10+B (min) 10+B (min) 10+B (min) 11+B (min) E 11+B (min) 11+B (min) 11+B (min) 11+B (min) 10+B (min) 10+B (min) E

Table content: <table><tr><td>Layer Index l</td><td>1</td><td>2</td><td>3</td><td>4</td><td>5</td><td>6</td></tr><tr><td>E W l →pm</td><td>2.41E+03</td><td>9.80E+02</td><td>1.22E+03</td><td>1.62E+03</td><td>1.52E+03</td><td>3.05E+03</td></tr><tr><td>B W l</td><td/><td/><td/><td/><td/><td/></tr><tr><td>A l →pm</td><td>7.32E-01</td><td>5.15E-01</td><td>1.29E-01</td><td>1.12E-01</td><td>7.31E-02</td><td>8.98E-02</td></tr><tr><td>B A l</td><td>5+B (min)</td><td>5+B (min)</td><td>4+B (min)</td><td>4+B (min)</td><td>3+B (min)</td><td>3+B (min)</td></tr><tr><td>Layer Index l</td><td>7</td><td>8</td><td>9</td><td>10</td><td>11</td><td>12</td></tr><tr><td>E W l →pm</td><td>1.47E+03</td><td>2.15E+03</td><td>2.74E+03</td><td>4.96E+03</td><td>4.23E+03</td><td>4.20E+03</td></tr><tr><td>B W l</td><td colspan="6">10+B (min) 11+B (min) 11+B (min) 11+B (min) 11+B (min) 12+B (min)</td></tr><tr><td>E A l →pm</td><td>7.70E-02</td><td>8.39E-02</td><td>6.38E-02</td><td>1.92E-01</td><td>1.54E-01</td><td>1.33E-01</td></tr><tr><td>B A l</td><td>3+B (min)</td><td>3+B (min)</td><td>3+B (min)</td><td>4+B (min)</td><td>4+B (min)</td><td>4+B (min)</td></tr><tr><td>Layer Index l</td><td>13</td><td>14</td><td>15</td><td>16</td><td>17</td><td>18</td></tr><tr><td>E W l →pm</td><td>7.25E+03</td><td>2.99E+03</td><td>2.86E+03</td><td>3.00E+03</td><td>5.02E+03</td><td>4.34E+03</td></tr><tr><td>B W l</td><td/><td/><td/><td/><td/><td/></tr><tr><td>A l →pm</td><td>1.13E-01</td><td>8.51E-02</td><td>6.57E-02</td><td>1.29E-01</td><td>6.51E-02</td><td>2.16E-02</td></tr><tr><td>B A l</td><td>4+B (min)</td><td>3+B (min)</td><td>3+B (min)</td><td>4+B (min)</td><td>3+B (min)</td><td>2+B (min)</td></tr><tr><td>Layer Index l</td><td>19</td><td>20</td><td>21</td><td>22</td><td/><td/></tr><tr><td>E W l →pm</td><td>1.41E+03</td><td>1.30E+03</td><td>1.08E+02</td><td>8.31E+00</td><td/><td/></tr><tr><td>B W l</td><td>9+B (min)</td><td>7+B (min)</td><td colspan="2">11+B (min) 11+B (min)</td><td/><td/></tr><tr><td>E A l →pm</td><td>4.80E-03</td><td>7.82E-04</td><td/><td/><td/><td/></tr><tr><td>B A l</td><td>1+B (min)</td><td>0+B (min)</td><td/><td/><td/><td/></tr></table>

Gradient Precisions: The spatial variance of the gradient tensors is used to determine the PDRs and the quantization step sizes of weight gradients. The singular values of the square-Jacobians are needed to determine the quantization step sizes of activation gradients. They were computed as follows:E.4 CIFAR-100 RESNET Feedforward Precisions: The quantization noise gains are used to obtain values for the precisions as a function of B (min) as summarized below: 10+B (min) 11+B (min) 11+B (min) 11+B (min) 11+B (min) 11+B(min)  E 11+B (min) 11+B (min) 11+B (min) 11+B (min) 11+B (min) 11+B (min) E

Table content: <table><tr><td>Layer Index l</td><td>1</td><td/><td>2</td><td/><td/><td>3</td><td/><td>4</td><td/><td/><td>5</td><td>6</td></tr><tr><td>E W l →pm</td><td>2.32E+03</td><td/><td colspan="2">8.23E+02</td><td colspan="2">1.18E+03</td><td/><td colspan="2">1.28E+03</td><td colspan="2">1.70E+03</td><td>2.78E+03</td></tr><tr><td>B W l</td><td>10+B (min)</td><td/><td colspan="2">9+B (min)</td><td colspan="8">10+B (min) 10+B (min) 10+B (min) 10+B (min)</td></tr><tr><td>E A l →pm</td><td>1.42E+00</td><td/><td colspan="2">7.84E-01</td><td colspan="2">2.52E-01</td><td/><td colspan="2">1.46E-01</td><td colspan="2">7.68E-02</td><td>7.40E-02</td></tr><tr><td>B A l</td><td>5+B (min)</td><td/><td colspan="2">4+B (min)</td><td colspan="2">4+B (min)</td><td/><td colspan="2">3+B (min)</td><td colspan="2">3+B (min)</td><td>3+B (min)</td></tr><tr><td>Layer Index l</td><td>7</td><td/><td>8</td><td/><td/><td>9</td><td/><td>10</td><td/><td/><td>11</td><td>12</td></tr><tr><td>E W l →pm</td><td>3.03E+03</td><td/><td colspan="2">5.80E+03</td><td colspan="2">7.29E+03</td><td/><td colspan="2">9.20E+03</td><td colspan="2">9.81E+03</td><td>1.41E+04</td></tr><tr><td>B W l</td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/></tr><tr><td>A l →pm</td><td>7.52E-02</td><td/><td colspan="2">8.70E-02</td><td colspan="2">1.38E-01</td><td/><td colspan="2">2.49E-01</td><td colspan="2">2.11E-01</td><td>1.51E-01</td></tr><tr><td>B A l</td><td>3+B (min)</td><td/><td colspan="2">3+B (min)</td><td colspan="2">3+B (min)</td><td/><td colspan="2">4+B (min)</td><td colspan="2">3+B (min)</td><td>3+B (min)</td></tr><tr><td>Layer Index l</td><td>13</td><td/><td>14</td><td/><td/><td>15</td><td/><td>16</td><td/><td/><td>17</td><td>18</td></tr><tr><td>E W l →pm</td><td>7.67E+03</td><td/><td colspan="2">1.40E+04</td><td colspan="2">1.13E+04</td><td/><td colspan="2">1.09E+04</td><td colspan="2">5.35E+03</td><td>3.97E+03</td></tr><tr><td>B W l</td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/></tr><tr><td>A l →pm</td><td>1.54E-01</td><td/><td colspan="2">1.09E-01</td><td colspan="2">1.93E-01</td><td/><td colspan="2">2.36E-01</td><td colspan="2">1.27E-01</td><td>3.01E-02</td></tr><tr><td>B A l</td><td>3+B (min)</td><td/><td colspan="2">3+B (min)</td><td colspan="2">3+B (min)</td><td/><td colspan="2">4+B (min)</td><td colspan="2">3+B (min)</td><td>2+B (min)</td></tr><tr><td>Layer Index l</td><td>19</td><td/><td>20</td><td/><td/><td>21</td><td/><td>22</td><td/><td/><td/></tr><tr><td>E W l →pm</td><td>8.35E+02</td><td/><td colspan="2">2.30E+01</td><td colspan="2">6.78E+03</td><td/><td colspan="2">6.03E+03</td><td/><td/></tr><tr><td>B W l</td><td>9+B (min)</td><td/><td colspan="2">7+B (min)</td><td colspan="5">11+B (min) 11+B (min)</td><td/><td/></tr><tr><td>E A l →pm</td><td>2.01E-02</td><td/><td colspan="2">1.80E-03</td><td/><td/><td/><td/><td/><td/><td/></tr><tr><td>B A l</td><td>2+B (min)</td><td/><td colspan="2">0+B (min)</td><td/><td/><td/><td/><td/><td/><td/></tr><tr><td/><td colspan="2">Index l 1</td><td>2</td><td>3</td><td>4</td><td>5</td><td>6</td><td>7</td><td>8</td><td colspan="3">9 10 11</td></tr><tr><td/><td>B W l</td><td colspan="11">14 13 13 13 13 14 13 14 14 14 14</td></tr><tr><td/><td>B A l</td><td>8</td><td>8</td><td>7</td><td>7</td><td>6</td><td>6</td><td>6</td><td>6</td><td>6</td><td>7</td><td>7</td></tr><tr><td colspan="13">Layer Index l 12 13 14 15 16 17 18 19 20 21 22</td></tr><tr><td/><td>B W l</td><td colspan="11">15 14 14 14 14 13 13 12 10 14 14</td></tr><tr><td/><td>B A l</td><td>7</td><td>7</td><td>6</td><td>6</td><td>7</td><td>6</td><td>5</td><td>4</td><td>3</td><td/></tr></table>

Peer review (iclr2019):
 This paper proposes a FX (fixed point) framework to calculate the reduced bit numbers, which can (I) use float numbers with the reduced bits to represent each NN layer*s weight values W_l, activation values A_l, gradients of weights G_l^W, gradients of nodes G_{l+1}^A, and the cumulated (updated) weight W_l^(acc); (II) with the reduced representations, the training loss and testing loss will not be sacrificed much comparing with the original FL framework, where each float number is 32-bit. Some positive points: (a) The proposed FX framework can reduce cost for both inference and training. (b) The experimental results looks promising. (c) The Criteria 1-5 seems systematic and the conditions in Claim 1 can be used to calculate the required bit numbers. The author proposed an implementation of Claim 1. Some negative points / questions: (a) The most important part of the paper is Criteria 1-5. Criteria 1 generalizes the idea from Sakr et al. 2017, to force the contributions of weight and activation almost at the same order, which seems reasonable to make the mismatch budget p_m smallest. But the other criteria (with their corresponding notions, e.g., Criterion 2 and the concept of clipping rate .beta) are introduced in a way which is not clear enough and make the audience confused. For example, why is clipping rate .beta and relative quantization bias .eta is needed here and what is their relationship with the usual weight gradient clip norm (5% target .beta and 1% .eta target correspond to what order of clip norm)? Criteria 4 & 5 are introduced in the same way with one sentence explained like heuristics. They seem to me are introduced just for reducing corresponding bit numbers. More motivation and explanation of introducing these criteria and notions are needed. (b) For W_l and A_l, the necessary bit numbers are calculated using Criteria 1 & EFQN condition. But for gradients, their PDRs and Deltas are calculated using other criteria & conditions. Then how to calculate their bit numbers from PDRs and Deltas? (c) The proof of Lemma 2 & 3 directly used CLT for mini-batch average gradient items. CLT is for asymptotic case and in finite sample case it is not true. So it is heuristic calculation rather than lemma with proof (If seeking proof then some finite-sample argument like Berry-Esseen theorem is needed to quantify the probability of the average is not Gaussian). And what is the mini-batch size used here? If it is too small then probably the error of taking the mini-batch SGD as Gaussian will be large. (d) The importance/minimality of each individual bit number in C_o is not investigated, and the claim of the near minimality of C_o cannot hold according to the current experiments. The experiments of C_{+1} and C_{-1} do not exclude the possibility that some items (not the whole C_o) are minimal. More experiments (changing one or more items while fixing the others) are needed to show every item in C_o is minimal (or not). And which one of them is the most important for training/testing (how sensitive the training/testing performance is to each bit number)? Also in C_{-1} and C_{+1}, are target .beta, .eta changed? What are these changed values? (e) From the Claim 1, it seems that these bit numbers are sufficient, and smaller than necessary to get the same training/testing performance (e.g., from the proof of Lemma 3, with the result .eta = 0.4% < 1%). So one question is what kind of .beta and .eta target are necessary for preserving the performance and what corresponding bit numbers are necessary to achieve the necessary .beta and .eta values? (f) The computational cost definition looks different with Eq (3) in Sakr et al. 2017. Why? Some typos: 1.Figure 2(a), B_{A_l} is 8 for Layer 1, but 9 in Appendix Table. 2. In the last third paragraph of Page 7, is it 2.6 = (148/56.5) instead of 2.6 * (148/56.5)? Same for the following numbers. ======================================== Revision: I have read the reply from the authors and it clarified several matters. I adjusted my rating of this paper (from 6 to 7).