- Research
- Open Access
- Published:

# A multi-block alternating direction method with parallel splitting for decentralized consensus optimization

*EURASIP Journal on Wireless Communications and Networking*
**volume 2012**, Article number: 338 (2012)

## Abstract

Decentralized optimization has attracted much research interest for resource-limited networked multi-agent systems in recent years. Decentralized _{T}*consensus* optimization, which is one of the decentralized optimization problems of great practical importance, minimizes an objective function that is the sum of the terms from individual agents over a set of variables on which all the agents should reach a consensus. This problem can be reformulated into an equivalent model with two blocks of variables, which can then be solved by the alternating direction method (ADM) with only communications between neighbor nodes. Motivated by a recently emerged class of so-called *multi-block* ADMs, this article demonstrates that it is more natural to reformulate a decentralized consensus optimization problem to one with multiple blocks of variables and solve it by a multi-block ADM. In particular, we focus on the multi-block ADM with parallel splitting, which has easy decentralized implementation. Convergence rate is analyzed in the setting of average consensus, and the relation between two-block and multi-block ADMs are studied. Numerical experiments demonstrate the effectiveness of the multi-block ADM with parallel splitting in terms of speed and communication cost and show that it has better network scalability.

## Introduction

In recent years, the communication, signal processing, control, and optimization communities have witnessed considerable research efforts on decentralized optimization for networked multi-agent systems [1–3]. A networked multi-agent system, such as a wireless sensor network (WSN) or a networked control system (NCS), is composed of multiple geographically distributed but interconnected agents which have sensing, computation, communication, and actuating abilities. This system generally has limited resources for communication, since battery power is limited and recharging is difficult, while communication between two agents is energy-consuming. Furthermore, the communication link is often vulnerable and bandwidth-limited. In this situation, decentralized optimization emerges as an effective approach to improve network scalability. In decentralized optimization, data and computation are decentralized. Each agent exchanges information with its neighbors and accomplishes an otherwise centralized optimization task.

This article focuses on the *decentralized consensus optimization* problem. We consider a network of *L* agents which cooperatively optimize a separable objective function [3–8]:

where ${f}_{i}\left(x\right):{\mathcal{R}}^{N}\to \mathcal{R}$ is a convex function known to agent *i* only. The goal is to minimize the objective subject to consensus on *x*.

### Related study

The decentralized consensus optimization formulation (1) arises in many practical applications, such as averaging [9–11], estimation [12–17], learning [18–21], etc. The form of *f*_{
i
}(*x*) can be least squares [11–13], *ℓ*_{1}-regularized least squares [14–17], or more general ones [18–21]. Note that this model can be extended to account for those with separable constraints, such as the network utility maximization (NUM) problem [22–24].

Existing approaches to solving (1) include: i) belief propagation based on graphical models and Markovian random fields [18–20]; ii) incremental optimization which minimizes the overall objective function along a predefined path on the network [7, 8]; iii) stochastic optimization with information exchange between neighboring agents [4–6]; and iv) optimization with explicit consensus constraints which can be handled with the alternating direction method (ADM) [3, 12–17]. The ADM approach is fully decentralized, does not make any assumptions on network infrastructure such as free of loop or with a predefined path, and generally has satisfactory convergence performance. In this article, we mainly discuss the application of ADMs in the decentralized consensus optimization problem.

Our research is along the line of information-driven signal processing and control of WSNs and NCSs [24–26]. Accompanied with the unprecedented data collection abilities offered by large-scale networked multi-agent systems, a new challenge also arises: *how should we process such a large amount of data to make estimates and produce control strategies given limited network resources?* Instead of processing the data in a fusion center, our solution is letting each agent autonomously make decisions aided by limited communication with its neighbors. From this perspective, each individual objective function *f*_{
i
}(*x*) in (1) is constructed from the data collected by agent *i*, and *x* is the global information common to all agents (e.g., estimates or control strategies) obtained based on the data collected by the whole network. Though this framework can be generalized to various signal processing and control problems, this article focuses on those can be formulated as (1). For problems such as dynamic control and Kalman filtering of networked multi-agent systems, interested readers are referred to [1, 2, 27, 28], respectively.

### Our contribution

Motivated by a series of recent articles on multi-block ADMs and their convergence analysis [29–31], this article describes their applications to the decentralized consensus optimization problem. The multi-block ADM with parallel spliting is reviewed in Section 3. Unlike the classical ADM (see textbooks [32, 33]), this multi-block ADM splits the optimization variables into *multiple* blocks and sequentially updates just one of them while fixing the others. The classical ADM, on the other hand, only has two blocks of variables. Hence in this article we refer to it by the two-block ADM. Our problem (1) does not naturally have two distinct blocks of variables, and to apply the two-block ADM one needs to introduce extra variables (see e.g., [15, 16, 32]). We review this in Section 2. On the other hand, it is simpler to apply the multi-block ADM to (1) and the resulting algorithm is readily decentralized.

In this article also analyzes the convergence rate of the multi-block ADM applied to the average consensus problem, which is a special case of (1) where ${f}_{i}\left(x\right)=\frac{1}{2}\parallel x-{b}_{i}{\parallel}_{2}^{2}$ for all *i*. In this setting, if the parameters of the multi-block ADM satisfy a certain formula, it is equivalent to the two-block ADM. Therefore, the two-block ADM can be considered as a special case of the multi-block ADM on average consensus problems. This relation also gives a guideline to select the parameters of the multi-block ADM so that it is not equivalent to and runs faster than the two-block ADM on all the tested decentralized consensus optimization problems, including the tested average consensus problems. The simulation results demonstrate that the multi-block ADM accelerates convergence, reduces communication cost, and thus improves network scalability.

### Paper organization

The rest of this article is organized as follows. Section 2 reviews a reformulation of the decentralized consensus optimization problem (1), to which the two-block ADM is applied. Section 3 reviews the multi-block ADM and applies a parallel-splitting version of it to (1). Section 4 elaborates on the convergence rate analysis on the average consensus problem, and shows that the two-block ADM is a special case of the multi-block ADM in this case. Section 5 presents numerical simulations of the two-block and multi-block ADMs. Finally, Section 6 concludes the article. Appendix Appendix 1 is placed in the last section.

## Problem formulation and the two-block ADM

In this section, we describe an equivalent formulation of the decentralized consensus optimization problem (1) and outline the algorithm design based on the two-block ADM.

### Problem formulation

We consider a networked multi-agent system described by an undirected connected communication graph $\mathcal{G}=(\mathcal{L},\mathcal{E})$, where $\mathcal{L}$ is the set of *L* vertexes (distributed agents) and $\mathcal{E}$ is the set of edges (communication links). There exists an edge $(i,j)\in \mathcal{E}$ between agents *i* and *j* if they can directly communicate with each other. The two agents are also called one-hop neighbors, or simply neighbors. The set of one-hop neighbors of agent *i* is denoted by ${\mathcal{N}}_{i}$, whose cardinality is denoted by $\left|{\mathcal{N}}_{i}\right|$.

