Distributed Non-Convex First-Order Optimization and Information Processing: Lower Complexity Bounds and Rate Optimal Algorithms

# Distributed Non-Convex First-Order Optimization and Information Processing: Lower Complexity Bounds and Rate Optimal Algorithms

Haoran Sun and Mingyi Hong H. Sun and M. Hong are with the Department of Electrical and Computer Engineering (ECE), University of Minnesota, Minneapolis, MN 55414, USA. Email: {sun00111,mhong}@umn.edu
###### Abstract

We consider a class of distributed non-convex optimization problems often arises in modern distributed signal and information processing, in which a number of agents connected by a network collectively optimize a sum of smooth (possibly non-convex) local objective functions. We address the following fundamental question: For a class of unconstrained non-convex problems with Lipschitz continuous gradient, by only utilizing local gradient information, what is the fastest rate that distributed algorithms can achieve, and how to achieve those rates.

We develop a lower bound analysis that identifies difficult problem instances for any first-order method. We show that in the worst-case it takes any first-order algorithm iterations to achieve certain -solution, where is the network diameter, and is the Lipschitz constant of the gradient. Further for a general problem class and a number of network classes, we propose optimal primal-dual gradient methods whose rates precisely match the lower bounds (up to a ploylog factor). To the best of our knowledge, this is the first time that lower rate bounds and optimal methods have been developed for distributed non-convex problems. Our results provide guidelines for future design of distributed optimization algorithms, convex and non-convex alike.

Keywords. Non-convex distributed optimization; Optimal methods; Lower Complexity Bounds; Primal-dual methods.

## 1 Introduction

### 1.1 Problem and motivation

In this work, we consider the following distributed optimization problem over a network

 miny∈RS¯f(y):=1MM∑i=1fi(y), (1)

where is a smooth and possibly non-convex function accessible to agent . There is no central controller, and the agents are connected by a network defined by an undirected and unweighted graph , with vertices and edges. Each agent can only communicate with its immediate neighbors, and it can access one component function .

A common way to reformulate problem (1) in the distributed setting is given below. Define the degree of node as ; Define the incidence matrix (IM) as follows: if and it connects vertex and with , then if , if and otherwise. Introduce local variables , and suppose the graph is connected, then the following formulation is equivalent to the global consensus problem

 minx∈RSMf(x):=1MM∑i=1fi(xi),s.t.(A⊗IS)x=0. (2)

The main benefit of the above formulation is that the objective function is now separable, and the linear constraint encodes the network connectivity pattern.

### 1.2 Distributed non-convex optimization

Non-convex distributed optimization has gained considerable attention recently, and has found applications in training neural networks , distributed information processing and machine learning [2, 3, 4], and distributed signal processing .

The problem (1) and (2) have been studied extensively in the literature when ’s are all convex; see for example [6, 7, 8, 9, 10, 11]. Primal methods such as distributed subgradient method , the EXTRA method , as well as primal-dual based methods such as distributed augmented Lagrangian method , Alternating Direction Method of Multipliers (ADMM) [12, 13, 14, 15] have been proposed.

On the contrary, only recently there have been some works addressing the more challenging problems without assuming convexity of ; see recent developments in [16, 17, 18, 4, 5, 19, 20, 21]. Reference  develops a non-convex ADMM based methods (with the first known global sublinear convergence rate) for solving the distributed consensus problem (1). However the network considered therein is a star network in which the local nodes are all connected to a central controller. Reference  proposes a primal-dual based method for unconstrained problem over a connected network, and derives the first global convergence rate for this setting. In  and follow up works [22, 23], the authors utilize certain gradient tracking idea to solve a constrained nonsmooth distributed problem over possibly time-varying networks. It is worth noting that the distributed algorithms proposed in all these works converge to first-order stationary solutions, which contain local maximum, local minimum and saddle points. Only recently, the authors of  developed first-order distributed algorithms that are capable of computing second-order stationary solutions (which under suitable conditions become local optimal solutions).

### 1.3 Lower and upper rate bounds analysis

Despite all the recent interests and contributions in this field, one major question remains open:

rgb]0.9,0.9,0.9

(Q)   What is the best convergence rate achievable by any distributed algorithms for the non-convex problem (1)?

