Decentralized Gossip-Based Stochastic Bilevel Optimization over Communication Networks

Shuoguang YangXuezhou ZhangMengdi Wang

article2022NeurIPS66 citations

Proposes a single-timescale gossip-based algorithm for decentralized stochastic bilevel optimization that achieves optimal per-agent sample complexity and a linear speedup with network size for both nonconvex and Polyak-Łojasiewicz objectives.

Listen

Modern distributed machine learning applications—such as hyperparameter tuning, meta-learning, and multi-agent reinforcement learning—increasingly rely on nested optimization structures known as bilevel optimization. In these problems, an outer decision depends directly on the outcome of an inner, lower-level task. In decentralized networks, data is spread across multiple devices or sensors that cannot share raw data due to privacy concerns and lack a central server to coordinate learning. Under these conditions, finding optimal solutions is exceptionally difficult because nodes cannot directly access the global calculations required to compute the necessary gradients.

The article develops and evaluates a decentralized, gossip-based stochastic approximation algorithm designed to solve these nested optimization problems collaboratively across communication networks. The framework enables participating nodes to solve both the inner and outer optimization levels simultaneously in a single timescale without relying on a central coordinator.

To establish mathematical and empirical credibility, the authors conducted theoretical convergence analyses across general nonconvex objectives and structured Polyak-Łojasiewicz conditions, which include strongly convex objectives. Nodes exchange only local parameter estimates with immediate network neighbors using a gossip protocol based on a doubly stochastic weighting matrix. The authors also evaluated the algorithm using simulated experiments on a ring network topology for two benchmark applications: tuning regularizers on a handwriting recognition dataset across up to 20 nodes and policy evaluation for multi-agent reinforcement learning across 100 states.

The analysis produced several key findings. First, the algorithm achieves optimal sample complexity, matching the theoretical performance of traditional single-server systems. It achieves a per-node sample complexity proportional to one over the number of nodes multiplied by the squared error tolerance for general nonconvex objectives, and one over the number of nodes multiplied by the error tolerance for structured objectives. Second, the algorithm demonstrates an exact linear speedup: as the number of network nodes increases, the per-node data samples required to reach a target accuracy decrease proportionally. Third, the mathematical analysis proves that network consensus errors diminish rapidly over time, meaning the specific communication topology does not degrade the long-term convergence rate. Finally, simulated experiments confirmed that the proposed approach converges faster and requires significantly fewer data samples to achieve target accuracy than baseline decentralized methods.

These findings demonstrate that organizations can deploy nested, multi-task machine learning over fully peer-to-peer networks without sacrificing computational efficiency or privacy. By communicating only parameter estimates rather than raw data, the method reduces privacy risks and data transfer costs. Furthermore, eliminating the central server prevents single-point-of-failure vulnerabilities, ensuring the system remains operational even if specific communication channels fail.

Decision-makers can consider implementing this peer-to-peer gossip framework for privacy-sensitive, multi-agent systems and edge-device networks. For future development, the article recommends investigating algorithmic variants with reduced iteration counts and lower per-round communication overhead, which would further optimize operational network bandwidth.

Confidence in these findings is supported by rigorous mathematical proofs and consistent numerical simulations. However, readers should note that the current experimental evaluations rely on simulated desktop environments and specific network configurations, such as ring topologies. Practical deployments across larger, highly irregular networks with real-world latency, packet drops, or asynchronous communication should be validated with targeted pilot testing.

arXiv: 2206.10870
Cover for Decentralized Gossip-Based Stochastic Bilevel Optimization over Communication Networks

Abstract

Bilevel optimization have gained growing interests, with numerous applications being found in meta learning, minimax games, reinforcement learning, and nested composition optimization. This paper studies the problem of decentralized distributed stochastic bilevel optimization over a network where each agent can only communicate with its neighbors, and gives examples from multi-task, multi-agent learning and federated learning. In this paper, we propose a gossip-based decentralized bilevel learning algorithm that allows networked agents to solve both the inner and outer optimization problems in a single timescale and share information through network propagation. We show that our algorithm enjoys the Õ(1/(Kϵ₂)) per-agent sample complexity for general nonconvex bilevel optimization and Õ(1/(Kϵ)) for Polyak-Łojasiewicz objectives, achieving a speedup that scales linearly with the network size K. The sample complexities are optimal in both ϵ and K. We test our algorithm on the examples of hyperparameter tuning and decentralized reinforcement learning. Simulated experiments confirmed that our algorithm achieves the state-of-the-art training efficiency and test accuracy.

Table of Contents

  • 1 Introduction
  • 1.1 Example applications of SBO
  • 1.2 Challenges with Distributed SBO
  • 2 Related Works
  • 3 Problem Setup
  • 4 Algorithm
  • 5 Theory
  • 5.1 Nonconvex Objectives
  • 5.2 PL Objectives
  • 6 Numerical Experiments
  • 7 Conclusion
  • Acknowledgments and Disclosure of Funding
  • References

