In [2]:
%display latex 
In [3]:
load("../specifier.py")

We consider the class $\mathrm{Av}(132)$. It contains no simple permutation. Hence we obtain its specification as follows:

In [4]:
S = Specification([],[[1,3,2]])
In [5]:
S
Out[5]:

In the specification algorithm, the families $\mathcal T_0, \ldots, \mathcal T_d$ are families of constrained permutations. We can see which families exactly as follows:

In [6]:
S.concrete()
Out[6]:
In [7]:
S.solve()
Out[7]:

We see that $T_0,T_1,T_3$ have radius of convergence $1/4$, while the others have strictly larger radius. Hence the critical families are $\mathcal T_0,\mathcal T_1,\mathcal T_3$. Hypothesis (SC) (strong connectivity among critical families) holds, and the specification is branching (see for instance the term $\ominus[\mathcal T_{3},\mathcal T_{0}]$).

Since the series are convergent, we may perform Boltzmann sampling right at the radius of convergence:

In [8]:
C = S.sampler(0.25, formal=True)
In [9]:
perm = C.run(10)
In [10]:
list_plot(perm,aspect_ratio=1)
Out[10]:

The limiting permuton seems to be the antidiagonal. This is confirmed by our results: indeed, we are in the essentially branching case, and the substitutions containing two critical series all have a $\ominus$ at the root (it is always $\ominus[\mathcal T_0,\mathcal T_3]$).

In [ ]:
 
In [ ]:
 
In [ ]:
 
In [ ]:
 
In [ ]: