Singularity Degree of PSD Matrix Completion Shin-ichi Tanigawa CWI - - PowerPoint PPT Presentation

singularity degree of psd matrix completion
SMART_READER_LITE
LIVE PREVIEW

Singularity Degree of PSD Matrix Completion Shin-ichi Tanigawa CWI - - PowerPoint PPT Presentation

Singularity Degree of PSD Matrix Completion Shin-ichi Tanigawa CWI and Kyoto July 29, 2016 1 / 13 Positive Semidefinite Matrix Completion PSD completion problem ( G , c ) Given G = ( V , E ) with V = { 1 , . . . , n } and edge weight c : E


slide-1
SLIDE 1

Singularity Degree of PSD Matrix Completion

Shin-ichi Tanigawa

CWI and Kyoto

July 29, 2016

1 / 13

slide-2
SLIDE 2

Positive Semidefinite Matrix Completion

PSD completion problem (G, c)

Given G = (V , E) with V = {1, . . . , n} and edge weight c : E → [−1, 1], find X ∈ Sn s.t. X[i, j] = c(ij) (ij ∈ E) X[i, i] = 1 (i ∈ V ) X 0

2 / 13

slide-3
SLIDE 3

Positive Semidefinite Matrix Completion

PSD completion problem (G, c)

Given G = (V , E) with V = {1, . . . , n} and edge weight c : E → [−1, 1], find X ∈ Sn s.t. X[i, j] = c(ij) (ij ∈ E) X[i, i] = 1 (i ∈ V ) X 0 min Ω, C s.t. Ω ∈ S+(G) where C[i, j] =      c(ij) (ij ∈ E) 1 (i = j) (otherwise) S(G) := {A ∈ Sn : A[i, j] = 0 ∀ij / ∈ V ∪ E} S+(G) := {A ∈ S(G) : A 0} S+(G) := {A ∈ S+(G) : A[i, j] = 0 ∀ij ∈ E}

2 / 13

slide-4
SLIDE 4

Geometric View

Given a completion problem (G, c), PSD completion X = PP⊤ with rank d ⇔ spherical embedding p : V → Sd−1 realizing c, i.e., pi · pj = c(ij) ∀ij ∈ E

◮ spherical (bar-joint) framework (G, p) 3 / 13

slide-5
SLIDE 5

Geometric View

Given a completion problem (G, c), PSD completion X = PP⊤ with rank d ⇔ spherical embedding p : V → Sd−1 realizing c, i.e., pi · pj = c(ij) ∀ij ∈ E

◮ spherical (bar-joint) framework (G, p)

dual optimal solution : Ω ∈ S+(G) with C, Ω = 0 C, Ω = 0 ⇔ X, Ω = 0 ⇔ ΩP = 0 ⇔ Ω[i, i]p(i) +

  • j∼i

Ω[i, j]p(j) = 0 (∀i ∈ V ) (1)

3 / 13

slide-6
SLIDE 6

Geometric View

Given a completion problem (G, c), PSD completion X = PP⊤ with rank d ⇔ spherical embedding p : V → Sd−1 realizing c, i.e., pi · pj = c(ij) ∀ij ∈ E

◮ spherical (bar-joint) framework (G, p)

dual optimal solution : Ω ∈ S+(G) with C, Ω = 0 C, Ω = 0 ⇔ X, Ω = 0 ⇔ ΩP = 0 ⇔ Ω[i, i]p(i) +

  • j∼i

Ω[i, j]p(j) = 0 (∀i ∈ V ) (1) Ω is called a stress (matrix) of (G, p) if Ω satisfies (1) Given (G, p), Ω ∈ S(G) is dual opt iff Ω is a PSD stress of (G, p).

3 / 13

slide-7
SLIDE 7

SDP Duality

For any primal and dual optimal pair (X, Ω), X, Ω = 0 ⇒ rank X + rank Ω ≤ n. high rank dual opt ⇒ low rank completion

4 / 13

slide-8
SLIDE 8

SDP Duality

For any primal and dual optimal pair (X, Ω), X, Ω = 0 ⇒ rank X + rank Ω ≤ n. high rank dual opt ⇒ low rank completion

Rank maximality certificate

