Semidefinite Farkas Lemma and Dimensional Rigidity of Bar Frameworks - - PowerPoint PPT Presentation

semidefinite farkas lemma and dimensional rigidity of bar
SMART_READER_LITE
LIVE PREVIEW

Semidefinite Farkas Lemma and Dimensional Rigidity of Bar Frameworks - - PowerPoint PPT Presentation

Semidefinite Farkas Lemma and Dimensional Rigidity of Bar Frameworks A. Y. Alfakih University of Windsor, ON Canada Workshop on Making Models: Stimulating research in rigidity Theory and Spatial-Visual Reasoning Fields Institute, Aug 2014


slide-1
SLIDE 1

Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks

  • A. Y. Alfakih

University of Windsor, ON Canada

Workshop on Making Models: Stimulating research in rigidity Theory and Spatial-Visual Reasoning Fields Institute, Aug 2014

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-2
SLIDE 2

Farkas Lemma

Given an m × n matrix A and m-vector b, Farkas Lemma states that exactly one of the following two statements holds:

there exists x ∈ Rn such that Ax = b, x ≥ 0, there exists y ∈ Rm such that ATy ≤ 0, bTy > 0.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-3
SLIDE 3

Farkas Lemma

Given an m × n matrix A and m-vector b, Farkas Lemma states that exactly one of the following two statements holds:

there exists x ∈ Rn such that Ax = b, x ≥ 0, there exists y ∈ Rm such that ATy ≤ 0, bTy > 0.

Farkas Lemma is at the heart of optimization theory.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-4
SLIDE 4

Farkas Lemma

Given an m × n matrix A and m-vector b, Farkas Lemma states that exactly one of the following two statements holds:

there exists x ∈ Rn such that Ax = b, x ≥ 0, there exists y ∈ Rm such that ATy ≤ 0, bTy > 0.

Farkas Lemma is at the heart of optimization theory. In one direction the proof is easy. In the other direction, the proof is based on the separation theorm. S a

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-5
SLIDE 5

Farkas Lemma

Given an m × n matrix A and m-vector b, Farkas Lemma states that exactly one of the following two statements holds:

there exists x ∈ Rn such that Ax = b, x ≥ 0, there exists y ∈ Rm such that ATy ≤ 0, bTy > 0.

Farkas Lemma is at the heart of optimization theory. In one direction the proof is easy. In the other direction, the proof is based on the separation theorm. S a

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-6
SLIDE 6

LP and SDP

The linear programming problem: p∗ = min cTx (P) : subject to Ax = b x ≥ 0

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-7
SLIDE 7

LP and SDP

The linear programming problem: p∗ = min cTx d∗ = max bTy (P) : subject to Ax = b (D) : subject to ATy ≤ c x ≥ 0

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-8
SLIDE 8

LP and SDP

The linear programming problem: p∗ = min cTx d∗ = max bTy (P) : subject to Ax = b (D) : subject to ATy ≤ c x ≥ 0 LP strong duality: If both (P) or (D) are feasible, then both (P) and (D) have optimal solutions x∗ and y ∗ respectively, and p∗ = cTx∗ = d∗ = bTy ∗.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-9
SLIDE 9

LP and SDP

The linear programming problem: p∗ = min cTx d∗ = max bTy (P) : subject to Ax = b (D) : subject to ATy ≤ c x ≥ 0 LP strong duality: If both (P) or (D) are feasible, then both (P) and (D) have optimal solutions x∗ and y ∗ respectively, and p∗ = cTx∗ = d∗ = bTy ∗. LP strict complementarity: If both (P) or (D) are feasible, then ∃ optimal pair x∗ ≥ 0 and s∗ = c − ATy ∗ ≥ 0 such that x∗

i s∗ i = 0 and x∗ i + s∗ i > 0 for all i.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-10
SLIDE 10

LP and SDP

The linear programming problem: p∗ = min cTx d∗ = max bTy (P) : subject to Ax = b (D) : subject to ATy ≤ c x ≥ 0 LP strong duality: If both (P) or (D) are feasible, then both (P) and (D) have optimal solutions x∗ and y ∗ respectively, and p∗ = cTx∗ = d∗ = bTy ∗. LP strict complementarity: If both (P) or (D) are feasible, then ∃ optimal pair x∗ ≥ 0 and s∗ = c − ATy ∗ ≥ 0 such that x∗

i s∗ i = 0 and x∗ i + s∗ i > 0 for all i.

Semidefinite Programming (SDP): p∗ = min trace(CX) d∗ = max

  • i biyi
  • s. t.

trace(AiX) = bi

  • s. t.
  • i Aiy C

X 0

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-11
SLIDE 11

Semidefinite Farkas Lemma

Farkas Lemma: Ax = b, x ≥ 0, ATy ≤ 0, bTy > 0.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-12
SLIDE 12

Semidefinite Farkas Lemma

Farkas Lemma: Ax = b, x ≥ 0, ATy ≤ 0, bTy > 0. Potential SD Farkas Lemma: trace(AiX) = bi, X 0,

  • i Aiyi 0, bTy > 0.
  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-13
SLIDE 13

Semidefinite Farkas Lemma

Farkas Lemma: Ax = b, x ≥ 0, ATy ≤ 0, bTy > 0. False SD Farkas Lemma: trace(AiX) = bi, X 0,

  • i Aiyi 0, bTy > 0.
  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-14
SLIDE 14

Semidefinite Farkas Lemma

Farkas Lemma: Ax = b, x ≥ 0, ATy ≤ 0, bTy > 0. False SD Farkas Lemma: trace(AiX) = bi, X 0,

  • i Aiyi 0, bTy > 0.

Correct SD Farkas Lemma: Given symmetric matrices A1, . . . , Am and m-vector b. Then exactly one of the following two statements holds:

∃ X ≻ 0: trace (AiX) = bi, for i = 1, . . . , m, ∃ Ω = y1A1 + · · · + ymAm: Ω 0, = 0,

i biyi ≤ 0.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-15
SLIDE 15

Semidefinite Farkas Lemma

Farkas Lemma: Ax = b, x ≥ 0, ATy ≤ 0, bTy > 0. False SD Farkas Lemma: trace(AiX) = bi, X 0,

  • i Aiyi 0, bTy > 0.

Correct SD Farkas Lemma: Given symmetric matrices A1, . . . , Am and m-vector b. Then exactly one of the following two statements holds:

∃ X ≻ 0: trace (AiX) = bi, for i = 1, . . . , m, ∃ Ω = y1A1 + · · · + ymAm: Ω 0, = 0,

i biyi ≤ 0.

This Farkas Lemma is used to prove SD strong duality under Slater Condition.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-16
SLIDE 16

Semidefinite Farkas Lemma Cont’d

SD Farkas Lemma: Given symmetric matrices A1, . . . , Am and m-vector b and assume that there exists X ∗ 0 such that trace (AiX ∗) = bi for i = 1, . . . , m. Then exactly one of the following two statements holds:

∃ X ≻ 0: trace (AiX) = bi, for i = 1, . . . , m, ∃ Ω = y1A1 + · · · + ymAm: Ω 0, = 0, trace (ΩX ∗) = 0.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-17
SLIDE 17

Semidefinite Farkas Lemma Cont’d

SD Farkas Lemma: Given symmetric matrices A1, . . . , Am and m-vector b and assume that there exists X ∗ 0 such that trace (AiX ∗) = bi for i = 1, . . . , m. Then exactly one of the following two statements holds:

∃ X ≻ 0: trace (AiX) = bi, for i = 1, . . . , m, ∃ Ω = y1A1 + · · · + ymAm: Ω 0, = 0, trace (ΩX ∗) = 0.

Let F = {X 0 :trace (AiX) = bi, for i = 1, . . . , m}. Then F has no interior iff F is contained in the hyperplane: {X : trace (ΩX) = 0}.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-18
SLIDE 18

A New SD Farkas Lemma

Given symmetric matrices A1, . . . , Am and m-vector b and assume that there exists X ∗ 0 such that trace (AiX ∗) = bi for i = 1, . . . , m and rank X ∗ = r. Then exactly one of the following two statements holds: ∃ X 0: trace (AiX) = bi, for i = 1, . . . , m, rank X ≥ r + 1, ∃ non-zero matrices Ω0, Ω1, . . . , Ωk, k ≤ n − r − 1, such that:

Ωj =

i yj i Ai,

Ω0 0, UT

1 Ω1U1 0, · · · , UT k ΩkUk 0,

trace (ΩjX ∗) = 0 for all j, rank Ω0 + rank (UT

1 Ω1U1) + · · · + rank (UT k ΩkUk) = n − r,

