Uniform Sampling of Subshifts of Finite Type Ir` ene Marcovici - - PowerPoint PPT Presentation

uniform sampling of subshifts of finite type
SMART_READER_LITE
LIVE PREVIEW

Uniform Sampling of Subshifts of Finite Type Ir` ene Marcovici - - PowerPoint PPT Presentation

Uniform Sampling of Subshifts of Finite Type Ir` ene Marcovici With the support of the European INTEGER project Institut Elie Cartan de Lorraine, Universit e de Lorraine, Nancy, France AofA15, Strobl Monday 8 June 2015 Ir` ene


slide-1
SLIDE 1

Uniform Sampling of Subshifts of Finite Type

Ir` ene Marcovici

With the support of the European INTEGER project

Institut ´ Elie Cartan de Lorraine, Universit´ e de Lorraine, Nancy, France

AofA’15, Strobl Monday 8 June 2015

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-2
SLIDE 2

Subshifts of finite type

Let us color the vertices of the lattice Zd using a finite number

  • f colors, with the constraint that some pairs of colors are not

allowed for adjacent sites.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-3
SLIDE 3

Subshifts of finite type

Let us color the vertices of the lattice Zd using a finite number

  • f colors, with the constraint that some pairs of colors are not

allowed for adjacent sites. Questions What do typical configurations look like?

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-4
SLIDE 4

Subshifts of finite type

Let us color the vertices of the lattice Zd using a finite number

  • f colors, with the constraint that some pairs of colors are not

allowed for adjacent sites. Questions What do typical configurations look like? How to define “the uniform distribution” on infinite allowed colorings?

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-5
SLIDE 5

Subshifts of finite type

Let us color the vertices of the lattice Zd using a finite number

  • f colors, with the constraint that some pairs of colors are not

allowed for adjacent sites. Questions What do typical configurations look like? How to define “the uniform distribution” on infinite allowed colorings? How to sample configurations uniformly?

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-6
SLIDE 6

Subshifts of finite type

Let us color the vertices of the lattice Zd using a finite number

  • f colors, with the constraint that some pairs of colors are not

allowed for adjacent sites. Questions What do typical configurations look like? How to define “the uniform distribution” on infinite allowed colorings? How to sample configurations uniformly? Given the set A of colors and the (finite) list of constraints, the set Σ of allowed configurations is called a: subshift of finite type (SFT).

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-7
SLIDE 7

Subshifts of finite type

Let us color the vertices of the lattice Zd using a finite number

  • f colors, with the constraint that some pairs of colors are not

allowed for adjacent sites. Questions What do typical configurations look like? How to define “the uniform distribution” on infinite allowed colorings? How to sample configurations uniformly? Given the set A of colors and the (finite) list of constraints, the set Σ of allowed configurations is called a: subshift of finite type (SFT). It is a subset of AZd, which is shift-invariant.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-8
SLIDE 8

Example

Fibonacci / golden mean / hard-core (or hard-square) subshift Set of configurations without two consecutive black squares.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-9
SLIDE 9

Example

Fibonacci / golden mean / hard-core (or hard-square) subshift Set of configurations without two consecutive black squares. A two-dimensional configuration

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-10
SLIDE 10

Example

Fibonacci / golden mean / hard-core (or hard-square) subshift Set of configurations without two consecutive black squares. A two-dimensional configuration A one-dimensional configuration

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-11
SLIDE 11

Example

Fibonacci / golden mean / hard-core (or hard-square) subshift Set of configurations without two consecutive black squares. A two-dimensional configuration A one-dimensional configuration

1

Graph of allowed transitions in one-dimension

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-12
SLIDE 12

One-dimensional SFT

Let A be an alphabet with n letters, and let A ∈ Mn({0, 1}).

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-13
SLIDE 13

One-dimensional SFT

Let A be an alphabet with n letters, and let A ∈ Mn({0, 1}). One-dimensional subshift of finite type The subshift of finite type associated to A is the set ΣA of words w ∈ AZ such that if Ai,j = 0, w does not contain the pattern ij. Ai,j = 1 if ij is an allowed pattern, 0 if ij is a forbidden pattern. ΣA = {w ∈ AZ; ∀k ∈ Z, Awk,wk+1 = 1}.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-14
SLIDE 14

One-dimensional SFT

Let A be an alphabet with n letters, and let A ∈ Mn({0, 1}). One-dimensional subshift of finite type The subshift of finite type associated to A is the set ΣA of words w ∈ AZ such that if Ai,j = 0, w does not contain the pattern ij. Ai,j = 1 if ij is an allowed pattern, 0 if ij is a forbidden pattern. ΣA = {w ∈ AZ; ∀k ∈ Z, Awk,wk+1 = 1}. In what follows, we assume that the matrix A is irreducible and aperiodic.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-15
SLIDE 15

The Parry measure

From Perron-Frobenius theory, the matrix A has a real eigenvalue λ > 0 such that |µ| ≤ λ for any other eigenvalue µ. Furthermore, there is a unique choice of r1, . . . , rn ≥ 0 such that n

i=1 ri = 1 and

A    r1 . . . rn    = λ    r1 . . . rn    .

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-16
SLIDE 16

The Parry measure

From Perron-Frobenius theory, the matrix A has a real eigenvalue λ > 0 such that |µ| ≤ λ for any other eigenvalue µ. Furthermore, there is a unique choice of r1, . . . , rn ≥ 0 such that n

i=1 ri = 1 and

A    r1 . . . rn    = λ    r1 . . . rn    . Definition of the Parry measure The Parry measure is the (shift-invariant) Markov measure π on AZ of transition matrix P defined, for any i, j ∈ A, by Pi,j = Ai,j rj λri .

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-17
SLIDE 17

The Parry measure

From Perron-Frobenius theory, the matrix A has a real eigenvalue λ > 0 such that |µ| ≤ λ for any other eigenvalue µ. Furthermore, there is a unique choice of r1, . . . , rn ≥ 0 such that n

i=1 ri = 1 and

A    r1 . . . rn    = λ    r1 . . . rn    . Definition of the Parry measure The Parry measure is the (shift-invariant) Markov measure π on AZ of transition matrix P defined, for any i, j ∈ A, by Pi,j = Ai,j rj λri . For a word a1 . . . ak ∈ Ak, π(a1 . . . ak) = π(a1)Pa1,a2 . . . Pak−1,ak.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-18
SLIDE 18

The Parry measure

The Parry measure π is “the uniform distribution” on ΣA.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-19
SLIDE 19

The Parry measure

The Parry measure π is “the uniform distribution” on ΣA. Proposition Let µk be the uniform measure on allowed patterns of length 2k + 1, centered at position 0. The sequence µk converges (weakly) to π on AZ .

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-20
SLIDE 20

Markov-uniform property

Proposition The Parry measure is Markov-uniform: for given k ≥ 1 and a, b ∈ A, the value π(awb) does not depend on the word w ∈ Ak such that awb is allowed.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-21
SLIDE 21

Markov-uniform property

Proposition The Parry measure is Markov-uniform: for given k ≥ 1 and a, b ∈ A, the value π(awb) does not depend on the word w ∈ Ak such that awb is allowed.

  • Proof. By definition, Pi,j = Ai,j

rj λri . If awb is allowed, then:

π(awb) = πaPa,w1Pw1,w2 . . . Pwk−1,wkPwk,b

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-22
SLIDE 22

Markov-uniform property

Proposition The Parry measure is Markov-uniform: for given k ≥ 1 and a, b ∈ A, the value π(awb) does not depend on the word w ∈ Ak such that awb is allowed.

  • Proof. By definition, Pi,j = Ai,j

rj λri . If awb is allowed, then:

π(awb) = πaPa,w1Pw1,w2 . . . Pwk−1,wkPwk,b = πa rw1 λra rw2 λrw1 · · · rwk λrwk−1 rb λrwk

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-23
SLIDE 23

Markov-uniform property

Proposition The Parry measure is Markov-uniform: for given k ≥ 1 and a, b ∈ A, the value π(awb) does not depend on the word w ∈ Ak such that awb is allowed.

  • Proof. By definition, Pi,j = Ai,j

rj λri . If awb is allowed, then:

π(awb) = πaPa,w1Pw1,w2 . . . Pwk−1,wkPwk,b = πa rw1 λra rw2 λrw1 · · · rwk λrwk−1 rb λrwk = πarb λk+1ra .

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-24
SLIDE 24

Measure of maximal entropy

Theorem Let MΣA be the set of translation invariant measures on the SFT ΣA, and let π ∈ MΣA. The following properties are equivalent. (i) π is the Parry measure associated to ΣA, (ii) π is a Markov-uniform measure on ΣA, (iii) π is the measure of maximal entropy of ΣA, (iv) the entropy of π is equal to the topological entropy h(ΣA).

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-25
SLIDE 25

Measure of maximal entropy

Theorem Let MΣA be the set of translation invariant measures on the SFT ΣA, and let π ∈ MΣA. The following properties are equivalent. (i) π is the Parry measure associated to ΣA, (ii) π is a Markov-uniform measure on ΣA, (iii) π is the measure of maximal entropy of ΣA, (iv) the entropy of π is equal to the topological entropy h(ΣA). On Zd, the equivalence between (ii), (iii), and (iv) can be extended to a class of multi-dimensional SFT.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-26
SLIDE 26

Measure of maximal entropy

Theorem Let MΣA be the set of translation invariant measures on the SFT ΣA, and let π ∈ MΣA. The following properties are equivalent. (i) π is the Parry measure associated to ΣA, (ii) π is a Markov-uniform measure on ΣA, (iii) π is the measure of maximal entropy of ΣA, (iv) the entropy of π is equal to the topological entropy h(ΣA). On Zd, the equivalence between (ii), (iii), and (iv) can be extended to a class of multi-dimensional SFT. But there can be several measures satisfying these properties...

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-27
SLIDE 27

One-dimensional Fibonacci SFT

Let A = {0, 1}. The one-dimensional Fibonacci SFT is the set

  • f words that do not contain two consecutive 1’s. It is given by:

A = 1 1 1

  • .

Its topological entropy is equal to log ϕ, where ϕ = 1+

√ 5 2

. The Parry measure is the Markov measure given by 1

1 ϕ 1 ϕ2

1 with π0 =

ϕ2 1+ϕ2 and π1 = 1 1+ϕ2 .

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-28
SLIDE 28

First rejection sampling

The Parry measure of the Fibonacci SFT can be generated by: choosing independently to write a 0 with probability r0 = 1

ϕ

and a 1 with probability r1 =

1 ϕ2 ,

rejecting the 1’s creating forbidden patterns.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-29
SLIDE 29

First rejection sampling

The Parry measure of the Fibonacci SFT can be generated by: choosing independently to write a 0 with probability r0 = 1

ϕ

and a 1 with probability r1 =

1 ϕ2 ,

rejecting the 1’s creating forbidden patterns. Lemma For any SFT, the Parry measure can be generated by independent draws of letters with probability (ri)i∈A, with reject of a letter if it creates a forbidden pattern.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-30
SLIDE 30

First rejection sampling

The Parry measure of the Fibonacci SFT can be generated by: choosing independently to write a 0 with probability r0 = 1

ϕ

and a 1 with probability r1 =

1 ϕ2 ,

rejecting the 1’s creating forbidden patterns. Lemma For any SFT, the Parry measure can be generated by independent draws of letters with probability (ri)i∈A, with reject of a letter if it creates a forbidden pattern. Proof. Pi,j = Ai,j rj λri = Ai,j rj

  • k∈A Ai,krk

= Ai,j rj

  • k∈S(i) rk

. where S(i) = {j ∈ A; Ai,j = 1} is the set of succesors of i.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-31
SLIDE 31

Second rejection sampling

The Parry measure of the Fibonacci SFT can be generated by: choosing independently to write a 0 with probability ˜ r0 =

1 ϕ2

and a 1 with probability ˜ r1 = 1

ϕ,

deleting pairs of consecutive 1’s.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-32
SLIDE 32

Second rejection sampling

The Parry measure of the Fibonacci SFT can be generated by: choosing independently to write a 0 with probability ˜ r0 =

1 ϕ2

and a 1 with probability ˜ r1 = 1

ϕ,

deleting pairs of consecutive 1’s. Confluent SFT A SFT is confluent if for any i, j, k ∈ A such that both ij and jk are forbidden, then i = k.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-33
SLIDE 33

Second rejection sampling

The Parry measure of the Fibonacci SFT can be generated by: choosing independently to write a 0 with probability ˜ r0 =

1 ϕ2

and a 1 with probability ˜ r1 = 1

ϕ,

deleting pairs of consecutive 1’s. Confluent SFT A SFT is confluent if for any i, j, k ∈ A such that both ij and jk are forbidden, then i = k. Proposition [J. Mairesse - I. Marcovici] For confluent SFT, the Parry measure can be generated by independent draws of letters and deletion of forbidden patterns.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-34
SLIDE 34

Coloring trees...

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-35
SLIDE 35

SFT defined on trees

Let A be a (symmetric) matrix defining allowed and forbidden patterns, and consider the corresponding SFT Σd

A on the infinite

regular tree of degree d + 1.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-36
SLIDE 36

SFT defined on trees

Let A be a (symmetric) matrix defining allowed and forbidden patterns, and consider the corresponding SFT Σd

A on the infinite

regular tree of degree d + 1. Question: how to construct Markov-uniform measures on Σd

A?

d = 2

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-37
SLIDE 37

SFT defined on trees

Let A be a (symmetric) matrix defining allowed and forbidden patterns, and consider the corresponding SFT Σd

A on the infinite

regular tree of degree d + 1. Question: how to construct Markov-uniform measures on Σd

A?

d = 2

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-38
SLIDE 38

SFT defined on trees

Idea 1: consider a (reversible) Markov chain P on the alphabet A,

  • f stationary distribution π.

Choose the letter at one given vertice according to π and then label the other vertices using P.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-39
SLIDE 39

SFT defined on trees

Idea 1: consider a (reversible) Markov chain P on the alphabet A,

  • f stationary distribution π.

Choose the letter at one given vertice according to π and then label the other vertices using P. Example:

i j k l

Probability of this pattern: π(i)Pi,jPj,kPj,l

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-40
SLIDE 40

SFT defined on trees

Idea 1: consider a (reversible) Markov chain P on the alphabet A,

  • f stationary distribution π.

Choose the letter at one given vertice according to π and then label the other vertices using P. Example:

i j k l

Probability of this pattern: π(i)Pi,jPj,kPj,l = π(j)Pj,iPj,kPj,l = . . .

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-41
SLIDE 41

SFT defined on trees

Idea 1: consider a (reversible) Markov chain P on the alphabet A,

  • f stationary distribution π.

Choose the letter at one given vertice according to π and then label the other vertices using P. Example:

i j k l

Probability of this pattern: π(i)Pi,jPj,kPj,l = π(j)Pj,iPj,kPj,l = . . . For given i, k, l, we want this value to be independent of the letter j such that the pattern is allowed.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-42
SLIDE 42

SFT defined on trees

Idea 2: like for the Parry measure, choose P under the form: Pi,j = Ai,j rj

  • s∈A Ai,srs

= Ai,j rj

  • s∈S(i) rs

, for some probability vector (ri)i∈A.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-43
SLIDE 43

SFT defined on trees

Idea 2: like for the Parry measure, choose P under the form: Pi,j = Ai,j rj

  • s∈A Ai,srs

= Ai,j rj

  • s∈S(i) rs

, for some probability vector (ri)i∈A. Then, π(i)Pi,jPj,kPj,l = π(i) rj

  • s∈A Ai,srs

rk

  • s∈A Aj,srs

rl

  • s∈A Aj,srs

i j k l

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-44
SLIDE 44

SFT defined on trees

Idea 2: like for the Parry measure, choose P under the form: Pi,j = Ai,j rj

  • s∈A Ai,srs

= Ai,j rj

  • s∈S(i) rs

, for some probability vector (ri)i∈A. Then, π(i)Pi,jPj,kPj,l = π(i) rj

  • s∈A Ai,srs

rk

  • s∈A Aj,srs

rl

  • s∈A Aj,srs

Let us try to choose (ri)i∈A such that:

  • s∈A

Aj,srs = λrj 1/2, for any j ∈ A !

i j k l

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-45
SLIDE 45

SFT defined on trees

For a tree of degree d + 1, the problem is to find a probability distribution (ri)i∈A such that for some λ > 0, A    r1 . . . rn    = λ    r1/d

1.

. . r1/d

n

   .

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-46
SLIDE 46

SFT defined on trees

For a tree of degree d + 1, the problem is to find a probability distribution (ri)i∈A such that for some λ > 0, A    r1 . . . rn    = λ    r1/d

1.

. . r1/d

n

   . Proposition Let A be an irreducible non-negative matrix, and let d ≥ 1. There exist λ > 0 and r1, . . . , rn > 0 satisfying n

i=1 ri = 1 and:

A    r1 . . . rn    = λ    r1/d

1.

. . r1/d

n

   .

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-47
SLIDE 47

SFT defined on trees

For a tree of degree d + 1, the problem is to find a probability distribution (ri)i∈A such that for some λ > 0, A    r1 . . . rn    = λ    r1/d

1.

. . r1/d

n

   . Proposition Let A be an irreducible non-negative matrix, and let d ≥ 1. There exist λ > 0 and r1, . . . , rn > 0 satisfying n

i=1 ri = 1 and:

A    r1 . . . rn    = λ    r1/d

1.

. . r1/d

n

   .

  • Proof. Fixed point theorem.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-48
SLIDE 48

SFT defined on trees

For a tree of degree d + 1, the problem is to find a probability distribution (ri)i∈A such that for some λ > 0, A    r1 . . . rn    = λ    r1/d

1.

. . r1/d

n

   . Proposition Let A be an irreducible non-negative matrix, and let d ≥ 1. There exist λ > 0 and r1, . . . , rn > 0 satisfying n

i=1 ri = 1 and:

A    r1 . . . rn    = λ    r1/d

1.

. . r1/d

n

   .

  • Proof. Fixed point theorem.
  • Remark. λ and (ri)i∈A may not be unique.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-49
SLIDE 49

SFT defined on trees

Proposition If the distribution of probability (ri)i∈A satisfies A    r1 . . . rn    = λ    r1/d

1.

. . r1/d

n

   for some λ > 0, then the Markov chain given by: Pi,j = Ai,j rj

  • s∈A Ai,srs

= Ai,j rj λr1/d

i

, defines a Markov-uniform measure on the SFT ΣA.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-50
SLIDE 50

Fibonacci SFT on trees

We search P = α 1 − α 1

  • , such that

1 1 1 α 1 − α

  • = λ
  • α1/d

(1 − α)1/d

  • .

For any d ≥ 1, there exists a unique solution, given by r0 = α and r1 = 1 − α, where α is the unique positive solution of the equation αd+1 = 1 − α. For d = 1, we recover r0 = 1

ϕ and r1 = 1 ϕ2 .

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-51
SLIDE 51

Fibonacci SFT on trees

We search P = α 1 − α 1

  • , such that

1 1 1 α 1 − α

  • = λ
  • α1/d

(1 − α)1/d

  • .

For any d ≥ 1, there exists a unique solution, given by r0 = α and r1 = 1 − α, where α is the unique positive solution of the equation αd+1 = 1 − α. For d = 1, we recover r0 = 1

ϕ and r1 = 1 ϕ2 .

But we also have examples of SFT for which there are several Markov-uniform measures...

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-52
SLIDE 52

Fibonacci SFT on trees

We search P = α 1 − α 1

  • , such that

1 1 1 α 1 − α

  • = λ
  • α1/d

(1 − α)1/d

  • .

For any d ≥ 1, there exists a unique solution, given by r0 = α and r1 = 1 − α, where α is the unique positive solution of the equation αd+1 = 1 − α. For d = 1, we recover r0 = 1

ϕ and r1 = 1 ϕ2 .

But we also have examples of SFT for which there are several Markov-uniform measures... Interpretation of these measures in terms of entropy?

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-53
SLIDE 53

Sampling using probabilistic cellular automata

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-54
SLIDE 54

One-dimensional Fibonacci SFT

Consider a configuration distributed according to the Parry measure π of the Fibonacci SFT.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-55
SLIDE 55

One-dimensional Fibonacci SFT

Consider a configuration distributed according to the Parry measure π of the Fibonacci SFT.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-56
SLIDE 56

One-dimensional Fibonacci SFT

Consider a configuration distributed according to the Parry measure π of the Fibonacci SFT.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-57
SLIDE 57

One-dimensional Fibonacci SFT

Consider a configuration distributed according to the Parry measure π of the Fibonacci SFT. By the Markov-uniform property, the new sequence is still distributed according to π.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-58
SLIDE 58

One-dimensional Fibonacci SFT

π X−2 X−1 X0 X1 X2 X3 X4 X5 X6 Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-59
SLIDE 59

One-dimensional Fibonacci SFT

π X−2 X−1 X0 X1 X2 X3 X4 X5 X6

For all i ∈ Z, if X2i = X2i+2 = 0, we flip the value of X2i+1 with probability 1/2.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-60
SLIDE 60

One-dimensional Fibonacci SFT

π X−2 X−1 X0 X1 X2 X3 X4 X5 X6

For all i ∈ Z, if X2i = X2i+2 = 0, we flip the value of X2i+1 with probability 1/2.

X−2 X−1 X0 X1 X2 X3 X4 X5 X6 FA FA π2 π2 Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-61
SLIDE 61

One-dimensional Fibonacci SFT

1 1 1 1 with probability 1/2 1 with probability 1/2 (with probability 1) X−2 X−1 X0 X1 X2 X3 X4 X5 X6 FA FA π2 π2

Proposition The projection π2 of the Parry measure on odd (resp. even) sites is an invariant measure of the probabilistic cellular automaton.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-62
SLIDE 62

General one-dimensional SFT

For a general one-dimensional SFT ΣA, we consider the PCA FA defined by:

k with proba 1/|{s ∈ A ; isj ∈ W3}| if ikj ∈ W3 (and with proba 0 otherwise) i j

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-63
SLIDE 63

General one-dimensional SFT

For a general one-dimensional SFT ΣA, we consider the PCA FA defined by:

k with proba 1/|{s ∈ A ; isj ∈ W3}| if ikj ∈ W3 (and with proba 0 otherwise) i j

Generalisation to Zd, d ≥ 2 and to infinite trees (bipartite graphs).

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-64
SLIDE 64

General one-dimensional SFT

For a general one-dimensional SFT ΣA, we consider the PCA FA defined by:

k with proba 1/|{s ∈ A ; isj ∈ W3}| if ikj ∈ W3 (and with proba 0 otherwise) i j

Generalisation to Zd, d ≥ 2 and to infinite trees (bipartite graphs). When the measure of maximal entropy is unique, if the PCA dynamics is ergodic (convergence to the measure of maximal entropy of the SFT from any initial configuration), we can use it for sampling.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-65
SLIDE 65

General one-dimensional SFT

For a general one-dimensional SFT ΣA, we consider the PCA FA defined by:

k with proba 1/|{s ∈ A ; isj ∈ W3}| if ikj ∈ W3 (and with proba 0 otherwise) i j

Generalisation to Zd, d ≥ 2 and to infinite trees (bipartite graphs). When the measure of maximal entropy is unique, if the PCA dynamics is ergodic (convergence to the measure of maximal entropy of the SFT from any initial configuration), we can use it for sampling. Monte Carlo method Perfect sampling via coupling from the past

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-66
SLIDE 66

PCA family of the Fibonacci SFT

1 1 1 1 with probability 1/2 1 with probability 1/2 (with probability 1) Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-67
SLIDE 67

PCA family of the Fibonacci SFT

1 1 1 1 with probability p 1 with probability 1 − p (with probability 1) Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-68
SLIDE 68

PCA family of the Fibonacci SFT

1 1 1 1 with probability p 1 with probability 1 − p (with probability 1)

Proposition [J. Martin - I. Marcovici] For any p ∈ (0, 1), the PCA above is ergodic.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-69
SLIDE 69

PCA family of the Fibonacci SFT

1 1 1 1 with probability p 1 with probability 1 − p (with probability 1)

Proposition [J. Martin - I. Marcovici] For any p ∈ (0, 1), the PCA above is ergodic. Furthermore, its envelope PCA is ergodic, meaning that we can sample its unique invariant measure perfectly by coupling from the past.

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-70
SLIDE 70

PCA family of the Fibonacci SFT

1 1 1 1 with probability p 1 with probability 1 − p (with probability 1)

Proposition [J. Martin - I. Marcovici] For any p ∈ (0, 1), the PCA above is ergodic. Furthermore, its envelope PCA is ergodic, meaning that we can sample its unique invariant measure perfectly by coupling from the past. General criterion ensuring the ergodicity of the PCA associated to a SFT?

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-71
SLIDE 71

PCA family of the Fibonacci SFT

1 1 1 1 with probability p 1 with probability 1 − p (with probability 1)

Proposition [J. Martin - I. Marcovici] For any p ∈ (0, 1), the PCA above is ergodic. Furthermore, its envelope PCA is ergodic, meaning that we can sample its unique invariant measure perfectly by coupling from the past. General criterion ensuring the ergodicity of the PCA associated to a SFT? And how to be sure that the coupling from the past algorithm will stop in finite time?

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-72
SLIDE 72

Conclusion

Different descriptions of the measures of maximal entropy using i.i.d. random variables

New results for confluent one-dimensional SFT Exploratory works for multi-dimensional SFT

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-73
SLIDE 73

Conclusion

Different descriptions of the measures of maximal entropy using i.i.d. random variables

New results for confluent one-dimensional SFT Exploratory works for multi-dimensional SFT

Introduction of a PCA dynamics

When is the PCA ergodic? In that case, can we always sample its invariant measure by coupling from the past?

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type

slide-74
SLIDE 74

Conclusion

Different descriptions of the measures of maximal entropy using i.i.d. random variables

New results for confluent one-dimensional SFT Exploratory works for multi-dimensional SFT

Introduction of a PCA dynamics

When is the PCA ergodic? In that case, can we always sample its invariant measure by coupling from the past?

Ir` ene Marcovici Uniform Sampling of Subshifts of Finite Type