A completion X for (G, c) attains the maximum rank if ∃ dual opt with rank n − rank X.

4 / 13

slide-9
SLIDE 9

Parameter ν and Unique Completability

Theorem (Connelly82, Laurent-Varvitsiotis14)

A completion X for (G, c) is unique if ∃ dual opt Ω with rank Ω = n − rank X and the SAP, i.e., ∄X ∈ Sn \ {0} with ΩX = 0 and X[i, j] = 0 for ij∈V∪E (G, p) is universally rigid in Sd−1 if (G, p) admits a PSD stress Ω with rank Ω = n − d and the SAP.

5 / 13

slide-10
SLIDE 10

Parameter ν and Unique Completability

Theorem (Connelly82, Laurent-Varvitsiotis14)

A completion X for (G, c) is unique if ∃ dual opt Ω with rank Ω = n − rank X and the SAP, i.e., ∄X ∈ Sn \ {0} with ΩX = 0 and X[i, j] = 0 for ij∈V∪E (G, p) is universally rigid in Sd−1 if (G, p) admits a PSD stress Ω with rank Ω = n − d and the SAP.

Colin de Verdi` ere Parameter ν

ν(G) := max{corank Ω : Ω ∈ S+(G) has the SAP}. ν(G) ≤ max{d : ∃ universally rigid (G, p) in Sd−1 }

5 / 13

slide-11
SLIDE 11

Strict Complementarity and Singularity Degree

Strict Complementarity

A primal and dual optimal pair (X, Ω) satisfies a strict complementarity condition if rank X + rank Ω = n For which problem the strict complementarity can be guaranteed? How far from the strict complementarity?

6 / 13

slide-12
SLIDE 12

Strict Complementarity and Singularity Degree

Strict Complementarity

A primal and dual optimal pair (X, Ω) satisfies a strict complementarity condition if rank X + rank Ω = n For which problem the strict complementarity can be guaranteed? How far from the strict complementarity? ⇒ singularity degree of SDP

6 / 13

slide-13
SLIDE 13

Strict Complementarity and Singularity Degree

Strict Complementarity

A primal and dual optimal pair (X, Ω) satisfies a strict complementarity condition if rank X + rank Ω = n For which problem the strict complementarity can be guaranteed? How far from the strict complementarity? ⇒ singularity degree of SDP ⇒ singularity degree of a graph G

6 / 13

slide-14
SLIDE 14

Strict Complementarity and Singularity Degree

Strict Complementarity

A primal and dual optimal pair (X, Ω) satisfies a strict complementarity condition if rank X + rank Ω = n For which problem the strict complementarity can be guaranteed? How far from the strict complementarity? ⇒ singularity degree of SDP ⇒ singularity degree of a graph G

Proposition

The following are equivalent for a graph G:

1

sd(G) = 1;

2

The strict complementarity holds for any PSD completion problem with underlying graph G;

6 / 13

slide-15
SLIDE 15

Strict Complementarity and Singularity Degree

Strict Complementarity

A primal and dual optimal pair (X, Ω) satisfies a strict complementarity condition if rank X + rank Ω = n For which problem the strict complementarity can be guaranteed? How far from the strict complementarity? ⇒ singularity degree of SDP ⇒ singularity degree of a graph G

Proposition

The following are equivalent for a graph G:

1

sd(G) = 1;

2

The strict complementarity holds for any PSD completion problem with underlying graph G;

3

The projection E(G) of the elliptope (the set of correlation matrices) onto RE is exposed (Druvyatskiy-Pataki-Wolkowicz15).

6 / 13

slide-16
SLIDE 16

Facial Reduction (Borwein-Wolkowitcz81)

A sequence {Ω1, . . . , Ωk} in Sn is iterated PSD if Ωi is positive semidefinite on Vi−1, where V0 = Rn and Vi = {x ∈ Rn : xx⊤, Ωj = 0 (j = 1, . . . , i − 1)}.

7 / 13

slide-17
SLIDE 17

Facial Reduction (Borwein-Wolkowitcz81)

A sequence {Ω1, . . . , Ωk} in Sn is iterated PSD if Ωi is positive semidefinite on Vi−1, where V0 = Rn and Vi = {x ∈ Rn : xx⊤, Ωj = 0 (j = 1, . . . , i − 1)}.

Theorem (Facial reduction)

For any feasible (G, c), ∃X and ∃Ω1, . . . , Ωk ∈ S(G) s.t.

1

the sequence is iterated PSD

2

C, Ωi = 0 for each i

3

rank X = dim Vk

7 / 13

slide-18
SLIDE 18

Facial Reduction (Borwein-Wolkowitcz81)

A sequence {Ω1, . . . , Ωk} in Sn is iterated PSD if Ωi is positive semidefinite on Vi−1, where V0 = Rn and Vi = {x ∈ Rn : xx⊤, Ωj = 0 (j = 1, . . . , i − 1)}.

Theorem (Facial reduction)

For any feasible (G, c), ∃X and ∃Ω1, . . . , Ωk ∈ S(G) s.t.

1

the sequence is iterated PSD

2

C, Ωi = 0 for each i

3

rank X = dim Vk the existence of a dual sequence characterizes the max rank of completions (Connelly-Gortler15) with the SAP, it characterize the unique completability (Connelly-Gortler15)

7 / 13

slide-19
SLIDE 19

Facial Reduction (Borwein-Wolkowitcz81)

A sequence {Ω1, . . . , Ωk} in Sn is iterated PSD if Ωi is positive semidefinite on Vi−1, where V0 = Rn and Vi = {x ∈ Rn : xx⊤, Ωj = 0 (j = 1, . . . , i − 1)}.

Theorem (Facial reduction)

For any feasible (G, c), ∃X and ∃Ω1, . . . , Ωk ∈ S(G) s.t.

1

the sequence is iterated PSD

2

C, Ωi = 0 for each i

3

rank X = dim Vk the existence of a dual sequence characterizes the max rank of completions (Connelly-Gortler15) with the SAP, it characterize the unique completability (Connelly-Gortler15)

Definition (Sturm 2000)

For a completion problem (G, c), the singularity degree sd(G, c) is the length of the shortest dual certificate sequence {Ω1, . . . , Ωk}.

7 / 13

slide-20
SLIDE 20

Singularity Degree of Graphs

Singularity degree of G

sd(G) = max

c

sd(G, c) Question (Druvyatskiy-Pataki-Wolkowicz15) Characterize G with sd(G) = 1 Question (So15) sd(G) = o(n)?

8 / 13

slide-21
SLIDE 21

Main Results

Theorem (T16)

sd(G) = 1 iff G is chordal. G is chordal if G has no Cn(n ≥ 4) as an induced subgraph

9 / 13

slide-22
SLIDE 22

Main Results

Theorem (T16)

sd(G) = 1 iff G is chordal. G is chordal if G has no Cn(n ≥ 4) as an induced subgraph

Theorem (T16)

If G has neither Wn (n ≥ 5) nor a proper splitting of Wn (n ≥ 4) as an induced subgraph, then sd(G) ≤ 2. If G has an induced subgraph which is a proper splitting of one of the above forbidden subgraphs, then sd(G) > 2. If tw(G) ≤ 2, then sd(G) ≤ 2.

9 / 13

slide-23
SLIDE 23

Main Results

Theorem (T16)

sd(G) = 1 iff G is chordal. G is chordal if G has no Cn(n ≥ 4) as an induced subgraph

Theorem (T16)

If G has neither Wn (n ≥ 5) nor a proper splitting of Wn (n ≥ 4) as an induced subgraph, then sd(G) ≤ 2. If G has an induced subgraph which is a proper splitting of one of the above forbidden subgraphs, then sd(G) > 2. If tw(G) ≤ 2, then sd(G) ≤ 2.

Theorem (T16)

For each n there is a graph G with n vertices and tw(G) = 3 whose singularity degree is ⌊ n−1

3 ⌋.

9 / 13

slide-24
SLIDE 24

Proof of the first theorem

Theorem (T16)

sd(G) = 1 iff G is chordal. ”⇐” (Druvyatskiy-Pataki-Wolkowicz15) ”⇒” Lemma sd(Cn) ≥ 2 if n ≥ 4.

  • Lemma. sd(G) ≥ sd(H) for any induced subgraph H of G.

10 / 13

slide-25
SLIDE 25

Proof of the first theorem

Theorem (T16)

sd(G) = 1 iff G is chordal. ”⇐” (Druvyatskiy-Pataki-Wolkowicz15) ”⇒” Lemma sd(Cn) ≥ 2 if n ≥ 4. Consider (G, p):

  • Lemma. sd(G) ≥ sd(H) for any induced subgraph H of G.

10 / 13

slide-26
SLIDE 26

Proof of the first theorem

Theorem (T16)

sd(G) = 1 iff G is chordal. ”⇐” (Druvyatskiy-Pataki-Wolkowicz15) ”⇒” Lemma sd(Cn) ≥ 2 if n ≥ 4. Consider (G, p): (G, p) is universally rigid

  • Lemma. sd(G) ≥ sd(H) for any induced subgraph H of G.

10 / 13

slide-27
SLIDE 27

Proof of the first theorem

Theorem (T16)

sd(G) = 1 iff G is chordal. ”⇐” (Druvyatskiy-Pataki-Wolkowicz15) ”⇒” Lemma sd(Cn) ≥ 2 if n ≥ 4. Consider (G, p): (G, p) is universally rigid there is a unique stress Ω with rank Ω = 1 < n − 2

  • Lemma. sd(G) ≥ sd(H) for any induced subgraph H of G.

10 / 13

slide-28
SLIDE 28

Nondegenerate Singularity Degree

A completion problem (G, c) is nondegenerate if c(ij) = ±1 for every ij ∈ E(G). Degenerate edges can easily be eliminated. Suppose c(ij) = 1 for ij ∈ E ... Any solution X of (G, c) satisfies X[i, k] = X[j, k] for every k Equivalently, any embedding p realizing c satisfies p(i) = p(j). The example in the last proof is degenerate...

11 / 13

slide-29
SLIDE 29

Nondegenerate Singularity Degree

A completion problem (G, c) is nondegenerate if c(ij) = ±1 for every ij ∈ E(G).

Nondegenerate Singularity Degree

sd∗(G) = max{sd(G, c) : nondegenerate (G, c)}.

11 / 13

slide-30
SLIDE 30

Nondegenerate Singularity Degree

A completion problem (G, c) is nondegenerate if c(ij) = ±1 for every ij ∈ E(G).

Nondegenerate Singularity Degree

sd∗(G) = max{sd(G, c) : nondegenerate (G, c)}. Theorem (T16) sd∗(G) = 1 iff G has neither Wn(n ≥ 5) nor a proper splitting of Wn (n ≥ 4) as an induced subgraph.

11 / 13

slide-31
SLIDE 31

Nondegenerate Singularity Degree

A completion problem (G, c) is nondegenerate if c(ij) = ±1 for every ij ∈ E(G).

Nondegenerate Singularity Degree

sd∗(G) = max{sd(G, c) : nondegenerate (G, c)}. Theorem (T16) sd∗(G) = 1 iff G has neither Wn(n ≥ 5) nor a proper splitting of Wn (n ≥ 4) as an induced subgraph.

11 / 13

slide-32
SLIDE 32

Nondegenerate Singularity Degree

A completion problem (G, c) is nondegenerate if c(ij) = ±1 for every ij ∈ E(G).

Nondegenerate Singularity Degree

sd∗(G) = max{sd(G, c) : nondegenerate (G, c)}. Theorem (T16) sd∗(G) = 1 iff G has neither Wn(n ≥ 5) nor a proper splitting of Wn (n ≥ 4) as an induced subgraph. If G has no forbidden induced subgraph listed above, then the hyperplane exposing the minimal face is determined by cliques and cycles satisfying the metric inequality with equality

11 / 13

slide-33
SLIDE 33

Nondegenerate Singularity Degree

A completion problem (G, c) is nondegenerate if c(ij) = ±1 for every ij ∈ E(G).

Nondegenerate Singularity Degree

sd∗(G) = max{sd(G, c) : nondegenerate (G, c)}. Theorem (T16) sd∗(G) = 1 iff G has neither Wn(n ≥ 5) nor a proper splitting of Wn (n ≥ 4) as an induced subgraph. If G has no forbidden induced subgraph listed above, then the hyperplane exposing the minimal face is determined by cliques and cycles satisfying the metric inequality with equality Lemma sd(G) ≤ sd∗(G) + 1. Corollary (T16) sd(G) ≤ 2 if G has no forbidden induced subgraph listed above.

11 / 13

slide-34
SLIDE 34

Example of Large Singularity Degree

Theorem (T16)

For each n there is a graph G with n vertices and tw(G) = 3 whose singularity degree is ⌈ n−1

3 ⌉.

v1 u2 v2 w2 u3 w3 v3 v4 v5 u4 u5 w4 w5

12 / 13

slide-35
SLIDE 35

Concluding Remarks

A similar result can be established for EDM

13 / 13

slide-36
SLIDE 36

Concluding Remarks

A similar result can be established for EDM Signed PSD matrix completion and the singularity degree of signed graphs:

◮ ”X[i, j] ≤ c(ij)” or ”X[i, j] ≥ c(ij)” instead of ”X[i, j] = c(ij)” 13 / 13

slide-37
SLIDE 37

Concluding Remarks

A similar result can be established for EDM Signed PSD matrix completion and the singularity degree of signed graphs:

◮ ”X[i, j] ≤ c(ij)” or ”X[i, j] ≥ c(ij)” instead of ”X[i, j] = c(ij)” ◮ primal — the theory of tensegrities by e.g., Connelly ◮ dual — signed Colin de Verdiere parameter by Arav et al. 13 / 13

slide-38
SLIDE 38

Concluding Remarks

A similar result can be established for EDM Signed PSD matrix completion and the singularity degree of signed graphs:

◮ ”X[i, j] ≤ c(ij)” or ”X[i, j] ≥ c(ij)” instead of ”X[i, j] = c(ij)” ◮ primal — the theory of tensegrities by e.g., Connelly ◮ dual — signed Colin de Verdiere parameter by Arav et al. ◮ (T16) sd(G, Σ) ≤ 2 if (G, Σ) is odd-K4-minor free 13 / 13

slide-39
SLIDE 39

Concluding Remarks

A similar result can be established for EDM Signed PSD matrix completion and the singularity degree of signed graphs:

◮ ”X[i, j] ≤ c(ij)” or ”X[i, j] ≥ c(ij)” instead of ”X[i, j] = c(ij)” ◮ primal — the theory of tensegrities by e.g., Connelly ◮ dual — signed Colin de Verdiere parameter by Arav et al. ◮ (T16) sd(G, Σ) ≤ 2 if (G, Σ) is odd-K4-minor free

  • Q. Characterize signed graphs (G, Σ) with sd(G, Σ) = 1.

13 / 13

slide-40
SLIDE 40

Concluding Remarks

A similar result can be established for EDM Signed PSD matrix completion and the singularity degree of signed graphs:

◮ ”X[i, j] ≤ c(ij)” or ”X[i, j] ≥ c(ij)” instead of ”X[i, j] = c(ij)” ◮ primal — the theory of tensegrities by e.g., Connelly ◮ dual — signed Colin de Verdiere parameter by Arav et al. ◮ (T16) sd(G, Σ) ≤ 2 if (G, Σ) is odd-K4-minor free

  • Q. Characterize signed graphs (G, Σ) with sd(G, Σ) = 1.
  • Q. Characterize graphs G with sd(G) ≤ 2.

13 / 13

slide-41
SLIDE 41

Concluding Remarks

A similar result can be established for EDM Signed PSD matrix completion and the singularity degree of signed graphs:

◮ ”X[i, j] ≤ c(ij)” or ”X[i, j] ≥ c(ij)” instead of ”X[i, j] = c(ij)” ◮ primal — the theory of tensegrities by e.g., Connelly ◮ dual — signed Colin de Verdiere parameter by Arav et al. ◮ (T16) sd(G, Σ) ≤ 2 if (G, Σ) is odd-K4-minor free

  • Q. Characterize signed graphs (G, Σ) with sd(G, Σ) = 1.
  • Q. Characterize graphs G with sd(G) ≤ 2.
  • Q. Bound sd(G) by other graph parameters.

13 / 13