Knowls

  1. Knowl 1 — Decentralized stochastic bilevel optimization problem

    definition

    The paper formulates decentralized stochastic bilevel optimization over KK agents, where agent kk has local outer objective fkf^k and inner objective gkg^k. The agents jointly solve

    min⁡x∈RdxF(x):=f(x,y⋆(x))=1K∑k=1Kfk(x,y⋆(x)),y⋆(x):=arg⁡min⁡y∈Rdyg(x,y),\min_{x\in\mathbb{R}^{d_x}} F(x):=f\bigl(x,y^\star(x)\bigr)=\frac{1}{K}\sum_{k=1}^{K}f^k\bigl(x,y^\star(x)\bigr), \qquad y^\star(x):=\arg\min_{y\in\mathbb{R}^{d_y}}g(x,y),

    where g(x,y):=1K∑k=1Kgk(x,y)g(x,y):=\frac{1}{K}\sum_{k=1}^{K}g^k(x,y), fk(x,y)=E[fk(x,y;ζk)]f^k(x,y)=\mathbb{E}[f^k(x,y;\zeta^k)], and gk(x,y)=E[gk(x,y;ξk)]g^k(x,y)=\mathbb{E}[g^k(x,y;\xi^k)]. The random variables ζk\zeta^k and ξk\xi^k may have heterogeneous distributions across agents, and each agent can communicate only with its network neighbors.

    When the inner objective is strongly convex in yy, the outer gradient is characterized by

    ∇F(x)=∇xf(x,y⋆(x))−∇xy2g(x,y⋆(x))[∇yy2g(x,y⋆(x))]−1∇yf(x,y⋆(x)),\nabla F(x)=\nabla_x f\bigl(x,y^\star(x)\bigr)-\nabla^2_{xy}g\bigl(x,y^\star(x)\bigr)\left[\nabla^2_{yy}g\bigl(x,y^\star(x)\bigr)\right]^{-1}\nabla_y f\bigl(x,y^\star(x)\bigr),

    where ∇xy2g∈Rdx×dy\nabla^2_{xy}g\in\mathbb{R}^{d_x\times d_y}, ∇yy2g∈Rdy×dy\nabla^2_{yy}g\in\mathbb{R}^{d_y\times d_y}, and ∇yf∈Rdy\nabla_y f\in\mathbb{R}^{d_y}. The decentralized difficulty is that the terms in this expression must be estimated from heterogeneous local data and aggregated through neighbor communication.

  2. Knowl 2 — Oracle, network, and regularity conditions

    assumption

    The theoretical guarantees assume the following conditions for the decentralized bilevel problem. Each agent can query an independent sampling oracle that returns unbiased stochastic estimates of ∇xfk\nabla_x f^k, ∇yfk\nabla_y f^k, ∇ygk\nabla_y g^k, ∇xy2gk\nabla^2_{xy}g^k, and ∇yy2gk\nabla^2_{yy}g^k at any (x,y)(x,y). The communication network is represented by a symmetric, nonnegative, doubly stochastic matrix W∈RK×KW\in\mathbb{R}^{K\times K}, and its disagreement factor satisfies

    ∥W−1K11T∥2=ρ<1,\left\|W-\frac{1}{K}\mathbf{1}\mathbf{1}^{\mathsf T}\right\|_2=\rho<1,

    where 1∈RK\mathbf{1}\in\mathbb{R}^{K} is the all-ones vector and ∥⋅∥2\|\cdot\|_2 is the spectral norm.

    The outer functions have an optimal solution, their gradients with respect to xx and yy are LfL_f-Lipschitz in (x,y)(x,y), and their stochastic gradient estimates have uniformly bounded second moments by Cf2C_f^2. The averaged inner objective g(x,y)g(x,y) is μg\mu_g-strongly convex in yy. Each gkg^k is twice continuously differentiable; ∇ygk\nabla_y g^k, ∇xy2gk\nabla^2_{xy}g^k, and ∇yy2gk\nabla^2_{yy}g^k are Lipschitz in (x,y)(x,y), with constants LgL_g and L~g\widetilde L_g as appropriate. The stochastic inner gradients and Hessians have bounded second moments, including bounds Cg2C_g^2 for ∇ygk\nabla_y g^k and Lg2L_g^2 for the Hessian estimates. Finally, for some 0<κg≤μg/Lg≤10<\kappa_g\leq \mu_g/L_g\leq 1, the Hessian samples satisfy

    E[∥I−1Lg∇yy2gk(x,y;ξk)∥22]≤(1−κg)2,\mathbb{E}\left[\left\|I-\frac{1}{L_g}\nabla^2_{yy}g^k(x,y;\xi^k)\right\|_2^2\right]\leq (1-\kappa_g)^2,

    where I∈Rdy×dyI\in\mathbb{R}^{d_y\times d_y} is the identity matrix. These conditions permit heterogeneity in both the outer and inner objectives.

  3. Knowl 3 — Gossip-based single-timescale DSBO algorithm

    algorithm

    The proposed method maintains, at every agent kk, primal variables (xtk,ytk)(x_t^k,y_t^k) and network-averaged estimators stk≈∇xfs_t^k\approx\nabla_x f, htk≈∇yfh_t^k\approx\nabla_y f, utk≈∇xy2gu_t^k\approx\nabla^2_{xy}g, and qtk≈[∇yy2g]−1q_t^k\approx[\nabla^2_{yy}g]^{-1}. Here xtk∈Rdxx_t^k\in\mathbb{R}^{d_x}, ytk∈Rdyy_t^k\in\mathbb{R}^{d_y}, stk∈Rdxs_t^k\in\mathbb{R}^{d_x}, htk∈Rdyh_t^k\in\mathbb{R}^{d_y}, utk∈Rdx×dyu_t^k\in\mathbb{R}^{d_x\times d_y}, and qtk∈Rdy×dyq_t^k\in\mathbb{R}^{d_y\times d_y}. Let Nk\mathcal{N}_k be agent kk's neighbors, wkjw_{kj} the entries of the gossip matrix, αt\alpha_t the outer stepsize, βt\beta_t the inner stepsize, γt\gamma_t the estimator stepsize, and bb the number of Hessian samples used for the inverse-Hessian approximation.

    At each iteration, every agent samples local first- and second-order information, gossips each estimator with its neighbors, updates the inner and outer variables, and constructs the inverse-Hessian estimator as follows:

    Input: Stepsizes {αt}\{\alpha_t\}, {βt}\{\beta_t\}, {γt}\{\gamma_t\}, iterations TT, sampling oracle, gossip matrix WW, Hessian bound LgL_g, and truncation length bb
    Initialize x0k=y0k=s0k=h0k=u0k=0x_0^k=y_0^k=s_0^k=h_0^k=u_0^k=0 and v0,ik=μgIv_{0,i}^k=\mu_g I for every agent kk and i=1,…,bi=1,\ldots,b
    For t=0,…,T−1t=0,\ldots,T-1
        At every agent kk, sample ∇xfk\nabla_x f^k, ∇yfk\nabla_y f^k, ∇ygk\nabla_y g^k, ∇xy2gk\nabla^2_{xy}g^k, and independent ∇yy2gk\nabla^2_{yy}g^k samples at (xtk,ytk)(x_t^k,y_t^k)
        xt+1k=∑j∈Nkwkjxtj−αt(stk−utkqtkhtk)x_{t+1}^k = \sum_{j\in\mathcal{N}_k}w_{kj}x_t^j-\alpha_t(s_t^k-u_t^kq_t^kh_t^k)
        yt+1k=∑j∈Nkwkjytj−βt∇ygk(xtk,ytk)y_{t+1}^k = \sum_{j\in\mathcal{N}_k}w_{kj}y_t^j-\beta_t\nabla_y g^k(x_t^k,y_t^k)
        st+1k=(1−γt)∑j∈Nkwkjstj+γt∇xfk(xtk,ytk)s_{t+1}^k=(1-\gamma_t)\sum_{j\in\mathcal{N}_k}w_{kj}s_t^j+\gamma_t\nabla_x f^k(x_t^k,y_t^k)
        ht+1k=(1−γt)∑j∈Nkwkjhtj+γt∇yfk(xtk,ytk)h_{t+1}^k=(1-\gamma_t)\sum_{j\in\mathcal{N}_k}w_{kj}h_t^j+\gamma_t\nabla_y f^k(x_t^k,y_t^k)
        ut+1k=(1−γt)∑j∈Nkwkjutj+γt∇xy2gk(xtk,ytk)u_{t+1}^k=(1-\gamma_t)\sum_{j\in\mathcal{N}_k}w_{kj}u_t^j+\gamma_t\nabla^2_{xy}g^k(x_t^k,y_t^k)
        Set Qt+1,0k=IQ_{t+1,0}^k=I
        For i=1,…,bi=1,\ldots,b
            vt+1,ik=(1−γt)∑j∈Nkwkjvt,ij+γt∇yy2gk(xtk,ytk)v_{t+1,i}^k=(1-\gamma_t)\sum_{j\in\mathcal{N}_k}w_{kj}v_{t,i}^j+\gamma_t\nabla^2_{yy}g^k(x_t^k,y_t^k)
            Qt+1,ik=I+(I−Lg−1vt+1,ik)Qt+1,i−1kQ_{t+1,i}^k=I+(I-L_g^{-1}v_{t+1,i}^k)Q_{t+1,i-1}^k
        End for
        qt+1k=Lg−1Qt+1,bkq_{t+1}^k=L_g^{-1}Q_{t+1,b}^k
    End for
    Output xˉt=K−1∑k=1Kxtk\bar{x}_t=K^{-1}\sum_{k=1}^{K}x_t^k

    The outer update uses the estimated bilevel gradient stk−utkqtkhtks_t^k-u_t^kq_t^kh_t^k, while the inner update is a stochastic gradient step on the local inner objective combined with gossip averaging. All agents perform the inner and outer updates on the same iteration timescale.

  4. Knowl 4 — Truncated stochastic inverse-Hessian estimator

    model/method

    The algorithm estimates the inverse of the averaged inner Hessian without requiring an unbiased inverse-Hessian sample. For a positive-definite Hessian H=∇yy2g(x,y)H=\nabla^2_{yy}g(x,y) and a scalar LgL_g satisfying the paper's Hessian bounds, the inverse is approximated by the truncated Neumann series

    H−1=1Lg(I−(I−1LgH))−1≈1Lg∑j=0b(I−1LgH)j,H^{-1}=\frac{1}{L_g}\left(I-\left(I-\frac{1}{L_g}H\right)\right)^{-1} \approx \frac{1}{L_g}\sum_{j=0}^{b}\left(I-\frac{1}{L_g}H\right)^j,

    where bb is the truncation length. Because the local inverse Hessians do not average to the inverse of the averaged Hessian, the method instead constructs bb independently gossiped Hessian estimators vt,1k,…,vt,bkv_{t,1}^k,\ldots,v_{t,b}^k. It recursively forms

    Qt,0k=I,Qt,ik=I+(I−1Lgvt,ik)Qt,i−1k,qtk=1LgQt,bk.Q_{t,0}^k=I, \qquad Q_{t,i}^k=I+\left(I-\frac{1}{L_g}v_{t,i}^k\right)Q_{t,i-1}^k, \qquad q_t^k=\frac{1}{L_g}Q_{t,b}^k.

    The resulting qtk∈Rdy×dyq_t^k\in\mathbb{R}^{d_y\times d_y} is used in utkqtkhtku_t^kq_t^kh_t^k to estimate the Hessian-vector correction in the bilevel gradient. The paper chooses b=Θ(log⁡T)b=\Theta(\log T) in both convergence regimes, making the truncation error decrease sufficiently quickly while requiring only logarithmically many additional Hessian samples per iteration.

  5. Knowl 5 — Nonconvex convergence and linear sample speedup

    theoretical result

    For a generally nonconvex outer objective FF, use constant stepsizes over a horizon of TT iterations:

    αt=C0KT,βt=γt=KT,b=Θ(log⁡T),\alpha_t=C_0\sqrt{\frac{K}{T}}, \qquad \beta_t=\gamma_t=\sqrt{\frac{K}{T}}, \qquad b=\Theta(\log T),

    where C0>0C_0>0 is sufficiently small. Under the sampling, network, smoothness, bounded-moment, and strong-convexity conditions stated for the decentralized bilevel problem, the network average xˉt=K−1∑k=1Kxtk\bar{x}_t=K^{-1}\sum_{k=1}^{K}x_t^k satisfies

    1T∑t=0T−1E[∥∇F(xˉt)∥2]≤O(1KT)+O(KT(1−ρ)2).\frac{1}{T}\sum_{t=0}^{T-1}\mathbb{E}\left[\left\|\nabla F(\bar{x}_t)\right\|^2\right] \leq O\left(\frac{1}{\sqrt{KT}}\right) +O\left(\frac{K}{T(1-\rho)^2}\right).

    The first term gives the main stochastic optimization rate, while the second is the network-consensus error. For a fixed connected network, the consensus term is lower order as TT grows, so the asymptotic rate is independent of the network topology. Ignoring logarithmic factors from bb, obtaining an ε\varepsilon-stationary point in the criterion T−1∑t=0T−1E[∥∇F(xˉt)∥2]≤εT^{-1}\sum_{t=0}^{T-1}\mathbb{E}[\|\nabla F(\bar{x}_t)\|^2]\leq\varepsilon requires O~(1/(Kε2))\widetilde O(1/(K\varepsilon^2)) samples per agent, exhibiting a linear speedup in KK.

  6. Knowl 6 — PL convergence and optimal sample complexity

    theoretical result

    Suppose the outer objective satisfies the Polyak–Łojasiewicz condition

    2μ(F(x)−F⋆)≤∥∇F(x)∥22\mu\bigl(F(x)-F^\star\bigr)\leq\|\nabla F(x)\|^2

    for some μ>0\mu>0, where F⋆=min⁡xF(x)F^\star=\min_xF(x). This class includes strongly convex objectives. With diminishing stepsizes

    αt=2μ(C1+t),βt=γt=C1C1+t,b=Θ(log⁡T),\alpha_t=\frac{2}{\mu(C_1+t)}, \qquad \beta_t=\gamma_t=\frac{C_1}{C_1+t}, \qquad b=\Theta(\log T),

    where C1C_1 is sufficiently large, the algorithm satisfies

    E[F(xˉT)]−F⋆≤O(1KT)+O(log⁡TT2(1−ρ)2).\mathbb{E}\left[F(\bar{x}_T)\right]-F^\star \leq O\left(\frac{1}{KT}\right) +O\left(\frac{\log T}{T^2(1-\rho)^2}\right).

    The network-consensus term decays faster than the leading O(1/(KT))O(1/(KT)) term for a fixed connected network. Consequently, finding an ε\varepsilon-optimal point with expected gap at most ε\varepsilon requires O(1/(Kε))O(1/(K\varepsilon)) iterations and O~(1/(Kε))\widetilde O(1/(K\varepsilon)) stochastic samples per agent. Thus the method achieves linear speedup with the number of agents and the optimal dependence on ε\varepsilon claimed for decentralized stochastic bilevel optimization.

  7. Knowl 7 — Communication properties and central-server specialization

    model/method

    Each agent communicates estimator states and model variables with only its neighbors rather than transmitting raw local data. In one iteration, agent kk performs O(∣Nk∣)O(|\mathcal{N}_k|) neighbor communications, where ∣Nk∣|\mathcal{N}_k| is its neighborhood size, instead of communicating with all KK agents. The algorithm therefore preserves local-data privacy at the communication level and can continue operating after an individual communication link fails, provided the remaining network stays connected.

    A central-server version is obtained by replacing neighbor gossip with synchronous aggregation by a server. In that setting all agents receive a common iterate, the consensus-error terms disappear, and the corresponding variant retains the nonconvex convergence rate O(1/KT)O(1/\sqrt{KT}).

  8. Knowl 8 — Federated hyperparameter-optimization experiment

    experimental setup

    The paper evaluates the method on federated hyperparameter optimization for binary handwriting recognition using the Australia handwriting dataset. Each example has a feature vector wi∈R14w_i\in\mathbb{R}^{14} and label zi∈{0,1}z_i\in\{0,1\}. The outer validation objective uses a sigmoid-based loss, while the inner training objective uses the training loss plus the strongly convex regularizer

    R(x,y)=∑i=1dxi22∥yi∥2,R(x,y)=\sum_{i=1}^{d}\frac{x_i^2}{2}\|y_i\|^2,

    where xx is the hyperparameter vector and yy is the model parameter vector. The data are randomly split into training and validation sets and then distributed across agents. The communication network is a ring with self-weight and neighbor weights wij=1/3w_{ij}=1/3 for j∈{i−1,i,i+1}j\in\{i-1,i,i+1\}.

    The proposed method uses b=200b=200, T=20,000T=20{,}000, αt=0.1K/T\alpha_t=0.1\sqrt{K/T}, and βt=γt=10K/T\beta_t=\gamma_t=10\sqrt{K/T}. Ten independent simulations compare the method with a decentralized extension of the double-loop BSA baseline. For K=5K=5, the proposed method reaches a given validation accuracy with fewer total samples than the baseline. Runs with K=5,10,20K=5,10,20 converge faster as KK increases, supporting the predicted decentralized speedup.

  9. Knowl 9 — Decentralized policy-evaluation experiment

    experimental setup

    The second experiment studies distributed policy evaluation in a finite-state Markov decision process. Let S\mathcal{S} be the state space, let ϕs∈Rm\phi_s\in\mathbb{R}^{m} be a state feature, and approximate the value function by V(s)=ϕsTxV(s)=\phi_s^{\mathsf T}x. For discount factor γ∈(0,1)\gamma\in(0,1) and regularization coefficient λ>0\lambda>0, the regularized Bellman objective is

    F(x)=12∣S∣∑s∈S(ϕsTx−Es′[r(s,s′)+γϕs′Tx∣s])2+λ2∥x∥2.F(x)=\frac{1}{2|\mathcal{S}|}\sum_{s\in\mathcal{S}} \left(\phi_s^{\mathsf T}x-\mathbb{E}_{s'}\left[r(s,s')+\gamma\phi_{s'}^{\mathsf T}x\mid s\right]\right)^2 +\frac{\lambda}{2}\|x\|^2.

    With heterogeneous reward functions rkr^k at KK agents, the inner response for state ss is

    ys⋆(x)=ϕsTx−Es′[1K∑k=1Krk(s,s′)+γϕs′Tx∣s], y_s^\star(x)=\phi_s^{\mathsf T}x-\mathbb{E}_{s'}\left[\frac{1}{K}\sum_{k=1}^{K}r^k(s,s')+\gamma\phi_{s'}^{\mathsf T}x\mid s\right],

    and the objective can be written as f(x,y⋆(x))=(2∣S∣)−1∑s∈S(ys⋆(x))2+(λ/2)∥x∥2f(x,y^\star(x))=(2|\mathcal{S}|)^{-1}\sum_{s\in\mathcal{S}}(y_s^\star(x))^2+(\lambda/2)\|x\|^2. The experiment uses ∣S∣=100|\mathcal{S}|=100, λ=1\lambda=1, ring networks with K=5,10,20K=5,10,20, and ten independent simulations per network size. The proposed method is compared with a double-loop DSGD baseline.

    The proposed method outperforms DSGD for K=5K=5. The log mean-squared error decreases with an empirical slope close to −1-1, consistent with the O(1/t)O(1/t) PL/strongly-convex rate. Achieving an error threshold of 10−610^{-6} requires approximately the same total number of samples for K=5,10,20K=5,10,20, so the samples required per agent decrease approximately linearly with KK.

  10. Knowl 10 — Remaining efficiency limitations

    limitation

    The paper establishes convergence and sample-complexity guarantees but does not provide a lower-iteration-complexity or lower-communication-cost algorithm. The authors identify reducing the number of iterations and the communication burden as future work. The reported experiments are also conducted primarily on simulated ring networks; larger network sizes and fully connected or randomly connected topologies are treated only as supplementary extensions rather than as a full characterization of topology-dependent practical costs.

Coverage note — Related-work comparisons, background applications, proof derivations, and supplementary topology experiments were omitted because they do not add separate load-bearing contributions beyond the extracted methods, guarantees, and main experiments.

References

  1. 1.Sanjeev Arora, Simon Du, Sham Kakade, Yuping Luo, and Nikunj Saunshi. Provable representation learning for imitation learning via bi-level optimization. In International Conference on Machine Learning, pages 367–376. PMLR, 2020.
  2. 2.Luca Bertinetto, Joao F Henriques, Philip HS Torr, and Andrea Vedaldi. Meta-learning with differentiable closed-form solvers. arXiv preprint arXiv:1805.08136, 2018.
  3. 3.Jerome Bracken and James T McGill. Mathematical programs with optimization problems in the constraints. Operations Research, 21(1):37–44, 1973.
  4. 4.Chih-Chung Chang and Chih-Jen Lin. Libsvm: a library for support vector machines. ACM transactions on intelligent systems and technology (TIST), 2(3):1–27, 2011. URL https://www.csie.ntu.edu.tw/~cjlin/libsvmtools/datasets/.
  5. 5.Tianyi Chen, Yuejiao Sun, and Wotao Yin. A single-timescale stochastic bilevel optimization method. arXiv preprint arXiv:2102.04671, 2021a.
  6. 6.Tianyi Chen, Yuejiao Sun, and Wotao Yin. Solving stochastic compositional optimization is nearly as easy as solving stochastic optimization. IEEE Transactions on Signal Processing, 69:4937–4948, 2021b.
  7. 7.Tianyi Chen, Yuejiao Sun, and Wotao Yin. Tighter analysis of alternating stochastic gradient method for stochastic nested problems. arXiv preprint arXiv:2106.13781, 2021c.
  8. 8.Nicolas Couellan and Wenjuan Wang. On the convergence of stochastic bi-level gradient methods. Optimization, 2016.
  9. 9.Evin J Cramer, John E Dennis, Jr, Paul D Frank, Robert Michael Lewis, and Gregory R Shubin. Problem formulation for multidisciplinary optimization. SIAM Journal on Optimization, 4(4):754–776, 1994.
  10. 10.Ofer Dekel, Ran Gilad-Bachrach, Ohad Shamir, and Lin Xiao. Optimal distributed online prediction using mini-batches. Journal of Machine Learning Research, 13(1), 2012.
  11. 11.Hamid Reza Feyzmahdavian, Arda Aytekin, and Mikael Johansson. An asynchronous mini-batch algorithm for regularized stochastic optimization. IEEE Transactions on Automatic Control, 61(12):3740–3754, 2016.
  12. 12.Luca Franceschi, Paolo Frasconi, Saverio Salzo, Riccardo Grazzi, and Massimiliano Pontil. Bilevel programming for hyperparameter optimization and meta-learning. In International Conference on Machine Learning, pages 1568–1577. PMLR, 2018.
  13. 13.Hongchang Gao and Heng Huang. Periodic stochastic gradient descent with momentum for decentralized training. arXiv preprint arXiv:2008.10435, 2020.
  14. 14.Jason Ge, Zhaoran Wang, Mengdi Wang, and Han Liu. Minimax-optimal privacy-preserving sparse pca in distributed systems. In International Conference on Artificial Intelligence and Statistics, pages 1589–1598. PMLR, 2018.
  15. 15.Saeed Ghadimi and Guanghui Lan. Accelerated gradient methods for nonconvex nonlinear and stochastic programming. Mathematical Programming, 156(1):59–99, 2016.
  16. 16.Saeed Ghadimi and Mengdi Wang. Approximation methods for bilevel programming. arXiv preprint arXiv:1802.02246, 2018.
  17. 17.Saeed Ghadimi, Andrzej Ruszczynski, and Mengdi Wang. A single timescale stochastic approximation method for nested stochastic optimization. SIAM Journal on Optimization, 30(1):960–979, 2020.
  18. 18.Zhishuai Guo, Quanqi Hu, Lijun Zhang, and Tianbao Yang. Randomized stochastic variance-reduced methods for multi-task stochastic bilevel optimization. arXiv preprint arXiv:2105.02266, 2021.
  19. 19.Pierre Hansen, Brigitte Jaumard, and Gilles Savard. New branch-and-bound rules for linear bilevel programming. SIAM Journal on scientific and Statistical Computing, 13(5):1194–1217, 1992.
  20. 20.Mingyi Hong, Hoi-To Wai, Zhaoran Wang, and Zhuoran Yang. A two-timescale framework for bilevel optimization: Complexity analysis and application to actor-critic. arXiv preprint arXiv:2007.05170, 2020.
  21. 21.Kaiyi Ji, Junjie Yang, and Yingbin Liang. Bilevel optimization: Convergence analysis and enhanced design. In International Conference on Machine Learning, pages 4882–4892. PMLR, 2021.
  22. 22.Prashant Khanduri, Siliang Zeng, Mingyi Hong, Hoi-To Wai, Zhaoran Wang, and Zhuoran Yang. A near-optimal algorithm for stochastic bilevel optimization via double-momentum. Advances in Neural Information Processing Systems, 34, 2021.
  23. 23.Anastasia Koloskova, Tao Lin, Sebastian U Stich, and Martin Jaggi. Decentralized deep learning with arbitrary communication compression. arXiv preprint arXiv:1907.09356, 2019.
  24. 24.Guanghui Lan and Yi Zhou. Random gradient extrapolation for distributed and stochastic optimization. SIAM Journal on Optimization, 28(4):2753–2782, 2018.
  25. 25.Guanghui Lan, Soomin Lee, and Yi Zhou. Communication-efficient algorithms for decentralized and stochastic optimization. Mathematical Programming, 180(1):237–284, 2020.
  26. 26.Xiangru Lian, Mengdi Wang, and Ji Liu. Finite-sum composition optimization via variance reduced gradient descent. In Artificial Intelligence and Statistics, pages 1159–1167. PMLR, 2017a.
  27. 27.Xiangru Lian, Ce Zhang, Huan Zhang, Cho-Jui Hsieh, Wei Zhang, and Ji Liu. Can decentralized algorithms outperform centralized algorithms? a case study for decentralized parallel stochastic gradient descent. Advances in Neural Information Processing Systems, 30, 2017b.
  28. 28.Xiangru Lian, Wei Zhang, Ce Zhang, and Ji Liu. Asynchronous decentralized parallel stochastic gradient descent. In International Conference on Machine Learning, pages 3043–3052. PMLR, 2018.
  29. 29.Jialin Liu, Yuantao Gu, and Mengdi Wang. Averaging random projection: A fast online solution for large-scale constrained stochastic optimization. In 2015 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), pages 3586–3590, 2015.
  30. 30.Eli Livne. Integrated aeroservoelastic optimization: status and direction. Journal of Aircraft, 36(1):122–145, 1999.
  31. 31.Brendan McMahan, Eider Moore, Daniel Ramage, Seth Hampson, and Blaise Aguera y Arcas. Communication-efficient learning of deep networks from decentralized data. In Artificial intelligence and statistics, pages 1273–1282. PMLR, 2017.
  32. 32.Takayuki Okuno, Akiko Takeda, Akihiro Kawana, and Motokazu Watanabe. On lp-hyperparameter learning via bilevel nonsmooth optimization. Journal of Machine Learning Research, 22(245):1–47, 2021.
  33. 33.Shi Pu and Angelia Nedic. Distributed stochastic gradient tracking methods. Mathematical Programming, 187(1):409–457, 2021.
  34. 34.Alexander Rakhlin, Ohad Shamir, and Karthik Sridharan. Making gradient descent optimal for strongly convex stochastic optimization. arXiv preprint arXiv:1109.5647, 2011.
  35. 35.A. Ruszczynski and A. Shapiro. Optimization of convex risk functions. Mathematics of Operations Research, 31(3):433–452, 2006.
  36. 36.Andrzej Ruszczynski. A stochastic subgradient method for nonsmooth nonconvex multilevel composition optimization. SIAM Journal on Control and Optimization, 59(3):2301–2320, 2021.
  37. 37.Chenggen Shi, Jie Lu, and Guangquan Zhang. An extended kuhn–tucker approach for linear bilevel programming. Applied Mathematics and Computation, 162(1):51–63, 2005.
  38. 38.Jake Snell, Kevin Swersky, and Richard Zemel. Prototypical networks for few-shot learning. Advances in neural information processing systems, 30, 2017.
  39. 39.Jaroslaw Sobieszczanski-Sobieski and Raphael T Haftka. Multidisciplinary aerospace design optimization: survey of recent developments. Structural optimization, 14(1):1–23, 1997.
  40. 40.Murtaza Taj and Andrea Cavallaro. Distributed and decentralized multicamera tracking. IEEE Signal Processing Magazine, 28(3):46–58, 2011.
  41. 41.Davoud Ataee Tarzanagh, Mingchen Li, Christos Thrampoulidis, and Samet Oymak. Fednest: Federated bilevel, minimax, and compositional optimization. arXiv preprint arXiv:2205.02215, 2022.
  42. 42.Shenyinying Tu. Two-Stage Decomposition Algorithms and Their Application to Optimal Power Flow Problems. PhD thesis, Northwestern University, 2021.
  43. 43.Shenyinying Tu, Andreas Wächter, and Ermin Wei. A two-stage decomposition approach for ac optimal power flow. IEEE Transactions on Power Systems, 36(1):303–312, 2020.
  44. 44.Heinrich Von Stackelberg and Stackelberg Heinrich Von. The theory of the market economy. Oxford University Press, 1952.
  45. 45.M. Wang and J. Liu. A stochastic compositional gradient method using markov samples. In Proceedings of the 2016 Winter Simulation Conference, pages 702–713. IEEE Press, 2016.
  46. 46.M. Wang, E. X. Fang, and H. Liu. Stochastic compositional gradient descent: Algorithms for minimizing compositions of expected-value functions. Mathematical Programming, 161(1-2):419–449, 2017a.
  47. 47.M. Wang, J. Liu, and E. X. Fang. Accelerating stochastic composition optimization. The Journal of Machine Learning Research, 18(1):3721–3743, 2017b.
  48. 48.Mengdi Wang and Dimitri P Bertsekas. Incremental constraint projection methods for variational inequalities. Mathematical Programming, 150(2):321–363, 2015.
  49. 49.Mengdi Wang and Dimitri P Bertsekas. Stochastic first-order methods with random constraint projection. SIAM Journal on Optimization, 26(1):681–717, 2016.
  50. 50.Mengdi Wang, Ji Liu, and Ethan Fang. Accelerating stochastic composition optimization. In Advances in Neural Information Processing Systems, pages 1714–1722, 2016.
  51. 51.Xiaohan Wang, Mengdi Wang, and Yuantao Gu. A distributed tracking algorithm for reconstruction of graph signals. IEEE Journal of Selected Topics in Signal Processing, 9(4):728–740, 2015.
  52. 52.Ran Xin, Usman K Khan, and Soummya Kar. Variance-reduced decentralized stochastic optimization with accelerated convergence. IEEE Transactions on Signal Processing, 68:6255–6271, 2020.
  53. 53.Ran Xin, Usman K Khan, and Soummya Kar. An improved convergence analysis for decentralized online stochastic non-convex optimization. IEEE Transactions on Signal Processing, 69:1842–1858, 2021.
  54. 54.Yue Xu, Zengde Deng, Mengdi Wang, Wenjun Xu, Anthony Man-Cho So, and Shuguang Cui. Voting-based multiagent reinforcement learning for intelligent iot. IEEE Internet of Things Journal, 8(4):2681–2693, 2020.
  55. 55.Junjie Yang, Kaiyi Ji, and Yingbin Liang. Provably faster algorithms for bilevel optimization. Advances in Neural Information Processing Systems, 34, 2021.
  56. 56.S. Yang, M. Wang, and E. X. Fang. Multilevel stochastic gradient methods for nested composition optimization. SIAM Journal on Optimization, 29(1):616–659, 2019.

Citation

MLA
Yang, S., et al. “Decentralized Gossip-Based Stochastic Bilevel Optimization over Communication Networks”. Advances in Neural Information Processing Systems, vol. 35, 2022, pp. 238–52, https://proceedings.neurips.cc/paper_files/paper/2022/file/01db36a646c07c64dd39a92b4eceb417-Paper-Conference.pdf.
APA
Yang, S., Zhang, X., & Wang, M. (2022). Decentralized Gossip-Based Stochastic Bilevel Optimization over Communication Networks. Advances in Neural Information Processing Systems, 35, 238–252. https://proceedings.neurips.cc/paper_files/paper/2022/file/01db36a646c07c64dd39a92b4eceb417-Paper-Conference.pdf
Chicago
Yang, S., X. Zhang, and M. Wang. 2022. “Decentralized Gossip-Based Stochastic Bilevel Optimization over Communication Networks”. Advances in Neural Information Processing Systems 35: 238–52. https://proceedings.neurips.cc/paper_files/paper/2022/file/01db36a646c07c64dd39a92b4eceb417-Paper-Conference.pdf.
Harvard
Yang, S., Zhang, X. and Wang, M. (2022) “Decentralized Gossip-Based Stochastic Bilevel Optimization over Communication Networks”, Advances in Neural Information Processing Systems. Curran Associates, Inc., pp. 238–252. Available at: https://proceedings.neurips.cc/paper_files/paper/2022/file/01db36a646c07c64dd39a92b4eceb417-Paper-Conference.pdf.
Vancouver
1. Yang S, Zhang X, Wang M (2022) Decentralized Gossip-Based Stochastic Bilevel Optimization over Communication Networks. In: Advances in Neural Information Processing Systems. Curran Associates, Inc., pp 238–252

BibTeX

@inproceedings{yang2022decentralized,
  title = {Decentralized Gossip-Based Stochastic Bilevel Optimization over Communication Networks},
  author = {Yang, Shuoguang and Zhang, Xuezhou and Wang, Mengdi},
  year = {2022},
  booktitle = {Advances in Neural Information Processing Systems},
  publisher = {Curran Associates, Inc.},
  volume = {35},
  pages = {238-252},
  url = {https://proceedings.neurips.cc/paper_files/paper/2022/file/01db36a646c07c64dd39a92b4eceb417-Paper-Conference.pdf}
}
Metadata:DOI registry

Access the Paper

This paper is available from its original source. Click below to access the PDF.

Open PDF
License: Authors