Our objective is to solve (1) with only information exchange between neighbors. To this end, define *x*^{(i)} as agent *i*’s local copy of *x* and impose consensus constraints *x*^{(i)}=*x*^{(j)}for all pairs of neighbors *i* and *j*. With these and given that the communication graph $\mathcal{G}$ is connected, we obtain the following equivalent formulation of (1) (see e.g., [13]):

### The two-block ADM

Let us consider the following convex program with separable equality constraints:

Here for *i*=1 and 2, ${g}_{i}:{\mathcal{R}}^{{N}_{i}^{\prime}}\to \mathcal{R}$ is convex, ${D}_{i}\in {\mathcal{R}}^{{M}^{\prime}\times {N}_{i}^{\prime}}$, $e\in {\mathcal{R}}^{{M}^{\prime}}$. The two-block ADM constructs the augmented Lagrangian function as:

Here $\lambda \in {\mathcal{R}}^{{M}^{\prime}}$ is a Lagrange multiplier and *c* is a positive constant. At the *t* th iteration, the two-block ADM updates the optimization variables *θ*_{1}(*t* + 1) and *θ*_{2}(*t* + 1) as:

and updates the Lagrange multiplier *λ*(*t* + 1) as:

The two-block ADM guarantees global convergence for any *c* > 0 [32]. More precisely, when each *g*_{
i
} is convex for *i* = 1 and 2, the dual sequence {*λ*(*t*)} converges to an optimal dual solution of (5); if further the primal sequence {*θ*_{1}(*t*)^{T}*θ*_{2}(*t*)^{T}^{T}} is bounded, the sequence converges to an optimal primal solution of (5).

### The two-block ADM for decentralized consensus optimization

The two-block ADM cannot be directly applied to problem (2) because its constraints interconnect all the variables pair by pair. There are no obvious two blocks. To overcome this, [32] describes a new block of auxiliary variables, and reformulates (2) as:

Here *z*_{
ij
} is an auxiliary variable attached to *x*^{(i)}and *x*^{(j)}.

Treating {*x*^{(i)}} and {*z*_{
ij
}} as two blocks of variables, the two-block ADM is applied to problem (4). This technique has been adopted in [15, 16] to solve the decentralized consensus optimization problem with neighboring consensus constraints. After eliminating {*z*_{
ij
}} from the iterative updates and further simplifications, the two-block ADM for (4) is given below as algorithm TB-ADM.

**Initialization:** Each agent *i* initializes *x*^{(i)}(0)=0 and *α*_{
i
}(0)=0.

**Step 1:** At time *t*, each agent *i* updates its local copy *x*^{(i)}as: ${x}^{\left(i\right)}(t+1)=arg\underset{{x}^{\left(i\right)}}{min}\phantom{\rule{1em}{0ex}}{f}_{i}\left({x}^{\left(i\right)}\right)+\underset{i}{\overset{T}{\alpha}}\left(t\right){x}^{\left(i\right)}+c\sum _{j\in {\mathcal{N}}_{i}}\left|\right|{x}^{\left(i\right)}-\frac{1}{2}\left({x}^{\left(i\right)}\left(t\right)+{x}^{\left(j\right)}\left(t\right)\right)|\underset{2}{\overset{2}{|}}$, where *α*_{
i
}is the Lagrange multiplier and *c* is a positive constant.

**Step 2:** At time *t*, each agent *i* updates its Lagrange multiplier *α*_{
i
}as: ${\alpha}_{i}(t+1)={\alpha}_{i}\left(t\right)+c\left|{\mathcal{N}}_{i}\right|{x}^{\left(i\right)}(t+1)-c\sum _{j\in {\mathcal{N}}_{i}}{x}^{\left(j\right)}(t+1)$.

**Step 3:** Repeat Step 1 and Step 2 until convergence.

TB-ADM is well suited for decentralized computation since the updates require only communication between agents *i* and *j*, who are one-hop neighbors. Detailed derivation of TB-ADM can be found in [15, 16, 32].

## The multi-block ADM

The fact that many practical optimization problems naturally have multiple blocks of variables motivates the development of a class of multi-block ADMs, such as the one with parallel splitting [29], with prediction-correction [30], and with Gaussian back substitution [31]. Due to the nature of the decentralized consensus optimization problem (2) and the need of parallelization, we choose the multi-block ADM with parallel splitting in [29].

### The multi-block ADM with parallel splitting

Consider an equality constrained convex program which can be separated to *L* parts:

Here for all *i*, ${g}_{i}:{\mathcal{R}}^{{N}_{i}^{\prime}}\to \mathcal{R}$ is convex, ${D}_{i}\in {\mathcal{R}}^{{M}^{\prime}\times {N}_{i}^{\prime}}$, $e\in {\mathcal{R}}^{{M}^{\prime}}$. At the *t* th iteration, the multi-block ADM with parallel splitting works as follows: *Step 1: Updating an auxiliary variable* *q*:

where *q* is an auxiliary variable, *λ* is a Lagrange multiplier, and *β* is a positive constant. *Step 2: Updating optimization variables*{*θ*_{
i
}}:

where *μ* is a positive constant. *Step 3: Updating the Lagrange multiplier* *λ*:

The multi-block ADM guarantees global convergence if the two positive constants *β* and *μ* are properly chosen. For the convergence proof and the settings of *β* and *μ*, the interested reader is referred to [29].

### The multi-block ADM for decentralized consensus optimization

Applying the multi-block ADM in (2) directly gets a decentralized algorithm, and does not need to introduce a new block of auxiliary variables and eliminate them, as we have done in the two-block ADM. We provide the algorithm to solve (2) based on the multi-block ADM with parallel splitting, denoted as MB-ADM. Detailed derivation of MB-ADM is given in Appendix Appendix 1.

**Initialization:** Each agent *i* initializes *q*_{
i
}(0)=0, *x*^{(i)}(0)=0, and *λ*_{
i
}(0)=0.

**Step 1:** At time *t*, each agent *i* updates its auxiliary variable *q*_{
i
}as: ${q}_{i}(t+1)={\lambda}_{i}\left(t\right)+\beta \left|{\mathcal{N}}_{i}\right|{x}^{\left(i\right)}\left(t\right)-\beta \sum _{j\in {\mathcal{N}}_{i}}{x}^{\left(j\right)}\left(t\right)$, where *β* is a positive constant.

