On Hadamards Maximal Determinant Problem Judy-anne Osborn MSI, ANU - - PowerPoint PPT Presentation

on hadamard s maximal determinant problem
SMART_READER_LITE
LIVE PREVIEW

On Hadamards Maximal Determinant Problem Judy-anne Osborn MSI, ANU - - PowerPoint PPT Presentation

On Hadamards Maximal Determinant Problem Judy-anne Osborn MSI, ANU April 2009 Judy-anne Osborn MSI, ANU On Hadamards Maximal Determinant Problem 0 1 0 1 0 0 0 1 1 . . . . . . . . . . . . . . .


slide-1
SLIDE 1

On Hadamard’s Maximal Determinant Problem

Judy-anne Osborn MSI, ANU April 2009

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-2
SLIDE 2

              . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .              

1 1 1 1 1 1 1

m m

max det =?

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-3
SLIDE 3

A Naive Computer Search

Order max det Time 1 1 fast 2 1 fast 3 2 fast 4 3 fast 5 5 fast 6 9

  • rder of days

7 32

  • rder of years

8 56

  • rder of the age of the Universe

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-4
SLIDE 4

As well by hand? I found nested Max Dets ...

        1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1         }1

1

}}

2} 3} 5} 9

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-5
SLIDE 5

As well by hand? I found nested Max Dets ...

        1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1         }1

1

}}

2} 3} 5} 9

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-6
SLIDE 6

As well by hand? I found nested Max Dets ...

        1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1         }1

1

}}

2} 3} 5} 9

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-7
SLIDE 7

As well by hand? I found nested Max Dets ...

        1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1         }1

1

}}

2} 3} 5} 9

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-8
SLIDE 8

As well by hand? I found nested Max Dets ...

        1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1         }1

1

}}

2} 3} 5} 9

◮ But no further!

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-9
SLIDE 9

The problem turns out to be famous

◮ Hadamard’s Maximal Determinant Problem was posed in 1893

Jacques Hadamard

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-10
SLIDE 10

◮ A little selected history on this century-old question ...

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-11
SLIDE 11

Observe:

◮ 2D:

1

  • 1
  • 1
  • r

1 1

  • 1

1

  • 1

1 1

  • Judy-anne Osborn MSI, ANU

On Hadamard’s Maximal Determinant Problem

slide-12
SLIDE 12

Geometry: max |det| = max (hyper-)Volume

◮ 3D:   1 1 1 1 1 1  

  1 1     1 1     1 1     1 1  

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-13
SLIDE 13

An equivalent problem: {+1, −1} matrices

1 1 1 1 1 1