where U1, . . . , Uk and W0, W1, . . . , Wk are full column rank matrices defined as follows: for i = 0, 1, . . . , k, R(Wi) = N(U T

i ΩiUi), and Ui+1 = UiWi, U0 = I.

back

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-19
SLIDE 19

The PSD Cone

Sn denotes the space of n × n symmetric matrices. Sn

+ denotes the cone of n × n symmetric PSD matrices. It is

closed and convex.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-20
SLIDE 20

The PSD Cone

Sn denotes the space of n × n symmetric matrices. Sn

+ denotes the cone of n × n symmetric PSD matrices. It is

closed and convex. A subset F ⊂ Sn

+ is a face of Sn + if:

∀X, Y ∈ Sn

+ : (X + Y ) ∈ F =

⇒ X ∈ F and Y ∈ F.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-21
SLIDE 21

The PSD Cone

Sn denotes the space of n × n symmetric matrices. Sn

+ denotes the cone of n × n symmetric PSD matrices. It is

closed and convex. A subset F ⊂ Sn

+ is a face of Sn + if:

∀X, Y ∈ Sn

+ : (X + Y ) ∈ F =

⇒ X ∈ F and Y ∈ F. The faces of Sn

+ are in one-to-one correspondence with the

subspaces of Rn. i.e., F is a face of Sn

+ iff:

F = {X ∈ Sn

+ : L ⊆ N(X)}, for some subspace L,

relint(F) = {X ∈ Sn

+ : L = N(X)},

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-22
SLIDE 22

The PSD Cone

Sn denotes the space of n × n symmetric matrices. Sn

+ denotes the cone of n × n symmetric PSD matrices. It is

closed and convex. A subset F ⊂ Sn

+ is a face of Sn + if:

∀X, Y ∈ Sn

+ : (X + Y ) ∈ F =

⇒ X ∈ F and Y ∈ F. The faces of Sn

+ are in one-to-one correspondence with the

subspaces of Rn. i.e., F is a face of Sn

+ iff:

F = {X ∈ Sn

+ : L ⊆ N(X)}, for some subspace L,

relint(F) = {X ∈ Sn

+ : L = N(X)},

The faces of Sn

+ are uniquely characterized by their relative

interior.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-23
SLIDE 23

Characterizing face(A)

Let A ∈ Sn

+, face(A) denotes the minimal face containing A.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-24
SLIDE 24

Characterizing face(A)

Let A ∈ Sn

+, face(A) denotes the minimal face containing A.

Theorem (Barker ’73, Barker and Carlson ’75)

Let A ∈ Sn

+ of rank r and let A = [W U]

Λ W T UT

  • be the

spectral decomposition of A. Then face (A) = {X ∈ Sn

+ : XU = 0},

= {X ∈ Sn

+ : X = WYW T for some Y ∈ Sr +}.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-25
SLIDE 25

Characterizing face(A)

Let A ∈ Sn

+, face(A) denotes the minimal face containing A.

Theorem (Barker ’73, Barker and Carlson ’75)

Let A ∈ Sn

+ of rank r and let A = [W U]

Λ W T UT

  • be the

spectral decomposition of A. Then face (A) = {X ∈ Sn

+ : XU = 0},

= {X ∈ Sn

+ : X = WYW T for some Y ∈ Sr +}.

Sn

+ = face (I).

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-26
SLIDE 26

Characterizing face(A)

Let A ∈ Sn

+, face(A) denotes the minimal face containing A.

Theorem (Barker ’73, Barker and Carlson ’75)

Let A ∈ Sn

+ of rank r and let A = [W U]

Λ W T UT

  • be the

spectral decomposition of A. Then face (A) = {X ∈ Sn

+ : XU = 0},

= {X ∈ Sn

+ : X = WYW T for some Y ∈ Sr +}.

Sn

+ = face (I).

Note that if A, B ∈ Sn

+ and R(B) ⊂ R(A), then face (B) ⊂

face (A).

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-27
SLIDE 27

Theorem (A ’14)

Given symmetric matrices A1, . . . , Am and m-vector b, let F = {X 0: trace (AiX) = bi, for i = 1, . . . , m}. Assume that X ∗ ∈ F: rank X ∗ = r. Then there does not exist X ∈ F: rank X ≥ r + 1 if and only if face (F) = face (X ∗).

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-28
SLIDE 28

Theorem (A ’14)

Given symmetric matrices A1, . . . , Am and m-vector b, let F = {X 0: trace (AiX) = bi, for i = 1, . . . , m}. Assume that X ∗ ∈ F: rank X ∗ = r. Then there does not exist X ∈ F: rank X ≥ r + 1 if and only if face (F) = face (X ∗). Recall that face (F) ⊂ Sn

