BEER: Fast O (1 /T ) Rate for Decentralized Nonconvex Optimization with Communication Compression
Abstract
Communication efficiency has been widely recognized as the bottleneck for largescale decentralized machine learning applications in multi-agent or federated environments.
To tackle the communication bottleneck, there have been many efforts to design communication-compressed algorithms for decentralized nonconvex optimization, where the clients are only allowed to communicate a small amount of quantized information (aka bits) with their neighbors over a predefined graph topology.
Despite significant efforts, the state-of-the-art algorithm in the nonconvex setting still suffers from a slower rate of convergence O (( G/T ) 2 / 3 ) compared with their uncompressed counterpart, where G measures the data heterogeneity across different clients, and T is the number of communication rounds.
This paper proposes BEER , which adopts communication compression with gradient tracking, and shows it converges at a faster rate of O (1 /T ) .
This significantly improves over the state-of-the-art rate, by matching the rate without compression even under arbitrary data heterogeneity.
Numerical experiments are also provided to corroborate our theory and confirm the practical superiority of BEER in the data heterogeneous regime.
Introduction
Decentralized machine learning is gaining attention in both academia and industry because of its emerging applications in multi-agent systems such as the internet-of-things (IoT) and networked autonomous systems [31, 43].
One of the key problems in decentralized machine learning is ondevice training, which aims to optimize a machine learning model using the datasets stored on (geographically) different clients and can be formulated as a decentralized optimization problem.
Decentralized optimization aims to solve the following optimization problem without sharing the local datasets with other clients:
[Equation]
where f ( x ; D i ) :
=
E ξ i ∼D i f ( x ; ξ i ) for i ∈ [ n ] , and n is the total number of clients.
Here, x ∈ R d is the machine learning model, f ( x ; D ) , f ( x ; D i ) , and f ( x ; ξ i ) denote the loss functions of the model x on the entire dataset D , the local dataset D i , and a random data sample ξ i , respectively.
Different from the widely studied distributed or federated learning setting where there is a central server to coordinate the parameter sharing across all clients, in the decentralized setting, each client can only communication with its neighbors over a communication network determined by a predefined network topology.
The main bottleneck of decentralized optimization—when it comes to large-scale machine learning applications—is communication efficiency, due to the large number of clients involved in the network [43] and the enormous size of machine learning models [4], exacerbated by resource constraints such as limited bandwidth availability and stringent delay requirements.
One way to reduce the communication cost is communication compression, which only transmits compressed messages (with fewer bits) between the clients using compression operators .
The compression operators come with many design choices and offer great flexibility in different trade-offs of communication and computation in practice.
Even though communication compression has been extensively applied to distributed or federated optimization with a central server [47, 12, 6, 26, 11, 25, 40, 41, 27], its use in the decentralized setting has been relatively sparse.
Most of the existing approaches only apply to the strongly convex setting [39, 17, 30, 19, 29, 23], and only a few consider the general nonconvex setting [16, 51, 45].
Our contributions This paper considers decentralized optimization with communication compression, focusing on the nonconvex setting due to its critical importance in modern machine learning, such as training deep neural networks [21], word embeddings, and other unsupervised learning models [42].
Unfortunately, existing algorithms [16, 51, 45] suffer from several important drawbacks in the nonconvex setting: they need strong bounded gradient or bounded dissimilarity assumptions to guarantee convergence, and the convergence rate is order-wise slower than their uncompressed counterpart in terms of the communication rounds (see Table 1).
In this paper, we introduce BEER , which is a decentralized optimization algorithm with communication compression using gradient tracking.
BEER not only removes the strong assumptions required in all prior works, but enjoys a faster convergence rate in the nonconvex setting.
Concretely, we have the following main contributions (see Tables 1 and 2).
We show that BEER converges at a fast rate of O (1 /T ) in the nonconvex setting, which improves over the state-of-the-art rate O (1 /T 2 / 3 ) of CHOCO-SGD [16] and Deepsqueeze [51], where T is the number of communication rounds.
This matches the rate without compression even under arbitrary data heterogeneity across the clients.
We also provide the analysis of BEER under the Polyak- Łojasiewicz (PL) condition (Assumption 4), and show that BEER converges at a linear rate (see Table 2).
Note that strong convexity implies the PL condition, and thus BEER also achieves linear convergence in the strongly convex setting.
We run numerical experiments on real-world datasets and show BEER achieves superior or competitive performance when the data are heterogeneous compared with state-of-the-art baselines with and without communication compression.
To the best of our knowledge, BEER is the first algorithm that achieves O (1 /T ) rate without the bounded gradient or bounded dissimilarity assumptions, supported by a strong empirical performance in the data heterogeneous setting.
Notation Throughout this paper, we use boldface letters to denote vectors, e.g., x ∈ R d .
Let [ n ] denote the set { 1 , 2 , · · · , n } , 1 be the all-one vector, I be the identity matrix, ∥·∥ denote the Euclidean norm of a vector, and ∥·∥ F denote the Frobenius norm of a matrix.
Let ⟨ u , v ⟩ denote the standard Euclidean inner product of two vectors u and v .
In addition, we use the standard order notation O ( · ) to hide absolute constants.
Due to space constraints, additional discussions of related work, further experiments, and the proof details can be found in the appendix.
Problem Setup
In this section, we formally define the decentralized optimization problem with communication compression, and introduce a few important quantities and assumptions that will be used in developing our algorithm and theory.
[Table]
[Table]
Decentralized optimization The goal of decentralized optimization is to solve
[Equation]
, where n is the number of clients, f ( x ) is the global objective function, and f i ( x ) :
= f ( x ; D i ) :=
E ξ i ∼D i f ( x ; ξ i ) is the local objective function, with x the parameter of interest, ξ i a random data sample drawn from the local dataset D i .
In the decentralized setting, the clients can only communicate with their local neighbors over a prescribed network topology, which is specified by an undirected weighted graph G ([ n ] , E ) .
Here, each node in [ n ] represents a client, and E is the set of possible communication links between different clients.
Information sharing across the clients is implemented mathematically by the use of a mixing matrix W =
[ w ij ] ∈ [0 , 1] n × n , which is defined in accordance with the network topology: we assign a positive weight w ij for any ( i, j ) ∈ E and w ij = 0 for all ( i, j ) / ∈ E .
We make the following standard assumption on the mixing matrix [36].
Assumption 1 (Mixing matrix)
The mixing matrix W = [ w ij ] ∈ [0 , 1] n × n is symmetric ( W ⊤ =
W ) and doubly stochastic ( W 1 = 1 , 1 ⊤ W = 1 ⊤ ).
Let its eigenvalues be 1 = | λ 1 ( W ) | > | λ 2 ( W ) | ≥ · · · ≥ | λ n ( W ) | .
The spectral gap is denoted by
[Equation]
The spectral gap of a mixing matrix is closely related to the network topology, see Nedi´c et al.
[36] for its scaling with respect to the network size (i.e. the number of clients n ) for representative network topologies.
Compression operators Compression, in the forms of quantization or sparsification, can be used to reduce the total communication cost.
We now introduce the notion of a randomized compression operator , which is widely used in the decentralized/federated optimization literature, e.g. Tang et al.
[49], Stich et al.
[47], Koloskova et al.
Definition 1 (Compression operator)
A randomized map C : R d (cid:55)→ R d is an α -compression operator if for all x ∈ R d , it satisfies
[Equation]
In particular, no compression ( C ( x ) ≡ x ) implies α = 1 .
Compared with the unbiased compression operator used in, e.g., Alistarh et al.
[1], Khirirat et al.
[15], Mishchenko et al.
[33], Li and Richtárik [24], the compression operator in Definition 1 does not impose the additional constraint on the expectation such that E [ C ( x )]
= x .
Besides, it is always possible to convert an unbiased compression operator into a biased one satisfying Definition 1, and thus Definition 1 is a generalization of the unbiased compression operator that allows biased compression.
Practical examples of the compression operators are provided in Appendix B.
Assumptions on functions We now state the assumptions on the functions { f i } and f .
Throughout this paper, we assume that f ∗ =
min x f ( x ) exists and f ∗ >
[Equation]
In the nonconvex setting, we assume that the functions { f i } i ∈ [ n ] are arbitrary functions that satisfy the following standard smoothness assumption.
Assumption 2 (Smoothness)
The function f is L -smooth if there exists L ≥ 0 such that
[Equation]
In addition, we allow local computation to be performed via stochastic gradient updates, where ˜ ∇ f i ( x ) :=
∇ f i ( x ; ξ i ) denotes a local stochastic gradient computed via a sample ξ i drawn i.i.d.
from
D i , and ˜ ∇ b f i ( x ) := 1 b (cid:80) bj =1 ∇ f i ( x ; ξ i,j ) denotes the stochastic gradient computed by a minibatch with size b drawn i.i.d.
from D i .
We assume ˜ ∇ f i ( x ) and ˜ ∇ b f i ( x ) have bounded variance, which is again standard in the decentralized/federated optimization literature [32, 13, 16].
Assumption 3 (Bounded variance)
There exists a constant σ ≥
such that for all i ∈ [ n ] and x ∈ R d ,
[Equation]
For a stochastic gradient with minibatch size b , we have
[Equation]
In addition, we consider the setting when the function f additionally satisfies the following PolyakŁojasiewicz (PL) condition [37], which can lead to fast linear convergence even when the function is nonconvex.
Assumption 4 (PL condition)
There exists some constant µ > 0 such that for any x ∈ R d ,
[Equation]
Note that the PL condition is a weaker assumption than strong convexity, which means that if the objective function f is µ -strongly convex, then the PL condition also holds with the parameter µ .
Proposed Algorithm
In this section, we introduce our proposed algorithm BEER for decentralized nonconvex optimization with compressed communication.
Before embarking on the description of BEER , we introduce some convenient matrix notation.
Since in a decentralized setting, the parameter estimates at different clients are typically different, we use X =
[ x 1 , x 2 , . . .
, x n ] to denote the collection of parameter estimates from all clients, where x i is from client i .
The average of { x i } i ∈ [ n ] is denoted by ¯ x := 1 n X 1 .
Other quantities are defined similarly.
With slight abuse of notation, we define
[Equation]
which collects the local gradients computed at the local parameters.
Similarly, the stochastic variant is defined as ˜ ∇ b F ( X ) :=
[ ˜ ∇ b f 1 ( x 1 ) , ˜ ∇ b f 2 ( x 2 ) , . .
. , ˜ ∇ b f n ( x n )] .
We also allow the compression operator to take vector values, which are applied in a column-wise fashion, i.e.,
C ( X ) :=
[ C ( x 1 ) , . . .
, C ( x n )]
∈ R d × n .
We now proceed to describe BEER , which is detailed in Algorithm 1 using the matrix notation introduced above.
At the t -th iteration, BEER maintains the current model estimates X t and the global gradient estimates V t across the clients.
At the crux of its design, BEER also tracks and maintains two control sequences H t and G t that serve as compressed surrogates of X t and V t , respectively.
In particular, these two control sequences are updated by aggregating the received compressed messages alone (cf.
Line 5 and Line 8).
It then boils down to how to carefully update these quantities in each iteration with communication compression, which we now explain in details.
To begin, note that for each client i , BEER not only maintains its own parameters { x ti , v ti , h ti , g ti } , but also the control variables from its neighbors, namely, { h tj } j ∈N ( i ) and { g tj } j ∈N ( i ) .
Update the model estimate Each client i first updates its model x t +1 i according to Line 3.
By thinking of { h tj } j ∈N ( i ) as a surrogate of { x tj } j ∈N ( i ) , the second term aims to achieve better consensus among the clients through mixing, while the last term performs a gradient descent update.
Update the global gradient estimate Each client i updates the global gradient estimate v t
+1 i according to Line 6, where the last correction term—based on the difference of the gradients at consecutive models—is known as a trick called gradient tracking [38, 7, 35].
The use of gradient tracking is critical: as shall be seen momentarily, it contributes to the key difference from CHOCOSGD that enables the fast rate of O (1 /T ) without any bounded dissimilarity or bounded gradient assumptions.
Indeed, if we remove the control sequence G t and substitute Lines 6-8 by V t +1 =
˜ ∇ b F ( X t +1 ) , we recover CHOCO-SGD from BEER .
Update the compressed surrogates with communication To update { h tj } j ∈N ( i ) , each client i first computes a compressed message q t +1 h,i that encodes the difference x t +1 i − h ti , and broadcasts to its neighbors (cf.
Line 4).
Then, each client i updates { h tj } j ∈N ( i ) by aggregating the received compressed messages { q t +1 h,j } j ∈N ( i ) following Line 5.
The updates of { g tj } j ∈N ( i ) can be performed similarly.
Moreover, all the compressed messages can be sent in a single communication round at one iteration, i.e., the communications in Lines 4 and 7 can be performed at once.
This leverages
EF21 [40] for communication compression, which is a better and simpler algorithm that deals with biased compression operators compared with the error feedback (or error compensation, EF/EC ) framework [12, 46].
Using the control sequence G t , BEER does not need to apply EF/EC explicitly and can deal with the error implicitly.
Convergence Guarantees
In this section, we show the convergence guarantees of BEER under different settings: the O (1 /T ) rate in the nonconvex setting in Section 4.1, and the improved linear rate under the PL condition (Assumption 4) in Section 4.2.
In Section 4.3, we briefly sketch the proof.
Our convergence guarantees are based on an appropriately designed Lyapunov function, given by
[Equation]
where the choice of constants { c i } 4 i =1 might be different from theorem to theorem, E f (¯ x t ) − f ∗ represents the sub-optimality gap, and the errors { Ω ti } 4 i =1 are defined by
[Equation]
Here, Ω t 1 and Ω t 2 denote the compression errors for X t and V t when approximated using the compressed surrogates H t and G t respectively, and Ω t 3 and Ω t 4 denote the consensus errors of X t and V t .
4.1 Convergence in the nonconvex setting
First, we present the following convergence result of BEER in the nonconvex setting when there is no local variance ( σ 2 = 0 ), i.e., we can use the local full gradient ∇ F ( X t ) instead of ˜ ∇ b F ( X t ) in
Line 6 of Algorithm 1.
Theorem 1 (Convergence in the nonconvex setting without local variance) Suppose
Assumptions 1, and 2 hold, and we can compute the local full gradient ∇ f i ( x ) for any x .
Then there exist absolute constants c 1 , c 2 , c 3 , c 4 , c γ , c η > 0 , such that if we set γ =
c γ αρ , η = c η γρ 2 /L , then for the
Lyapunov function Φ t in (4), it holds
[Equation]
Theorem 1 shows that BEER converges at a rate of O (1 /T )
when there is no local variance ( σ 2 = 0 ), which is faster than the O (1 /T 2 / 3 ) rate by CHOCO-SGD [16] and DeepSqueeze [51], and
the
O (1 / √ T ) rate by SQuARM-SGD [45]; see also Table 1.
More specifically, to achieve 1 T (cid:80) T − 1 t =
E ∥∇ f
(¯ x t ) ∥ 2 ≤ ϵ 2 , BEER needs O (cid:16) 1 ρ 3 αϵ 2 (cid:17) iterations or communication rounds, where ρ and α are the spectral gap (cf.
and the compression parameter (cf.
(3)), respectively.
In comparison, the state-of-the-art algorithm CHOCO-SGD [16] converges at a rate of O (( G/ρ 2 αT ) 2 / 3 ) , which translates to an iteration complexity of O (cid:16) G ρ 2 αϵ 3 (cid:17) , with G being the bounded gradient parameter, namely, E ξ i ∼D i ∥∇ f
( x , ξ i ) ∥ 2 ≤ G 2 .
Therefore, BEER improves over CHOCO-SGD not only in terms of a better dependency on ϵ , but also removing the bounded gradient assumption, which is significant since in practice, G can be excessively large due to data heterogeneity across the clients.
The dependency on α of BEER is consistent with other compression schemes, such as CHOCO-SGD ,
DeepSqueeze and SQuARM-SGD for the nonconvex setting, as well as LEAD [30] and EF-C-GT [29] for the strongly convex setting.
As for the dependency on ρ , BEER is slightly worse than CHOCO-SGD , where CHOCO-SGD has a dependency of O (1 /ρ 2 ) whereas BEER has a dependency of O (1 /ρ 3 ) .
This degeneration is also seen in the analysis of uncompressed decentralized algorithms using gradient tracking [48, 54], where the rate O (1 /ρ 2 ) is worse than the rate of O (1 /ρ ) for basic decentralized SGD algorithms [14, 28] by a factor of ρ .
In addition, both BEER and CHOCO-SGD use small mixing step size γ to guarantee convergence, which makes the dependency on ρ worse than their uncompressed counterparts.
Stochastic gradient oracles BEER also supports the use of stochastic gradient oracles with bounded local variance (Assumption 3).
More specifically, we have the following theorem, which generalizes
Theorem 1.
Theorem 2 (Convergence in the nonconvex setting) Suppose Assumptions 1, 2 and 3 hold.
Then there exist absolute constants c 1 , c 2 , c 3 , c 4 , c γ , c η > 0 , such that if we set γ = c γ αρ , η = c
η γρ 2 /L ,
[Table]
then for the Lyapunov function Φ t in (4), it holds
[Equation]
In the presence of local variance, the squared gradient norm of BEER has an additional term that scales on the order of O (cid:16) σ 2 αb (cid:17) (ignoring other parameters).
By choosing a sufficiently large minibatch size b , i.e. b ≥ O (cid:16) σ 2 αϵ 2 (cid:17) , BEER maintains the iteration complexity O (cid:16) 1 ρ 3 αϵ 2 (cid:17) to reach
[Equation]
T
[Equation]
E ∥∇ f
(¯ x t ) ∥ 2 ≤ ϵ 2 , without the bounded gradient assumption, thus inheriting similar advantages over CHOCO-SGD as discussed earlier.
While our focus is on communication efficiency, to gain more insights, Table 3 summarizes both the communication rounds and the gradient complexity for different decentralized stochastic methods.
While BEER does not require the bounded gradient assumption, it may lead to a worse gradient complexity in the data homogeneous setting due to the use of large minibatch size.
Fortunately, this only impacts the local computation cost, and does not exacerbate the communication complexity, which is often the bottleneck.
It is of great interest to further refine the design and analysis of BEER in terms of the gradient complexity.
4.2 Linear convergence with PL condition
Now, we show that the convergence of BEER can be strengthened to a linear rate with the addition of the PL condition (Assumption 4).
Similar to the nonconvex setting, we first show the convergence result without local gradient variance ( σ 2 = 0 ).
Theorem 3 (Convergence under the PL condition without local variance) Suppose
Assumptions 1, 2, and 4 hold, and we can compute the local full gradient ∇ f i ( x ) for any x .
Then there exist constants c 1 , c 2 , c 3 , c 4 , c γ , c η > 0 , such that if we set γ =
c γ αρ , η = c η γρ 2 /L , then for the
Lyapunov function Φ t in (4), it holds
[Equation]
Theorem 3 demonstrates that under the PL condition, BEER converges linearly to the global optimum f ∗ , where it finds an ϵ -optimal solution in O (cid:16) κρ 3 α log (cid:0)
ϵ
cid:1)(cid:17) iterations, with κ :=
L/µ the condition number.
Note that in the strongly convex case, Liao et al.
[29] proposed an algorithm that also uses error feedback compression and gradient tracking simultaneously, which achieves a linear rate of convergence with unclear dependencies with salient problem parameters.
In comparison, BEER achieves an explicit linear rate of convergence in the strongly convex case as well, given strong convexity implies the PL condition.
In fact, our analysis for the nonconvex setting—as will be made evident in our proof—almost implies immediately the linear convergence under the PL condition, thus provides a major step forward compared with prior analyses that only considered the strongly convex setting.
Stochastic gradient oracles Under the PL condition, BEER also supports the use of stochastic gradient oracles with bounded local variance (Assumption 3).
The following theorem shows that
BEER linearly converges to a neighborhood of size O (cid:16) σ 2 αb (cid:17) around the global optimum.
Theorem 4 (Convergence under PL condition) Suppose Assumptions 1, 2, 3, and 4 hold.
Then there exist absolute constants c 1 , c 2 , c 3 , c 4 , c γ , c η > 0 , such that if we set γ = c γ αρ , η = c
η γρ 2 /L , then for the Lyapunov function Φ t in (4), it holds
[Equation]
4.3 Proof sketch
We now provide a proof sketch of Theorem 1, which establishes the O (1 /T ) rate of BEER in the nonconvex setting using full gradient, highlighting the technical reason of the rate improvement of
BEER over CHOCO-SGD .
The full proofs of our theorems are delegated to Appendix D.
Recalling the quantities Ω t 1 to Ω t 4 from (5), which capture the approximation errors using compression and the consensus errors of X t and V t , we would like to control these errors by obtaining inequalities of the form:
[Equation]
where 0 < a i < 1 denotes the size of the contraction, and b i > 0 wraps together other terms which may be dependent on Ω tj for j ̸
= i as well as the expected squared gradient norm of ¯ v t , i.e.,
[Equation]
Then, by choosing the Lyapunov function properly (cf.
(4)), we can show that the Lyapunov function actually descends, and small manipulations lead to the claimed convergence results in Theorem 1.
We now explain briefly how gradient tracking helps in BEER .
Note that CHOCO-SGD also has the control variable H t for the model X t , therefore in its analysis, it also deals with the quantities Ω tand Ω t 3 .
However, CHOCO-SGD also needs to bound the term ∥ V t ∥ 2F , where V t = ∇ F ( X t ) for
CHOCO-SGD when using full gradients.
Thus, CHOCO-SGD needs to assume the bounded gradient assumption and only obtain a slower O (1 /T 2 / 3 ) convergence rate.
In contrast, BEER deals with the term ∥ V t ∥ 2F by decomposing it using Young’s inequality, leading to
[Equation]
for some β > 0 .
Here, Ω t 4 can be controlled via the gradient tracking technique (see Line 6 in
Algorithm 1) without the bounded gradient assumption, and Ω t 5 can be handled using the smoothness assumption (Assumption 2).
Numerical Experiments
This section presents numerical experiments on real-world datasets to showcase BEER ’s superior ability to handle data heterogeneity across the clients, by running each experiment on unshuffled datasets and comparing the performances with the state-of-the-art baseline algorithms both with and without communication compression.
The code can be accessed at:
[Equation]
Experiment setup We run experiments on two nonconvex problems to compare with the baseline algorithms both with and without communication compression: logistic regression with a nonconvex regularizer [52] on the a9a dataset [5], and training a 1-hidden layer neural network on the MNIST dataset [20].
For logistic regression with a nonconvex regularizer, following Wang et al.
[52], the objective function over a datum ( a , b ) is defined as
[Equation]
Communication rounds
T r a i n i n g g r a d
i e n t n o r m 500 1 , 000 1 , 500 2 , 000 0 .0 .0 .0 .
Communication rounds
T e s t i n g a cc u r a c y 0 0 .
Communication bits
T r a i n i n g g r a d
i e n t n o r m 0 0 .
Communication bits
T e s t i n g a cc u r a c y
BEERCHOCO-SGDDSGDD where the last term is the nonconvex regularizer and the regularization parameter is set to α = 0 .
For 1-hidden layer neural network training, we use 32 hidden neurons, sigmoid activation functions and cross-entropy loss.
The objective function over a datum ( a , b ) is defined as f ( x ; ( a , b ))
=
ℓ ( softmax ( W 2 sigmoid ( W 1 a + c 1 ) + c 2 ) , b ) , where ℓ ( · , · ) denotes the cross-entropy loss, the optimization variable is collectively denoted by x = vec ( W 1 , c 1 , W 2 , c 2 ) , where the dimensions of the network parameters W 1 , c 1 , W 2 , c 2 are 64 × 784 , 64 × 1 , 10 × 64 , and 10 × 1 , respectively.
For both experiments, we split the unshuffled datasets evenly to 10 clients that are connected by a ring topology.
By using unshuffled data, we can simulate the scenario with high data heterogeneity across clients.
Approximately, for the a9a dataset, 5 clients receive data with label 1 and others receive data with label 0 ; for the MNIST dataset, client i receives data with label i .
We use the FDLA matrix [53] as the mixing matrix to perform weighted information aggregation to accelerate convergence.
Results We compare BEER with 1) CHOCO-SGD [17], which is the state-of-the-art nonconvex decentralized optimizing algorithm using communication compression, and 2) DSGD [28] and
D 2 [50], which are decentralized optimization algorithms without compression.
All algorithms are initialized in the same experiment by the same initial point.
Moreover, we use the same best-tuned learning rate η = 0 .
1 , batch size b = 100 , and biased compression operator (gsgd b ) [1] for BEER and CHOCO-SGD on both experiments.
Figure 1 and Figure 2 plot the training gradient norm and testing accuracy against communication rounds and communication bits for logistic regression with nonconvex regularization and 1-hiddenlayer neural network training, respectively.
In the nonconvex logistic regression experiment (cf.
Figure 1), the algorithms with communication compression ( BEER and CHOCO-SGD ) converge faster than the uncompressed algorithms ( DSGD and D 2 ) in terms of the communication bits.
However, CHOCO-SGD fails to converge to a small gradient norm solution since it cannot tolerate a high level of data dissimilarity across different clients, and its performance is not comparable to D 2 .
In contrast, BEER can converge to a point with a relatively smaller gradient norm, which is comparable to D 2 .
The performance of testing accuracy is similar to that of the training gradient norm, where BEER achieves the best testing accuracy and also learns faster than the uncompressed algorithms.
Turning to the neural network experiment (cf.
Figure 2), in terms of the final training gradient norm, BEER converges to a solution comparable to D 2 , but at a lower communication cost, while
CHOCO-SGD and DSGD cannot converge due to the data heterogeneity.
In terms of testing accuracy,
BEER and D 2 have very similar performance, and outperform CHOCO-SGD and DSGD .
Convolutional neural network training We further compare the performance of BEER and
CHOCO-SGD on training a convolutional neural network using the unshuffled MNIST dataset.
Specifically, the network is consist of three modules: the first module is a 2-d convolution layer (1 input channel, 16 output channels, kernel size 5, stride 1 and padding 2) followed by 2-d batch normalization, ReLU activation and 2-d max pooling (kernel size 2 and stride 2); the second module is the same as the first module, except the convolution layer has 16 input channels and 32 output channels; the last module is a fully-connected layer with 1568 inputs and 10 outputs.
We adopt the standard
[Figure]
i
[Figure]
s t
i
[Figure]
u
[Figure] [Figure]
T r a i n i n g g r a d
i e n t n o
[Figure]
m
[Figure] [Figure]
[Figure] [Figure] [Figure]
[Figure] [Figure] [Figure]
[Figure] [Figure]
e s t
i n g
[Figure]
cc
[Figure]
a c y
[Figure]
cross-entropy loss, and simply average each agent’s model with its neighbors.
Figure 3 shows the testing gradient norm and accuracy against the communication bits.
It can be seen that BEER outperforms CHOCO-SGD in terms of both testing gradient norm and testing accuracy.
Both algorithms converge fast initially, however, due to to the extreme data heterogeneity, their convergence speeds significantly degenerate after a short time.
BEER keeps improving the objective when CHOCO-SGD hits its error floor
[Figure]
which highlights BEER ’s advantage to deal with data heterogeneity
[Figure]
In summary, BEER has much better performance in terms of communication efficiency than CHOCOSGD in heterogeneous data scenario, which corroborates our theory.
BEER also performs similarly or even better than the uncompressed baseline algorithm D 2 , and much better than DSGD .
In addition, by leveraging different communication compression schemes, BEER allows more flexible trade-offs between communication and computation, making it an appealing choice in practice.
Conclusion
This paper presents BEER , which achieves a faster O (1 /T ) convergence rate for decentralized nonconvex optimization with communication compression, without imposing the bounded dissimilarity or bounded gradient assumptions.
In addition, a faster linear rate of convergence is established for
BEER under the PL condition.
Numerical experiments are provided to corroborate our theory on the advantage of BEER in the data heterogeneous scenario.
An interesting direction of future work is to investigate the lower bounds for decentralized (nonconvex) optimization with communication compression.
In addition, improving the dependency of BEER with the network topology parameter ρ , possibly leveraging the analysis in Koloskova et al.
[18], is of interest.
Acknowledgement
The work of H. Zhao is supported in part by NSF, ONR, Simons Foundation, DARPA and SRC through awards to S. Arora.
The work of B. Li, Z. Li and Y. Chi is supported in part by ONR N0001419-1-2404, by AFRL under FA8750-20-2-0504, and by NSF under CCF-1901199, CCFand CNS-2148212.
B. Li is also gratefully supported by Wei Shen and Xuehong Zhang Presidential
Fellowship at Carnegie Mellon University.
The work of P. Richtárik is supported by KAUST Baseline
[-Lemma 
-]
Research Fund.