( (

2 2 2 2 2 2

− − − − − −

( (

1 1 1 1 2 2 2 2 2 2

− − − − − −

( (

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

− − − − − −

( (

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

− − − − − −

( (

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

− − − − − −

( (

− −

×(−2) border add row 1 column ops row ops m × m matrix Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-14
SLIDE 14

An equivalent problem: {+1, −1} matrices

1 1 1 1 1 1

( (

2 2 2 2 2 2

− − − − − −

( (

1 1 1 1 2 2 2 2 2 2

− − − − − −

( (

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

− − − − − −

( (

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

− − − − − −

( (

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

− − − − − −

( (

− −

×(−2) border add row 1 column ops row ops

|detnew| = 2m|detold|

m × m matrix Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-15
SLIDE 15

Volume interpretation ⇒ upper bound on |max det|

max

  • ±1

· · · ±1 . . . ... . . . ±1 · · · ±1

  • }

}

n n

◮ What is the upper bound?

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-16
SLIDE 16

Volume interpretation ⇒ upper bound on |max det|

max

  • ±1

· · · ±1 . . . ... . . . ±1 · · · ±1

  • }

}

n n

◮ What is the upper bound?

  • (±1)2 + · · · + (±1)2

n = nn/2

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-17
SLIDE 17

Volume interpretation ⇒ upper bound on |max det|

max

  • ±1

· · · ±1 . . . ... . . . ±1 · · · ±1

  • }

}

n n

◮ What is the upper bound?

  • (±1)2 + · · · + (±1)2

n = nn/2

◮ Why?

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-18
SLIDE 18

Volume interpretation ⇒ upper bound on |max det|

max

  • ±1

· · · ±1 . . . ... . . . ±1 · · · ±1

  • }

}

n n

◮ What is the upper bound?

  • (±1)2 + · · · + (±1)2

n = nn/2

◮ Why? (Columns/rows orthogonal)

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-19
SLIDE 19

When is the bound tight?

◮ Tight when {+1, −1} square matrix H of order n satisfies

HHT = nI

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-20
SLIDE 20

When is the bound tight?

◮ Tight when {+1, −1} square matrix H of order n satisfies

HHT = nI

◮ H is called a Hadamard Matrix.

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-21
SLIDE 21

When is the bound tight?

◮ Tight when {+1, −1} square matrix H of order n satisfies

HHT = nI

◮ H is called a Hadamard Matrix. ◮ A necessary condition on existence of H is:

n = 1, 2

  • r

n ≡ 0 (mod 4)

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-22
SLIDE 22

When is the bound tight?

◮ Tight when {+1, −1} square matrix H of order n satisfies

HHT = nI

◮ H is called a Hadamard Matrix. ◮ A necessary condition on existence of H is:

n = 1, 2

  • r

n ≡ 0 (mod 4)

◮ Hadamard Conjecture (Paley, 1933):

this is also sufficient.

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-23
SLIDE 23

Evidence for Hadamard Conjecture

◮ Many constructions for infinite families, including

◮ Sylvester, ∀ 2r ◮ First Paley, using finite fields, ∀ pr + 1, p prime ◮ Second Paley, using finite fields, ∀ 2pr + 2, p prime Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-24
SLIDE 24

Evidence for Hadamard Conjecture

◮ Many constructions for infinite families, including

◮ Sylvester, ∀ 2r ◮ First Paley, using finite fields, ∀ pr + 1, p prime ◮ Second Paley, using finite fields, ∀ 2pr + 2, p prime

◮ Other ‘constructions’ and ‘ad hoc’ examples due to people

including

◮ Williamson ◮ Jenny Seberry Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-25
SLIDE 25

Evidence for Hadamard Conjecture

◮ Many constructions for infinite families, including

◮ Sylvester, ∀ 2r ◮ First Paley, using finite fields, ∀ pr + 1, p prime ◮ Second Paley, using finite fields, ∀ 2pr + 2, p prime

◮ Other ‘constructions’ and ‘ad hoc’ examples due to people

including

◮ Williamson ◮ Jenny Seberry

◮ Smallest n ≡ 0 (mod 4) currently undecided:

n = 668.

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-26
SLIDE 26

More evidence

Order Number of inequivalent Hadamard matrices – see Sloan’s sequence A007299 1 1 2 1 4 1 8 1 12 1 16 5 20 3 24 60 28 487 32 ≥ 3 578 006 36 ≥ 18 292 717

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-27
SLIDE 27

Max Dets for non-Hadamard orders?

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-28
SLIDE 28

Max Dets for non-Hadamard orders?

n ≡ 1 1 5 9 13 17 21 25

|max det| 2n−1

1 3 × 11 7 × 23 15 × 35 20 × 47 29 × 59 42 × 611 n ≡ 2 2 6 10 14 18

|max det| 2n−1

1 5 × 11 18 × 23 39 × 35 68 × 47 n ≡ 3 3 7 11 15

|max det| 2n−1

1 9 × 11 40 × 23 105 × 35

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-29
SLIDE 29

Max Dets for non-Hadamard orders?

n ≡ 1 1 5 9 13 17 21 25

|max det| 2n−1

1 3 × 11 7 × 23 15 × 35 20 × 47 29 × 59 42 × 611 n ≡ 2 2 6 10 14 18

|max det| 2n−1

1 5 × 11 18 × 23 39 × 35 68 × 47 n ≡ 3 3 7 11 15

|max det| 2n−1

1 9 × 11 40 × 23 105 × 35

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-30
SLIDE 30

Max Dets for non-Hadamard orders?

n ≡ 1 1 5 9 13 17 21 25

|max det| 2n−1

1 3 × 11 7 × 23 15 × 35 20 × 47 29 × 59 42 × 611 n ≡ 2 2 6 10 14 18

|max det| 2n−1

1 5 × 11 18 × 23 39 × 35 68 × 47 n ≡ 3 3 7 11 15

|max det| 2n−1

1 9 × 11 40 × 23 105 × 35

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-31
SLIDE 31

Tighter upper bounds?

◮ The Hadamard bound of nn/2 holds for all orders but is never

tight for n ≡ 0 (mod 4) when n > 2.

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-32
SLIDE 32

Tighter upper bounds?

◮ The Hadamard bound of nn/2 holds for all orders but is never

tight for n ≡ 0 (mod 4) when n > 2.

◮ Better upper bounds for non-Hadamard orders were proved by

◮ Barba in 1933 ◮ Ehlich in 1962, 64 ◮ Wojtas in 1964 ◮ Cohn (proved a new bound tightness result) in 2000 Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-33
SLIDE 33

The best known upper bounds are:

◮ The Barba-Ehlich bound holds for n ≡ 1 (mod 4):

√ 2n − 1(n − 1)(n−1)/2

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-34
SLIDE 34

The best known upper bounds are:

◮ The Barba-Ehlich bound holds for n ≡ 1 (mod 4):

√ 2n − 1(n − 1)(n−1)/2

◮ The Ehlich-Wojtas bound holds for n ≡ 2 (mod 4):

(2n − 2)(n − 2)(n/2)−1

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-35
SLIDE 35

The best known upper bounds are:

◮ The Ehlich bound holds for n ≡ 3 (mod 4):

(n−3)

n−s 2 (n−3+4r) u 2 (n+1+4r) v 2

  • 1 −

ur n − 3 + 4r − v(r + 1) n + 1 + 4r where s = 3 for n = 3, s = 5 for n = 7, s = 5 or 6 for n = 11, s = 6 for n = 15, 19, ..., 59, and s = 7 for n ≥ 63, r = ⌊ n

s ⌋,

n = rs + v and u = s − v.

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-36
SLIDE 36

Percentages of bounds met: summary from Will Orrick’s:

◮ www.indiana.edu/∼maxdet

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-37
SLIDE 37

The idea behind the bounds

◮ Let R be a maximal determinant square ±1 matrix of order n.

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-38
SLIDE 38

The idea behind the bounds

◮ Let R be a maximal determinant square ±1 matrix of order n. ◮ Consider ‘gram matrix’

G := RRT

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-39
SLIDE 39

The idea behind the bounds

◮ Let R be a maximal determinant square ±1 matrix of order n. ◮ Consider ‘gram matrix’

G := RRT

  • 1. G has all n’s on the diagonal

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-40
SLIDE 40

The idea behind the bounds

◮ Let R be a maximal determinant square ±1 matrix of order n. ◮ Consider ‘gram matrix’

G := RRT

  • 1. G has all n’s on the diagonal
  • 2. G is positive definite ⇒ off-diagonal entries have size < n

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-41
SLIDE 41

The idea behind the bounds

◮ Let R be a maximal determinant square ±1 matrix of order n. ◮ Consider ‘gram matrix’

G := RRT

  • 1. G has all n’s on the diagonal
  • 2. G is positive definite ⇒ off-diagonal entries have size < n
  • 3. G has all off-diagonal entries ≡ n (mod 2)

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-42
SLIDE 42

The idea behind the bounds

◮ Let R be a maximal determinant square ±1 matrix of order n. ◮ Consider ‘gram matrix’

G := RRT

  • 1. G has all n’s on the diagonal
  • 2. G is positive definite ⇒ off-diagonal entries have size < n
  • 3. G has all off-diagonal entries ≡ n (mod 2)
  • 4. G has all off-diagonal entries ≡ n (mod 4) with R normalized

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-43
SLIDE 43

The idea behind the bounds

◮ Let R be a maximal determinant square ±1 matrix of order n. ◮ Consider ‘gram matrix’

G := RRT

  • 1. G has all n’s on the diagonal
  • 2. G is positive definite ⇒ off-diagonal entries have size < n
  • 3. G has all off-diagonal entries ≡ n (mod 2)
  • 4. G has all off-diagonal entries ≡ n (mod 4) with R normalized
  • 5. G is symmetric

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-44
SLIDE 44

The idea behind the bounds

◮ Let R be a maximal determinant square ±1 matrix of order n. ◮ Consider ‘gram matrix’

G := RRT

  • 1. G has all n’s on the diagonal
  • 2. G is positive definite ⇒ off-diagonal entries have size < n
  • 3. G has all off-diagonal entries ≡ n (mod 2)
  • 4. G has all off-diagonal entries ≡ n (mod 4) with R normalized
  • 5. G is symmetric

◮ Let Gn be the set of all gram matrices, G, of order n

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-45
SLIDE 45

The idea behind the bounds, continued

◮ Let Gn be the set of matrices for which properties 1 – 5 hold.

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-46
SLIDE 46

The idea behind the bounds, continued

◮ Let Gn be the set of matrices for which properties 1 – 5 hold. ◮ Then

Gn ⊇ Gn

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-47
SLIDE 47

The idea behind the bounds, continued

◮ Let Gn be the set of matrices for which properties 1 – 5 hold. ◮ Then

Gn ⊇ Gn

◮ Hence

  • max det
  • Gn
  • ≥ |max det (Gn)|

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-48
SLIDE 48

Case: n ≡ 1 (mod 4)

◮ The matrix which was proven by Barba and Ehlich to have

largest determinant in Gn is    n 1 ... 1 n   

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-49
SLIDE 49

Case: n ≡ 2 (mod 4)

◮ The matrix which was proven by Ehlich and Wojtas to have

largest determinant in Gn is F F

  • , where F =

   n 2 ... 2 n   

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-50
SLIDE 50

Case: n ≡ 3 (mod 4)

◮ We expect a matrix with largest determinant in Gn to be:

   n −1 ... −1 n   

◮ In general, this is wrong ◮ Ehlich proved: the correct best determinant matrix in Gn is a

block form with off-diagonal entries from the set {−1, +3}.

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-51
SLIDE 51

Structure in the n ≡ 1 (mod 4) case

◮ What is a necessary condition for tightness of: Barba-Ehlich:

√ 2n − 1(n − 1)(n−1)/2?

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-52
SLIDE 52

Structure in the n ≡ 1 (mod 4) case

◮ What is a necessary condition for tightness of: Barba-Ehlich:

√ 2n − 1(n − 1)(n−1)/2?

◮ Lemma: The number 2n − 1 is a perfect square iff ∃q ∈ N

such that n = q2 + (q + 1)2

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-53
SLIDE 53

In the literature, a conjecture on tightness:

◮ Conjecture:

The Barba-Ehlich bound is tight whenever n is a sum of two consecutive squares: n = q2 + (q + 1)2

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-54
SLIDE 54

In the literature, a conjecture on tightness:

◮ Conjecture:

The Barba-Ehlich bound is tight whenever n is a sum of two consecutive squares: n = q2 + (q + 1)2

◮ Evidence:

True for q = 2, 4

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-55
SLIDE 55

In the literature, a conjecture on tightness:

◮ Conjecture:

The Barba-Ehlich bound is tight whenever n is a sum of two consecutive squares: n = q2 + (q + 1)2

◮ Evidence:

True for q = 2, 4 and when q = pr for p an odd prime – proved by Brouwer’s Construction.

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-56
SLIDE 56

What chance an exact maxdet formula ∀ n ≡ 1 (mod 4)?

n ≡ 1 1 5 9 13 17 21 25 29

|max det| 2n−1

1 3 × 11 7 × 23 15 × 35 20 × 47 29 × 59 42 × 611 48 × 713 +4 +8 +5 +9 +13 +6 Guess

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-57
SLIDE 57

What chance an exact maxdet formula ∀ n ≡ 1 (mod 4)?

n ≡ 1 1 5 9 13 17 21 25 29

|max det| 2n−1

1 3 × 11 7 × 23 15 × 35 20 × 47 29 × 59 42 × 611 48 × 713 +4 +8 +5 +9 +13 +6 Guess ◮ A guess/conjecture due to Will Orrick is that |max det| 2n−1

for n = 4k + 1 is always divisible by k2k−1,

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-58
SLIDE 58

What chance an exact maxdet formula ∀ n ≡ 1 (mod 4)?

n ≡ 1 1 5 9 13 17 21 25 29

|max det| 2n−1

1 3 × 11 7 × 23 15 × 35 20 × 47 29 × 59 42 × 611 48 × 713 +4 +8 +5 +9 +13 +6 Guess ◮ A guess/conjecture due to Will Orrick is that |max det| 2n−1

for n = 4k + 1 is always divisible by k2k−1, with coefficients growing quadratically between n’s for which n = q2 + (q + 1)2.

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-59
SLIDE 59

What can we hope to compute?

◮ The ideas of Ehlich and Wojtas were reused by Chadjipantelis,

Kounias and Moissiadis in the 1980’s to find max det matrices for n = 17 and n = 21.

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-60
SLIDE 60

What can we hope to compute?

◮ The ideas of Ehlich and Wojtas were reused by Chadjipantelis,

Kounias and Moissiadis in the 1980’s to find max det matrices for n = 17 and n = 21.

◮ Will Orrick used similar ideas in the 2000’s to prove maximal

an n = 15 matrix that had previously been found by Cohn; as well as filling in some gaps in CKM’s published proofs.

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-61
SLIDE 61

What can we hope to compute?

◮ The ideas of Ehlich and Wojtas were reused by Chadjipantelis,

Kounias and Moissiadis in the 1980’s to find max det matrices for n = 17 and n = 21.

◮ Will Orrick used similar ideas in the 2000’s to prove maximal

an n = 15 matrix that had previously been found by Cohn; as well as filling in some gaps in CKM’s published proofs.

◮ Are we within reach of n = 29?

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-62
SLIDE 62

Basic Idea

Two steps:

  • 1. Find candidate gram matrices in Gn.
  • 2. Check if candidates decompose in the form

RRT.

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-63
SLIDE 63

Essentials of Step 1.

◮ There exist theorems which bound the determinants of

candidate gram matrices in terms of their sub-matrices.

◮ So we can set a target determinant and build candidates:

n n n n n pruning too-small sub-matrices as we go.

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-64
SLIDE 64

Computational considerations for Step 1.

◮ We must have efficient ways to calculate determinants! ◮ ‘Rank-One Update’ Theorems:

O(size3) → O(size2), at the expense of some book-keeping.

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-65
SLIDE 65

Computational considerations for Step 1.

◮ Need to prune equivalent gram matrices:

7 1 1 3 1 1 1 1 7 3 1 1 1 1 1 3 7 1 1 1 1 3 1 1 7 1 1 1 1 1 1 1 7 3 1 1 1 1 1 3 7 1 1 1 1 1 1 1 7

− − − − − − − − − − − − − − − − − − − − − − − − − − − − − − − − − − − −

( (

7 3 1 1 1 1 1 3 7 1 1 1 1 1 1 1 7 3 1 1 1 1 1 3 7 1 1 1 1 1 1 1 7 3 1 1 1 1 1 3 7 1 1 1 1 1 1 1 7

− − − − − − − − − − − − − − − − − − − − − − − − − − − − − − − − − − − −

( (

eg. under simultaneous row and column permutation

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-66
SLIDE 66

Computational considerations for Step 1.

◮ So need a (partial ordering) which

◮ prunes heavily enough, and ◮ is practical to compute on-the-fly

?

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-67
SLIDE 67

Computational considerations for Step 1.

◮ So need a (partial ordering) which

◮ prunes heavily enough, and ◮ is practical to compute on-the-fly

?

◮ Make sure we don’t miss any valid candidates!

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-68
SLIDE 68

Essentials of Step 2.

◮ Consider gram candidates A and B whose determinants agree.

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-69
SLIDE 69

Essentials of Step 2.

◮ Consider gram candidates A and B whose determinants agree. ◮ These are candidates for

RRT and RTR

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-70
SLIDE 70

Essentials of Step 2.

◮ Consider gram candidates A and B whose determinants agree. ◮ These are candidates for

RRT and RTR

◮ Implement two kinds of constraints:

◮ Linear, ◮ Quadratic – which use A and B simultaneously Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-71
SLIDE 71

Essentials of Step 2. – the linear constraints

◮ Let A = (aij) and B = (bij).

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-72
SLIDE 72

Essentials of Step 2. – the linear constraints

◮ Let A = (aij) and B = (bij). ◮ Assume

R =    r0 . . . rn−1    =

  • c0

· · · cn−1

  • Judy-anne Osborn MSI, ANU

On Hadamard’s Maximal Determinant Problem

slide-73
SLIDE 73

Essentials of Step 2. – the linear constraints

◮ Let A = (aij) and B = (bij). ◮ Assume

R =    r0 . . . rn−1    =

  • c0

· · · cn−1

  • ◮ Then

aij = ri.rj and bij = ci.cj

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-74
SLIDE 74

Essentials of Step 2. – the linear constraints

◮ Let A = (aij) and B = (bij). ◮ Assume

R =    r0 . . . rn−1    =

  • c0

· · · cn−1

  • ◮ Then

aij = ri.rj and bij = ci.cj

◮ We need to implement an ordering to prune duplicates

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-75
SLIDE 75

Considerations for the linear constraints of Step 2.

A =      17 −3 1 · · · −3 17 1 · · · 1 1 17 · · · . . . . . . . . . ...      eg. If we have already built r0 = (1, −, −, 1, 1, −, −, −, −, −, 1, 1, 1, 1, 1, 1, 1) r1 = (−, −, 1, −, −, −, −, 1, 1, 1, −, −, −, 1, 1, 1, 1) and then r2 breaks into blocks: r2 = (a; b; c; d, e; f, g; h, i, j; k, l, m; n, o, p, q)

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-76
SLIDE 76

Considerations for the linear constraints of Step 2.

A =      17 −3 1 · · · −3 17 1 · · · 1 1 17 · · · . . . . . . . . . ...      eg. If we have already built r0 = (1, −, −, 1, 1, −, −, −, −, −, 1, 1, 1, 1, 1, 1, 1) r1 = (−, −, 1, −, −, −, −, 1, 1, 1, −, −, −, 1, 1, 1, 1) and then r2 breaks into blocks: r2 = (a; b; c; d, e; f, g; h, i, j; k, l, m; n, o, p, q)

◮ Adding more rows is a process of successive block refinement

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-77
SLIDE 77

Considerations for the linear constraints of Step 2.

◮ Because we work with both rows and columns, we need a way

  • f refining row-blocks and column blocks simultaneously!

1 − − 1 1 − − − − 1 − − − − − 1 1 1 − − 1 1

?

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-78
SLIDE 78

Essentials of Step 2. – the quadratic constraints

◮ To derive the quadratic constraints, write key ingredients in

block form:

R =

  • 1

yT x R′

  • , RT =
  • 1

xT y R′T

  • , A =
  • n

aT a A′

  • , B =
  • n

bT b B′

  • Judy-anne Osborn MSI, ANU

On Hadamard’s Maximal Determinant Problem

slide-79
SLIDE 79

Essentials of Step 2. – the quadratic constraints

◮ To derive the quadratic constraints, write key ingredients in

block form:

R =

  • 1

yT x R′

  • , RT =
  • 1

xT y R′T

  • , A =
  • n

aT a A′

  • , B =
  • n

bT b B′

  • ◮ This allow several quadratic constraints to be found, eg.

det(A′ − xxT) = a perfect square = det(B′ − yyT)

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-80
SLIDE 80

Computational Considerations for Step 2.

◮ We need to decide in which order to implement the various

quadratic and linear constraints.

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-81
SLIDE 81

Computational Considerations for Step 2.

◮ We need to decide in which order to implement the various

quadratic and linear constraints.

◮ How to decide?

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem

slide-82
SLIDE 82

THE END

Judy-anne Osborn MSI, ANU On Hadamard’s Maximal Determinant Problem