+ = face (I).

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-29
SLIDE 29

Theorem (A ’14)

Given symmetric matrices A1, . . . , Am and m-vector b, let F = {X 0: trace (AiX) = bi, for i = 1, . . . , m}. Assume that X ∗ ∈ F: rank X ∗ = r. Then there does not exist X ∈ F: rank X ≥ r + 1 if and only if face (F) = face (X ∗). Recall that face (F) ⊂ Sn

+ = face (I).

Borwein and Wolkowicz ’81 devised a facial reduction algorithm for finding face (F).

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-30
SLIDE 30

Theorem (A ’14)

Given symmetric matrices A1, . . . , Am and m-vector b, let F = {X 0: trace (AiX) = bi, for i = 1, . . . , m}. Assume that X ∗ ∈ F: rank X ∗ = r. Then there does not exist X ∈ F: rank X ≥ r + 1 if and only if face (F) = face (X ∗). Recall that face (F) ⊂ Sn

+ = face (I).

Borwein and Wolkowicz ’81 devised a facial reduction algorithm for finding face (F). At each step, a smaller dimensional face of Sn

+ is found. Thus,

this algorithm will find matrices U1, . . . , Uk such that: face (F) = face (Uk+1U T

k+1) ⊂ · · · ⊂ face (U1U T 1 ) ⊂ face (I).

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-31
SLIDE 31

Gram Matrices

A bar framework (G, p) is a simple graph G whose nodes are points p1, . . . , pn; and whose edges are line segments.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-32
SLIDE 32

Gram Matrices

A bar framework (G, p) is a simple graph G whose nodes are points p1, . . . , pn; and whose edges are line segments. (G, p) is r-dimensional if the dimension of the affine span of p1, . . . , pn is equal to r.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-33
SLIDE 33

Gram Matrices

A bar framework (G, p) is a simple graph G whose nodes are points p1, . . . , pn; and whose edges are line segments. (G, p) is r-dimensional if the dimension of the affine span of p1, . . . , pn is equal to r. The configuration p = (p1, . . . , pn) will be represented by the n × n Gram matrix B = (bij) where bij = (pi)Tpj.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-34
SLIDE 34

Gram Matrices

A bar framework (G, p) is a simple graph G whose nodes are points p1, . . . , pn; and whose edges are line segments. (G, p) is r-dimensional if the dimension of the affine span of p1, . . . , pn is equal to r. The configuration p = (p1, . . . , pn) will be represented by the n × n Gram matrix B = (bij) where bij = (pi)Tpj. B is PSD of rank r.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-35
SLIDE 35

Gram Matrices

A bar framework (G, p) is a simple graph G whose nodes are points p1, . . . , pn; and whose edges are line segments. (G, p) is r-dimensional if the dimension of the affine span of p1, . . . , pn is equal to r. The configuration p = (p1, . . . , pn) will be represented by the n × n Gram matrix B = (bij) where bij = (pi)Tpj. B is PSD of rank r. We fix the centroid of p at the origin, i.e., we require Be = 0 where e is the vector of all 1’s.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-36
SLIDE 36

Gram Matrices

A bar framework (G, p) is a simple graph G whose nodes are points p1, . . . , pn; and whose edges are line segments. (G, p) is r-dimensional if the dimension of the affine span of p1, . . . , pn is equal to r. The configuration p = (p1, . . . , pn) will be represented by the n × n Gram matrix B = (bij) where bij = (pi)Tpj. B is PSD of rank r. We fix the centroid of p at the origin, i.e., we require Be = 0 where e is the vector of all 1’s. B ⊂ face (VV T), where V Te = 0 and V has full column rank.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-37
SLIDE 37

Dimensional Rigidity

(G, p) and (G, p′) are equivalent if the Euclidean length of the corresponding edges in both frameworks are the same.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-38
SLIDE 38

Dimensional Rigidity

(G, p) and (G, p′) are equivalent if the Euclidean length of the corresponding edges in both frameworks are the same. An r-dimensional bar framework (G, p) is dimensionally rigid if there does not exist an s-dimensional bar framework (G, p′), that is equivalent to (G, p), for any s ≥ r + 1.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-39
SLIDE 39

Dimensional Rigidity