Question seeks to find a “best convergence rate”, which is a characterization of the smallest number of iterations required to achieve certain high quality solutions, among all distributed algorithms. Clearly, understanding provides fundamental insights to distributed optimization and information processing. For example, the answer to can provide meaningful optimal estimates on the total amount of communication effort required to achieve a given level of accuracy. Further, the identified optimal strategies capable of attaining the best convergence rates will also help guide the practical design of distributed information processing algorithms.

Question is easy to state, but formulating it rigorously is quite involved and a number of delicate issues have to be clarified. Below we provide a high level discussion on some of the issues.

(1) Fix Problem and Network Classes. A class of problems and networks of interest should be fixed. Roughly speaking, in this work is the family of smooth unconstrained problem (1), and is defined over the set of connected and unweighted graphs with finite number of nodes.

(2) Characterize High-Quality Solutions. For a properly defined error constant , one needs to define a high-quality solution in distributed and non-convex setting. Differently from the centralized case, the following questions have to be addressed: Should the solution quality be evaluated based on the averaged iterates among all the agents, or on the individual iterates? Shall we include some consensus measure in the solution characterization? Different solution notion could potentially lead to different lower and upper rate bounds.

(3) Fix Algorithm Classes. A class of algorithms has to be fixed. In the classical complexity analysis in (centralized) optimization, it is common to define the class of algorithms by the information structures that they utilize . In the distributed and non-convex setting, it is necessary to specify both function information that can be used by individual nodes, as well as the communication protocols allowed.

(4) Develop Sharp Upper Bounds. It is necessary to develop algorithms within class , which possess provable and sharp global convergence rate for problem/network class . These algorithms provide achievable upper bounds on the global convergence rates.

(5) Identify Lower Bounds. It is important to characterize the worst rates achievable by any algorithm in class for problem/network class . This task involves identifying instances in that are difficult for algorithm class .

(6) Match Lower and Upper Bounds. The key task is to investigate whether the developed algorithms are rate optimal, in the sense that rate upper bounds match the worst-case lower bounds. Roughly speaking, matching two bounds requires that for the class of problem and networks , the following quantities should be matched between the lower and upper bounds: i) the order of the error constants ; ii) the order of problem parameters such as , or that of network parameters such as the spectral gap, diameter, etc.

Convergence rate analysis (aka iteration complexity analysis) for convex problems dates back to Nesterov, Nemirovsky and Yudin [26, 27], in which lower bounds and optimal first-order algorithms have been developed; also see . In recent years, many accelerated first-order algorithms achieving those lower bounds for different kinds of convex problems have been derived; see e.g., [29, 30, 31], including those developed for distributed convex optimization . In those works, the optimality measure used is , and the lower bound can be expressed as [28, Theorem 2.2.2]

 f(xr)−f(x∗)≤∥x0−x∗∥L(r+2)2, (3)

where is the Lipschitz constant for ; (resp. ) is the global optimal solution (resp. the initial solution); is the iteration index. Therefore to achieve -optimal solution in which , one needs iterations. Recently the above approach has been extended to distributed strongly convex optimization in . In particular, the authors consider problem (1) in which each is strongly convex, and they provide lower and upper rate bounds for a class of algorithms in which the local agents can utilize both and its Fenchel conjugate . Unfortunately this result is not directly related to the class of “first-order” method of interest, since beyond the first-order gradient information, the Fenchel conjugate is also needed, but computing this quantity requires performing certain exact minimization, which itself involves solving a strongly convex optimization problem (with unknown complexity). Recently related works in this direction also include . Therefore the optimal first-order distributed algorithm for strongly convex problems are still left open, not to mention for general convex and non-convex distributed problems.

When the problem becomes non-convex, the size of the gradient function can be used as a measure of solution quality. In particular, let , then it has been shown that the classical gradient descent (GD) method achieves the following rate [28, page 28]

 h∗T≤c0L(f(x0)−f(x∗))T+1,where% c0>0 is some constant.

It has been shown in  that the above rate is (almost) tight for GD. Recently,  has further shown that the above rate is optimal for any first-order methods that only utilizes the gradient information, when applied to problems with Lipschitz gradient. However, no lower bound analysis has been developed for distributed non-convex problem (14); there are even not many algorithms that provide achievable upper rate bounds (except for the recent works [4, 20, 37, 38]), not to mention any analysis on the tightness/sharpness of these upper bounds.

### 1.4 Contribution of this work

In this work, we address various issues arise in answering . Our main contributions are given below:

1) We identify a class of non-convex problems and networks , a class of distributed first-order algorithms , and rigorously define the -optimality gap that measures the progress of the algorithms;