**Step 2:** At time *t*, each agent *i* updates its local copy *x*^{(i)}as: ${x}^{\left(i\right)}(t+1)=arg\underset{{x}^{\left(i\right)}}{min}\phantom{\rule{2.77695pt}{0ex}}\phantom{\rule{2.77695pt}{0ex}}{f}_{i}\left({x}^{\left(i\right)}\right)+2\underset{i}{\overset{T}{q}}(t+1){x}^{\left(i\right)}+\mu \left|{\mathcal{N}}_{i}\right|\left|\right|{x}^{\left(i\right)}-{x}^{\left(i\right)}\left(t\right)|\underset{2}{\overset{2}{|}}$, where *μ* is a positive constant.

**Step 3:** At time *t*, each agent *i* updates its Lagrange multiplier *λ*_{
i
}as: ${\lambda}_{i}(t+1)={\lambda}_{i}\left(t\right)+\beta \left|{\mathcal{N}}_{i}\right|{x}^{\left(i\right)}(t+1)-\beta \sum _{j\in {\mathcal{N}}_{i}}{x}^{\left(j\right)}(t+1)$.

**Step 4:** Repeat Step 1 to Step 3 until convergence.

In each iteration, to update *q*_{
i
}(*t* + 1) and *λ*_{
i
}(*t*), agent *i* needs *x*^{(j)}(*t*) with the size of *N*×1 from all neighbors $j\in {\mathcal{N}}_{i}$; to optimize *x*^{(i)}(*t* + 1), agent *i* only needs local information *q*_{
i
}(*t* + 1) and *x*^{(i)}(*t*). In all, each agent only needs to broadcast an *N*×1 vector of its local copy (i.e., *x*^{(i)}(*t*)) to its neighbors per iteration. MB-ADM and TB-ADM have the same per-iteration communication cost. At the *t* th iteration, agent *i* needs to update *x*^{(i)}(*t*), *q*_{
i
}(*t*), and *λ*_{
i
}(*t*) in its memory for MB-ADM. Hence the memory requirement is slightly higher than that of TB-ADM, for which only *x*^{(i)}(*t*) and *α*_{
i
}(*t*) need to be updated.

## Convergence rate analysis

Convergence rate is an significant issue for decentralized algorithms, since it directly influences the overall communication cost. With respect to general separable convex programs, [29, 34] proves the sublinear convergence rates of $\sim \frac{1}{t}$ for the multi-block and two-block ADMs, respectively. However, when they are applied to the average consensus problems, much faster convergence can be observed. For this reason, we improve the convergence rate in this section.

The average consensus problem gives rise to problem (2) with ${f}_{i}\left({x}^{\left(i\right)}\right)=\frac{1}{2}\left|\right|{x}^{\left(i\right)}-{b}_{i}|{|}_{2}^{2},\forall i$[9–11]; namely, agents aims at averaging their original measurements {*b*_{
i
}} via one-hop communication. Without loss of generality, we assume that *x*^{(i)} and *b*_{
i
} are both scalars since their dimensions have no effect on the convergence rate.

### Convergence rate of MB-ADM

In analyzing the convergence rate of MB-ADM for the average consensus problem, we first rewrite MB-ADM as a state transition equation form and then use the spectral analysis tools to provide a bound of convergence rate. Our train of thought is similar to that in [35] for the two-block ADM.

According to the derivation in Appendix Appendix 1, we can rewrite MB-ADM in a state transition equation form. Let us define a state vector *s*_{
M
}(*t* + 1)=[*x*^{(1)}(*t* + 1),…,*x*^{(L)}(*t* + 1),*x*^{(1)}(*t*),…,*x*^{(L)}(*t*)]^{T}and the corresponding state transition equation of MB-ADM is:

Here the state transition matrix *Φ*_{
M
}is defined as:

with *Γ*_{
M
} being an *L*×*L* matrix whose (*i*,*i*)th entry is $\frac{1-4\beta \left|{\mathcal{N}}_{i}\right|+4\mu \left|{\mathcal{N}}_{i}\right|}{1+2\mu \left|{\mathcal{N}}_{i}\right|}$ and (*i*,*j*)th entry is $\frac{4\beta}{1+2\mu \left|{\mathcal{N}}_{i}\right|}$ if *i* and *j* are neighbors, and *Ω*_{
T
}being an *L*×*L* matrix whose (*i*,*i*)th entry is $\frac{2\beta \left|{\mathcal{N}}_{i}\right|-2\mu \left|{\mathcal{N}}_{i}\right|}{1+2\mu \left|{\mathcal{N}}_{i}\right|}$ and (*i*,*j*)th entry is $-\frac{2\beta}{1+2\mu \left|{\mathcal{N}}_{i}\right|}$ if *i* and *j* are neighbors. We can see that summation of each row of *Φ*_{
M
}is 1. The initial state is ${s}_{M}\left(1\right)=[\frac{{b}_{1}}{1+2\mu \left|{\mathcal{N}}_{1}\right|},\dots ,\frac{{b}_{L}}{1+2\mu \left|{\mathcal{N}}_{L}\right|},0,\dots ,0]$ when each agent *i* initializes *q*_{
i
}(0)=0, *x*^{(i)}(0)=0, and *λ*_{
i
}(0)=0.

#### Proposition 1

(**convergence and convergence rate of MB-ADM on average consensus**) The state transition equation (6) defined above has the following properties:

#### Property 1

The matrix *Φ*_{
M
} has an eigenvalue *ρ*_{M 1}=1 with multiplicity 1, and its corresponding left and right eigenvectors are:

and

respectively. Note that *l*_{M 1} and *r*_{M 1} are chosen subject to *l*_{M 1}*r*_{M 1}=1.

#### Property 2

Define:

where *ρ*_{
Mi
} is the *i* th eigenvalue of *Φ*_{
M
}. If *ρ*_{
M
}<1, then the limit property of *s*_{
M
}(*t*) is:

Further, denoting that *κ*_{
M
} is the size of the largest Jordan block of *Φ*_{
M
}, the convergence rate is:

Proof of Property 1 is given in Appendix Appendix 1. Property 2 comes from the classical convergence rate analysis of state transition equations. If *ρ*_{M 1}=1 and *ρ*_{
M
}<1, then there exists a unique ${s}_{M}\left(\infty \right)$ and the convergence rate is $\left|\right|{s}_{M}(t+1)-{s}_{M}\left(\infty \right)|{|}_{2}\sim {t}^{({\kappa}_{M}-1)}{\rho}_{M}^{t}$ (see [36], Fact 3). Next we try to find one possible (and hence unique) ${s}_{M}\left(\infty \right)$. By definition, *Φ*_{
M
}*r*_{M 1}=*ρ*_{M 1}*r*_{M 1}=*r*_{M 1}. Hence $\underset{t\to \infty}{lim}{\Phi}_{M}^{t}{r}_{M1}={r}_{M1}$. Similarly, $\underset{t\to \infty}{lim}{l}_{M1}{\Phi}_{M}={l}_{M1}$. These two facts mean that *r*_{M 1}*l*_{M 1} is a possible limit point of $\underset{t\to \infty}{lim}{\Phi}_{M}^{t}$. Therefore, *r*_{M 1}*l*_{M 1}*s*_{
M
}(1) is a possible (and hence unique) limit point of ${s}_{M}\left(\infty \right)$.