(G, p) and (G, p′) are equivalent if the Euclidean length of the corresponding edges in both frameworks are the same. An r-dimensional bar framework (G, p) is dimensionally rigid if there does not exist an s-dimensional bar framework (G, p′), that is equivalent to (G, p), for any s ≥ r + 1. 1 2 3 4 (G1, p1) 1 2 3 4 1 2 3 4

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-40
SLIDE 40

Dimensional Rigidity

(G, p) and (G, p′) are equivalent if the Euclidean length of the corresponding edges in both frameworks are the same. An r-dimensional bar framework (G, p) is dimensionally rigid if there does not exist an s-dimensional bar framework (G, p′), that is equivalent to (G, p), for any s ≥ r + 1. 1 2 3 4 (G1, p1) 1 2 3 4 1 2 3 4 1 2 3 (G2, p2) 1 2 3

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-41
SLIDE 41

Stress Matrices

Given r-dimensional framework (G, p) in Rr The matrix P =    (p1)T . . . (pn)T    is called the configuration matrix

  • f (G, p). Thus the Gram matrix B = PPT.
  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-42
SLIDE 42

Stress Matrices

Given r-dimensional framework (G, p) in Rr The matrix P =    (p1)T . . . (pn)T    is called the configuration matrix

  • f (G, p). Thus the Gram matrix B = PPT.

a symmetric matrix Ω is called a stress matrix of (G, p) if:

PTΩ = 0, eTΩ = 0, Ωij = 0 for all {i, j} ∈ E(G).

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-43
SLIDE 43

Stress Matrices

Given r-dimensional framework (G, p) in Rr The matrix P =    (p1)T . . . (pn)T    is called the configuration matrix

  • f (G, p). Thus the Gram matrix B = PPT.

a symmetric matrix Ω is called a stress matrix of (G, p) if:

PTΩ = 0, eTΩ = 0, Ωij = 0 for all {i, j} ∈ E(G).

a symmetric matrix Ω is quasi-stress matrix of (G, p) if:

PTΩP = 0, eTΩ = 0, Ωij = 0 for all {i, j} ∈ E(G).

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-44
SLIDE 44

Gale Matrices

A Gale matrix of an r-dimensional framework (G, p) in Rr is the matrix Z whose columns form a basis of the null space of : PT eT

  • .
  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-45
SLIDE 45

Gale Matrices

A Gale matrix of an r-dimensional framework (G, p) in Rr is the matrix Z whose columns form a basis of the null space of : PT eT

  • .

In Polytope theory, the rows of Z are called Gale transforms of p1, . . . , pn.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-46
SLIDE 46

Gale Matrices

A Gale matrix of an r-dimensional framework (G, p) in Rr is the matrix Z whose columns form a basis of the null space of : PT eT

  • .

In Polytope theory, the rows of Z are called Gale transforms of p1, . . . , pn. The Gale matrix Z encodes the affine dependency among the points p1, . . . , pn.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-47
SLIDE 47

Gale Matrix Z and Stress Matrix S

Theorem (Alfakih ’07)

Let Ω and Z be, respectively, a stress matrix and a Gale matrix of (G, p). Then Ω = ZΨZ T for some symmetric matrix Ψ. On the other hand, let Ψ′ be any symmetric matrix such that zi TΨ′zj = 0 for all {i, j} ∈ E(G), where zi is the ith row of Z. Then ZΨ′Z T is a stress matrix of (G, p).

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-48
SLIDE 48

Gale Matrix Z and Stress Matrix S

Theorem (Alfakih ’07)

Let Ω and Z be, respectively, a stress matrix and a Gale matrix of (G, p). Then Ω = ZΨZ T for some symmetric matrix Ψ. On the other hand, let Ψ′ be any symmetric matrix such that zi TΨ′zj = 0 for all {i, j} ∈ E(G), where zi is the ith row of Z. Then ZΨ′Z T is a stress matrix of (G, p). Thus results in terms of stress matrices an be equivalently stated in terms of Gale transform.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-49
SLIDE 49

Dimensional Rigidity

Theorem (A ’07)

Let (G, p) be an r-dimensional framework (G, p) on n nodes in Rr. Then (G, p) is dimensionally rigid if (G, p) admits a positive semidefinite stress matrix Ω with rank n − r − 1.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-50
SLIDE 50

Dimensional Rigidity

Theorem (A ’07)