2) We develop the first lower rate bound for class to solve class : to achieve -optimality, it is necessary for any to perform rounds of communication among all the nodes;

3) We design a class of algorithms belonging to [named distributed gradient primal-dual algorithm (D-GPDA), and its variants], which achieves -optimality condition with provable global rates [in the order of ];

4) We show that the D-GPDA class are optimal methods in for problem class as well as a number of its refinements, in that they precisely achieve the lower complexity bounds under all these situations (up to a ploylog factor).

The main results of the paper are outlined in Table 1.

Notations. For a given symmetric matrix , we use , and to denote the maximum, the minimum and the minimum nonzero eigenvalues for a symmetric matrix ; We use to denote an identity matrix with size . We use to denote the set . For a vector we use to denote its th element. We use to denote where is the problem dimension. We use to denote two connected nodes and .

## 2 Preliminaries

### 2.1 The class P, N, A

We present the classes of problems, networks and algorithms to be studied, as well as some useful results. We parameterize these classes using a few key parameters so that we can specify their subclasses when needed.

Problem Class. A problem is in class if it satisfies the following conditions.

• The objective is a sum of functions; see (1).

• Each component function ’s has Lipschitz gradient:

 ∥∇fi(xi)−∇fi(zi)∥≤Li∥xi−zi∥,∀ xi,zi∈RS,∀ i. (4)

Define the matrix of Lipschitz constants as:

 L=diag([L1,⋯,LM])⊗IS∈RMS×MS. (5)
• The function is lower bounded over , i.e.,

 f––:=infxf(x)>−∞. (6)

These assumptions are rather mild. For example an satisfies [A2-A3] is not required to be second-order differentiable. Below we provide a few non-convex functions that satisfy Assumption [A2-A3], and each of those can be the component function ’s. Note that the first four functions are of particular interest in learning neural networks, as they are commonly used as activation functions.

(1) The sigmoid function is given by We have , , therefore [A2-A3] are true with .

(2) The function satisfies , . So [A2-A3] hold with .

(3) The function satisfies so [A2-A3] hold with .

(4) The logit function is related to the function as follows

 2logit(x)=2exex+1=1+tanh(x/2),

then Assumptions [A2-A3] are again satisfied.

(5) The function has applications in structured matrix factorization . Clearly it is lower bounded. Its second order derivative is and it is also bounded.

(6) Other functions like , , are easy to verify. Consider where . This function is interesting because it is not second-order differentiable; nonetheless we can verify that [A2-A3] are satisfied with .

Network Class. Consider the following descriptions.

• The network is represented by an undirected and unweighted graph , with vertices and edges, and edge weights all being .

• The graph is connected, with graph Laplacian matrix , and the normalized Lapacian (NL) matrix 

 LN:=P−1/2ATAP−1/2,withP=% diag[d1,⋯,dM]. (7)
• The graph has diameter , which is the largest distance path among any pair of nodes:

 D:=maxu,v∈Vdist(u,v). (8)

We use to denote a class of network described in with nodes and being its diameter. Based on different definitions of Laplacian matrices, we define the spectral gap of the graph as

 ξ(G)=λ––min(L)λmax(L),ξN(G)=λ––min(LN)λmax(LN). (9)

Algorithm Class. Define the neighbor set for node as

 Ni:={i∣i∼i,j≠i}. (10)

We say that a distributed, first-order algorithm is in class if it satisfies the following conditions.

• At iteration , each node can obtain some network related constants, such as , , , , etc.

• At iteration , the output of node is a linear combination of the historical output of its neighboring set, as well as historical local gradients, i.e.,

 xr+1i∈span{{(xtj)rt=1}j∈Ni∪i,{∇fi(xti)}rt=1}. (11)

The linear combination coefficients can be dependent on those constants obtained at iteration .

Clearly is a class of first-order methods because only historical gradient information is used in the computation. It is also a class of distributed algorithms because at each iteration the nodes are only allowed to communicate with its neighbors. Note that we do not specify the per-iteration communication patterns, therefore each can receive either one or multiple vectors of size , as long as they are coming from the set specified in (11). This adds flexibility to the practical algorithm design, for example one can allow asynchronous communication in which at a given iteration, delayed messages from the past are also received (although we do not investigate these possibilities in the current work).

### 2.2 Solution Quality Measure