#### Remark 1

Note that the $\sim {t}^{({\kappa}_{M}-1)}{\rho}_{M}^{t}$ rate, though still loose, is tighter than the $\sim \frac{1}{t}$ rate of the multi-block ADM for general separable convex programs [29]. Indeed, from numerical experiments, we find that *κ*_{
M
}, the size of the largest Jordan block of *Φ*_{
M
}, is often equal to 1 (it means that *Φ*_{
M
}is diagonalizable). In this case, the convergence rate can be as fast as $\sim {\rho}_{M}^{t}$.

In Property 2, there is a condition that *ρ*_{
M
}<1. It is not necessarily for true any choices of *μ* and *β*. Next we show two nontrivial special cases where the condition in Property 2 satisfy. The first special case connects MB-ADM with TB-ADM. Analysis of these two special cases as well as numerical simulations provide guidelines for parameter selection in MB-ADM.

#### Proposition 2

(**two nontrivial special cases**) We have *ρ*_{
M
}<1 in either one of the following two cases: Case 1: The parameters *μ* and *β* are chosen such that *μ* = 2*β* > 0; further, $\frac{2\beta \left|{\mathcal{N}}_{j}\right|}{1+2\mu \left|{\mathcal{N}}_{j}\right|}<\frac{1}{4}$ and $\frac{2\mu \left|{\mathcal{N}}_{j}\right|}{1+2\mu \left|{\mathcal{N}}_{j}\right|}<\frac{1}{2}$ for all *j*=1,2,…,*L*. Case 2: The parameters *μ* and *β* are chosen such that *μ*=*β*>0; further, $\frac{2\beta \left|{\mathcal{N}}_{j}\right|}{1+2\mu \left|{\mathcal{N}}_{j}\right|}<\frac{1}{2}$ and $\frac{2\mu \left|{\mathcal{N}}_{j}\right|}{1+2\mu \left|{\mathcal{N}}_{j}\right|}<\frac{1}{2}$ for all *j*=1,2,…,*L*.

#### Remark 2

The proof of Proposition 2 is given in Appendix Appendix 1. In case 1, we set *μ*=2*β*>0, which indeed leads to the equivalence between MB-ADM and TB-ADM, as we will show in the next subsection. In case 2, we set *μ*=*β*>0, which brings faster convergence for the average consensus problem according to numerical simulations (see Section 5.2). Hence we recommend to set *β*=*τμ* with a fixed ratio $\frac{1}{2}\le \tau \le 1$, and just tune the value of *μ*. This setting also works well for the general decentralized consensus optimization problem (1). Tuning *μ* for MB-ADM is similar to tuning *c* for TB-ADM; both algorithms have 1 parameter subject to the user choice. Note that the conditions in Proposition 2 are merely sufficient; $\frac{2\beta \left|{\mathcal{N}}_{j}\right|}{1+2\mu \left|{\mathcal{N}}_{j}\right|}$ and $\frac{2\mu \left|{\mathcal{N}}_{j}\right|}{1+2\mu \left|{\mathcal{N}}_{j}\right|}$ can be larger than their upper bounds given above.

### Connection between MB-ADM and TB-ADM

To show the connection between MB-ADM and TB-ADM, we also write TB-ADM as a state transition equation form. Note that [35] considers another kind of two-block ADM for the average consensus problem, where consensus constraints are quadratically penalized by different weights in the augmented Lagrangian function. In TB-ADM, the consensus constraints are quadratically penalized by the same weight *c*.

We define a state vector *s*_{
T
}(*t* + 1)=[*x*^{(1)}(*t* + 1),…,*x*^{(L)}(*t* + 1),*x*^{(1)}(*t*),…,*x*^{(L)}(*t*)]^{T}and the corresponding state transition equation, according to the derivation in Appendix 1:

Here the state transition matrix *Φ*_{
T
}is defined as:

with *I*_{L×L} being the *L*×*L* identity matrix, 0_{L×L} being the *L*×*L* zero matrix, *Γ*_{
T
} being an *L*×*L* matrix whose (*i*,*i*)th entry is 1 and (*i*,*j*)th entry is $\frac{2c}{1+2c\left|{\mathcal{N}}_{i}\right|}$ if *i* and *j* are neighbors, and *Ω*_{
T
}being an *L*×*L* matrix whose (*i*,*i*)th entry is $-\frac{c\left|{\mathcal{N}}_{i}\right|}{1+2c\left|{\mathcal{N}}_{i}\right|}$ and (*i*,*j*)th entry is $-\frac{c}{1+2c\left|{\mathcal{N}}_{i}\right|}$ if *i* and *j* are neighbors. The initial state is ${s}_{T}\left(1\right)=[\frac{{b}_{1}}{1+2c\left|{\mathcal{N}}_{1}\right|},\dots ,\frac{{b}_{L}}{1+2c\left|{\mathcal{N}}_{L}\right|},0,\dots ,0]$ when each agent *i* initializes *x*^{(i)}(0)=0 and *α*_{
i
}(0)=0.

Comparing the state transition equations of MB-ADM and TB-ADM, we can find that TB-ADM is indeed a special case of MB-ADM when *c*=*μ*=2*β*>0. In this sense, MB-ADM provides more flexibility in parameter selection than TB-ADM. According to our simulations in Section 5.2, setting *β*=*τμ* with $\frac{1}{2}\le \tau \le 1$ makes MB-ADM faster than TB-ADM.

Let *ρ*_{
Ti
} be the *i* th eigenvalue of *Φ*_{
T
}. Apparently *ρ*_{T 1}=1. Defining:

and denoting *κ*_{
T
} as the size of the largest Jordan block of *Φ*_{
T
}, we can prove that TB-ADM has a similar $\sim {t}^{({\kappa}_{T}-1)}{\rho}_{T}^{t}$ convergence rate to the optimal solution given the conditions in Case 1 of Proposition 2. Interestingly, the upper bounds of $\frac{2c\left|{\mathcal{N}}_{j}\right|}{1+2c\left|{\mathcal{N}}_{j}\right|}<\frac{1}{2}$ for all *j*=1,2,…,*L* are no longer needed since TB-ADM guarantees global convergence for any *c*>0.

## Numerical Experiments