Let (G, p) be an r-dimensional framework (G, p) on n nodes in Rr. Then (G, p) is dimensionally rigid if (G, p) admits a positive semidefinite stress matrix Ω with rank n − r − 1. This sufficient condition is not necessary.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-51
SLIDE 51

Dimensional Rigidity

Theorem (A ’07)

Let (G, p) be an r-dimensional framework (G, p) on n nodes in Rr. Then (G, p) is dimensionally rigid if (G, p) admits a positive semidefinite stress matrix Ω with rank n − r − 1. This sufficient condition is not necessary. A framework (G, p) is generic if all coordinates of p1, . . . , pn are algebraically independent

  • ver the integers.
  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-52
SLIDE 52

Dimensional Rigidity

Theorem (A ’07)

Let (G, p) be an r-dimensional framework (G, p) on n nodes in Rr. Then (G, p) is dimensionally rigid if (G, p) admits a positive semidefinite stress matrix Ω with rank n − r − 1. This sufficient condition is not necessary. A framework (G, p) is generic if all coordinates of p1, . . . , pn are algebraically independent

  • ver the integers.

Theorem (Gortler and Thurston ’10)

Let (G, p) be a generic r-dimensional framework (G, p) on n nodes in Rr. Assume that (G, p) is dimensionally rigid. Then (G, p) admits a positive semidefinite stress matrix Ω with rank n − r − 1.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-53
SLIDE 53

Dimensional Rigidity Cont’d

Let F = {B 0 : Be = 0, tr(BF ij) = ||pi − pj||2∀{i, j} ∈ E(G)}, where F ij = (ei − ej)(ei − ej)T and ei is the ith unit vector.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-54
SLIDE 54

Dimensional Rigidity Cont’d

Let F = {B 0 : Be = 0, tr(BF ij) = ||pi − pj||2∀{i, j} ∈ E(G)}, where F ij = (ei − ej)(ei − ej)T and ei is the ith unit vector. A framework (G, p′), with Gram matrix B′, is equivalent to (G, p) iff B′ ∈ F.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-55
SLIDE 55

Dimensional Rigidity Cont’d

Let F = {B 0 : Be = 0, tr(BF ij) = ||pi − pj||2∀{i, j} ∈ E(G)}, where F ij = (ei − ej)(ei − ej)T and ei is the ith unit vector. A framework (G, p′), with Gram matrix B′, is equivalent to (G, p) iff B′ ∈ F. An r-dimensional framework (G, p) is dimensionally rigid iff there does not exist B′ ∈ F : rank(B′) ≥ r + 1. (∗)

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-56
SLIDE 56

Dimensional Rigidity Cont’d

Let F = {B 0 : Be = 0, tr(BF ij) = ||pi − pj||2∀{i, j} ∈ E(G)}, where F ij = (ei − ej)(ei − ej)T and ei is the ith unit vector. A framework (G, p′), with Gram matrix B′, is equivalent to (G, p) iff B′ ∈ F. An r-dimensional framework (G, p) is dimensionally rigid iff there does not exist B′ ∈ F : rank(B′) ≥ r + 1. (∗) Condition (*) can be easily put in the form of Statement 1 in the new SD Farkas Lemma.

to Lemma

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-57
SLIDE 57

A Characterization of Dim Rigidity

This Theorem is a refined version of a result by Bob and Gortler ’14

Theorem

Let (G, p) be an r-dimensional bar framework on n vertices in Rr, r ≤ (n − 2). Then (G, p) is dimensionally rigid iff ∃ non-zero quasi-stress matrices Ω0, Ω1, . . . , Ωk, k ≤ n − r − 2, such that: Ω0 0, U T

1 Ω1U1 0, · · · , U T k ΩkUk 0,

rank Ω0 + rank (U T

1 Ω1U1) + · · · + rank (U T k ΩkUk) = n − r − 1,

PTΩ1ρ1 = 0, . . . , PTΩkρk = 0, where ρ1, U1, . . . , Uk and ξ1, . . . , ξk are full column rank matrices defined as follows: R(ρ1) = N(   Ω0 PT eT  ), R(ρi) = N(ρT

i Ωiρi), and

Ui = [P ρi] for i = 1, . . . , k and ρi+1 = ρiξi for i = 1, . . . , k − 1.

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18

slide-58
SLIDE 58

Thank You! Happy Birthday Walter

  • A. Y. Alfakih

(University of Windsor, ON Canada) Semidefinite Farkas’ Lemma and Dimensional Rigidity of Bar Frameworks Workshop on Making Models: Stimulating / 18