Next we provide definition for the quality of the solution. Since we consider using first-order methods to solve non-convex problems, it is expected that in the end some first-order stationary solution will be computed, in which the size of the gradient of is small.

Our first definition is related to a global variable . We say that is a global--solution if the following holds

 yr∈span{xri}Mi=1,minr∈[T]∥∇g(yr)∥2≤ϵ. (12)

This definition is conceptually simple and it is identical to the centralized criteria in Section 1.3. However it has the following issues. First, no global variable will be formed in the entire network, so criteria (12) is difficult to evaluate. Second, there is no characterization of how close the local variables ’s are. To see the second point, consider the following toy example arise in the non-convex setting.

Example 1: Consider a network with and and . Suppose local variable and . Then if we pick , we have

 ∇g(yT)=12(∇f1(yT)+∇f2(yT))=0.

This suggests that at iteration , there exists one linear combination that makes measure (12) precisely zero. However one can hardly say that the current solution is a good solution for problem (2).

To address the above issues, we provide a second definition which is directly related to local variables . We say that is a local -solution if the following holds

 minr∈[T]∥∥∥M∑i=1∇fi(xi)M∥∥∥2+1M2∑i,j√LiLj∥xi−xj∥2≤ϵ. (13)

Clearly this definition takes into consideration the consensus error as well as the size of the local gradients. When applied to Example 1, this measure will be large. We will use short-handed notation to denote the left hand side of (13).

In our work we will focus on providing answers to the following question: For any given , how many iterations (as a function of ) does it need for any algorithm in class to achieve the local -solution.

### 2.3 Some Useful Facts

Below we provide a few facts about the above classes.

On Lipschitz constants. Assume that each has Lipschitz continuous gradient with constant in (4). Then we have

 ∥∇¯f(y1)−∇¯f(y2)∥≤M∑i=11MLi∥y1−y2∥:=¯L∥y1−y2∥. (14)

Then we have the following [the matrix is defined in (5)]

 ∥∇f(x)−∇f(z)∥2 =1M2M∑i=1∥∇fi(xi)−∇fi(zi)∥2,

which implies

 ∥∇f(x)−∇f(z)∥≤1M∥L(x−z)∥,∀ x,z∈RMS. (15)

On Quantities for Graph . Define the following matrices:

 Σ:=diag[σ1,⋯,σE]≻0,Υ:=diag([β1,⋯,βM])≻0. (16)

Define where the absolute value is taken component-wise. Then we have

 12(ATA+BTB)=P=diag[d1,⋯,dM]∈RM×M (17) 12(ATΣ2A+BTΣ2B)=% diag({∑j:i∼jσ2ij}j∈N):=Δ,

where is the degree matrix defined in (7). The NL matrix defined in (7) can be expressed as

 [LN]ij=⎧⎪ ⎪ ⎪⎨⎪ ⎪ ⎪⎩1if i=j−1√didjif i∼j,i≠j0otherwise.

For two diagonal matrices and of appropriate sizes, the normalized generalized Laplacian (NGL) matrix is defined as

 ˜L=Υ−1ATΣ2AΥ−1, (18)

and in particular

 [˜L]ij=⎧⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎨⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎩∑q:i∼qσ2iqβ2iif i=j−σ2ijβi×βjif (ij)∈E,i≠j0otherwise.

Below we list some useful results about NL  . First, all eigenvalues of lie in the interval . The eigenvalues of NL for a number of special graphs are given below:

1) Complete Graph: The eigenvalues are and (with multiplicity ), so ;

2) Star Graph: The eigenvalues are and (with multiplicity ), and , so ;

3) Path Graph: The eigenvalues are for , and .

4) Cycle Graph: The eigenvalues are for , and .

5) -Path+Star Graph: Consider the following graph. First construct a path graph with nodes, then divide the remaining nodes into groups, each with either or nodes. Starting from the two ends, connecting each group with one of the nodes using a star topology. It is easy to verify that the resulting graph has the following property:

 maxidi≤MD−1+2,minidi=1,∑idi≤2M. (19)

Then for this new graph we have . This fact is due to [40, Lemma 1.9], that

 λ––min(LN)≥1D∑idi≥12DM. (20)

Since for NL , we get the desired bound.

Below we show an extension of the result [40, Section 1.4], on the lower bound for the second smallest eigenvalue of .

###### Lemma 2.1

Suppose that is a connected graph. We have

 λ––min(˜L)≥min(i,j)∈Eσ2ijD∑iβ2i. (21)