In this section, we present numerical simulations and demonstrate the performance of MB-ADM on the decentralized consensus optimization problems. Particularly, we are interested in how the communication cost scales to the network size.

### Simulation Settings

In the numerical experiments, we consider the case that the agents cooperatively solve a least-squares problem. Each agent *i* has a measurement matrix ${A}_{i}\in {\mathcal{R}}^{M\times N}$ and a measurement vector ${b}_{i}\in {\mathcal{R}}^{M}$. The objective function in (1) is thus $f\left(x\right)=\sum _{i=1}^{L}{f}_{i}\left(x\right)=\frac{1}{2}\sum _{i=1}^{L}\left|\right|{A}_{i}x-{b}_{i}|{|}_{2}^{2}$. The elements of the true signal vector *x*_{0} and the entries of the measurement matrices {*A*_{
i
}} follow the normal distribution $\mathcal{N}(0,1)$. The measurement vector *b*_{
i
}=*A*_{
i
}*x*_{0} + *η*_{
i
}; the elements of the noise vector *η*_{
i
}follow the normal distribution $\mathcal{N}(0,0.1)$. In the tests of average consensus, {*A*_{
i
}} reduce to identity matrices and are no longer random.

In the simulation, we assume that *L* agents are uniformly randomly deployed in a 100×100 area. All agents have a common communication range *r*_{
C
}, which is chosen such that the networked multi-agent system is connected. Given *r*_{
C
}, the average node degree *d* can be calculated. We consider the following three scenarios: *#* 1) *L*=50, *M*=1, *N*=1, {*A*_{
i
}=1}, *r*_{
C
}=30, *d*≃12; *#* 2) *L*=50, *M*=10, *N*=5, *r*_{
C
}=30, *d*≃12; *#* 3) *L*=200, *M*=10, *N*=5, *r*_{
C
}=15, *d*≃12. Scenario *#* 1 is the average consensus test. Throughout the simulations, we set *β*=*τμ* in MB-ADM with *τ*=0*.* 9.

### Convergence rate for average consensus

Under different choices of *c*, *μ*, and *β*, the values of *ρ*_{
T
} for TB-ADM and the values of *ρ*_{
M
} for MB-ADM with respect to scenario *#* 1 are shown in Figure 1. For TB-ADM, *ρ*_{
T
}sharply reduces when *c* increases from 0; after a certain turning point (at *c*^{∗}≃0*.* 17) which corresponds to the fastest convergence rate, *ρ*_{
T
}steadily increases. The curve of *ρ*_{
M
}for MB-ADM shows to be more complicated due to the existence of two parameters, *μ* and *β*. For each *μ*, *ρ*_{
M
} steadily reduces when *β* increases from 0, then sharply goes to be larger than 1 which corresponds to divergence. The larger *μ*, the wider convergence range for *β*; but the side-effect is the relatively slower convergence rate. The curve of particular interest to us is *μ*=*c*^{∗}≃0*.* 17. In this curve, 2*β*=*c*^{∗}≃0*.* 17 corresponds to Case 1 in Proposition 2; namely, when MB-ADM reduces to TB-ADM. Increasing *β* from *c*^{∗}, *ρ*_{
M
} still decreases until reaching a turning point 2*β*=2*c*^{∗}≃0*.* 34, which corresponds to Case 2 in Proposition 2. This simulation validates our analysis in Section 5.2, as well as the proposed parameter selection rule (namely, setting a ratio *τ*, $\frac{1}{2}\le \tau \le 1$, such that *β*=*τμ*).

Simulation results about the actual convergence properties are shown in Figure 2. By absolute error we denote the ℓ_{2}-norm of the distance between the current solution and the centralized optimal solution. Though the convergence rates of MB-ADM and TB-ADM are at the same magnitude, MB-ADM shows to be slightly superior to TB-ADM.

According to the theoretical analysis in Sections 4.1 and 4.2, the estimated convergence rates of MB-ADM and TB-ADM are $\sim {t}^{({\kappa}_{M}-1)}{\rho}_{M}^{t}$ and $\sim {t}^{({\kappa}_{T}-1)}{\rho}_{T}^{t}$, respectively. However, numerical simulations show that they are loose bounds; the actual convergence rate, as we can observe from Figure 2, are linear.

### Performance Comparison

Figures 3 and 4 depict the convergence properties of the two decentralized consensus optimization algorithms for scenarios #2 and #3, respectively. The parameters μ and c are tune to be near the best ones with. Here we still have β=τμ with τ=0.9. For either the medium network in scenario #2 or the large network in scenario #3, both algorithms linearly converge to the optimal solution. Comparing the two decentralized algorithms, MB-ADM outperforms TB-ADM in each scenario regarding convergence rate.

What of particular interest to us is whether the decentralized algorithms are scalable to network size. Observing Figure 3 with *L*=50 agents and *d*≃12, and Figure 4 with *L*=200 agents and *d*≃12, we can find that the convergence rates of the two algorithm are more dependent on the average node degree other than on the network size. These numerical experiments verify the well-recognized claim that decentralized optimization may improve the performance of a networked multi-agent system with respect to network scalability.

### Communication cost

Communication cost, in terms of energy consumption and bandwidth, is the major design consideration of a resource-limited networked multi-agent system, and can be approximately evaluated by the volume of information exchange during the decentralized consensus optimization process. Ignoring the extra burden of coordinating the network, for each agent, the communication cost is proportional to the number of iterations multiplied by the volume of information exchange per iteration. Therefore, reducing the information exchange per iteration is of critical importance to the design of lightweight algorithms.

Comparing the two decentralized consensus optimization algorithms, the information exchange per iteration is decided by the communication mode of agents, namely, broadcast or unicast. In the broadcast mode, one agent can send one piece of information to all of its neighbors with one transmission; contrarily, in the unicast mode, the agent needs multiple transmissions to do so. The two modes both have their pros and cons. The broadcast mode utilizes the characteristic of wireless communication, but may brings difficulties in coordinating the network, such as avoiding collisions. Though the unicast mode consumes much more transmissions, the randomized-gossip-like scheme is very useful in communication for the sake of robustness [37]. The average volume of information exchange per iteration of the four decentralized consensus optimization algorithms are outlined in Table 1, for both the broadcast and unicast modes.

In summary, the decentralized consensus optimization algorithms, no matter with the broadcast or unicast mode, are scalable to the network size. Since the number of iterations is proportional to the average node degree *d*, the overall average volume of information exchange is ∼*Nd* for the broadcast mode and ∼*N* *d*^{2} for the unicast mode. As a comparison, consider a centralized networked multi-agent system uniformly randomly deployed in a two-dimensional area with a fusion center which collects measurement vectors from all agents. The average volume of information exchange is $\sim M\sqrt{L}$ while the worst one is ∼*ML* for agents near the fusion center. When the network size *L* increases, the communication cost caused by the centralized network infrastructure is unaffordable and the decentralized network infrastructure is hence superior.

