This paper studies distributed convex optimization with both affine equality and nonlinear inequality couplings through the duality analysis. We first formulate the dual of the coupling-constraint problem and reformulate it as a consensus optimization problem over a connected network. To efficiently solve this dual problem and hence the primal problem, we design an accelerated linearized algorithm that, at each round, a look-ahead linearization of the separable objective is combined with a quadratic penalty on the Laplacian constraint, a proximal step, and an aggregation of iterations. On the theory side, we prove non-ergodic rates for both the primal optimality error and the feasibility error. On the other hand, numerical experiments show a faster decrease of optimality error and feasibility residual than augmented-Lagrangian tracking and distributed subgradient baselines under the same communication budget.
পেপার আইডি : 2511.19708শিরোনাম : সমতা এবং অসমতা সংযোগ সীমাবদ্ধতা সহ একটি ত্বরান্বিত বিতরণকৃত অ্যালগরিদমলেখক : চেনইয়াং কিউ, ইয়াংইয়াং কিয়ান, জোংলি লিন, ইয়াকভ এ. শামাশলেখকদের প্রতিষ্ঠান : ভার্জিনিয়া বিশ্ববিদ্যালয় (কিউ, কিয়ান, লিন), স্টোনি ব্রুক বিশ্ববিদ্যালয় (শামাশ)শ্রেণীবিভাগ : math.OC (অপ্টিমাইজেশন এবং নিয়ন্ত্রণ), cs.SY (সিস্টেম এবং নিয়ন্ত্রণ), eess.SY (সিস্টেম এবং নিয়ন্ত্রণ)জমা দেওয়ার সময় : ২০২৫ সালের নভেম্বর ২৪পেপার লিঙ্ক : https://arxiv.org/abs/2511.19708 এই পেপারটি অ্যাফাইন সমতা সীমাবদ্ধতা এবং অরৈখিক অসমতা সীমাবদ্ধতা সহ বিতরণকৃত উত্তল অপ্টিমাইজেশন সমস্যা অধ্যয়ন করে। দ্বৈত বিশ্লেষণের মাধ্যমে, সংযোগ সীমাবদ্ধতা সমস্যাটি সংযুক্ত নেটওয়ার্কে একটি সামঞ্জস্য অপ্টিমাইজেশন সমস্যায় রূপান্তরিত হয়। দ্বৈত সমস্যা দক্ষতার সাথে সমাধান করার জন্য এবং মূল সমস্যা সমাধান করার জন্য, একটি ত্বরান্বিত রৈখিকীকরণ অ্যালগরিদম ডিজাইন করা হয়েছে যা প্রতিটি পুনরাবৃত্তিতে বিভাজ্য উদ্দেশ্য ফাংশনের এগিয়ে যাওয়া রৈখিকীকরণ, ল্যাপ্লাসিয়ান সীমাবদ্ধতার দ্বিঘাত শাস্তি পদ, প্রক্সিমাল পদক্ষেপ এবং পুনরাবৃত্তিমূলক সমন্বয় একত্রিত করে। তাত্ত্বিকভাবে মূল সমস্যার সর্বোত্তমতা ত্রুটি এবং সম্ভাব্যতা ত্রুটির অ-ঐতিহ্যবাহী সংগ্রহ হার প্রমাণ করা হয়েছে। সংখ্যাগত পরীক্ষা-নিরীক্ষা দেখায় যে একই যোগাযোগ বাজেটের অধীনে, অ্যালগরিদম সর্বোত্তমতা ত্রুটি এবং সম্ভাব্যতা অবশিষ্টাংশের হ্রাসের গতিতে বিদ্যমান অত্যাধুনিক অ্যালগরিদমকে ছাড়িয়ে যায়।
বিতরণকৃত অপ্টিমাইজেশন স্থানীয় গণনা এবং যোগাযোগের মাধ্যমে বহু-এজেন্ট সিস্টেমে বৈশ্বিক উদ্দেশ্য ফাংশন কমানোর লক্ষ্য রাখে। এই পেপারটি যে সংযোগ সীমাবদ্ধতা সমস্যা (Coupling-Constraint Problem, CCP) এর উপর দৃষ্টি নিবদ্ধ করে তা বিশেষভাবে চ্যালেঞ্জিং কারণ এজেন্টদের বৈশ্বিক সংযোগ সীমাবদ্ধতা পূরণ করার সময় স্থানীয় সিদ্ধান্ত সমন্বয় করতে হবে।
এই ধরনের সমস্যা বাস্তব প্রয়োগে ব্যাপকভাবে বিদ্যমান:
স্মার্ট গ্রিড : অর্থনৈতিক প্রেরণ সমস্যায়, বৈশ্বিক অ্যাফাইন সমতা সীমাবদ্ধতা শক্তি ভারসাম্য অবস্থা প্রতিনিধিত্ব করে (মোট বিদ্যুৎ উৎপাদন মোট চাহিদা পূরণ করে)সম্পদ বরাদ্দ : ব্যক্তিগত সীমা এবং বৈশ্বিক সীমা উভয়ই পূরণ করার প্রয়োজননির্গমন সীমাবদ্ধতা : নেটওয়ার্ক ক্ষমতা সীমাবদ্ধতা এবং অন্যান্য শারীরিক সীমাবদ্ধতা সংযোগ অসমতা সীমাবদ্ধতা হিসাবে মডেল করা হয়সমতা সীমাবদ্ধতা পরিচালনা : ADMM, মিরর পদ্ধতি, গ্রেডিয়েন্ট ট্র্যাকিং ইত্যাদির মতো বিদ্যমান পদ্ধতিগুলি প্রধানত সমতা সীমাবদ্ধতার জন্যঅসমতা সীমাবদ্ধতা পরিচালনা : অ্যাফাইন অসমতা সীমাবদ্ধতার জন্য পদ্ধতিগুলি অরৈখিক সীমাবদ্ধতার জন্য প্রযোজ্য নয়সংগ্রহ হার সমস্যা : বৈশ্বিক সংযোগ অরৈখিক অসমতা সীমাবদ্ধতা পরিচালনা করে এমন বিদ্যমান অ্যালগরিদমগুলির নিম্নলিখিত সীমাবদ্ধতা রয়েছে:
অ্যাসিম্পটোটিক সংগ্রহ 13,17,18 ঐতিহ্যবাহী সংগ্রহ হার: O(ln N/√N)14 、O(1/√N)15 、O(1/N)16 ত্বরণ এবং অ-ঐতিহ্যবাহী সংগ্রহ গ্যারান্টির অভাব বেশিরভাগ বিদ্যমান বিতরণকৃত অ্যালগরিদম ত্বরান্বিত সংগ্রহ বিবেচনা করে না, সংগ্রহ হার তুলনামূলকভাবে ধীর। এই পেপারটি প্রমাণযোগ্য ত্বরান্বিত অ-ঐতিহ্যবাহী সংগ্রহ হার সহ একটি বিতরণকৃত অ্যালগরিদম বিকাশের লক্ষ্য রাখে, শাস্ত্রীয় প্রথম-ক্রম পদ্ধতির তাত্ত্বিক গ্যারান্টি সাধারণ (সম্ভবত অ-মসৃণ) খরচ ফাংশন সহ CCP ফ্রেমওয়ার্কে প্রসারিত করে।
অ্যালগরিদম উদ্ভাবন : একটি নতুন ত্বরান্বিত বিতরণকৃত অপ্টিমাইজেশন অ্যালগরিদম প্রস্তাব করা হয়েছে যা অ্যাফাইন সমতা সীমাবদ্ধতা এবং অরৈখিক অসমতা সংযোগ সীমাবদ্ধতা উভয়ই পরিচালনা করতে পারেতাত্ত্বিক অগ্রগতি : অ-ঐতিহ্যবাহী সংগ্রহ হার প্রতিষ্ঠা করা হয়েছে:মূল সমস্যার সর্বোত্তমতা ত্রুটি: O(1/N²) + O(1/N) সীমাবদ্ধতা লঙ্ঘন ত্রুটি: O(1/N²) + O(1/N) বিদ্যমান কাজের ঐতিহ্যবাহী বা অ্যাসিম্পটোটিক সংগ্রহ গ্যারান্টি উল্লেখযোগ্যভাবে উন্নত করা হয়েছে দ্বৈত পুনর্গঠন : CCP কে দ্বৈত সমস্যায় রূপান্তরিত করা হয়েছে, বিভাজ্যতা ব্যবহার করে এটিকে একটি সামঞ্জস্য অপ্টিমাইজেশন সমস্যা হিসাবে ব্যাখ্যা করা হয়েছেপরীক্ষামূলক যাচাইকরণ : সংখ্যাগত পরীক্ষা-নিরীক্ষা দেখায় যে একই পুনরাবৃত্তি বাজেটের অধীনে, অ্যালগরিদম সর্বোত্তমতা ত্রুটি এবং সম্ভাব্যতা অবশিষ্টাংশের হ্রাসের গতিতে ALT এবং বিতরণকৃত সাব-গ্রেডিয়েন্টের মতো অত্যাধুনিক অ্যালগরিদমকে ছাড়িয়ে যায়মূল সমস্যা (সমস্যা 1) :
min x ∈ X f ( x ) = ∑ i = 1 n f i ( x i ) \min_{x \in X} f(x) = \sum_{i=1}^{n} f_i(x_i) min x ∈ X f ( x ) = ∑ i = 1 n f i ( x i )
সীমাবদ্ধতা শর্ত:
সমতা সংযোগ সীমাবদ্ধতা: ∑ i = 1 n B i x i = ∑ i = 1 n b i \sum_{i=1}^{n} B_i x_i = \sum_{i=1}^{n} b_i ∑ i = 1 n B i x i = ∑ i = 1 n b i অসমতা সংযোগ সীমাবদ্ধতা: ∑ i = 1 n h i ( x i ) ≤ 0 \sum_{i=1}^{n} h_i(x_i) \leq 0 ∑ i = 1 n h i ( x i ) ≤ 0 স্থানীয় সীমাবদ্ধতা: x i ∈ X i ⊆ R p x_i \in X_i \subseteq \mathbb{R}^p x i ∈ X i ⊆ R p যেখানে:
x = [ x 1 T , x 2 T , … , x n T ] T ∈ R n p x = [x_1^T, x_2^T, \ldots, x_n^T]^T \in \mathbb{R}^{np} x = [ x 1 T , x 2 T , … , x n T ] T ∈ R n p B i ∈ R d × p B_i \in \mathbb{R}^{d \times p} B i ∈ R d × p , b i ∈ R d b_i \in \mathbb{R}^d b i ∈ R d h i : R p → R m h_i: \mathbb{R}^p \to \mathbb{R}^m h i : R p → R m সম্ভবত অরৈখিক ফাংশনমূল অনুমান :
অনুমান 1 : f i f_i f i উপযুক্ত μ f \mu_f μ f -দৃঢ়ভাবে উত্তল ফাংশন; h i h_i h i উত্তল এবং l h l_h l h -লিপশিৎজ ক্রমাগতঅনুমান 2 : X i X_i X i সংক্ষিপ্ত উত্তল সেট; স্লেটার শর্ত সন্তুষ্ট (কঠোর সম্ভাব্য বিন্দু বিদ্যমান)ল্যাগ্রেঞ্জ গুণক μ ∈ R d \mu \in \mathbb{R}^d μ ∈ R d (সমতা সীমাবদ্ধতা) এবং δ ∈ R + m \delta \in \mathbb{R}_+^m δ ∈ R + m (অসমতা সীমাবদ্ধতা) প্রবর্তন করা হয়েছে, ল্যাগ্রেঞ্জ ফাংশন হল:
L ( x , μ , δ ) = ∑ i = 1 n ( F i ( x i ) + ⟨ μ , B i x i − b i ⟩ + ⟨ δ , h i ( x i ) ⟩ ) L(x, \mu, \delta) = \sum_{i=1}^{n} \left( F_i(x_i) + \langle \mu, B_i x_i - b_i \rangle + \langle \delta, h_i(x_i) \rangle \right) L ( x , μ , δ ) = ∑ i = 1 n ( F i ( x i ) + ⟨ μ , B i x i − b i ⟩ + ⟨ δ , h i ( x i )⟩ )
যেখানে F i = f i + 1 X i F_i = f_i + \mathbb{1}_{X_i} F i = f i + 1 X i (1 X i \mathbb{1}_{X_i} 1 X i সূচক ফাংশন)।
দ্বৈত সমস্যা :
min μ ∈ R d , δ ∈ R + m ∑ i = 1 n g i ( μ , δ ) \min_{\mu \in \mathbb{R}^d, \delta \in \mathbb{R}_+^m} \sum_{i=1}^{n} g_i(\mu, \delta) min μ ∈ R d , δ ∈ R + m ∑ i = 1 n g i ( μ , δ )
যেখানে g i ( μ , δ ) = − min x i L i ( x i , μ , δ ) g_i(\mu, \delta) = -\min_{x_i} L_i(x_i, \mu, \delta) g i ( μ , δ ) = − min x i L i ( x i , μ , δ ) ।
প্রতিটি এজেন্ট i i i দ্বৈত পরিবর্তনশীল অনুলিপি y i = [ μ i T , δ i T ] T ∈ Y = R d × R + m y_i = [\mu_i^T, \delta_i^T]^T \in Y = \mathbb{R}^d \times \mathbb{R}_+^m y i = [ μ i T , δ i T ] T ∈ Y = R d × R + m বজায় রাখে, দ্বৈত সমস্যাটি পুনর্গঠন করা হয়েছে:
min y ∈ Y G ( y ) = ∑ i = 1 n g i ( y i ) \min_{y \in \mathcal{Y}} G(y) = \sum_{i=1}^{n} g_i(y_i) min y ∈ Y G ( y ) = ∑ i = 1 n g i ( y i ) s.t. y 1 = y 2 = ⋯ = y n \text{s.t. } y_1 = y_2 = \cdots = y_n s.t. y 1 = y 2 = ⋯ = y n
ল্যাপ্লাসিয়ান ম্যাট্রিক্স H H H এবং W = H ⊗ I d + m W = H \otimes I_{d+m} W = H ⊗ I d + m ব্যবহার করে, সামঞ্জস্য সীমাবদ্ধতা W 1 / 2 y = 0 W^{1/2}y = 0 W 1/2 y = 0 এর সমতুল্য, সংক্ষিপ্ত ফর্ম (সমস্যা 4) পাওয়া যায়:
min y ∈ Y G ( y ) s.t. W 1 / 2 y = 0 \min_{y \in \mathcal{Y}} G(y) \quad \text{s.t. } W^{1/2}y = 0 min y ∈ Y G ( y ) s.t. W 1/2 y = 0
বর্ধিত ল্যাগ্রেঞ্জ ফাংশন :
L ρ ( y , v ) = G ( y ) − ⟨ v , W 1 / 2 y ⟩ + ρ 2 ∥ W 1 / 2 y ∥ 2 \mathcal{L}_\rho(y, v) = G(y) - \langle v, W^{1/2}y \rangle + \frac{\rho}{2} \|W^{1/2}y\|^2 L ρ ( y , v ) = G ( y ) − ⟨ v , W 1/2 y ⟩ + 2 ρ ∥ W 1/2 y ∥ 2
অ্যালগরিদম পুনরাবৃত্তি (অ্যালগরিদম 1) :
প্রাথমিকীকরণ: ŷ_{i,1} = y_{i,1} ∈ Y, λ_{i,1} = 0
k = 1, 2, ..., N এর জন্য:
1. বহিরাগমন পদক্ষেপ:
ỹ_{i,k} = (1 - α_k)ŷ_{i,k} + α_k y_{i,k}
2. স্থানীয় অপ্টিমাইজেশন (গ্রেডিয়েন্ট গণনা):
x_{i,k} = argmin_x {F_i(x) + ⟨[B_i x - b_i; h_i(x)], ỹ_{i,k}⟩}
∇g_i(ỹ_{i,k}) = -[B_i x_{i,k} - b_i; h_i(x_{i,k})]
3. তথ্য বিনিময়:
t_{i,k} = Σ_{j∈N_i} H_{ij}(y_{i,k} - y_{j,k})
4. প্রক্সিমাল আপডেট:
y_{i,k+1} = P_Y{y_{i,k} - 1/η_k(∇g_i(ỹ_{i,k}) - λ_{i,k} - θ_k t_{i,k})}
5. সমন্বয় পদক্ষেপ:
ŷ_{i,k+1} = (1 - α_k)ŷ_{i,k} + α_k y_{i,k+1}
6. দ্বৈত পরিবর্তনশীল আপডেট:
λ_{i,k+1} = λ_{i,k} - β_k t_{i,k}
পরামিতি সেটিং :
α k = 2 k + 1 \alpha_k = \frac{2}{k+1} α k = k + 1 2 (নেস্টেরভ ত্বরণ পরামিতি)θ k = ρ N k \theta_k = \frac{\rho N}{k} θ k = k ρN (স্ব-অভিযোজিত ল্যাপ্লাসিয়ান শাস্তি)β k = ρ k N \beta_k = \frac{\rho k}{N} β k = N ρ k (দ্বৈত পদক্ষেপ দৈর্ঘ্য)η k = 2 l g + ρ N ∥ W ∥ k \eta_k = \frac{2l_g + \rho N \|W\|}{k} η k = k 2 l g + ρN ∥ W ∥ (প্রক্সিমাল পরামিতি)যেখানে l g = 2 μ f 2 ( ∥ B i ∥ 2 + l h 2 ) ⋅ max { ∥ B i ∥ 2 , l h 2 } l_g = \sqrt{\frac{2}{\mu_f^2}(\|B_i\|^2 + l_h^2)} \cdot \max\{\|B_i\|^2, l_h^2\} l g = μ f 2 2 ( ∥ B i ∥ 2 + l h 2 ) ⋅ max { ∥ B i ∥ 2 , l h 2 } হল g i g_i g i এর লিপশিৎজ ধ্রুবক।
তিন-পরিবর্তনশীল সহযোগিতা প্রক্রিয়া :y ~ k \tilde{y}_k y ~ k : বহিরাগমন পূর্বাভাস বিন্দু, গ্রেডিয়েন্ট মূল্যায়নের জন্য ব্যবহৃত, গতিশীলতা প্রভাব প্রবর্তন করেy k y_k y k : প্রক্সিমাল সংশোধন বিন্দু, স্থিতিশীলতা নিশ্চিত করেy ^ k \hat{y}_k y ^ k : ট্র্যাজেক্টরি মসৃণকরণ বিন্দু, সর্বোত্তম সংগ্রহ বিশ্লেষণ বাস্তবায়ন করেস্ব-অভিযোজিত পরামিতি সময়সূচী :θ k \theta_k θ k এবং β k \beta_k β k পুনরাবৃত্তি সংখ্যার সাথে স্ব-অভিযোজিতভাবে সামঞ্জস্য করা হয়, সংগ্রহ গতি এবং স্থিতিশীলতার ভারসাম্য রাখেপরামিতি ডিজাইন অ-ঐতিহ্যবাহী O(1/N²) ত্বরান্বিত হার নিশ্চিত করে রৈখিকীকরণ কৌশল :অ-বিভাজ্য দ্বিঘাত পদ ρ 2 ∥ W 1 / 2 y ∥ 2 \frac{\rho}{2}\|W^{1/2}y\|^2 2 ρ ∥ W 1/2 y ∥ 2 রৈখিকীকৃত করা হয়েছে বর্তমান বিন্দু গ্রেডিয়েন্টের পরিবর্তে এগিয়ে যাওয়া গ্রেডিয়েন্ট ∇ G ( y ~ k ) \nabla G(\tilde{y}_k) ∇ G ( y ~ k ) সহ একত্রিত করা হয়েছে বিতরণকৃত বাস্তবায়ন :প্রতিটি নোড শুধুমাত্র স্থানীয় সাব-সমস্যা সমাধান করতে হবে (সমীকরণ 14) শুধুমাত্র একটি রাউন্ড প্রতিবেশী তথ্য বিনিময় প্রয়োজন (ধাপ 6 এ t i , k t_{i,k} t i , k ) কোন বৈশ্বিক সমন্বয়কারী প্রয়োজন নেই সিন্থেটিক অপ্টিমাইজেশন সমস্যা :
min x i ∈ X i ∑ i = 1 n ( x i T A i x i + b i T x i + ∥ x i ∥ 1 ) \min_{x_i \in X_i} \sum_{i=1}^{n} \left( x_i^T A_i x_i + b_i^T x_i + \|x_i\|_1 \right) min x i ∈ X i ∑ i = 1 n ( x i T A i x i + b i T x i + ∥ x i ∥ 1 )
সীমাবদ্ধতা শর্ত:
সমতা: ∑ i = 1 n C i x i = 0 p \sum_{i=1}^{n} C_i x_i = 0_p ∑ i = 1 n C i x i = 0 p অসমতা: ∑ i = 1 n ∥ x i − r i ∥ 1 ≤ ∑ i = 1 n d i \sum_{i=1}^{n} \|x_i - r_i\|_1 \leq \sum_{i=1}^{n} d_i ∑ i = 1 n ∥ x i − r i ∥ 1 ≤ ∑ i = 1 n d i পরামিতি সেটিং :
এজেন্ট সংখ্যা: n = 20 n = 20 n = 20 স্থানীয় মাত্রা: p = 5 p = 5 p = 5 বক্স সীমাবদ্ধতা: x i ∈ X i = { x ∈ R p ∣ x ‾ i ≤ x ≤ x ˉ i } x_i \in X_i = \{x \in \mathbb{R}^p | \underline{x}_i \leq x \leq \bar{x}_i\} x i ∈ X i = { x ∈ R p ∣ x i ≤ x ≤ x ˉ i } x ‾ i ∼ U [ − 10 , − 9 ] \underline{x}_i \sim U[-10, -9] x i ∼ U [ − 10 , − 9 ] , x ˉ i ∼ U [ 9 , 10 ] \bar{x}_i \sim U[9, 10] x ˉ i ∼ U [ 9 , 10 ] ম্যাট্রিক্স A i = U i Λ i U i T A_i = U_i \Lambda_i U_i^T A i = U i Λ i U i T :
U i U_i U i র্যান্ডম অর্থোগোনাল ম্যাট্রিক্সΛ i \Lambda_i Λ i এর আইজেনভ্যালু [ 1 , 100 ] [1, 100] [ 1 , 100 ] এ রৈখিকভাবে বিতরণ করা হয় (শর্ত সংখ্যা κ = 100 \kappa = 100 κ = 100 ) C i , b i ∼ N ( 0 , I p ) C_i, b_i \sim \mathcal{N}(0, I_p) C i , b i ∼ N ( 0 , I p ) d i ∼ U ( 1 , 6 ) d_i \sim U(1, 6) d i ∼ U ( 1 , 6 ) যোগাযোগ নেটওয়ার্ক :
সংযুক্ত অনির্দেশিত গ্রাফ: প্রতিটি নোড নিকটতম এবং দ্বিতীয় নিকটতম প্রতিবেশীর সাথে সংযুক্ত প্রান্ত সেট: ( i , i + 1 ) (i, i+1) ( i , i + 1 ) 1 ≤ i ≤ 19 1 \leq i \leq 19 1 ≤ i ≤ 19 এর জন্য, যোগ করুন ( 1 , 20 ) (1, 20) ( 1 , 20 ) মূল সমস্যার সর্বোত্তমতা ত্রুটি :
∣ f ( x k ) − f ( x ∗ ) ∣ 2 ∣ f ( x 1 ) − f ( x ∗ ) ∣ 2 \frac{|f(x_k) - f(x^*)|^2}{|f(x_1) - f(x^*)|^2} ∣ f ( x 1 ) − f ( x ∗ ) ∣ 2 ∣ f ( x k ) − f ( x ∗ ) ∣ 2 সীমাবদ্ধতা লঙ্ঘন পরম ত্রুটি :
∥ ∑ i = 1 n C i x i , k ∥ + [ ∑ i = 1 n ( ∥ x i , k − r i ∥ 1 − d i ) ] + \left\| \sum_{i=1}^{n} C_i x_{i,k} \right\| + \left[ \sum_{i=1}^{n} (\|x_{i,k} - r_i\|_1 - d_i) \right]_+ ∥ ∑ i = 1 n C i x i , k ∥ + [ ∑ i = 1 n ( ∥ x i , k − r i ∥ 1 − d i ) ] + বিতরণকৃত সাব-গ্রেডিয়েন্ট 14 : বিতরণকৃত সাব-গ্রেডিয়েন্ট অ্যালগরিদমALT (বর্ধিত ল্যাগ্রেঞ্জ ট্র্যাকিং) 17 : বর্ধিত ল্যাগ্রেঞ্জ ট্র্যাকিং অ্যালগরিদমIPLUX (একীভূত প্রাইমাল-ডুয়াল প্রক্সিমাল) 16 : একীভূত প্রাইমাল-ডুয়াল প্রক্সিমাল অ্যালগরিদমবেঞ্চমার্ক সমাধান : YALMIP সহ MOSEK সমাধক ব্যবহার করে সর্বোত্তম সমাধান x ∗ x^* x ∗ পাওয়া যায়
সমস্ত অ্যালগরিদম একই প্রাথমিকীকরণ ব্যবহার করে পুনরাবৃত্তি সংখ্যা: N = 1200 N = 1200 N = 1200 এই পেপারের অ্যালগরিদম পরামিতি উপপাদ্য 1 অনুযায়ী সেট করা হয়েছে চিত্র 1: মূল সমস্যার সর্বোত্তমতা ত্রুটি
এই পেপারের অ্যালগরিদম : k = 1200 k=1200 k = 1200 এ 10 − 6 10^{-6} 1 0 − 6 নির্ভুলতা অর্জন করেALT : একঘেয়ে হ্রাস কিন্তু ধীর, শেষে প্রায় 10 − 2 10^{-2} 1 0 − 2 বিতরণকৃত সাব-গ্রেডিয়েন্ট : সবচেয়ে ধীর হ্রাস, 10 − 1 10^{-1} 1 0 − 1 -10 0 10^0 1 0 0 পরিসরে থাকেIPLUX : ALT এবং এই পেপারের অ্যালগরিদমের মধ্যে কর্মক্ষমতাচিত্র 2: সীমাবদ্ধতা লঙ্ঘন পরম ত্রুটি
এই পেপারের অ্যালগরিদম : সবচেয়ে তাড়াতাড়ি 10 − 4 10^{-4} 1 0 − 4 এর নিচে পৌঁছায়অন্যান্য অ্যালগরিদম : স্পষ্টভাবে ধীর সংগ্রহসংগ্রহ গতি : এই পেপারের অ্যালগরিদম একই পুনরাবৃত্তি সংখ্যায় সমস্ত তুলনা পদ্ধতির চেয়ে উল্লেখযোগ্যভাবে দ্রুত সংগ্রহ করেনির্ভুলতা সুবিধা :সর্বোত্তমতা ত্রুটি প্রায় 4 অর্ডার মাত্রা হ্রাস (10 − 2 10^{-2} 1 0 − 2 থেকে 10 − 6 10^{-6} 1 0 − 6 ) সম্ভাব্যতা ত্রুটি প্রায় 2 অর্ডার মাত্রা হ্রাস ত্বরণ প্রভাব স্পষ্ট : অ-ঐতিহ্যবাহী সংগ্রহ হারের তাত্ত্বিক সুবিধা পরীক্ষায় যাচাই করা হয়েছেশক্তিশালীতা : অ্যালগরিদম অ-মসৃণ উদ্দেশ্য ফাংশন (ধারণ করে ℓ 1 \ell_1 ℓ 1 নর্ম) এবং অরৈখিক সীমাবদ্ধতার অধীনে স্থিতিশীল কর্মক্ষমতা প্রদর্শন করেADMM পদ্ধতি 6,7 : বিকল্প দিকনির্দেশনা গুণক পদ্ধতিমিরর পদ্ধতি 8 : মিরর ডিসেন্টের উপর ভিত্তি করে বিতরণকৃত অ্যালগরিদমগ্রেডিয়েন্ট ট্র্যাকিং 9 : দ্বৈত সমস্যার গ্রেডিয়েন্ট ট্র্যাকিংঅ্যাফাইন অসমতা 10-12 : বিতরণকৃত প্রক্সিমাল অ্যালগরিদম, সমন্বয় অপ্টিমাইজেশনঅরৈখিক অসমতা 13-18 :
দ্বৈত সাব-গ্রেডিয়েন্ট পদ্ধতি 13 অপারেটর বিভাজন প্রাইমাল-ডুয়াল ফ্রেমওয়ার্ক 14 গতিশীল গড় সামঞ্জস্য 15 বিরল/ঘন সীমাবদ্ধতা পরিচালনা 16 ALT অ্যালগরিদম 17 নেস্টেরভ ত্বরণ 19 : সীমাবদ্ধতা-মুক্ত উত্তল অপ্টিমাইজেশনে O(1/N²) হারFISTA 20 : দ্রুত পুনরাবৃত্তিমূলক সংকোচন থ্রেশহোল্ড অ্যালগরিদমদ্রুত ল্যাগ্রেঞ্জ পদ্ধতি 21,22 : উত্তল অপ্টিমাইজেশনে ত্বরান্বিত ল্যাগ্রেঞ্জ পদ্ধতিবিতরণকৃত ত্বরণ 23,24 : DCatalyst, শক্তি সংরক্ষণ নীতিপ্রথমবার নেস্টেরভ ত্বরণ সমতা এবং অরৈখিক অসমতা সংযোগ সীমাবদ্ধতা উভয়ই সহ বিতরণকৃত CCP তে প্রসারিত করা হয়েছেঅ-ঐতিহ্যবাহী সংগ্রহ গ্যারান্টি প্রদান করে (গড়ের উপর নির্ভর করে না), বিদ্যমান ঐতিহ্যবাহী বা অ্যাসিম্পটোটিক ফলাফল উন্নত করেঅ-মসৃণ উদ্দেশ্য ফাংশনের জন্য প্রযোজ্যদ্বৈত ফাংশনের লিপশিৎজ মসৃণতা :
∥ ∇ g i ( z 1 ) − ∇ g i ( z 2 ) ∥ ≤ l g ∥ z 1 − z 2 ∥ \|\nabla g_i(z_1) - \nabla g_i(z_2)\| \leq l_g \|z_1 - z_2\| ∥∇ g i ( z 1 ) − ∇ g i ( z 2 ) ∥ ≤ l g ∥ z 1 − z 2 ∥
যেখানে l g = 2 μ f 2 ( ∥ B i ∥ 2 + l h 2 ) ⋅ max { ∥ B i ∥ 2 , l h 2 } l_g = \sqrt{\frac{2}{\mu_f^2}(\|B_i\|^2 + l_h^2)} \cdot \max\{\|B_i\|^2, l_h^2\} l g = μ f 2 2 ( ∥ B i ∥ 2 + l h 2 ) ⋅ max { ∥ B i ∥ 2 , l h 2 }
প্রমাণ চিন্তাধারা :
F i F_i F i এর দৃঢ় উত্তলতা এবং h i h_i h i এর উত্তলতা ব্যবহার করা হয়েছেড্যানস্কিন উপপাদ্য দ্বারা গ্রেডিয়েন্ট অভিব্যক্তি পাওয়া যায় দৃঢ় উত্তলতা এবং লিপশিৎজ ক্রমাগততা একত্রিত করে অসমতা প্রতিষ্ঠা করা হয়েছে সংগ্রহ হার :
সম্ভাব্যতা ত্রুটি :
∥ ∑ i = 1 n B i x i , N + 1 − b i ∥ + ∥ [ ∑ i = 1 n h i ( x i , N + 1 ) ] + ∥ ≤ ε c \left\| \sum_{i=1}^{n} B_i x_{i,N+1} - b_i \right\| + \left\| \left[ \sum_{i=1}^{n} h_i(x_{i,N+1}) \right]_+ \right\| \leq \varepsilon_c ∥ ∑ i = 1 n B i x i , N + 1 − b i ∥ + [ ∑ i = 1 n h i ( x i , N + 1 ) ] + ≤ ε c যেখানে:
ε c = ( 2 l g N ( N + 1 ) + ρ N + 1 ∥ W ∥ ) ∥ y 1 − y ∗ ∥ 2 + 1 ρ ( N + 1 ) λ 2 ( W ) \varepsilon_c = \left( \frac{2l_g}{N(N+1)} + \frac{\rho}{N+1}\|W\| \right) \|y_1 - y^*\|^2 + \frac{1}{\rho(N+1)\lambda_2(W)} ε c = ( N ( N + 1 ) 2 l g + N + 1 ρ ∥ W ∥ ) ∥ y 1 − y ∗ ∥ 2 + ρ ( N + 1 ) λ 2 ( W ) 1
সর্বোত্তমতা ত্রুটি :
− ε p ≤ f ( x N + 1 ) − f ( x ∗ ) ≤ ε ˉ p -\varepsilon_p \leq f(x_{N+1}) - f(x^*) \leq \bar{\varepsilon}_p − ε p ≤ f ( x N + 1 ) − f ( x ∗ ) ≤ ε ˉ p যেখানে ε p \varepsilon_p ε p এবং ε ˉ p \bar{\varepsilon}_p ε ˉ p একই রকম O(1/N²) + O(1/N) ফর্ম রয়েছে।
প্রমাণ মূল পদক্ষেপ :
শক্তি ফাংশন নির্মাণ :
Φ k = G ( y ^ k ) − G ( y ∗ ) − ⟨ λ , y ^ k − y ∗ ⟩ \Phi_k = G(\hat{y}_k) - G(y^*) - \langle \lambda, \hat{y}_k - y^* \rangle Φ k = G ( y ^ k ) − G ( y ∗ ) − ⟨ λ , y ^ k − y ∗ ⟩ পুনরাবৃত্তিমূলক অসমতা : উত্তলতা এবং মসৃণতা ব্যবহার করে:
k ( k + 1 ) Φ k + 1 − k ( k − 1 ) Φ k ≤ 2 k [ telescoping terms ] k(k+1)\Phi_{k+1} - k(k-1)\Phi_k \leq 2k[\text{telescoping terms}] k ( k + 1 ) Φ k + 1 − k ( k − 1 ) Φ k ≤ 2 k [ telescoping terms ] সমষ্টি কৌশল : k = 1 k=1 k = 1 থেকে N N N পর্যন্ত সমষ্টি, ইনস্ট্রুমেন্টাল সম্পত্তি ব্যবহার করা হয়েছেপরামিতি নির্বাচন : সাবধানে ডিজাইন করা α k , θ k , β k , η k \alpha_k, \theta_k, \beta_k, \eta_k α k , θ k , β k , η k দ্বারা ত্বরণ বাস্তবায়ন করা হয়েছেঅ্যালগরিদম অবদান : সমতা এবং অরৈখিক অসমতা সংযোগ সীমাবদ্ধতা উভয়ই সহ প্রথম ত্বরান্বিত বিতরণকৃত অ্যালগরিদম প্রস্তাব করা হয়েছেতাত্ত্বিক গ্যারান্টি : অ-ঐতিহ্যবাহী O(1/N²) + O(1/N) সংগ্রহ হার প্রতিষ্ঠা করা হয়েছে, বিদ্যমান ফলাফল উল্লেখযোগ্যভাবে উন্নত করেব্যবহারিকতা : প্রতিটি পুনরাবৃত্তি গণনা সহজ (একটি স্থানীয় সাব-সমস্যা + একটি রাউন্ড প্রতিবেশী যোগাযোগ), বড় আকারের স্থাপনার জন্য উপযুক্তপরীক্ষামূলক যাচাইকরণ : প্রতিনিধিত্বমূলক পরীক্ষা সেটে, অ্যালগরিদম একই পুনরাবৃত্তি বাজেটে উচ্চতর সম্ভাব্যতা এবং কম ত্রুটি অর্জন করেদৃঢ় উত্তলতা অনুমান : অ্যালগরিদম এবং তাত্ত্বিক বিশ্লেষণ উদ্দেশ্য ফাংশনের দৃঢ় উত্তলতার উপর নির্ভর করে (অনুমান 1), প্রযোজ্য পরিসীমা সীমিত করেস্লেটার শর্ত : কঠোর সম্ভাব্য বিন্দু বিদ্যমান প্রয়োজন (অনুমান 2), কিছু বাস্তব সমস্যায় পূরণ না হতে পারেসংক্ষিপ্ত সেট অনুমান : অনুমান 2 স্থানীয় সীমাবদ্ধতা সেট X i X_i X i সংক্ষিপ্ত প্রয়োজন, সীমাহীন সীমাবদ্ধতা বাদ দেয়পরামিতি সমন্বয় : যদিও তাত্ত্বিক পরামিতি সেটিং প্রদান করা হয়েছে, বাস্তব প্রয়োগে নির্দিষ্ট সমস্যার জন্য সূক্ষ্ম সমন্বয় প্রয়োজন হতে পারেযোগাযোগ জটিলতা : যোগাযোগ জটিলতা স্পষ্টভাবে বিশ্লেষণ করা হয়নি, শুধুমাত্র পুনরাবৃত্তি জটিলতায় ফোকাস করা হয়েছেঅ-উত্তল সম্প্রসারণ : তাত্ত্বিক এবং অ্যালগরিদম ফ্রেমওয়ার্ক অ-উত্তল অপ্টিমাইজেশন সমস্যা অন্তর্ভুক্ত করে নাদৃঢ় উত্তলতা অনুমান শিথিল করা : সাধারণ উত্তল এমনকি অ-উত্তল সমস্যায় সম্প্রসারণর্যান্ডম/অনলাইন সংস্করণ : বড় আকারের ডেটা পরিচালনার জন্য র্যান্ডম গ্রেডিয়েন্ট সংস্করণ বিকাশঅ্যাসিঙ্ক্রোনাস যোগাযোগ : অ্যাসিঙ্ক্রোনাস যোগাযোগ প্রোটোকলের অধীনে সংগ্রহ গবেষণাসময়-পরিবর্তনশীল নেটওয়ার্ক : গতিশীল পরিবর্তনশীল যোগাযোগ টপোলজিতে সম্প্রসারণবাস্তব প্রয়োগ : স্মার্ট গ্রিড, অনুসরণকারী ড্রোন সমন্বয় ইত্যাদি বাস্তব সিস্টেমে যাচাইকরণতাত্ত্বিক উদ্ভাবনী শক্তি :সমতা এবং অরৈখিক অসমতা সংযোগ সীমাবদ্ধতা সহ বিতরণকৃত অপ্টিমাইজেশনে প্রথমবার O(1/N²) ত্বরণ অর্জন করা হয়েছে অ-ঐতিহ্যবাহী সংগ্রহ গ্যারান্টি বিদ্যমান ঐতিহ্যবাহী বা অ্যাসিম্পটোটিক ফলাফলের চেয়ে উন্নত কঠোর গাণিতিক প্রমাণ, যুক্তি স্পষ্ট অ্যালগরিদম ডিজাইন সূক্ষ্ম :তিন-পরিবর্তনশীল সহযোগিতা প্রক্রিয়া (y ~ k , y k , y ^ k \tilde{y}_k, y_k, \hat{y}_k y ~ k , y k , y ^ k ) কার্যকরভাবে ত্বরণ বাস্তবায়ন করে স্ব-অভিযোজিত পরামিতি সময়সূচী সংগ্রহ গতি এবং স্থিতিশীলতার ভারসাম্য রাখে রৈখিকীকরণ কৌশল গণনা বিভাজ্যতা বজায় রাখে পরীক্ষা পর্যাপ্ত :তিনটি অত্যাধুনিক অ্যালগরিদমের সাথে তুলনা পরীক্ষার ফলাফল স্পষ্টভাবে ত্বরণ প্রভাব প্রদর্শন করে চার্ট গুণমান উচ্চ, উপসংহার স্পষ্ট ব্যবহারিক মূল্য উচ্চ :অ্যালগরিদম সম্পূর্ণ বিতরণকৃত, বড় আকারের স্থাপনার জন্য উপযুক্ত প্রতিটি পুনরাবৃত্তি গণনা পরিমাণ যুক্তিসঙ্গত অ-মসৃণ উদ্দেশ্য ফাংশনের জন্য প্রযোজ্য লেখা স্পষ্ট :কাঠামো যুক্তিসঙ্গত, যুক্তি কঠোর প্রতীক সংজ্ঞা স্পষ্ট প্রমাণ বিস্তারিত এবং বোঝা সহজ অনুমান শক্তিশালী :দৃঢ় উত্তলতা অনুমান প্রযোজ্য পরিসীমা সীমিত করে (অনেক বাস্তব সমস্যা শুধুমাত্র উত্তল বা অ-উত্তল) সংক্ষিপ্ত সেট এবং স্লেটার শর্ত কিছু প্রয়োগে যাচাই করা কঠিন পরীক্ষা সীমাবদ্ধতা :শুধুমাত্র সিন্থেটিক ডেটায় পরীক্ষা, বাস্তব প্রয়োগ পরিস্থিতি যাচাইকরণের অভাব বড় আকারের নেটওয়ার্ক পরীক্ষা করা হয়নি (n=20 ছোট) যোগাযোগ ওভারহেড এবং গণনা সময় বিশ্লেষণ করা হয়নি পরামিতি নির্ভরতা :অ্যালগরিদম কর্মক্ষমতা সমস্যা পরামিতির উপর নির্ভর করে (μ f , l h , ∥ B i ∥ \mu_f, l_h, \|B_i\| μ f , l h , ∥ B i ∥ ইত্যাদি) বাস্তব প্রয়োগে এই পরামিতিগুলি অজানা বা অনুমান করা কঠিন হতে পারে সংগ্রহ ধ্রুবক :তাত্ত্বিক সংগ্রহ হারে ধ্রুবক পদ বড় হতে পারে সংগ্রহ হারের নিম্ন সীমা বা সর্বোত্তমতা বিশ্লেষণ প্রদান করা হয়নি বিশ্লেষণ অনুপস্থিত :প্রাথমিকীকরণের প্রতি অ্যালগরিদমের সংবেদনশীলতা আলোচনা করা হয়নি পরামিতি নির্বাচনের সংগ্রহে প্রভাব বিশ্লেষণ করা হয়নি ব্যর্থতার ক্ষেত্রে বা কঠিন পরিস্থিতির আলোচনা অনুপস্থিত একাডেমিক মূল্য :বিতরণকৃত সীমাবদ্ধতা অপ্টিমাইজেশনের জন্য নতুন তাত্ত্বিক সরঞ্জাম প্রদান করে ত্বরণ কৌশল অন্যান্য বিতরণকৃত অ্যালগরিদম ডিজাইনকে অনুপ্রাণিত করতে পারে অপ্টিমাইজেশন নিয়ন্ত্রণ ক্ষেত্রে উচ্চ উদ্ধৃতি প্রত্যাশিত ব্যবহারিক মূল্য :স্মার্ট গ্রিড অর্থনৈতিক প্রেরণে সরাসরি প্রযোজ্য মাল্টি-রোবট সমন্বয়, সেন্সর নেটওয়ার্ক ইত্যাদিতে সম্প্রসারণযোগ্য অ্যালগরিদম 1 স্পষ্ট বাস্তবায়ন নির্দেশিকা প্রদান করে পুনরুৎপাদনযোগ্যতা :অ্যালগরিদম বর্ণনা বিস্তারিত, বাস্তবায়ন সহজ পরীক্ষা সেটআপ স্পষ্ট লেখকদের কোড ওপেন সোর্স করার সুপারিশ করা হয় প্রয়োগ প্রচার করতে দৃঢ় সুপারিশ পরিস্থিতি :
স্মার্ট গ্রিড অর্থনৈতিক প্রেরণ (দৃঢ় উত্তলতা এবং সংক্ষিপ্ত সেট অনুমান পূরণ করে) সম্পদ বরাদ্দ সমস্যা (উত্তল খরচ ফাংশন) বিতরণকৃত মেশিন লার্নিং (দৃঢ় উত্তল নিয়মিতকরণ) সাবধানে ব্যবহার পরিস্থিতি :
অ-উত্তল অপ্টিমাইজেশন সমস্যা (তত্ত্ব প্রযোজ্য নয়) সীমাহীন সীমাবদ্ধতা সেট (সংক্ষিপ্ত সেট অনুমান লঙ্ঘন) রিয়েল-টাইম সিস্টেম (পুনরাবৃত্তি সংখ্যা অনেক হতে পারে) উন্নতির প্রয়োজন পরিস্থিতি :
বড় আকারের নেটওয়ার্ক (স্কেলেবিলিটি যাচাই প্রয়োজন) সময়-পরিবর্তনশীল পরিবেশ (অ্যালগরিদম সম্প্রসারণ প্রয়োজন) যোগাযোগ সীমিত (যোগাযোগ দক্ষতা বিবেচনা প্রয়োজন) 6 T.-H. Chang et al., "Multi-agent distributed optimization via inexact consensus ADMM," IEEE Trans. Signal Process., 2014.
14 S. Liang and G. Yin, "Distributed dual subgradient algorithms with iterate-averaging feedback," IEEE Trans. Cybernetics, 2019.
16 X. Wu et al., "Distributed optimization with coupling constraints," IEEE Trans. Automatic Control, 2022.
17 A. Falsone and M. Prandini, "Augmented Lagrangian tracking for distributed optimization," Automatica, 2023.
19 Y. Nesterov, "A method for unconstrained convex minimization problem with the rate of convergence O(1/k²)," Dokl. Akad. Nauk. SSSR, 1983.
সামগ্রিক মূল্যায়ন : এটি একটি উচ্চ মানের তাত্ত্বিক পেপার যা বিতরণকৃত অপ্টিমাইজেশন ক্ষেত্রে গুরুত্বপূর্ণ অবদান রাখে। অ্যালগরিদম ডিজাইন সূক্ষ্ম, তাত্ত্বিক বিশ্লেষণ কঠোর, পরীক্ষার ফলাফল প্রভাবশালী। যদিও কিছু অনুমান সীমাবদ্ধতা রয়েছে, তবে এর প্রযোজ্য পরিসরে উল্লেখযোগ্য সুবিধা রয়েছে। বাস্তব সিস্টেমে আরও যাচাইকরণ এবং দৃঢ় উত্তলতা অনুমান শিথিল করার সম্ভাবনা অন্বেষণের সুপারিশ করা হয়।