## 3 Lower Complexity Bounds

In this section we develop the lower complexity bounds for algorithms in class to solve problems over network . Our proof is for the case where ’s have uniform Lipschitz constants, i.e., , , but our analysis approach can be extended to non-uniform cases as well. Our proof combines ideas from the classical proof in Nesterov , as well as two recent constructions  (for centralized non-convex problems) and  (for strongly convex distributed problems). Our construction differs from the previous works in a number of ways, in particular, the constructed functions are only first-order differentiable, but not second-order differentiable. Further, we use the local- solution (13) to measure the quality of the solution, which makes the analysis more involved compared with the existing global error measures in [25, 36, 33].

To begin with, we construct the following two non-convex functions

 h(x):=1MM∑i=1hi(xi),f(x):=1MM∑i=1fi(xi), (22)

as well as their average version

 ¯h(y):=1MM∑i=1hi(y),¯f(y):=1MM∑i=1fi(y). (23)

Here we have , for all , , and . Later we make our construction so that functions and are easy to analyze, while and will be in the desired function class in . Without loss of generality, in the construction we will assume will be Lipschitz with constant , for all .

### 3.1 Path Graph (D=m−1)

First we consider the extreme case in which the nodes form a path graph with nodes and each node has its own local function , shown in Figure 1.

For notational simplicity assume that is a multiple of , that is for some integer . Also assume that is an odd number without loss of generality.

Let us define the component functions ’s in (22) as follows.

 hi(xi)=⎧⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎨⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎩Θ(xi,1)+3⌊T/2⌋∑j=1Θ(xi,2j),i∈[1,M3]Θ(xi,1),i∈[M3+1,2M3]Θ(xi,1)+3⌊T/2⌋∑j=1Θ(xi,2j+1),i∈[2M3+1,M]] (24)

where we have defined the following functions

 Θ(xi,j) :=Ψ(−xi[j−1])Φ(−xi[j])−Ψ(xi[j−1])Φ(xi[j]),∀ j≥2 (25a) Θ(xi,1) :=−Ψ(1)Φ(xi). (25b)

The component functions are given as below

 Ψ(w):={0w≤01−e−w2w>0,andΦ(w):=4arctanw+2π.

Suppose , then the average function becomes:

 ¯h(y) :=1MM∑j=1hi(y)=Θ(y,1)+T∑i=2Θ(y,i) =−Ψ(1)Φ(y)+T∑i=2[Ψ(−y[i−1])Φ(−y[i])−Ψ(y[i−1])Φ(y[i])].

Further for a given error constant and a given averaged Lipschitz constant , let us define

 fi(xi):=150πϵUhi(xiU75π√2ϵ). (26)

Therefore we also have, if , then

 ¯f(y):=1MM∑i=1fi(y)=150πϵU¯h(yU75π√2ϵ). (27)

First we present some properties of the component functions ’s.

###### Lemma 3.1

The functions and satisfy the following.

1. For all , , .

2. We have the following bounds for the functions and their first and second-order derivatives:

 0≤Ψ(w)<1,  0≤Ψ′(w)≤√2e,  −4e32≤Ψ′′(w)≤2,  ∀w>0
 0<Φ(w)<4π,  0<Φ′(w)≤4,  −3√32≤Φ′′(w)≤3√32,  ∀w∈R
3. We have the following key property:

 Ψ(w)Φ′(v)>1,∀ w≥1,|v|<1. (28)
4. The function is lower bounded as follows:

 hi(0)−infxihi(xi)≤10πT,h(0)−infxh(x)≤10πT.
5. The first order derivative of (resp. ) is Lipschitz continuous with some constant (resp. , ).

Proof. Property 1) is obviously true.

To prove Property 2), note that following holds for (note, is not second order differentiable over ).

 Ψ(w)=1−e−w2,  Ψ′(w)=2e−w2w,  Ψ′′(w)=2e−w2−4e−w2w2. (29)

Obviously, is an increasing function (), therefore the lower and upper bounds are ; is increasing on and decreasing on , where , therefore the lower and upper bounds are ; is decreasing on and increasing on , therefore the lower and upper bounds are , i.e.,

 0≤Ψ(w)<1,  0≤Ψ′(w)≤√2e,  −4e32≤Ψ′′(w)≤2,  ∀w>0