## Conclusion

This article considers solving the decentralized consensus optimization problem with the parallel version multi-block ADM in a networked multi-agent system. The traditional ADM can be used but it requires the introduction of a second block of auxiliary variables whereas our method takes advantages of the problem’s nature of having multiple blocks of variables. We analyze the rate of convergence of our method applied to the average consensus problem. Analysis results that the two-block ADM is a special case of the multi-block ADM on average consensus. With extensive numerical experiments, we demonstrate the effectiveness of the proposed algorithm.

In the implementation of a networked multi-agent system, practical issues such as packet loss, asynchronization, and quantization are inevitable. This article assumes that the communication links are reliable, the network time is slotted and well synchronized, and the exchanged information is not quantized. We would like to address these issues in future research.

## Appendix 1

This section provides some theoretical results in the article.

### Development of MB-ADM

The decentralized consensus optimization problem (2) with neighboring consensus constraints can be rewritten as the form of (5). Apparently, *g*_{
i
}=*f*_{
i
}, *θ*_{
i
}=*x*^{(i)}, and *e* is an *L*^{2}*N*×1 zero vector. Each *D*_{
i
} is an *L*^{2}*N*×*N* matrix with *L*^{2} blocks of *N*×*N* matrices. Each block of *D*_{
i
} can be defined as follows. Consider an *L*×*L* matrix *U*^{(i)}, whose $({i}^{\prime},{j}^{\prime})$th entry ${U}_{{i}^{\prime}{j}^{\prime}}^{\left(i\right)}=1$ if ${j}^{\prime}=i$ and ${i}^{\prime}\in {\mathcal{N}}_{i}$, ${U}_{{i}^{\prime}{j}^{\prime}}^{\left(i\right)}=-1$ if ${i}^{\prime}=i$ and ${j}^{\prime}\in {\mathcal{N}}_{i}$, and ${U}_{{i}^{\prime}{j}^{\prime}}^{\left(i\right)}=0$ otherwise. The $({i}^{\prime}+L{j}^{\prime}-L)$th block of *D*_{
i
}is ${U}_{{i}^{\prime}{j}^{\prime}}^{\left(i\right)}{I}_{N}$, where *I*_{
N
} is an *N*×*N* identity matrix. Substituting them to the multi-block ADM, we can find that in optimizing *x*^{(i)}, agent *i* only needs its local information as well as part of *q*; on the other hand, to update its corresponding part of *q* and *λ*, each agent only needs based on the information from itself and its neighbors. The resulting algorithm is hence fully decentralized due to the nice structure of {*D*_{
i
}}.

At time *t*, the multi-block ADM works as follows: *Step 1: Updating the auxiliary variables* {*q*_{
ij
}}:

*Step 2: Optimizing the local copies*{*x*^{(i)}}:

*Step 3: Updating the Lagrange multipliers* {*λ*_{
ij
}}:

Note that *β* and *μ* are positive constant parameters used by the multi-block ADM.

The updating rules (8), (9), and (10) can also be further simplified. Since we often set {*λ*_{
ij
}(0)} as 0, (8) and (10) imply that *q*_{
ij
}(*t* + 1)=−*q*_{
ji
}(*t* + 1) and *λ*_{
ij
}(*t* + 1)=−*λ*_{
ji
}(*t* + 1). Summing up the two sides of (8) and (10) and defining a new auxiliary variable ${q}_{i}=\sum _{j\in {\mathcal{N}}_{i}}{q}_{\mathrm{ij}}$ as well as a new Lagrange multiplier ${\lambda}_{i}=\sum _{j\in {\mathcal{N}}_{i}}{\lambda}_{\mathrm{ij}}$, their updating rules are:

Hence (9) simplifies to:

### State transition equation of MB-ADM

Combining the updating rules of *q*_{
i
}(*t* + 1) and *q*_{
i
}(*t*) in (11) and the updating rule of *λ*_{
i
}(*t*) in (12), we get:

Substituting ${f}_{i}\left({x}^{\left(i\right)}\right)=\frac{1}{2}\left|\right|{x}^{\left(i\right)}-{b}_{i}|{|}_{2}^{2}$ into (13), the optimality condition for *x*^{(i)}(*t* + 1) is:

Combining the optimality conditions of *x*^{(i)}(*t* + 1) and *x*^{(i)}(*t*) and (14) leads to:

The initial state is ${s}_{M}\left(1\right)=\left[\frac{{b}_{1}}{1+2\mu \left|{\mathcal{N}}_{1}\right|},\dots ,\frac{{b}_{L}}{1+2\mu \left|{\mathcal{N}}_{L}\right|},\right.$$\left(\right)close="]">0,\dots ,0\phantom{\rule{-8.0pt}{0ex}}$ when each agent *i* initializes *q*_{
i
}(0)=0, *x*^{(i)}(0)=0, and *λ*_{
i
}(0)=0.

### Proof of Property 1 in Proposition 1

It is straightforward to show that *ρ*_{M 1}=1 is an eigenvalue of *Φ*_{
M
}, as well as *l*_{M 1} and *r*_{M 1} are its corresponding left and right eigenvectors. Next we prove that *ρ*_{M 1} = 1 is with multiplicity 1 by contradiction. If *ρ*_{M 1} = 1 belongs to a larger Jordan block, there exists a vector ${[{w}^{T},{\stackrel{\u0304}{w}}^{T}]}^{T}$, such that ${\Phi}_{M}{[{w}^{T},{\stackrel{\u0304}{w}}^{T}]}^{T}={[{w}^{T},{\stackrel{\u0304}{w}}^{T}]}^{T}+{[1,\dots ,1,1,\dots ,1]}^{T}$. Here *w* and $\stackrel{\u0304}{w}$ are both *L*×1 vectors (see [36], Fact 2). Observing the lower half of *Φ*_{
M
}, apparently $\stackrel{\u0304}{w}=w-{[1,\dots ,1]}^{T}$. Suppose that *w*_{
k
}has the largest real part among all elements of *w*. Picking up the *k* th row of ${\Phi}_{M}{[{w}^{T},{\stackrel{\u0304}{w}}^{T}]}^{T}$$={[{w}^{T},{\stackrel{\u0304}{w}}^{T}]}^{T}+{[1,\dots ,1,1,\dots ,1]}^{T}$, we have:

or equivalent to:

Denote the real part of *w*_{
k
}and *w*_{
j
} as Re(*w*_{
k
}) and Re(*w*_{
j
}), respectively. Recalling that Re(*w*_{
k
})≥Re(*w*_{
j
}) and picking up the real part of (18), we have:

This leads to contradiction. Hence *ρ*_{M 1}=1 is an eigenvalue of *Φ*_{
M
}with multiplicity 1.

### Proof of Proposition 2

Denote the *i* th eigenvalue of *Φ*_{
i
}as *ρ*_{
Mi
}. Apparently, its eigenvectors should have the form of [*ρ*_{
Mi
}*v*^{T},*v*^{T}]^{T}where *v*^{T}=[*v*_{1},…,*v*_{
L
}]^{T} is a nonzero vector, since the lower half of *Φ*_{
M
}is [*I*_{L×L},0_{L×L}]. Suppose that *v*_{
k
} has the largest norm (here we use |·| to denote the norm of a complex number) among all elements of *v*. Then picking up the *k* th row of *Φ*_{
M
}[*ρ*_{
Mi
}*v*^{T},*v*^{T}]^{T} = *ρ*_{
Mi
}[*ρ*_{
Mi
}*v*^{T},*v*^{T}]^{T}, we have:

or equivalently:

Since *v*_{
k
}has the largest norm among all elements of *v*, taking norms for the both sides of (21) leads to:

Notice that the inequalities turn to equalities when and only when ${v}_{j}={v}_{k},\forall j\in {\mathcal{N}}_{k}$. As *v*_{
k
} has the largest norm among all elements of *v*, any *v*_{
j
}with $j\in {\mathcal{N}}_{k}$ also has such inequalities, and the inequalities turn to equalities when and only when ${v}_{{j}^{\prime}}={v}_{j},\forall {j}^{\prime}\in {\mathcal{N}}_{j}$. Because the network is connected, we can deduce that these inequalities turn to equalities when and only when {*v*_{
i
}} are all equal. This corresponds to the eigenvalue *ρ*_{M 1}=1. Canceling |*v*_{
k
}| from the both sides and defining ${d}_{1}=\frac{2\beta \left|{\mathcal{N}}_{k}\right|}{1+2\mu \left|{\mathcal{N}}_{k}\right|}$ and ${d}_{2}=\frac{2\mu \left|{\mathcal{N}}_{k}\right|}{1+2\mu \left|{\mathcal{N}}_{k}\right|}$, (22) is equivalent to:

Let us consider the two nontrivial special cases.

#### Case 1

The parameters *μ* and *β* are chosen such that *μ*=2*β*>0; further, $\frac{2\beta \left|{\mathcal{N}}_{j}\right|}{1+2\mu \left|{\mathcal{N}}_{j}\right|}<\frac{1}{4}$ and $\frac{2\mu \left|{\mathcal{N}}_{j}\right|}{1+2\mu \left|{\mathcal{N}}_{j}\right|}<\frac{1}{2}$ for all *j*=1,2,…,*L*.

In this case, ${d}_{1}=\frac{2\beta \left|{\mathcal{N}}_{k}\right|}{1+2\mu \left|{\mathcal{N}}_{k}\right|}<\frac{1}{4}$ and ${d}_{2}=\frac{2\mu \left|{\mathcal{N}}_{k}\right|}{1+2\mu \left|{\mathcal{N}}_{k}\right|}<\frac{1}{2}$. Let us choose $d={d}_{1}=\frac{{d}_{2}}{2}$, $\frac{1}{4}>d>0$. Hence, (23) simplifies to:

Define $w={\rho}_{\mathrm{Mi}}-\frac{1}{2}$, we have:

Recall that $\left|w\right|=\frac{1}{2}$ only for *ρ*_{M 1}=1. For any other eigenvalues, $\left|w\right|<\frac{1}{2}$, and hence |*ρ*_{
Mi
}|<1 for *i*≠1.

#### Case 2

The parameters *μ* and *β* are chosen such that *μ*=*β*>0; further, $\frac{2\beta \left|{\mathcal{N}}_{j}\right|}{1+2\mu \left|{\mathcal{N}}_{j}\right|}<\frac{1}{2}$ and $\frac{2\mu \left|{\mathcal{N}}_{j}\right|}{1+2\mu \left|{\mathcal{N}}_{j}\right|}<\frac{1}{2}$ for all *j*=1,2,…,*L*.

In this case, ${d}_{1}=\frac{2\beta \left|{\mathcal{N}}_{k}\right|}{1+2\mu \left|{\mathcal{N}}_{k}\right|}<\frac{1}{2}$ and ${d}_{2}=\frac{2\mu \left|{\mathcal{N}}_{k}\right|}{1+2\mu \left|{\mathcal{N}}_{k}\right|}<\frac{1}{2}$. Let us choose *d*=*d*_{1}=*d*_{2}, $\frac{1}{2}>d>0$. Hence, (23) simplifies to:

Let us prove the conclusion by contradiction. Suppose that there exists a *ρ*_{
Mi
} with |*ρ*_{
Mi
}|≥1 satisfies (26), then:

Again, the inequalities turns to equalities only for *ρ*_{M 1}=1. For any other eigenvalue *ρ*_{
Mi
}, we have |*ρ*_{
Mi
}|<1 which contradicts with |*ρ*_{
Mi
}|≥1. Therefore, |*ρ*_{
Mi
}|<1 for *i*≠1.

### State transition equation of TB-ADM

Substituting ${f}_{i}\left({x}^{\left(i\right)}\right)=\frac{1}{2}\left|\right|{x}^{\left(i\right)}-{b}_{i}|{|}_{2}^{2}$ into:

the optimality condition for *x*^{(i)}(*t* + 1) is:

Considering *x*^{(i)}(*t*), the optimality condition is correspondingly:

Combining (28) and (29) with:

the state transition equation for agent *i* is:

## References

- 1.
Ren W, Beard R, Atkins E: Information consensus in multivehicle cooperative control: collective group behavier through local interaction.

*IEEE Control Systs. Mag*2007, 27: 71-82. - 2.
Olfati-Saber R: Kalman-consensus filter: optimality, stability, and performance. In

*Proceedings of CDC*. Shanghai, China; 2009:7036-7042. - 3.
Boyd S, Parikh N, Chu E, Peleato B, Eckstein J: Distributed optimization and statistical learning via the alternating direction method of multipliers.

*Foundation Trends Mach. Learn*2010, 3: 1-122. 10.1561/2200000016 - 4.
Tsitsiklis J:

*Problems in decentralized decision making and computation*. MIT, Ph.D Thesis; 1984. - 5.
Nedic A, Ozdaglar A: Distributed subgradient methods for multi-agent optimization.

*IEEE Trans. Autom. Control*2009, 54: 48-61. - 6.
Srivastava K, Nedic A: Distrbited asynchronous constrained stochastic optimization.

*IEEE J. Sel. Topics Signal Process*2011, 5: 772-790. - 7.
Rabbat M, Nowak R: Distributed optimization in sensor networks. In

*Proceedings of IPSN*. Berkeley, USA; 2004:20-27. - 8.
Rabbat M, Nowak R: Quantized incremental algorithms for distributed optimization.

*IEEE J. Sel. Areas Commun*2006, 23: 798-808. - 9.
Xiao L, Boyd S, Kim S: Distributed average consensus with least-mean-square deviation.

*J. Parallel Distrib. Comput*2007, 67: 33-46. 10.1016/j.jpdc.2006.08.010 - 10.
Kar S, Moura J: Distributed consensus algorithms in sensor networks: quantized data and random link failures.

*IEEE Trans. Signal Process*2010, 58: 1383-1400. - 11.
Olshevsky A:

*Efficient Information Aggregation Strategies for Distributed Control and Signal Processing*. Ph.D Thesis, MIT; 2010. - 12.
Schizas I, Ribeiro A, Giannakis G: Consensus in ad hoc WSNs with noisy links - Part I: distributed estimation of deterministic signals.

*IEEE Trans. Signal Process*2008, 56: 350-364. - 13.
Mateos G, Schizas I, Giannakis G: Distributed recursive least-squares for consensus-based in-network adaptive estimation.

*IEEE Trans. Signal Process*2009, 57: 4583-4588. - 14.
Ling Q, Tian Z: Decentralized sparse signal recovery for compressive sleeping wireless sensor networks.

*IEEE Trans. Signal Process*2010, 58: 3816-3827. - 15.
Bazerque J, Giannakis G: Distributed spectrum sensing for cognitive radio networks by exploiting sparsity.

*IEEE Trans. Signal Process*2010, 58: 1847-1862. - 16.
Mateos G, Bazerque J, Giannakis G: Distributed sparse linear regression.

*IEEE Trans. Signal Process*2010, 58: 5262-5276. - 17.
Jakovetic D, Xavier J, Moura J: Cooperative convex optimization in networked systems: augmented Lagrangian algorithms with direct gossip communication.

*IEEE Trans. Signal Process*2011, 59: 3889-3902. - 18.
Cetin M, Chen L, Fisher I I I J, Ihler A, Moss R, Wainwright M, Willsky A: Distributed fusion in sensor networks.

*IEEE Signal Process. Mag*2006, 23: 42-55. - 19.
Predd J, Kulkarni S, Poor V: Distributed learning in wireless sensor networks.

*IEEE Signal Process. Mag*2007, 24: 56-69. - 20.
Predd J, Kulkarni S, Poor H: A collaborative training algorithm for distributed learning.

*IEEE Trans. Inf. Theory*2009, 55: 1856-1871. - 21.
Khan U, Kar S, Moura J: Higher dimensional consensus: learning in large-scale networks.

*IEEE Trans. Signal Process*2010, 58: 2836-2849. - 22.
Jadbabaie A, Ozdaglar A, Zargham M: A distributed Newton method for network optimization. In

*Proceedings of CDC*. Shanghai, China; 2009:2736-2741. - 23.
Koshal J, Nedic A, Shanbhag U: Multiuser optimization: distributed algorithms and error analysis.

*SIAM J. Optimiz*2011, 21: 1046-1081. 10.1137/090770102 - 24.
Wan P, Lemmon M: Distributed network utility maximization using event-triggered augmented Lagrangian methods. In

*Proceedings of ACC*. St. Louis, USA; 2009:3298-3303. - 25.
Zhao F, Shin J, Reich J: Information-driven dynamic sensor collaboration.

*IEEE Signal Process. Mag*2002, 19: 61-72. 10.1109/79.985685 - 26.
Zhao F, Guibas L:

*Wireless Sensor Networks: an Information Processing Approach*. Morgan Kaufmann, Burlington, USA; 2004. - 27.
Schizas I, Giannakis G, Roumeliotis S, Ribeiro A: Consensus in ad hoc WSNs with noisy links – part II: distributed estimation and smoothing of random signals.

*IEEE Trans. Signal Process*2008, 56(4):1650-1666. - 28.
Ribeiro A, Schizas I, Roumeliotis S, Giannakis G: Kalman filtering in wireless sensor networks: reducing communication cost in state estimation problems.

*IEEE Control Systs. Mag*2010, 30: 66-86. - 29.
Tao M: Some parallel splitting methods for separable convex programming with O(1/t) convergence rate. in press

- 30.
He B, Tao M, Xu M, Yuan X: Alternating directions based contraction method for generally separable linearly constrained convex programming problems. in press

- 31.
He B, Tao M, Yuan X: Alternating direction method with Gaussian back substitution for separable convex programming.

*SIAM J. Optim*2012, 22: 313-340. 10.1137/110822347 - 32.
Bertsekas D, Tsitsiklis J:

*Parallel and Distributed Computation: Numerical Methods*. Athena Scientific, Nashua, USA; 1997. - 33.
Bertsekas D:

*Numerical Optimization*. Athena Scientific, Nashua, USA; 1999. - 34.
He B, Yuan X: On the O(1/n) convergence rate of Douglas-Rachford alternating direction method.

*SIAM J. Num. Anal*2012, 50: 700-709. 10.1137/110836936 - 35.
Erseghe T, Zennaro D, Dall’Anese E, Vangelista L: Fast consensus by the alternating direction multipliers method.

*IEEE Trans. Signal Process*2011, 59: 5523-5537. - 36.
Rosenthal J: Convergence rates for Markov chains.

*SIAM Rev*1995, 37: 387-405. 10.1137/1037083 - 37.
Boyd S, Ghosh A, Prabhakar B, Shah D: Randomized gossip algorithms.

*IEEE Trans. Inf. Theory*2006, 52: 2508-2530.

## Acknowledgements

The work of Qing Ling is supported in part by NSFC grant 61004137 and Fundamental Research Funds for the Central Universities. The work of Wotao Yin is supported in part by ARL and ARO grant W911NF-09-1-0383 and NSF grants DMS-0748839 and ECCS-1028790. The work of Xiaoming Yuan is supported in part by the General Research Fund No. 203311 from Hong Kong Research Grants Council.

## Author information

### Affiliations

### Corresponding author

## Additional information

### Competing interests

The authors declare that they have no competing interests.

## Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

## Rights and permissions

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## About this article

### Cite this article

Ling, Q., Tao, M., Yin, W. *et al.* A multi-block alternating direction method with parallel splitting for decentralized consensus optimization.
*J Wireless Com Network* **2012, **338 (2012). https://doi.org/10.1186/1687-1499-2012-338

Received:

Accepted:

Published:

### Keywords

- Convergence Rate
- Communication Cost
- Fusion Center
- Alternate Direction Method
- Average Node Degree