Further, for all , the following holds

 Φ(w)=4arctanw+2π,  Φ′(w)=4w2+1,  Φ′′(w)=−8w(w2+1)2. (30)

Similarly, we have

 0<Φ(w)<4π,  0<Φ′(w)≤4,  −3√32≤Φ′′(w)≤3√32,  ∀w∈R.

To show Property 3), note that for all and ,

 Ψ(w)Φ′(v)>Ψ(1)Φ′(1)=2(1−e−1)>1

where the first inequality is true because is strictly increasing and is strictly decreasing for all , and that .

Next we show Property 4). Note that and . Therefore we have and using the construction in (24)

 infxihi(xi) (31) ≥−4π−6πT≥−10πT (32)

where the first inequality follows and second follows , we reach the conclusion.

Finally we show Property 5), using the fact that a function is Lipschitz if it is piecewise smooth with bounded derivative. From construction (24), the first order partial derivative of can be expressed below.

Case I) If is even, we have

 (33)

Case II) If is odd but not 1, we have

 (34)

Case III) If , we have

 ∂hq∂y=⎧⎨⎩−Ψ(1)Φ′(y)+3(−Ψ′(−y)Φ(−y)−Ψ′(y)Φ(y)),q∈[1,M3]−Ψ(1)Φ′(y),q∈[M3+1,M]. (35)

Obviously, is a piecewise smooth function for any , either equals zero or separated at the non-differentiable point because of the function .

Further, fix a point and a unit vector where . Define

 gq(θ;y,v):=hq(y+θv)

to be the directional projection of on to the direction at point . We will show that there exists such that for all .

 g′′q(0;y,v) =T∑i1,i2=1∂2∂y[i1]∂y[i2]hq(y)v[i1]v[i2]=∑δ∈{0,1,−1}T∑i=1∂2∂y[i]∂y[i+δ]hq(y)v[i]v[i+δ],

where we take and .

The second order partial derivative of () is given as follows when is even:

 ∂2hq∂y[i]∂y[i]=⎧⎪ ⎪ ⎪⎨⎪ ⎪ ⎪⎩3(Ψ(−y[i−1])Φ′′(−y[i])−Ψ(y[i−1])Φ′′(y[i])),q∈[1,M3]0,q∈[M3+1,2M3]3(Ψ′′(−y[i])Φ(−y[i+1])−Ψ′′(y[i])Φ(y[i+1])),q∈[2M3+1,M]
 ∂2hq∂y[i]∂y[i+1]=⎧⎨⎩0,q∈[1,2M3]3(Ψ′(−y[i])Φ′(−y[i+1])−Ψ′(y[i])Φ′(y[i+1])),q∈[2M3+1,M]
 ∂2hq∂y[i]∂y[i−1]=⎧⎨⎩3(Ψ′(−y[i−1])Φ′(−y[i])−Ψ′(y[i−1])Φ′(y[i])),q∈[1,M3]0,q∈[M3+1,M]

Apply Lemma 3.1 i) and ii), we have:

 ∣∣∣∂2hq∂y[i1]∂y[i2]∣∣∣ ≤3max{sup|Ψ′′(w)|sup|Φ(v)|,sup|Ψ(w)|sup|Φ′′(v)|,sup|Ψ′(w)|sup|Φ′(v)|} =3max{8π,3√32,4√2e}<25π.

Note the bound here is exactly the same for is odd but not , while for we have following

 ∂2hq∂y∂y=⎧⎨⎩−Ψ(1)Φ′′(y)+3(−Ψ′′(−y)Φ(−y)−Ψ′′(y)Φ(y)),q∈[1,M3]−Ψ(1)Φ′′(y),q∈[M3+1,M]
 ∂2hq∂y∂y=⎧⎨⎩3(−Ψ′(−y)Φ′(−y)−Ψ′(y)Φ′(y)),q∈[1,M3]0,q∈[M3+1,M]
 ∣∣∣∂2hq∂y[i1]∂y[i2]∣∣∣ ≤max{sup|Ψ(1)Φ′′(w)|+3sup|Ψ′′(w)|sup|Φ(v)|,3sup|Ψ′(w)|sup|Φ′(v)|} =max{3√32(1−e−1)+24π,12√2e}<25π.

Therefore we have

 |g′′q(0;y,v)| =∑δ∈{0,1,−1}T∑i=1∂2∂y[i]∂y[i+δ]hq(y)v[i]v[i+δ] ≤25π∑δ∈{0,1,−1}