1. Introduction
An iterated function system (IFS) is a finite collection of contracting maps on a complete metric space. In this article, we consider IFSs 
$\{A_u+\boldsymbol{v}_u\}_{u\in\mathcal U}$ with 
$\mathcal U$ a finite set that are given by a collection 
$(A_u)_{u \in \mathcal U} \in GL_2(\mathbb R)^{\# \mathcal U}$ of matrices of the form
\begin{equation}
A_u=\left[\begin{array}{cc}
b_u & 0 \\
0 & c_u
\end{array}\right], \quad \text{with }\, 0 \lt |b_u|, |c_u| \lt 1, \, \, u \in \mathcal U,
\end{equation} and a collection 
$(\boldsymbol{v}_u)_{u \in \mathcal U}$ of vectors of the form
\begin{align*} \boldsymbol{v}_u = \left[\begin{array}{cc}
\beta_u\\
\gamma_u
\end{array}\right] \in \mathbb R^2, \quad u \in \mathcal U.\end{align*} The unique non-empty compact subset 
$\Lambda \subseteq [0,1]^2$ that satisfies
\begin{equation}
\Lambda = \bigcup_{u \in \mathcal U} (A_u + \boldsymbol{v}_u)(\Lambda),
\end{equation} which exists by Hutchinson’s theorem [Reference John33], is called the attractor or self-affine set of the diagonally affine IFS 
$\{A_u + \boldsymbol{v}_u\}_{u \in \mathcal U}$.
We are interested in multifractal properties of certain types of self-affine sets. Let 
$\sigma: \mathcal U^{\mathbb N} \to \mathcal U^{\mathbb N}$ denote the left shift, so 
$\sigma (\xi_n)_{n \geqslant 1} = (\xi_{n+1})_{n \geqslant 1}$. For any continuous potential 
$\Phi: \mathcal U^{\mathbb N} \rightarrow \mathbb R^d$, 
$d \geqslant 1$, and for a given vector 
$\boldsymbol{\alpha} = (\alpha_1, \ldots, \alpha_d) \in \mathbb R^d$, the symbolic level set is given by
\begin{equation} E_{\Phi}(\boldsymbol{\alpha}):=\left\{\xi \in \mathcal U^{\mathbb N} \, :\, \lim _{n \rightarrow \infty} \frac{1}{n} \sum_{i=0}^{n-1} \Phi (\sigma^i(\xi))=\boldsymbol{\alpha}\right\}.
\end{equation} For 
$\xi \in \mathcal U^{\mathbb N}$ and 
$n \geqslant 1$, let 
$A_{\xi|_{n-1}} = A_{\xi_1} \cdots A_{\xi_{n-1}}$ and define 
$\pi: \mathcal U^{\mathbb N} \to \Lambda $ by
\begin{equation}
\pi(\xi) = \sum_{n \geqslant 1} A_{\xi|_{n-1}} \boldsymbol{v}_{\xi_n}.
\end{equation} The α-level set for Φ on Λ is the set 
$\pi(E_{\Phi}(\boldsymbol{\alpha}))$. There are various known results on the size of level sets, both in terms of Hausdorff dimension and topological entropy (in the sense of [Reference Bowen17]), see for instance [Reference Barral and Mensi7, Reference Barreira and Saussol10–Reference Barreira, Saussol and Schmeling12, Reference Fan, Feng and Jun25, Reference Feng and Huang27, Reference Käenmäki and Reeve36, Reference Mohammadpour43, Reference Mohammadpour45, Reference Reeve52]. Multifractal results for self-affine sets in terms of Lyapunov dimensions were obtained in [Reference Bárány, Jordan, Käenmäki and Rams2] for collections of matrices 
$(A_u)_{u \in \mathcal U}$ under certain strong-irreducibility and proximality conditions. One can find more information about the multifractal formalism in [Reference Barreira and Gelfert8, Reference Climenhaga21, Reference Mohammadpour44]. Here, we consider the Hausdorff dimension of level sets for self-affine sets.
Let 
$\mathcal F = \{f_u: \mathbb R \to \mathbb R \}_{u \in \mathcal U}$ be an IFS of real-valued affine maps of the form 
$f_u(y) = d_u y + \delta_u$, so 
$|d_u| \lt 1$ for each 
$u \in \mathcal U$. For a sequence 
$\boldsymbol{u}=u_1\cdots u_n \in \mathcal U^n$, 
$n \geqslant 1$, we write
\begin{align*} f_{\boldsymbol{u}} (y) := f_{u_1} \circ \cdots \circ f_{u_n} (y) = d_{\boldsymbol{u}}y + \delta_{\boldsymbol{u}}.\end{align*} We say that 
$\mathcal F$ has exact overlaps if there are 
$\boldsymbol{u}, \boldsymbol{u}' \in \mathcal U^n$ for some 
$n \geqslant 1$ such that 
$f_{\boldsymbol{u}}= f_{\boldsymbol{u}'}$. For 
$\boldsymbol{u},\boldsymbol{u}' \in \mathcal U^n$, we define the distance
\begin{align*}\text{dist}(f_{\boldsymbol{u}}, f_{\boldsymbol{u}'}):=\left\{
\begin{array}{@{}ll@{}}
|\delta_{\boldsymbol{u}} -\delta_{\boldsymbol{u}'}|, & \text{if}\ d_{\boldsymbol{u}}= d_{\boldsymbol{u}'}; \\
\infty, & \text{otherwise}.
\end{array}\right.\end{align*} In his breakthrough result, Hochman [Reference Hochman29] introduced the exponential separation condition (ESC) to calculate the dimension of self-similar measures. We say that 
$\mathcal F$ satisfies the ESC if there exists a constant c > 0 and infinitely many integers 
$n \geqslant 1$ such that for all 
$\boldsymbol{u}, \boldsymbol{u}' \in \mathcal U^n$,
\begin{align*} \text{dist}(f_{\boldsymbol{u}}, f_{\boldsymbol{u}'}) \geqslant c^n.\end{align*} We say that a diagonally affine IFS 
$\{A_u+\boldsymbol{v}_u\}_{u\in\mathcal U}$ satisfies the strong open set condition (SOSC) if there is an open set 
$V \subseteq \mathbb R^2$ such that all the sets 
$(A_u + \boldsymbol{v}_u)(V)$ are pairwise disjoint, 
$\bigcup_{u \in \mathcal U} (A_u + \boldsymbol{v}_u) (V) \subseteq V$ and 
$\Lambda \cap V \neq \emptyset$, where Λ is as in (1.2).
Our first result is for the following class of diagonally affine IFSs. Let 
$\mathcal D$ be the collection of all IFSs 
$\{A_u+ \boldsymbol{v}_u\}_{u\in \mathcal U}$ with matrices as in (1.1) that satisfy the SOSC together with either (D) or (Dʹ):
(D)
$|b_u| \gt |c_u|$ for all $u \in \mathcal U$ and(a) the IFS obtained from projecting to the first coordinate
$\mathcal G_1:= \{g_{1,u}(y)= b_u y +\beta_u\}_{u \in \mathcal U}$ satisfies the ESC, or,(b) bu is algebraic for all
$u \in \mathcal U$ and the IFS obtained from projecting to the first coordinate $\mathcal G_1:= \{g_{1,u}(y)= b_u y +\beta_u\}_{u \in \mathcal U}$ has no exact overlaps;
(Dʹ)
$|b_u| \lt |c_u|$ for all $u \in \mathcal U$ and(a) the IFS obtained from projecting to the second coordinate
$\mathcal G_2:= \{g_{2,u}(y)= c_u y +\gamma_u\}_{u \in \mathcal U}$ satisfies the ESC, or,(b) cu is algebraic for all
$u \in \mathcal U$ and the IFS obtained from projecting to the second coordinate $\mathcal G_2:= \{g_{2,u}(y)= c_u y +\gamma_u\}_{u \in \mathcal U}$ has no exact overlaps.
An IFS satisfying condition (Dʹ)(a) is shown in figure 1(a).
Let 
$\mathcal{M}(\mathcal U^{\mathbb N}, \sigma)$ denote the set of all σ-invariant Borel probability measures on 
$\mathcal U^{\mathbb N}$ and let
\begin{equation}
\begin{aligned}
L_{\Phi} & :=\left\{\boldsymbol{\alpha} \in \mathbb{R}^d\, :\, \exists\, \xi \in \mathcal U^{\mathbb N} \text {with } \lim _{n \rightarrow \infty} \frac{1}{n} \sum_{i=0}^{n-1} \Phi (\sigma^i \xi)=\boldsymbol{\alpha} \right\} \\
& =\left\{\boldsymbol{\alpha} \in \mathbb{R}^d\, :\, \exists \, \mu \in \mathcal{M}(\mathcal U^{\mathbb N}, \sigma) \text {with } \int_{\mathcal U^{\mathbb N}} \Phi \, \mathrm{d} \mu=\boldsymbol{\alpha} \right\}
\end{aligned}
\end{equation}Two examples of IFSs. The coloured rectangles indicate the images of the unit square under the maps in the IFS.
be the collection of all vectors α for which the corresponding level set is non-empty, known as the spectrum of Φ. We use 
$\mathring{L}_{\Phi}$ to denote the interior of 
$L_\Phi$. For any 
$\mu \in \mathcal{M}(\mathcal U^{\mathbb N}, \sigma)$, let 
$\operatorname{dim}_{\mathrm{L}}(\mu)$ denote the Lyapunov dimension of µ. Let P denote the topological pressure. Our first result is as follows.
Theorem 1.1 Let 
$\left\{A_u+ \boldsymbol{v}_u\right\}_{u \in \mathcal U} \in \mathcal D$ and let 
$\Phi: \mathcal U^{\mathbb N} \rightarrow \mathbb{R}^d$, 
$d \geqslant 1$, be a continuous potential. Then for each 
$\boldsymbol{\alpha} \in \mathring{L}_{\Phi}$,
\begin{align*} \begin{aligned}
\operatorname{dim}_{\mathrm{H}}\left(\pi(E_{\Phi}(\boldsymbol{\alpha}))\right) & =\sup \bigg\{\operatorname{dim}_{\mathrm{L}}(\mu)\, :\, \mu \in \mathcal{M}(\mathcal U^{\mathbb N}, \sigma) \text {and } \int_{\mathcal U^{\mathbb N}} \Phi\, \mathrm{d} \mu=\boldsymbol{\alpha} \bigg\} \\
& =\sup \bigg\{s \geqslant 0\, :\, \inf _{q \in \mathbb{R}^d} P\left(\log \varphi^s+\langle q, \Phi-\boldsymbol{\alpha}\rangle\right) \geqslant 0\bigg\}.
\end{aligned}
\end{align*} The second family of IFS sets that we consider in this article is motivated by a specific type of representations of real numbers called generalized Lüroth series (GLS) expansions as described in [Reference Barrionuevo, Burton, Dajani and Kraaikamp13]. As the name suggests, GLS expansions are generalizations of Lüroth expansions, which were introduced in 1883 by Lüroth [Reference Lüroth42] and for 
$x \in [0,1]$ have the form
\begin{align*} x = \sum_{n \geqslant 1} \frac{a_n}{\prod_{i=1}^n a_i(a_i+1)}, \quad a_n \in \mathbb N, \, n \geqslant 1.\end{align*} Lüroth expansions can be obtained from the IFS 
$\{l_k:[0,1]\to [0,1]\}_{k \in \mathbb N}$ where 
$l_k(x) = \frac{k+x}{k(k+1)}$, see e.g. [Reference Jager and de Vroedt31]. If for 
$x \in [0,1]$, there is a sequence 
$(a_n)_{n \geqslant 1} \in \mathbb N^{\mathbb N}$ such that
\begin{align*} x = \lim_{n \to \infty} l_{a_1} \circ \cdots \circ l_{a_n}(0),\end{align*} then x has a Lüroth expansion with digits given by 
$(a_n)_{n \geqslant 1} \in \mathbb N^{\mathbb N}$. While Lüroth expansions take their digits from the infinite digit set 
$\mathbb N$ and all terms in the expansion are positive, a GLS number system can have either finite or infinite digit sets and the corresponding GLS expansions can have both positive and negative terms. Given a finite or countably infinite digit set 
$\mathcal I$, a partition 
$\{[\ell_k,r_k] \}_{k\in \mathcal I}$ of 
$[0,1]$ into closed intervals and a vector 
$ (\varepsilon_k)_{k\in \mathcal I} \in \{0,1\}^{\# \mathcal I}$, one can consider the IFS
\begin{equation}
\{g_k:[0,1]\to [0,1]\}_{k\in \mathcal I},
\end{equation} where gk maps the interval 
$[0,1]$ affinely onto 
$[\ell_k,r_k]$ in an orientation-preserving manner if 
$\varepsilon_k=0$ and in an orientation-reversing manner if 
$\varepsilon_k=1$. In other words, if we write 
$K_k = (r_k-\ell_k)^{-1}$, then 
$g_k(x) = \ell_k + \frac{\varepsilon_k + (-1)^{\varepsilon_k}x }{K_k}$. Since 
$g_k([0,1]) = [\ell_k,r_k]$ for each k, it follows that for each 
$x \in [0,1]$, there is a sequence 
$(a_n)_{n \geqslant 1}$ such that
\begin{align*} x = \lim_{n \to \infty} g_{a_1} \circ \cdots \circ g_{a_n}(0).\end{align*}Thus, x can be expressed as
\begin{equation}
x = \sum_{n \geqslant 1} (-1)^{\sum_{i=1}^{n-1}\varepsilon_{a_i}} \frac{\ell_{a_n}K_{a_n}+\varepsilon_{a_n}}{\prod_{i=1}^n K_{a_i}},
\end{equation} which is called a GLS expansion of x with digit set 
$\mathcal I$. Here, we let 
$\sum_{i=1}^0 \varepsilon_{a_i}=0$ and 
$\prod_{i=1}^0K_{a_i}=1$. One recovers the Lüroth expansions by taking 
$\mathcal I=\mathbb N$, 
$[\ell_k,r_k] = [\frac1{k+1}, \frac1k]$ and 
$\varepsilon_k=0$ for each 
$k \geqslant 1$ and one obtains integer base N-expansions by setting 
$\mathcal I = \{0,1, \ldots, N-1\}$ and taking 
$[\ell_k,r_k] = [\frac{k}{N}, \frac{k+1}{N}]$ and 
$\varepsilon_k=0$. The expansions from (1.7) can also be seen as signed versions of Cantor base expansions, as introduced by Cantor in [Reference Cantor18]. GLS expansions have been considered previously in [Reference Arroyo1, Reference Boonstra and Kalle16, Reference Huang and Kalle30, Reference Jordan, Munday and Sahlsten34, Reference Kautzsch, Kesseböhmer, Samuel and Stratmann40, Reference Munday47] and recently also in relation to neural networks [Reference Balakrishnan, Kathpalia, Saha and Nagaraj6, Reference Remya Ajai, Harikrishnan and Nagaraj53]. Level sets for Lüroth expansions and more generally GLS expansions have been considered in particular with respect to digit frequencies, see [Reference Barreira and Iommi9]. Such level sets are known as Besicovitch–Eggleston sets due to the results from [Reference Besicovitch14] by Besicovitch and [Reference Eggleston24] by Eggleston on the Hausdorff dimension on digit frequency level sets for integer base expansions.
In the above setting, for any given GLS number system, all but countably many numbers in 
$[0,1]$ have a unique GLS expansion in that system and the numbers that do not have a unique expansion have exactly two expansions. In this article, we consider IFSs that correspond to GLS number systems with redundancy, that is, in which all numbers have uncountably many different representations in the system. Number systems with redundancy have proven interesting in several settings, including signed binary expansions where they are used to find so-called minimal weight expansions, i.e. expansions that maximize the number of digits 0, see e.g. [Reference Dajani and Kalle22, Reference Koyama and Tsuruoka41, Reference Morain and Olivos46], and in non-integer base expansions in relation to applications in analogue-to-digital converters and random number generation, see e.g. [Reference Daubechies, DeVore, Güntürk and Vaishampayan23, Reference Jitsumatsu and Matsumura32]. Number systems with redundancy have also been considered in [Reference Kalle, Kempton and Verbitskiy37–Reference Kalle, Matache, Tsujii and Verbitskiy39] for continued fraction expansions and Lüroth expansions. To obtain a GLS number system with redundancy, we let 
$J \in \mathbb N_{\geqslant 2}$ and start with J IFSs that correspond to J different GLS number systems with finite digit sets. We combine these into one diagonally affine IFS on 
$\mathbb R^2$, which we call a GLS IFS, by using a positive probability vector 
$ (p_j)_{0 \leqslant j \lt J}$, so 
$p_j \gt 0$ for all 
$0 \leqslant j \lt J$ and 
$\sum_{0 \leqslant j \lt J} p_j=1$. This vector 
$(p_j)_{0 \leqslant j \lt J}$ can be thought of as the probabilities with which the jth GLS number system is chosen to generate the nth digit in the expansions for any 
$n \geqslant 1$. Therefore, a GLS IFS is given by the following data:
(i) an integer
$J \in \mathbb N_{\geqslant 2}$ and a positive probability vector $(p_j)_{0 \leqslant j \lt J}$;(ii) for each
$0 \leqslant j \lt J$, a number $B_j \in \mathbb N_{\geqslant 2}$, a partition $0=r_{(j,0)} \lt r_{(j,1)} \lt \cdots \lt r_{(j,B_j)}=1$ and a vector $(\varepsilon_{(j,k)})_{0 \leqslant k \lt B_j} \in \{0,1\}^{B_j}$.
If we set 
$\mathcal E = \{(j,k) \, : \, 0 \leqslant j \lt J, \, 0 \leqslant k \lt B_j\}$ and for each 
$(j,k) \in \mathcal E$ let 
$q_{(j,k)}=r_{(j,k+1)}-r_{(j,k)}$ and
\begin{equation}
A_{(j,k)} = \left[\begin{array}{cc}
\displaystyle (-1)^{\varepsilon_{(j,k)}}q_{(j,k)} & 0 \\
0 & p_j
\end{array}\right], \quad \boldsymbol{v}_{(j,k)} = \left[\begin{array}{c}
\displaystyle r_{(j,k)} + \varepsilon_{(j,k)}q_{(j,k)} \\
\sum_{i=0}^{j-1}p_i
\end{array} \right],
\end{equation} then we call the IFS 
$\{A_e + \boldsymbol{v}_e\}_{e \in \mathcal E}$ a GLS IFS. See figure 1(b) for an example.
We mention a few particular properties of GLS IFSs. Each GLS IFS satisfies the SOSC and has 
$\Lambda = [0,1]^2$. For the projection onto the first coordinate, we use 
$\mathcal G_1=\{h_e:[0,1] \to [0,1]\}_{e \in \mathcal E}$, where for each 
$e \in \mathcal E$, we set
\begin{align*} h_e(x) = r_e + q_e(\varepsilon_e + (-1)^{\varepsilon_e} x).\end{align*} Without additional assumptions, the GLS IFS need not fall into one of the categories (D)(a) or (D)(b). The projection on the second coordinate 
$\mathcal G_2 = \{g_{2,(j,k)}(y)= p_j y + \sum_{i=0}^{j-1} p_i \}_{(j,k) \in \mathcal E}$ of a GLS IFS contains several duplicates of each map. Therefore, GLS IFSs do not fall into the class of diagonally affine IFSs that satisfy (Dʹ), but by removing these duplicates they can potentially contain a subcollection of contractions that satisfies (Dʹ) as shown in figure 1.
We can obtain number expansions from a GLS IFS in the following way. For each 
$x \in [0,1]$, there are sequences 
$(e_m)_{m \in \mathbb N} \in \mathcal E^{\mathbb N}$ such that x can be written as
\begin{equation}
x = \lim_{m \to \infty} h_{e_1} \circ \cdots \circ h_{e_m} (0).
\end{equation} If for each 
$m \in N$, we write
\begin{equation} s_m = \varepsilon_{e_m}, \quad K_m = q_{e_m}^{-1}, \quad t_m = r_{e_m} + \varepsilon_{e_m}K_m^{-1},
\end{equation}then it follows from (1.9) that
\begin{equation}
x = \sum_{m \geqslant 1} (-1)^{\sum_{i=1}^{m-1}s_i} \frac{t_m}{\prod_{i=1}^{m}K_i}
\end{equation} and we can see the resemblance with Lüroth expansions. It is shown in proposition 4.1 that, under the additional assumption on the GLS IFS that 
$h_e \neq h_{e'}$ whenever 
$e \neq e'$, indeed all numbers 
$x \in [0,1]$ have uncountably many different representations of the form (1.9). We give several examples of GLS IFSs and the associated number expansions at the end of the article.
For GLS IFSs, we consider the potential that captures digit frequencies. For 
$e \in \mathcal E$, let 
$[e] \subseteq \mathcal E^{\mathbb N}$ denote the cylinder set of those sequences that have 
$e \in \mathcal E$ as their first term and 
$\mathbb 1_{[e]}: \mathcal E^{\mathbb N} \to \{0,1\}$ the indicator function on 
$[e]$. Define the continuous potential 
$\mathbb 1: \mathcal E^{\mathbb N} \to \{0,1\}^{\#\mathcal E}$ by 
$\mathbb 1(\omega) = (\mathbb 1_{[e]}(\omega))_{e \in \mathcal E}$. For each 
$e \in \mathcal E$ and 
$\omega \in \mathcal E^{\mathbb N}$, write
\begin{align*} \tau_e (\omega) = \lim_{n \to \infty} \frac{\# \{1 \leqslant i \leqslant n \, : \, \omega_i = e\}}{n}\end{align*} for the frequency of the digit e in ω. Consider a frequency vector 
$\boldsymbol{\alpha} = (\alpha_e)_{e \in\mathcal E} \in [0,1]^{\# \mathcal E}$, i.e. that satisfies 
$\sum_{e \in \mathcal E} \alpha_e =1$, and let 
$F(\boldsymbol{\alpha}) = \pi (E_{\mathbb 1}(\boldsymbol{\alpha}))$. Then
\begin{equation}
F(\boldsymbol{\alpha}) = \{(x,y) \in [0,1]^2 \, : \, \exists\, \omega\in\pi^{-1}\{(x,y)\}\, \textrm{s.t. } \tau_e(\omega) = \alpha_e \, \text{for all } \, e \in \mathcal E \}
\end{equation} is the GLS digit frequency level set or Besicovitch–Eggleston set for α. Results on 
$\dim_{\mathrm H} (F(\boldsymbol{\alpha}))$ have been obtained in [Reference Nielsen49, theorem 1] in the specific case of Bedford–McMullen carpets, that is with 
$p_j = \frac1{J}$ for each 
$0 \leqslant j \lt J$ and 
$(-1)^{\varepsilon_e}=1$ and 
$q_e=\frac1N$ for some fixed 
$N \in \mathbb N_{\geqslant 2}$ and all 
$e \in \mathcal E$. This result was extended in [Reference Reeve52, corollary 1] for Lalley–Gatzouras carpets, which are similar to our setting but have the additional requirements that 
$\varepsilon_{(j,k)} =0$ and 
$q_{(j,k)} \leqslant p_j$ for all 
$(j,k) \in \mathcal E$. In the current setting, a lower bound for 
$\dim_{\mathrm H}(F(\boldsymbol{\alpha}))$ in terms of the Ledrappier–Young formula for the Bernoulli measure 
$\mu_{\boldsymbol{\alpha}}$ can be deduced from [Reference Bárány and Käenmäki3, theorem 2.3 and corollary 2.8] in case the two Lyapunov exponents of the Bernoulli measure of the system differ and the frequency vector α is strictly positive. In case the two Lyapunov exponents of the Bernoulli measure 
$\mu_{\boldsymbol{\alpha}}$ are equal, one can apply [Reference Feng and Hu28] to obtain a similar lower bound in terms of the Ledrappier–Young formula for 
$\mu_{\boldsymbol{\alpha}}$.
Here, we will instead, for fixed 
$y \in [0,1]$, focus on the fibre level sets
\begin{equation}
F_y({\boldsymbol{\alpha}}):= \left\{x\in[0,1]\, :\, \exists\,\omega\in\pi^{-1}\{(x,y)\}\, \textrm{s.t. }\tau_e(\omega)=\alpha_e\,\, \text{for all } e\in\mathcal E\right\}.
\end{equation} We only consider frequency vectors α with 
$\alpha_j:= \sum_{k=0}^{B_j-1} \alpha_{(j,k)} \gt 0$ for all 
$0 \leqslant j \lt J$ (otherwise we could just as well have considered a smaller GLS IFS). Let
\begin{equation}
W(\boldsymbol{\alpha}) := \{y \in [0,1] \, : \, F_y(\boldsymbol{\alpha}) \neq \emptyset \}.
\end{equation} Let 
$\mu_{\boldsymbol{\alpha}}$ be the α-Bernoulli measure on 
$\mathcal E^{\mathbb N}$. For each 
$0 \leqslant j \lt J$, let 
$f_j:[0,1]\to [0,1]$ be the map given by
\begin{equation}
f_j(y) = p_jy + \sum_{i=0}^{j-1}p_i
\end{equation}and define the map
\begin{equation} \pi_2: \mathcal E^{\mathbb N} \to [0,1]; \, (j_m,k_m)_{m \geqslant 1} \mapsto \lim_{m \to \infty} f_{j_1} \circ \cdots \circ f_{j_m}(0).
\end{equation} Set 
$\nu_{\boldsymbol{\alpha}} = \mu_{\boldsymbol{\alpha}} \circ \pi_2^{-1}$. As we will see later, 
$\nu_{\boldsymbol{\alpha}}(W(\boldsymbol{\alpha}))=1$. We have the following results on the Hausdorff dimension of the fibre Besicovitch–Eggleston sets.
Theorem 1.2 Let 
$\{A_e + \boldsymbol{v}_e\}_{e \in \mathcal E}$ be a GLS IFS and 
${\boldsymbol{\alpha}}=(\alpha_e)_{e\in\mathcal E} \in [0,1]^{\# \mathcal E}$ a frequency vector. Then
\begin{equation*}\dim_{\mathrm{H}}(F_y({\boldsymbol{\alpha}}))
\leq\frac{\sum_{e \in \mathcal E} \alpha_e \log \alpha_e - \sum_{0 \leqslant j \lt J} \alpha_j \log \alpha_j}{\sum_{e\in\mathcal E} \alpha_e\log q_e}\end{equation*} for all 
$y\in W({\boldsymbol{\alpha}})$. Furthermore, if α satisfies that for each 
$0 \leqslant j \lt J$ there are 
$k, \ell \in B_j$ with 
$k \neq \ell$, 
$\alpha_{(j,k)} \gt 0$ and 
$\alpha_{(j, \ell)} \gt 0$, then
\begin{equation*}\dim_{\mathrm{H}}(F_y({\boldsymbol{\alpha}}))
\geq\frac{\sum_{e \in \mathcal E} \alpha_e \log \alpha_e - \sum_{0 \leqslant j \lt J} \alpha_j \log \alpha_j}{\sum_{e\in\mathcal E} \alpha_e\log q_e}\end{equation*} for 
$\nu_{\boldsymbol{\alpha}}$-a.e. 
$y\in W(\boldsymbol{\alpha})$.
Fibrewise results similar in spirit to theorem 1.2 were obtained in [Reference Nakajima and Takahasi48], where the authors study real numbers with a semi-regular continued fraction expansion that satisfies a certain growth condition on its digits.
The article is outlined as follows. In §2, we provide the necessary preliminaries. We prove theorem 1.1 in §3. Section 4 is devoted to GLS IFSs. Here, we show that given a GLS IFS that has 
$h_e \neq h_{e'}$ whenever 
$e \neq e'$, all 
$x \in [0,1]$ have uncountably many expansions of the form (1.9). We then continue with some results on the spectrum of the Besicovitch–Eggleston sets 
$F(\boldsymbol{\alpha})$ and on the sets 
$W(\boldsymbol{\alpha})$, which will be used in the proof of theorem 1.2. This section also contains the proof of theorem 1.2. Finally, §5 contains some examples.
2. Preliminaries
In this section, we introduce notation and collect several bits of information that are used for the results in the later sections.
2.1. Strings and sequences
Let 
$\mathcal U$ be a finite set of symbols and denote by 
$\mathcal U^{\mathbb N}$ the set of one-sided infinite sequences of elements in 
$\mathcal U$. For each 
$n \geqslant 0$, the set 
$\mathcal U^n$ is the set of words of length n, where we let 
$\mathcal U^0 = \{\varnothing \}$ be the set containing only the empty word, which we denote by 
$\varnothing$. Let 
$\mathcal U^* = \bigcup_{n \geqslant 0}\mathcal U^n$ be the set of all words. For a word 
$\boldsymbol{u} \in \mathcal U$, we use the notation 
$|\boldsymbol{u}|$ for its length, so 
$|\boldsymbol{u}| = n$ if 
$\boldsymbol{u} \in \mathcal U^n$. If 
$\boldsymbol{u} = u_1 \cdots u_n \in \mathcal U^*$, then for each 
$k \leqslant n$, we use the notation 
$\boldsymbol{u}|_k = u_1 \cdots u_k$. Similarly for a sequence 
$\xi \in \mathcal U^{\mathbb N}$ and any 
$n \geqslant 1$, we set 
$\xi|_n = \xi_1 \cdots \xi_n$. The cylinder set corresponding to a word 
$\boldsymbol{u} \in \mathcal U^n$, 
$n \geqslant 0$, is denoted by
\begin{align*} [\boldsymbol{u}] = \{\xi \in \mathcal U^{\mathbb N} \, : \, \xi|_n = \boldsymbol{u} \}.\end{align*} For any sequence 
$\xi \in \mathcal U^{\mathbb N}$, any symbol 
$u \in \mathcal U$ and any 
$n \in \mathbb N$, we use the notation
\begin{align*} \tau_u(\xi,n) = \# \{1 \leqslant m \leqslant n \, : \, \xi_m = u\} \end{align*}for the number of times the symbol u occurs in the first n elements of ξ and
\begin{align*} \tau_u(\xi) = \lim_{n \to \infty} \frac{\tau_u(\xi,n)}{n}\end{align*}for the frequency of the digit u in ξ if it exists. We use this notation in §4.
We can equip 
$\mathcal U^{\mathbb N}$ with a metric η to obtain a compact metric space by setting
\begin{align*} \eta : \mathcal U^{\mathbb N} \to [0,1]; \, (\xi,\upsilon) \mapsto \begin{cases}
2^{-\min \{n \geqslant 1\, : \, \xi_n \neq \upsilon_n \}}, & \text{if } \xi \neq \upsilon,\\
0, & \text{if } \xi=\upsilon.
\end{cases}\end{align*} The left shift is denoted by 
$\sigma: \mathcal U^{\mathbb N} \to \mathcal U^{\mathbb N}$, i.e. 
$\sigma(\xi)_n = \xi_{n+1}$ for each 
$n \geqslant 1$. With a slight abuse of notation, we will use σ to denote the left shift on any sequence space without specifying the alphabet as a subscript whenever no confusion can arise. Cylinder sets are both open and closed and generate the Borel σ-algebra on 
$\mathcal U^{\mathbb N}$. Let 
$\mathcal M (\mathcal U^{\mathbb N}, \sigma)$ denote the set of all shift-invariant Borel probability measures on 
$\mathcal U^{\mathbb N}$. For 
$\mu \in \mathcal M (\mathcal U^{\mathbb N}, \sigma)$, we use 
$h_{\mu}(\sigma)$ to denote the measure-theoretic entropy of µ with respect to σ, which is defined by
\begin{align*}
h_{\mu}(\sigma):=-\lim_{n \to \infty}\frac{1}{n}\sum_{\boldsymbol{u} \in \mathcal{U}_n} \mu([\boldsymbol{u}])\log \mu([\boldsymbol{u}]),
\end{align*} where 
$0\log 0=0.$ Given a probability vector 
$\boldsymbol{p} = (p_u)_{u \in \mathcal U}$, the p-Bernoulli measure 
$\mu_{\boldsymbol{p}}$ is the probability measure on 
$(\mathcal U^{\mathbb N}, \sigma)$ that is defined on the cylinder 
$[\boldsymbol{u}]=[u_1 \cdots u_n]$ by
\begin{align*} \mu_{\boldsymbol{p}}([\boldsymbol{u}])=p_{u_1} \cdots p_{u_n}.\end{align*} Moreover, the measure-theoretic entropy of 
$\mu_{\boldsymbol{p}}$ with respect to σ is given by
\begin{equation}
h_{\mu_{\boldsymbol{p}}}(\sigma) = - \sum_{u \in \mathcal U} \mu_{\boldsymbol{p}}([u]) \log (\mu_{\boldsymbol{p}}([u])).
\end{equation}2.2. Matrix products
Let 
$(A_u)_{u \in \mathcal U} \in GL_2(\mathbb R)^{\# \mathcal U}$ be a collection of matrices as in (1.1). Recall that for a sequence 
$\xi=(\xi_n)_{n \geqslant 1} \in \mathcal U^{\mathbb N}$ and 
$n \in \mathbb N$, we set 
$A_{\xi|_n} = A_{\xi_1} \cdots A_{\xi_n}$. For the entries on the diagonal of 
$A_{\xi|_n}$, write 
$b_{\xi|_n} = b_{\xi_1} \cdots b_{\xi_n}$ and 
$c_{\xi|_n}=c_{\xi_1} \cdots c_{\xi_n}$. For 
$\boldsymbol{u} = u_1 \cdots u_n \in \mathcal U^n$, we similarly write 
$A_{\boldsymbol{u}} = A_{u_1} \cdots A_{u_n}$ with 
$b_{\boldsymbol{u}} = b_{u_1} \cdots b_{u_n}$ and 
$c_{\boldsymbol{u}}=c_{u_1} \cdots c_{u_n}$ for the diagonal entries.
Let 
$\mathbb{P}^1$ be the real projective line, which is the set of all lines through the origin in 
$\mathbb R^2$. We say that a proper subset 
$\mathcal{C} \subset \mathbb{P}^1$ is a cone if it is a closed projective interval and a multicone if it is a finite union of cones. The collection 
$ (A_u)_{u \in \mathcal U}$ of diagonal matrices as in (1.1) is called dominated if there exists a multicone 
$\mathcal{C} \subset \mathbb{P}^1$ such that 
$\bigcup_{u \in \mathcal U} A_u \mathcal{C} \subset \mathring{\mathcal{C}}$. It was shown in [Reference Bochi and Gourmelon15, theorem B] that 
$ (A_u)_{u \in \mathcal U}$ is dominated if and only if there exist constants C > 0 and 
$0 \lt \tau \lt 1$ such that
\begin{equation}
\frac{b_{\boldsymbol{u}}\cdot c_{\boldsymbol{u}}}{\max\{b_{\boldsymbol{u}},c_{\boldsymbol{u}}\}^2} \leqslant C \tau^n
\end{equation} for all 
$n \in \mathbb N$ and 
$\boldsymbol{u} \in \mathcal U^n$.
For each diagonal matrix Au as in (1.1), the singular value function is given by
\begin{align*} \varphi^s (A_u) := \begin{cases}
\max \{b_u, c_u \}^s , & \text{if } 0 \leqslant s \lt 1,\\
\max \{b_u, c_u\} \min \{b_u, c_u \}^{s-1}, & \text{if } 1 \leqslant s \lt 2.
\end{cases}\end{align*} The Lyapunov exponents of the collection 
$(A_u)_{u \in \mathcal U}$ with respect to a measure 
$\mu \in \mathcal{M}(\mathcal U^{\mathbb N}, \sigma)$ are defined as
\begin{align*}
\begin{aligned}
&\chi_1(\mu):= -\lim _{n\to\infty}\frac1n\int_{\mathcal U^{\mathbb N}}\log\max\left\{b_{\xi|_n},c_{\xi|_n}\right\}\,\mathrm{d}\mu(\xi),\\
&\chi_2(\mu):= -\lim _{n\to\infty}\frac1n\int_{\mathcal U^{\mathbb N}}\log\min\left\{b_{\xi|_n},c_{\xi|_n}\right\}\,\mathrm{d}\mu(\xi).
\end{aligned}
\end{align*} The Lyapunov dimension of 
$\mu \in \mathcal{M}(\mathcal U^{\mathbb N}, \sigma)$ is defined to be
\begin{align*} \operatorname{dim}_{\mathrm{L}}(\mu):= \min\left\{\frac{h_{\mu}(\sigma)}{\chi_1(\mu)}, 1+\frac{h_{\mu}(\sigma)-\chi_1(\mu)}{\chi_2(\mu)}\right\}.\end{align*} For a continuous potential 
$\Phi: \mathcal U^{\mathbb N} \rightarrow \mathbb R^d$, 
$d \geqslant 1$, write 
$S_n \Phi=\sum_{k=0}^{n-1} \Phi \circ \sigma^k$ for its Birkhoff sum. The topological pressure of Φ and 
$ (A_u)_{u \in \mathcal U}$ is given by
\begin{align*} P\left(\log \varphi^s+\Phi\right)=\lim _{n \rightarrow \infty} \frac{1}{n} \log \sum_{\boldsymbol{u} \in \mathcal U^n} \varphi^s(A_{\boldsymbol{u}}) \sup _{\xi \in[\boldsymbol{u}]} \exp \left(S_n \Phi(\xi)\right),\end{align*}where the existence of the limit is guaranteed by the sub-additivity of the potential.
2.3. Hausdorff dimension
For a subset 
$F \subseteq \mathbb R^n$, 
$n \geqslant 1$, and δ > 0, a δ-cover of F is a collection 
$\{U_i\}$ of subsets of 
$\mathbb R^n$ that each have diameter 
$\textrm{diam}(U_i)$ at most δ and satisfy 
$F \subseteq \bigcup_{i} U_i$. For s > 0, the s-dimensional Hausdorff outer measure is defined as
\begin{align*} \mathcal H^s(F) = \lim_{d \downarrow 0} \inf \left\{\sum_{i} \textrm{diam}(U_i)^s \, : \, \{U_i\} \text{is a $\delta$-cover of $F$} \right\}.\end{align*}The Hausdorff dimension of the set F is
\begin{align*} \operatorname{dim}_{\mathrm{H}}(F) = \inf\{s \gt 0 \, : \, \mathcal H^s(F) =0\}.\end{align*}Let µ be a finite Borel measure on F. The Hausdorff dimension of µ is
\begin{align*} \operatorname{dim}_{\mathrm{H}}(\mu) = \inf\{\operatorname{dim}_{\mathrm{H}}(Z) \, : \, \mu(Z) =1\}.\end{align*} The lower pointwise dimension of µ at a point 
$x \in F$ is defined by
\begin{align*} \underline{d}_{\mu}(x) = \liminf_{r \to 0} \frac{\log \mu(B(x,r))}{\log r},\end{align*} where 
$B(x,r)$ denotes the open ball in 
$\mathbb R^n$ with radius r centred at x. The following result can be found in, e.g., [Reference Pesin50, theorem 7.1 and theorem 7.2].
Lemma 2.1. Let 
$F \subseteq \mathbb R^n$ be a Borel set and µ a finite Borel measure on 
$\mathbb R^n$. The following statements hold.
(i) If
$\underline{d}_{\mu}(x) \leqslant c$ for some c > 0 and every $x \in F$, then $\operatorname{dim}_{\mathrm{H}}(F) \leqslant c$.(ii) If
$\underline{d}_{\mu}(x) \geqslant c$ for some c > 0 and µ-a.e. $x \in F$, then $\operatorname{dim}_{\mathrm{H}}(\mu) \geqslant c$.
3. Dominated diagonally affine IFSs
In this section, we prove theorem 1.1. Recall the definition of the natural projection 
$\pi: \mathcal U^{\mathbb N} \to \Lambda $ from (1.4). Also recall the definitions of the sets 
$E_{\Phi}(\boldsymbol{\alpha})$ and 
$L_\Phi$ from (1.3) and (1.5), respectively. We have the following upper bound for the Hausdorff dimension of 
$\pi(E_\Phi(\boldsymbol{\alpha}))$. In the proof, we make use of [Reference Bárány, Jordan, Käenmäki and Rams2, proposition 3.2], which holds for general affine IFSs (including the diagonally affine case) on 
$\mathbb R^2$.
Lemma 3.1. Let 
$\left\{A_u+ \boldsymbol{v}_u\right\}_{u \in \mathcal U}$ be a diagonally IFS on 
$\mathbb R^2$ such that the collection 
$ (A_u)_{u \in \mathcal U}$ is dominated. Let 
$\Phi: \mathcal U^{\mathbb N} \rightarrow \mathbb{R}^d$, 
$d \geqslant 1$, be a continuous potential. Then for each 
$\boldsymbol{\alpha} \in \mathring{L}_{\Phi}$,
\begin{align*} \begin{aligned}
\operatorname{dim}_{\mathrm{H}}\left(\pi(E_{\Phi}(\boldsymbol{\alpha}))\right) & \leqslant \sup \bigg\{\operatorname{dim}_{\mathrm{L}}(\mu)\, :\, \mu \in \mathcal{M}(\mathcal U^{\mathbb N}, \sigma) \text {and } \int_{\mathcal U^{\mathbb N}} \Phi\, \mathrm{d} \mu=\boldsymbol{\alpha} \bigg\} \\
& =\sup \bigg\{s \geqslant 0\, :\, \inf _{q \in \mathbb{R}^d} P\left(\log \varphi^s+\langle q, \Phi-\boldsymbol{\alpha}\rangle\right) \geqslant 0\bigg\}.
\end{aligned}
\end{align*}Proof. It follows directly from [Reference Bárány, Jordan, Käenmäki and Rams2, lemma 3.1 and proposition 3.2] that for any diagonally affine IFS 
$\{A_u+ \boldsymbol{v}_u\}_{u\in \mathcal U}$ on 
$\mathbb R^2$ and continuous potential 
$\Phi: \mathcal U^{\mathbb N} \to \mathbb R^d$, 
$d\geqslant 1$, and any 
$\boldsymbol{\alpha} \in \mathring{L}_{\Phi}$,
\begin{equation}
\dim_{\mathrm{H}} (\pi(E_{\Phi}(\boldsymbol{\alpha}))) \leqslant \sup\left\{s \geqslant 0 \, : \, \inf_{q \in \mathbb R^{d}} P (\log \varphi^s + \langle q,\Phi-\boldsymbol{\alpha} \rangle) \geqslant 0 \right\}.
\end{equation} A measure 
$\nu \in \mathcal{M}(\mathcal U^{\mathbb N}, \sigma^{n})$ is called an n-step Bernoulli measure if it is a Bernoulli measure on 
$\left(\mathcal U^{\mathbb N}, \sigma^{n} \right)$. For n-step Bernoulli measures 
$\nu \in \mathcal M (\mathcal U^{\mathbb N}, \sigma^{n})$, write
\begin{equation}
\tilde \nu = \frac1n \sum_{k=0}^{n-1} \nu \circ \sigma^{-k}.
\end{equation} Then 
$\tilde \nu \in \mathcal M (\mathcal U^{\mathbb N}, \sigma)$ and 
$\tilde \nu$ is ergodic. Since 
$( A_u)_{u \in \mathcal U}$ is dominated, it follows by [Reference Bárány, Jordan, Käenmäki and Rams2, proposition 4.3] and (3.1) that for any 
$\boldsymbol{\alpha} \in \mathring{L}_{\Phi}$,
\begin{align*}
\begin{split}
\sup & \left\{s \geqslant 0: \inf _{q \in \mathbb{R}^d} P\left(\log \varphi^s+\langle q, \Phi-\boldsymbol{\alpha}\rangle\right) \geqslant 0\right\} \\
\leqslant \ & \sup \bigg\{\dim_{\mathrm{L}}(\tilde{\nu}): \nu\,\text {fully supported } n\text{-step Bernoulli and } \int_{\mathcal U^{\mathbb N}} \Phi \, \mathrm{d} \tilde{\nu}=\boldsymbol{\alpha} \bigg\} \\
\leqslant \ & \sup \left\{\dim_{\mathrm{L}}(\mu): \mu \in \mathcal{M}(\mathcal U^{\mathbb N}, \sigma)\text {and } \int_{\mathcal U^{\mathbb N}} \Phi \, \mathrm{d} \mu=\boldsymbol{\alpha}\right\}.
\end{split}
\end{align*} On the other hand, if we let 
$\mu \in \mathcal{M}(\mathcal U^{\mathbb N}, \sigma)$ be such that 
$\int_{\mathcal U^{\mathbb N}} \Phi\, \mathrm{d} \mu=\boldsymbol{\alpha}$, then for any 
$0 \leqslant t \lt \operatorname{dim}_{\mathrm{L}}(\mu)$, it holds by the sub-additive variational principle (see [Reference Cao, Feng and Huang19]) that for all 
$q \in \mathbb{R}^d$,
\begin{align*} P\left(\log \varphi^t+\langle q, \Phi-\boldsymbol{\alpha}\rangle\right) \geqslant h_{\mu}(\sigma)+\lim _{n \rightarrow \infty} \frac{1}{n} \int_{\mathcal U^{\mathbb N}} \log \varphi^t\left(A_{\left.\xi\right|_n}\right)\, \mathrm{d} \mu(\xi) \gt 0.\end{align*} Hence, 
$t \leqslant \sup \left\{s \geqslant 0: \inf _{q \in \mathbb{R}^d} P\left(\log \varphi^s+\langle q, \Phi-\boldsymbol{\alpha}\rangle\right) \geq\right.$ 
$0\}$ and thus,
\begin{align*} \begin{split}
\sup \bigg\{\dim_{\mathrm{L}}(\mu) \,:\, \mu \in \mathcal{M}(\mathcal U^{\mathbb N}, \sigma) & \text{and } \int_{\mathcal U^{\mathbb N}} \Phi \, \mathrm{d} \mu=\boldsymbol{\alpha}\bigg\} \\
\leqslant \ & \sup \left\{s \geqslant 0: \inf _{q \in \mathbb{R}^d} P\left(\log \varphi^s+\langle q, \Phi-\boldsymbol{\alpha}\rangle\right) \geqslant 0\right\} .
\end{split}\end{align*}This gives the result.
Remark 3.1. Note that the proof of lemma 3.1 shows that in fact
\begin{align*} \begin{split}
\sup & \left\{s \geqslant 0: \inf _{q \in \mathbb{R}^d} P\left(\log \varphi^s+\langle q, \Phi-\boldsymbol{\alpha}\rangle\right) \geqslant 0\right\} \\
= \ & \sup \bigg\{\dim_{\mathrm{L}}(\tilde{\nu}): \nu\,\text {fully supported } n\text{-step Bernoulli and } \int_{\mathcal U^{\mathbb N}} \Phi \, \mathrm{d} \tilde{\nu}=\boldsymbol{\alpha} \bigg\} \\
= \ & \sup \left\{\dim_{\mathrm{L}}(\mu): \mu \in \mathcal{M}(\mathcal U^{\mathbb N}, \sigma)\text {and } \int_{\mathcal U^{\mathbb N}} \Phi \, \mathrm{d} \mu=\boldsymbol{\alpha}\right\}.
\end{split}\end{align*} Under the additional conditions mentioned in the statement of theorem 1.1, we can prove that this upper bound in fact equals the Hausdorff dimension of the level set. Note that it would be possible to combine Hochman [Reference Hochman29] and Jordan and Simon [Reference Jordan and Simon35] to obtain a similar result for almost all vectors vu but our theorem 1.1 is proved for all vectors vu. Bárány et al. [Reference Bárány, Jordan, Käenmäki and Rams2] proved a similar result for affine IFSs satisfying the SOSC under the assumption that the set of matrices 
$ (A_u)_{u \in \mathcal U}$ is strongly irreducible such that the generated subgroup of the normalized matrices is not relatively compact. Theorem 1.1 is inspired by their result.
Proof. Proof of theorem 1.1
Let 
$\{A_u + \boldsymbol{v}_u\}_{u \in \mathcal U} \in \mathcal D$. For each 
$u \in \mathcal U$, it holds that
\begin{align*} \frac{b_u \cdot c_u}{\max\{b_u, c_u\}^2} = \frac{\min \{b_u, c_u\} }{\max\{b_u, c_u\}} \lt 1,\end{align*} since either 
$|b_u| \gt |c_u|$ for all 
$u \in \mathcal U$ or 
$|b_u| \lt |c_u|$ for all 
$u \in \mathcal U$. Take 
$\tau = \max_{u \in \mathcal U} \left\{\frac{\min \{b_u, c_u\} }{\max\{b_u, c_u \}} \right\}$. Then τ < 1 and, since each Au is a diagonal matrix, we get (2.2) with C = 1. Hence, 
$\{A_u\}_{u \in \mathcal U}$ is dominated and therefore the desired upper bound for the Hausdorff dimension of 
$\pi(E_{\Phi}(\boldsymbol{\alpha}))$ is given by lemma 3.1.
For the lower bound, suppose that 
$\nu \in \mathcal M(\mathcal U^{\mathbb N}, \sigma^{n})$ is a fully supported n-step Bernoulli measure with 
$\int_{\mathcal U^{\mathbb N}} \Phi\, \mathrm{d} \tilde{\nu}=\boldsymbol{\alpha}$ with 
$\tilde \nu$ as defined in (3.2). The existence of the measure ν is guaranteed by [Reference Bárány, Jordan, Käenmäki and Rams2, proposition 4.3]. Then 
$\tilde \nu \in \mathcal M (\mathcal U^{\mathbb N}, \sigma)$ and 
$\tilde \nu$ is ergodic and therefore from 
$\int_{\mathcal U^{\mathbb N}} \Phi\, \mathrm{d} \tilde{\nu}=\boldsymbol{\alpha}$, we get that
\begin{equation}
\tilde{\nu}\left(\left\{\xi \in \mathcal U^{\mathbb N}: \lim _{n \rightarrow \infty} \frac{1}{n} S_n \Phi (\xi)=\boldsymbol{\alpha}\right\} \right)=1.
\end{equation} Let 
$\hat{\nu}=\tilde{\nu} \circ \pi^{-1} $. Assume that 
$|b_u| \gt |c_u|$ for all 
$u \in \mathcal U$ so that we are in the situation of condition (D) (the proof for the case (Dʹ) goes similarly). Then the strong stable direction of the collection 
$(A_u)_{u \in \mathcal U}$ is equal to the subspace parallel to the y-axis (see [Reference Bochi and Gourmelon15]). Let 
$P_x \hat{\nu}$ be the measure on 
$[0,1]$ given by the canonical projection onto the x-coordinate of 
$\hat \nu$. Since the matrices Au are diagonal, 
$P_x \hat{\nu}$ is a self-similar measure for the IFS 
$\mathcal G_1$, i.e. there is a probability vector 
$\hat{\boldsymbol{p}} = (\hat p_u)_{u \in \mathcal U}$ such that
\begin{align*}P_x \hat{\nu} (B)=\sum_{u \in \mathcal U} \hat p_u P_x \hat{\nu} (g_{1,u}^{-1}(B))\end{align*} for each Borel set 
$B \subseteq [0,1]$. Then condition (D)(a) together with [Reference Hochman29, theorem 1.1] or condition (D)(b) together with [Reference Rapaport51, theorem 1.2] yields
\begin{equation}
\dim_{\mathrm{H}} (P_x \hat{\nu})=\min\bigg\{1, \frac{h_{\tilde{\nu}}(\sigma)}{\chi_1(\tilde{\nu})}\bigg\}.
\end{equation}It then follows from [Reference Bárány and Käenmäki3, corollaries 2.7 and 2.8] and (3.4) that
\begin{align*} \dim_{\mathrm{H}} (\hat{\nu})=\operatorname{dim}_{\mathrm{L}}(\tilde{\nu}).\end{align*} This and (3.3) yield 
$\dim_{\mathrm{H}}\left(\pi(E_{\Phi}(\boldsymbol{\alpha}))\right) \geqslant \operatorname{dim}_{\mathrm{L}}(\tilde{\nu})$. Since this holds for arbitrary fully supported n-step Bernoulli measures ν with 
$\int_{\mathcal U^{\mathbb N}} \Phi\, \mathrm{d} \tilde{\nu}=\boldsymbol{\alpha}$, the result follows from remark 3.1.
Remark 3.2. We make a small remark on the conditions (D) and (Dʹ). It was shown by Hochman in [Reference Hochman29, proof of theorem 1.5] that an IFS satisfies the ESC if it does not have exact overlaps and all parameters 
$b_u, c_u,\beta_u, \gamma_u$ are algebraic numbers over 
$\mathbb Q$. In [Reference Bárány and Käenmäki4, Reference Baker5, Reference Chen20], it was shown that there exist IFSs that do not contain exact overlaps while there are cylinders which are super-exponentially close at all small scales, i.e. the ESC does not hold. What is needed in the proof of theorem 1.1 is (3.4), which is also guaranteed by [Reference Rapaport51] under the assumption of having algebraic 
$b_u,c_u$ and no exact overlaps.
4. Digit frequencies for finite GLS expansions
We now move to the second type of IFS we consider. Fix a GLS IFS 
$\{A_e +\boldsymbol{v}_e\}_{e \in \mathcal E}$. We start by proving some properties of the expansions from (1.11).
4.1. Multiple representations
First, consider the representations of the points 
$y \in [0,1]$. Recall the definition of the maps fj from (1.15). The IFS 
$\{f_j \}_{0 \leqslant j \lt J}$ satisfies the SOSC and has the interval 
$[0,1]$ as its attractor. Let 
$\pi_J: \{0,1, \ldots, J-1\}^{\mathbb N} \to [0,1]$ be the map given by
\begin{align*} \pi_J((j_m)_{m \geqslant 1}) = \lim_{m \to \infty}f_{j_1} \circ f_{j_2} \circ \cdots \circ f_{j_m}(0).\end{align*} One easily sees that to all but countably many 
$y \in [0,1]$, there corresponds a unique sequence 
$\zeta \in \{0,1, \ldots, J-1\}^{\mathbb N}$ such that 
$y = \pi_J(\zeta)$ and otherwise 
$\#\pi^{-1}_J\{y\} =2$ and there is one sequence ending in an infinite string of 0’s and one ending in an infinite string of 
$(J-1)$’s. We make the following observation, which we will use later on. Recall the definition of the set 
$W(\boldsymbol{\alpha})$ from (1.14).
Lemma 4.1. Let 
$\boldsymbol{\alpha} = (\alpha_e)_{e \in \mathcal E} \in [0,1]^{\# \mathcal E}$ be a frequency vector with 
$\alpha_j \gt 0$ for each 
$0 \leqslant j \lt J$. Then 
$\#\pi_J^{-1}\{y\}=1$ for any 
$y \in W(\boldsymbol{\alpha})$.
Proof. Let 
$y \in [0,1]$ be such that 
$\#\pi_J^{-1}\{y\}=2$ and let 
$x \in [0,1]$. Then any 
$\omega = (j_m,k_m)_{m \geqslant 1} \in \pi^{-1}\{(x,y)\}$ either has 
$j_m=0$ for all m large enough or 
$j_m = J-1$ for all m large enough. In the first case, 
$\sum_{k=0}^{B_j-1}\tau_{(j,k)}(\omega) =0 \neq \alpha_j$ for all j ≠ 0 and in the second case, 
$\sum_{k=0}^{B_j-1} \tau_{(j,k)}(\omega) =0 \neq \alpha_j$ for all 
$j \neq J-1$. Hence, 
$(x,y) \not \in F(\boldsymbol{\alpha})$ and thus 
$W(\boldsymbol{\alpha}) = \emptyset$.
For a fixed 
$y \in [0,1]$, we can consider the expansions one obtains from the GLS IFS for 
$x \in [0,1]$. We define the fibre fundamental intervals corresponding to y by setting for each 
$m \geqslant 1$ and 
$k_1, \ldots, k_m$ satisfying 
$0 \leqslant k_i \lt B_{j_i}$ for all 
$1 \leqslant i \leqslant m$,
\begin{equation}
\Delta_y(k_1, \ldots, k_m) := h_{(j_1,k_1)} \circ \cdots \circ h_{(j_m,k_m)} ([0,1]),
\end{equation} where we let 
$(j_m)_{m \geqslant 1}$ be the lexicographically smallest sequence in 
$ \pi_J^{-1}\{y\}$. For 
$y \in W(\boldsymbol{\alpha})$, this means that 
$(j_m)_{m \geqslant 1} \in \pi_J^{-1}\{y\}$ is the unique sequence that satisfies 
$\tau_j((j_m)_{m \geqslant 1})=\alpha_j$ for each 
$0 \leqslant j \lt J$ and for 
$y \in [0,1]\setminus W(\boldsymbol{\alpha})$ the sequence 
$(j_m)_{m \geqslant 1}$ is the one ending in an infinite string of 
$(J-1)$’s.
If we fix 
$y\in W(\boldsymbol{\alpha})$ and take 
$x\in[0,1]$ such that 
$\#\pi^{-1}\{(x,y)\} \gt 1$, then we know by lemma 4.1 that 
$\#\pi_J^{-1}\{y\}=1$, say 
$y=\pi_J((j_m)_{m\geq1})$. Consequently, x must have multiple expansions along the fibre y and so must lie on the boundary of a fibre fundamental interval 
$\Delta_y(k_1,\dotsc,k_m)$ for some 
$0\leqslant k_i \lt B_{j_i}$ (
$1\leqslant i\leqslant m$) and some 
$m\in\mathbb{N}$. Since each fibre fundamental interval has two boundary points and there are only countably many fibre fundamental intervals, the set of such points x must be countable.
Now, fix an 
$x \in [0,1]$. Since the GLS IFS 
$\{A_e + \boldsymbol{v}_e\}_{e \in \mathcal E}$ has 
$[0,1]^2$ as its attractor, to any 
$y \in [0,1]$, there corresponds a sequence 
$\omega \in \mathcal E^{\mathbb N}$ such that 
$\pi(\omega)=(x,y)$. Therefore, to any sequence 
$(j_m)_{m \geqslant 1} \in \{0,1, \ldots, J-1\}^{\mathbb N}$, there corresponds a sequence 
$(k_m)_{m \geqslant 1}$ with 
$0 \leqslant k_m \lt B_{j_m}-1$ for each 
$m \in \mathbb N$ such that
\begin{align*} x= \lim_{m \to \infty} h_{(j_1,k_1)} \circ \cdots \circ h_{(j_m,k_m)}(0).\end{align*} We show that if 
$h_e \neq h_e'$ whenever 
$e \neq e'$, then each of the sequences 
$(j_m)_{m \geqslant 1} \in \{0,1, \ldots, J-1\}^{\mathbb N}$ yields a different GLS expansion for x as in (1.11).
The GLS expansions from (1.11) are given by the triples of digits 
$(s_m,K_m,t_m)$, 
$m \in \mathbb N$, from (1.10). Therefore, if we set
\begin{align*} \mathcal A = \left\{\left(\varepsilon_e, q_e^{-1}, r_e+ \varepsilon_e q_e \right) \, : \, e \in \mathcal E \right\},\end{align*} then we can think of 
$\mathcal A$ as the GLS digit set corresponding to 
$\{A_e + \boldsymbol{v}_e\}_{e \in \mathcal E}$ and we can map sequences 
$ ((j_m,k_m))_{m \geqslant 1} \in \mathcal E^{\mathbb N}$ to sequences 
$(s_m,K_m,t_m)_{m \geqslant 1} \in \mathcal A^{\mathbb N}$ through the identification given in (1.10). Let 
$(j_m)_{m \geqslant 1}, (j_m')_{m \geqslant 1} \in \{0,1, \ldots, J-1\}^{\mathbb N}$ be two different sequences, so there is an 
$m \in \mathbb N$ such that 
$j_m \neq j_m'$. Let 
$(j_m,k_m)_{m \geqslant 1}, (j_m',k_m')_{m \geqslant 1} \in \mathcal E^{\mathbb N}$ be two sequences that both project to x in the second coordinate under π. Since 
$j_m \neq j_m'$, it holds that 
$(j_m,k_m) \neq (j_m',k_m')$. If we assume that 
$h_e \neq h_e'$ whenever 
$e \neq e'$, then it would follow that 
$\varepsilon_{(j_m,k_m)} \neq \varepsilon_{(j_m',k_m')}$ or 
$q_{(j_m,k_m)} \neq q_{(j_m',k_m')}$ and thus that the digits from 
$\mathcal A$ corresponding to 
$(j_m,k_m)$ and 
$(j_m',k_m')$ differ. Therefore, we immediately find the following result.
Proposition 4.1. Let 
$\{A_e +\boldsymbol{v}_e\}_{e \in \mathcal E}$ be a GLS IFS with the additional assumption that 
$h_e \neq h_{e'}$ whenever 
$e \neq e'$. Then for each 
$x \in [0,1]$, there are uncountably many different digit sequences 
$(s_m,K_m,t_m)_{m \geqslant 1} \in \mathcal A^{\mathbb N}$ with
\begin{align*} x =\sum_{m\geqslant 1} (-1)^{\sum_{i=1}^{m-1}s_i} \, \frac{t_m}{\prod_{i=1}^m{K_i}}.\end{align*} The above also shows that there is a one-to-one correspondence between the sequences in 
$\mathcal E^{\mathbb N}$ and in 
$ \mathcal A^{\mathbb N}$, which justifies considering the elements of 
$\mathcal E$ as digits in the GLS expansions.
4.2. Non-empty level sets
In this section, we determine for which frequency vectors α the level set 
$F(\boldsymbol{\alpha})$ from (1.12) and the set 
$W(\boldsymbol{\alpha})$ from (1.14) are non-empty. We first consider the level sets 
$F(\boldsymbol{\alpha})$.
Proposition 4.2. The set 
$F(\boldsymbol{\alpha})$ is non-empty for any frequency vector 
$\boldsymbol{\alpha} = (\alpha_e)_{e \in \mathcal E} \in [0,1]^{\# \mathcal E}$.
Proof. Fix a frequency vector 
$\boldsymbol{\alpha} =(\alpha_e)_{e \in \mathcal E} \in [0,1]^{\# \mathcal E}$. It is sufficient to construct a sequence 
$\omega = (\omega_n)_{n \geqslant 1} \in \mathcal E^{\mathbb N}$ such that 
$\tau_e(\omega)=\alpha_e$ for each 
$e \in \mathcal E$ since then 
$\pi(\omega) \in F(\boldsymbol{\alpha})$. Denote by 
$\lfloor\cdot\rceil$ the nearest-integer function. Order the elements in 
$\mathcal E$ by setting 
$(j,k) \prec (j',k')$ if either 
$j \lt j'$ or if 
$j=j'$ and 
$k \lt k'$. For each 
$n \geqslant 1$, set
\begin{align*} E_n = \left\{e \in \mathcal E \, : \, \lfloor n \alpha_e \rceil = \lfloor (n-1) \alpha_e \rceil +1 \right\} = \left\{e_{n,1} \prec \cdots \prec e_{n,m_n} \right\},\end{align*} where En can be empty and thus 
$m_n=0$ for some n. Define 
$\omega \in \mathcal E^{\mathbb N}$ by setting for each 
$n \geqslant 1$ and 
$1 \leqslant m \leqslant m_n$,
\begin{align*} \omega_{m+ \sum_{i=1}^{n-1} m_i} = e_{n,m},\end{align*} where we let 
$\sum_{i=1}^0 m_i=0$. Clearly, there are infinitely many n for which 
$E_n \neq \emptyset$ so ω is well-defined.
Now observe that for each n the number of terms of ω we have defined using 
$\bigcup_{i=1}^n E_i$ is
\begin{align*} \sum_{i=1}^n m_i = \sum_{e\in \mathcal E}\lfloor n\alpha_e \rceil.\end{align*} Since 
$\sum_{e\in \mathcal E}n\alpha_e=n$ and 
$\mathfrak m:=\# \mathcal E \lt \infty$, we must have
\begin{align*} n- \mathfrak m \leqslant \sum_{i=1}^n m_i \leqslant n+\mathfrak m.\end{align*} Thus for each 
$e \in \mathcal E$,
\begin{align*} \begin{split}
\frac{\#\{0\leqslant m \lt n\, :\, \omega_m=e\}}{n} \leqslant \ & \frac{\#\{0\leqslant m \lt \mathfrak m + \sum_{i=1}^n m_i \, :\, \omega_m=e\}}{n}\\
\leqslant \ & \frac{\#\{0\leqslant m \lt \sum_{i=1}^n m_i \, : \, \omega_m=e\} + \mathfrak m }{n}\\
\leqslant \ & \frac{n \alpha_e +1 + \mathfrak m }{n}.
\end{split}\end{align*}Similarly, it holds that
\begin{align*} \frac{\#\{0\leqslant m \lt n:\omega_m=e\}}{n} \geqslant \frac{n \alpha_e -1-\mathfrak m}{n}.\end{align*} Taking the limit as 
$n \to \infty$ yields 
$\tau_e(\omega) = \alpha_e$ for all 
$e \in \mathcal E$.
For a fixed frequency vector α, we would like to determine the set 
$W(\boldsymbol{\alpha})$ of points 
$y\in [0,1]$ for which there exists an 
$x\in [0,1]$ such that the point (x, y) has digit frequencies given by α, see (1.14). Recall the definition of the Borel measure 
$\nu_{\boldsymbol{\alpha}}=\mu_{\boldsymbol{\alpha}}\circ\pi_2^{-1}$ from §1. We have the following result.
Proposition 4.3. Let 
$\boldsymbol{\alpha} = (\alpha_e)_{e \in \mathcal E} \in [0,1]^{\# \mathcal E}$. Then
\begin{equation*}W({\boldsymbol{\alpha}})=\left\{y\in[0,1]\, :\, \exists \, \omega \in \pi_2^{-1} \{y\} \, \text{s.t. } \, \sum_{0 \leqslant k \lt B_j} \tau_{(j,k)}(\omega)= \alpha_j \,\text{for all } 0 \leqslant j \lt J\right\}.\end{equation*} In particular, 
$\nu_{\boldsymbol{\alpha}}(W({\boldsymbol{\alpha}}))=1$.
Proof. 
$(\subset)$ Set
\begin{align*} W = \left\{y\in[0,1]\, :\, \exists \, \omega \in \pi_2^{-1} \{y\} \, \text{s.t. } \, \sum_{0 \leqslant k \lt B_j} \tau_{(j,k)}(\omega)= \alpha_j \,\text{for all } 0 \leqslant j \lt J\right\}.\end{align*} First, let 
$y \in W(\boldsymbol{\alpha})$. This means that 
$F_y({\boldsymbol{\alpha}})\neq\emptyset$ so there exists an 
$\omega\in\mathcal E^{\mathbb N}$ such that 
$\pi(\omega)=(x,y)$ for some 
$x\in[0,1]$ and 
$\tau_e(\omega)=\alpha_e$ for all 
$e\in\mathcal E$. Consequently, we find that
\begin{equation*}\sum_{0\leqslant k \lt B_j}\tau_{(j,k)}(\omega)=\sum_{0\leqslant k \lt B_j}\alpha_{(j,k)} = \alpha_j\end{equation*} for all 
$0 \leqslant j \lt J$ and so 
$y \in W$.
$(\supset)$ Conversely, let 
$y\in W$. Then there exists an 
$\omega' = ( j_\ell', k_\ell')_{\ell \geqslant 1} \in \pi_1^{-1} \{y\}$ for which it holds that 
$\sum_{0 \leqslant k \lt B_j} \tau_{(j,k)}(\omega')= \alpha_j$ for all 
$0 \leqslant j \lt J$. Write 
$\zeta=(j_\ell')_{\ell \geqslant 1} \in \{0, \ldots, J-1\}^{\mathbb N}$. For each 
$0 \leqslant j \lt J$ and 
$n \geqslant 1$, set
\begin{align*} E_n^{(j)}:= \left\{(j,k)\in\mathcal E \, :\, \left\lfloor \frac{n\alpha_{(j,k)}}{\alpha_j}\right\rceil=\left\lfloor\frac{(n-1)\alpha_{(j,k)}}{\alpha_j}\right\rceil+1\right\}=\left\{e_{n,1}^{(j)}\prec\dotsb\prec e_{n,m_n^{(j)}}^{(j)}\right\}\end{align*} and let 
$\omega^{(j)} \in \mathcal E^{\mathbb N}$ be the sequence obtained from concatenating all elements from the sets 
$E_n^{(j)}$ as in proposition 4.2, so
\begin{align*} \omega^{(j)} = e_{1,1}^{(j)} \cdots e_{1,m_1^{(j)}}^{(j)} e_{2,1}^{(j)} \cdots e_{2,m_2^{(j)}}^{(j)} e_{3,1}^{(j)} \cdots.\end{align*} Then as in proposition 4.2, we obtain that for each 
$0 \leqslant k \lt B_j$,
\begin{align*} \lim_{n \to \infty} \frac{\# \{1 \leqslant m \leqslant n \, : \, \omega^{(j)}_m = (j,k) \} }{n} = \frac{\alpha_{(j,k)}}{\alpha_j}.\end{align*} We now weave the sequences 
$\omega^{(j)}$ together to construct a sequence 
$\omega = (j_\ell', k_\ell)_{\ell \geqslant 1} \in \mathcal E^{\mathbb N}$ that satisfies 
$\tau_e(\omega) = \alpha_e$ for all 
$e \in \mathcal E$. Then 
$\pi(\omega) = (x,y)$ for some 
$x \in F_y(\boldsymbol{\alpha})$, which shows that 
$F_y(\boldsymbol{\alpha}) \neq \emptyset$. For each 
$\ell \geqslant 1$, let 
$\omega_\ell$ be the 
$\tau_{j_\ell'}(\zeta,\ell)$th element of the sequence 
$\omega^{(j_\ell')}$. So, 
$\omega_1 = \omega^{(j'_1)}_1 = e_{1,1}^{(j'_1)}$, ω 2 either equals 
$\omega^{(j_1')}_2$ if 
$j'_1=j'_2$ or 
$\omega^{(j_2')}_1$ if 
$j'_1 \neq j'_2$, etc. As the sequences in the first coordinates of 
$\omega = (j_\ell', k_\ell)_{\ell \geqslant 1}$ and 
$\omega' = (j_\ell', k_\ell')_{\ell \geqslant 1}$ coincide, we have for each 
$0 \leqslant j \lt J$,
\begin{align*} \sum_{0 \leqslant k \lt B_j} \tau_{(j,k)}(\omega) = \sum_{0 \leqslant k \lt B_j} \tau_{(j,k)}(\omega')= \alpha_j. \end{align*} Moreover, for each 
$e = (j,k) \in \mathcal E$ and 
$n \geqslant 1$,
\begin{equation}
\# \{1 \leqslant m \leqslant n \, : \, \omega_m =e \} = \# \{1 \leqslant m \leqslant \tau_j(\zeta,n) \, : \, \omega^{(j)}_m = e\}.
\end{equation} If 
$\alpha_j \gt 0$, then 
$\tau_j(\zeta,n) \gt 0$ for all n large enough and for any 
$e = (j,k) \in \mathcal E$, we obtain
\begin{align*} \tau_e(\omega) & = \lim_{n \to \infty} \frac{\# \{1 \leqslant m \leqslant \tau_j(\zeta,n) \, : \, \omega^{(j)}_m = e\} }{\tau_j(\zeta,n)} \cdot \frac{\tau_j(\zeta,n)}{n}\\& = \frac{\alpha_{(j,k)}}{\alpha_j} \sum_{0 \leqslant k \lt B_j} \tau_{(j,k)}(\omega) = \alpha_{(j,k)}.\end{align*} If 
$ \alpha_j=0$, then 
$\alpha_{(j,k)}=0$ for each 
$0 \leqslant k \lt B_j$ and by (4.2),
\begin{align*} 0 \leqslant \tau_{(j,k)}(\omega) \leqslant \lim_{n \to \infty} \frac{\tau_j(\zeta,n)}{n} = \alpha_j =0.\end{align*}This gives the first part of the statement.
As 
$\left\{\omega \in \mathcal E^{\mathbb N} \, : \, \tau_e(\omega) = \alpha_e \,\text{for all } e \in \mathcal E \right\} \subseteq \pi_1^{-1}(W(\boldsymbol{\alpha}))$, it follows from the definition of 
$\mu_{\boldsymbol{\alpha}}$ that 
$\nu_{\boldsymbol{\alpha}}(W({\boldsymbol{\alpha}}))=1$.
4.3. The Hausdorff dimension of the Besicovitch–Eggleston sets
In this section, we prove theorem 1.2. The proof is similar to [Reference Barreira and Iommi9, theorem 3.1] and [Reference Fan, Liao, Ma and Wang26, theorem 1.1], which both treat digit frequencies for expansions with infinite digit sets that can be generated by an IFS on 
$\mathbb R$ as in (1.6). Their results do not apply to our setting because the IFS 
$\{h_e:[0,1] \to [0,1]\}_{e \in \mathcal E}$ on 
$\mathbb R$ is not of this type. Nevertheless, since we have a finite digit set, we can adapt the method of proof from [Reference Barreira and Iommi9, theorem 3.1].
Fix a 
$y \in [0,1]$. Recall the definition of the fibre fundamental intervals 
$\Delta_y (k_1, \ldots, k_m)$ from (4.1). Note that we obtain a semi-algebra of sets generating the Borel σ-algebra 
$\mathcal B([0,1])$ on 
$[0,1]$ by taking the collection of all intervals (open, closed, and half-open) that can be formed by the endpoints of the fibre fundamental intervals. Suppose that the frequency vector α satisfies the following additional property: for each 
$0 \leqslant j \lt J$, there are 
$k, \ell \in B_j$ with 
$k \neq \ell$ and 
$\alpha_{(j,k)} \gt 0$ and 
$\alpha_{(j, \ell)} \gt 0$. Let 
$m_{y,{\boldsymbol{\alpha}}}$ be the measure on 
$([0,1], \mathcal B([0,1]))$ determined by
\begin{equation*}m_{y,{\boldsymbol{\alpha}}}(\Delta_y(k_1, \ldots, k_m))=\prod_{i=1}^m \frac{\alpha_{(j_i,k_i)}}{\alpha_{j_i}}, \quad 0 \leqslant k_i \lt B_{j_i}, \, 1 \leqslant i \leqslant m, \, m \geqslant 1,\end{equation*} and by the same quantity for any interval determined by the same endpoints. This immediately implies that any endpoint x of a fibre fundamental interval has 
$m_{y,{\boldsymbol{\alpha}}}(\{x\})=0$. If 
$x\in [0,1]$ is not an endpoint of a fibre fundamental interval, then there is a sequence 
$(k_m)_{m \geqslant 1}$ such that
\begin{align*} \bigcap_{m \in \mathbb N} \Delta_y(k_1, \ldots, k_m) = \{x\}.\end{align*}This implies that
\begin{align*}m_{y,{\boldsymbol{\alpha}}}(\{x\}) = \lim_{m \to \infty} \prod_{i=1}^m \frac{\alpha_{(j_i,k_i)}}{\alpha_{j_i}}.\end{align*} By the additional assumption on α, there is a constant 
$0 \lt c \lt 1$ such that 
$\frac{\alpha_{j,k}}{\alpha_j} \lt c \lt 1$ for all 
$(j,k) \in \mathcal E$. Therefore, 
$\mu_{y,\boldsymbol{\alpha}}(\{x\})=0.$ We will need the following property of 
$m_{y,{\boldsymbol{\alpha}}}$.
Lemma 4.2. For 
$\nu_{\boldsymbol{\alpha}}$-a.e. 
$y\in W({\boldsymbol{\alpha}})$, it holds that 
$m_{y,{\boldsymbol{\alpha}}}(F_y({\boldsymbol{\alpha}}))=1$.
Proof. Observe that for each 
$(j,k)\in\mathcal E$ by proposition 4.3,
\begin{equation} \begin{split}
\int_{[0,1]}\int_{[0,1]} & \mathbb 1_{\pi([(j,k)])}(x,y)\,\mathrm{d}m_{y,{\boldsymbol{\alpha}}}(x)\,\mathrm{d}\nu_{\boldsymbol{\alpha}}(y)\\
&=\int_{[\sum_{i=0}^{j-1}p_i, \sum_{i=0}^j p_i]}
\int_{[0,1]}\mathbb 1_{\Delta_y(k)}(x)\,\mathrm{d}m_{y,{\boldsymbol{\alpha}}}(x)\,\mathrm{d}\nu_{\boldsymbol{\alpha}}(y)\\
&= \frac{\alpha_{(j,k)}}{\alpha_j}
\nu_{\boldsymbol{\alpha}}\left(\left[\sum_{i=0}^{j-1}p_i, \sum_{i=0}^j p_i\right] \right)\\
&=\alpha_{(j,k)} =\int_{[0,1]^2}\mathbb 1_{\pi([(j,k)])}\,\mathrm{d}\mu_{\boldsymbol{\alpha}}\circ\pi^{-1}.
\end{split}\end{equation} Since the collection 
$\{\pi([e_1, \ldots , e_n]): e_i\in\mathcal E, \,\, 1 \leqslant i \leqslant n \}$ generates the Borel σ-algebra 
$\mathcal B([0,1]^2)$, we can conclude from (4.3) that
\begin{equation}
\int_{[0,1]} \int_{[0,1]} f \, \mathrm{d}m_{y,{\boldsymbol{\alpha}}}\,\mathrm{d}\nu_{\boldsymbol{\alpha}} = \int_{[0,1]^2} f \, \mathrm{d}\mu_{\boldsymbol{\alpha}}\circ\pi^{-1}
\end{equation} for all 
$f\in L^1([0,1]^2,\mathcal B([0,1]^2),\mu_{\boldsymbol{\alpha}}\circ\pi^{-1})$.
Let 
$E := \{y \in W(\boldsymbol{\alpha}) \, : \, m_{y, \boldsymbol{\alpha}} (F_y(\boldsymbol{\alpha})) \lt 1 \}$ and suppose that 
$\nu_{\boldsymbol{\alpha}} (E) \gt 0$. From (4.4) with 
$f=\mathbb 1_{F({\boldsymbol{\alpha}})}$ together with proposition 4.3, we then find that
\begin{align*}
\mu_{\boldsymbol{\alpha}}\circ\pi^{-1}(F({\boldsymbol{\alpha}}))
&=\int_{W({\boldsymbol{\alpha}})}\int_{[0,1]}\mathbb 1_{F_y({\boldsymbol{\alpha}})}(x)\,\mathrm{d}m_{y,{\boldsymbol{\alpha}}}(x)\,\mathrm{d}\nu_{\boldsymbol{\alpha}}(y)\\
&=\int_{W({\boldsymbol{\alpha}})\setminus E}1\,\mathrm{d}\nu_{\boldsymbol{\alpha}}
+\int_E\int_{[0,1]}\mathbb 1_{F_y({\boldsymbol{\alpha}})}(x)\,\mathrm{d}m_{y,{\boldsymbol{\alpha}}}(x)\,\mathrm{d}\nu_{\boldsymbol{\alpha}}(y)\\
& \lt \int_{W({\boldsymbol{\alpha}})}1\,\mathrm{d}\nu_{\boldsymbol{\alpha}}\\
&=1.
\end{align*} On the other hand, recall that 
$E_{\mathbb 1}(\boldsymbol{\alpha})$ is the symbolic Besicovitch–Eggleston set containing all sequences 
$\omega \in \mathcal E^{\mathbb N}$ with 
$\tau_e (\omega)=\alpha_e$ for each 
$e \in \mathcal E$. Therefore, by the definition of 
$\mu_{\boldsymbol{\alpha}}$,
\begin{align*} 1 = \mu_{\boldsymbol{\alpha}}(E_{\mathbb 1}(\boldsymbol{\alpha})) \leqslant \mu_{\boldsymbol{\alpha}} \circ \pi^{-1} (F(\boldsymbol{\alpha})) \lt 1.\end{align*} This gives a contradiction. It follows that 
$m_{y,{\boldsymbol{\alpha}}}(F_y({\boldsymbol{\alpha}}))=1$ for 
$\nu_{\boldsymbol{\alpha}}$-a.e. 
$y\in W({\boldsymbol{\alpha}})$.
Before we move to the proof of theorem 1.2, to simplify notation, we put 
$p_{(j,k)}=p_j$ for all 
$(j,k)\in\mathcal E$. Also, let 
$\mathbb P_{\boldsymbol{\alpha}}$ be the 
$(\alpha_j)$-Bernoulli measure on 
$\{0, \ldots, J-1\}^{\mathbb N}$.
Proof. Proof of theorem 1.2
Fix a 
$y \in W(\boldsymbol{\alpha})$. Recall that the lower pointwise dimension of 
$m_{y,{\boldsymbol{\alpha}}}$ at the point 
$x \in [0,1]$ is defined by
\begin{equation*}\underline{d}_{m_{y,{\boldsymbol{\alpha}}}}(x) =\liminf_{r\to0}\frac{\log m_{y,{\boldsymbol{\alpha}}}(B(x,r))}{\log r},\end{equation*} where 
$B(x,r)$ is the open interval of length 2r centred at x. One can verify that the collection 
$\{\Delta_y(k_1 \cdots k_n) \, : \, n \in \mathbb N\}$ satisfies conditions (CB1)–(CB3) of the Moran-type construction from [Reference Pesin50, Section 15]. Moreover, for any intervals 
$\Delta_y(k_1 \cdots k_n)$, 
$\Delta_y(k_1 \cdots k_n, k_{n+1})$, it holds that
\begin{align*} \textrm{diam}(\Delta_y(k_1 \cdots k_n)) &\leqslant (\max_{e \in \mathcal E} q_e )^n, \quad (\min_{e \in \mathcal E} q_e) \cdot\textrm{diam}(\Delta_y(k_1 \cdots k_n))\\& \leqslant \textrm{diam}(\Delta_y(k_1 \cdots k_n k_{n+1})). \end{align*} Therefore, by, e.g., [Reference Pesin50, theorem 15.3(1)], we can replace the balls 
$B(x,r)$ in the definition of 
$\underline{d}_{m_{y,{\boldsymbol{\alpha}}}}$ with the fibre fundamental intervals 
$ \Delta_y (k_1, \ldots, k_n)$ to obtain an upper bound for 
$\underline{d}_{m_{y,{\boldsymbol{\alpha}}}}(x)$ for all 
$x\in F_y(\boldsymbol{\alpha})$ and a lower bound for 
$m_{y,{\boldsymbol{\alpha}}}$-a.e. 
$x\in[0,1]$ in the case that 
$m_{y,\boldsymbol{\alpha}}(F_y(\boldsymbol{\alpha}))=1$. To be more precise, for 
$x\in F_y(\boldsymbol{\alpha})$ with 
$\omega = ((j_\ell,k_\ell))_{\ell\geqslant 1} \in \pi^{-1}\{(x,y)\}$ that have 
$\tau_e(\omega)= \alpha_e$ for each 
$e \in \mathcal E$, we find that
\begin{align*}
\underline{d}_{m_{y,{\boldsymbol{\alpha}}}}(x) &\leq\lim_{n\to\infty}\frac{\log m_{y,{\boldsymbol{\alpha}}}(\Delta_y(k_1\dotsb k_n))}{\log \textrm{diam}(\Delta_y(k_1\dotsb k_n))}\\
&=\lim_{n\to\infty}\frac{\log\prod_{1\leq\ell\leqslant n}\frac{\alpha_{(j_\ell,k_\ell)}}{\alpha_{j_\ell}}}{\log\prod_{1\leq\ell\leqslant n}q_{(j_\ell,k_\ell)}}\\
&=\lim_{n\to\infty}\frac{\frac1n\sum_{1\leq\ell\leqslant n}\log\alpha_{(j_\ell,k_\ell)}-\frac1n\sum_{1\leq\ell\leqslant n}\log\alpha_{j_\ell}}
{\frac1n\sum_{1\leq\ell\leqslant n}\log q_{(j_\ell,k_\ell)}}.
\end{align*}By collecting like terms, we find that
\begin{align*}
\sum_{1\leq\ell\leqslant n}\log\alpha_{(j_\ell,k_\ell)} &=\sum_{e\in\mathcal E}\#\{1 \leqslant i \leqslant n \, : \, \omega_i=e\}\log\alpha_e,\\
\sum_{1\leq\ell\leqslant n}\log\alpha_{j_\ell} &=\sum_{0 \leqslant j \lt J}\#\{1 \leqslant i \leqslant n \, : \, j_i=j\}\log\alpha_j,\\
\sum_{1\leq\ell\leqslant n}\log q_{(j_\ell,k_\ell)} &=\sum_{e\in\mathcal E}\#\{1 \leqslant i \leqslant n \, : \, \omega_i=e\}\log q_e.
\end{align*} Since 
$x\in F_y({\boldsymbol{\alpha}})$, we have for each 
$e \in \mathcal E$ and 
$0 \leqslant j \lt J$ that
\begin{align*} \begin{split}
\lim_{n\to\infty}\frac{\#\{1 \leqslant i \leqslant n \, : \, \omega_i=e\}}n=\ & \tau_e(\omega)=\alpha_e,\\
\lim_{n\to\infty}\frac{\#\{1 \leqslant i \leqslant n \, : \, j_i=j\}}n=\ & \lim_{n\to\infty}\sum_{0\leqslant k \lt B_j}\frac{\#\{1\leqslant i\leqslant n:\omega_{(j_i,k)}=(j,k)\}}n=\alpha_j.
\end{split}\end{align*}Thus, recalling the definition of measure-theoretic entropy from (2.1), we find that
\begin{equation*}\underline{d}_{m_{y,{\boldsymbol{\alpha}}}}(x)\leq\frac{\sum_{e\in\mathcal E}\alpha_e\log\alpha_e-\sum_{0 \leqslant j \lt J}\alpha_j\log\alpha_j}{\sum_{e\in\mathcal E}\alpha_e\log q_e}=\frac{h_{\mu_{\boldsymbol{\alpha}}}(\sigma)-h_{\mathbb P_{\boldsymbol{\alpha}}}(\sigma_J)}{-\sum_{e \in \mathcal E} \alpha_e \log q_e}\end{equation*} for all 
$x\in F_y({\boldsymbol{\alpha}})$. Therefore, it follows from lemma 2.1(i) that
\begin{equation*}\dim_{\mathrm{H}}(F_y({\boldsymbol{\alpha}}))\leq\frac{h_{\mu_{\boldsymbol{\alpha}}}(\sigma)-h_{\mathbb P_{\boldsymbol{\alpha}}}(\sigma_J)}{-\sum_{e \in \mathcal E} \alpha_e \log q_e}.\end{equation*} To prove the second statement, fix 
$y\in W({\boldsymbol{\alpha}})$ such that 
$m_{y,\boldsymbol{\alpha}}(F_y(\boldsymbol{\alpha}))=1$, which holds for 
$\nu_{\boldsymbol{\alpha}}$-a.e. 
$y \in W(\boldsymbol{\alpha})$ by lemma 4.2. Let 
$(j_m)_{m \geqslant 1} \in \pi_J^{-1}\{y\}$ be the unique sequence with 
$\tau_j((j_m)_{m \geqslant 1}) = \alpha_j$ for each 
$0 \leqslant j \lt J$, see lemma 4.1. By the above computations for the upper bound of 
$\dim_{\mathrm{H}}(F_y({\boldsymbol{\alpha}}))$ together with [Reference Pesin50, theorem 15.3(2)], we have for 
$m_{y,{\boldsymbol{\alpha}}}$-a.e. 
$x\in [0,1]$ that
\begin{equation}
\underline{d}_{m_{y,{\boldsymbol{\alpha}}}}(x) \geqslant \inf \lim_{n\to\infty}\frac{\log m_{y,{\boldsymbol{\alpha}}}(\Delta_y(k_1\dotsb k_n))}{\log\textrm{diam}(\Delta_y(k_1\dotsb k_n))},
\end{equation} where the infimum is taken over all sequences 
$(k_m)_{m \geqslant 1}$ such that 
$(j_m,k_m)_{k \geqslant 1} \in \pi^{-1}\{(x,y)\}$. We have seen that the set of x for which 
$\#\pi^{-1}\{(x,y)\} \gt 1$ is countable so is therefore a 
$m_{y,{\boldsymbol{\alpha}}}$-null set. Consequently, the infimum on the right-hand side of (4.5) is over a single sequence for 
$m_{y,{\boldsymbol{\alpha}}}$-a.e. 
$x\in[0,1]$. Fix 
$x\in F_y(\boldsymbol{\alpha})$ such that (4.5) and 
$\#\pi^{-1}\{(x,y)\}=1$ both hold. Then
\begin{align*}
\underline{d}_{m_{y,{\boldsymbol{\alpha}}}}(x) &\geqslant \lim_{n\to\infty}\frac{\log m_{y,{\boldsymbol{\alpha}}}(\Delta_y(k_1\dotsb k_n))}{\log\textrm{diam}(\Delta_y(k_1\dotsb k_n))}\\
&=\lim_{n\to\infty}\frac{\sum_{e\in\mathcal E}\frac{\#\{1 \leqslant i \leqslant n \, : \, \omega_i=e\}}n\log\alpha_e-\sum_{0 \leqslant j \lt J}\frac{\#\{1 \leqslant i \leqslant n \, : \, j_i=j\}}n \log\alpha_j}{\sum_{e\in\mathcal E}\frac{\#\{1 \leqslant i \leqslant n \, : \, \omega_i=e\}}n\log q_e}.
\end{align*} Since 
$x \in F_y(\boldsymbol{\alpha})$ for each 
$e \in \mathcal E$ and 
$0 \leqslant j \lt J$, we find
\begin{align*}\begin{split}
\tau_e(\omega)=\ &\lim_{n\to\infty}\frac{\#\{1 \leqslant i \leqslant n \, : \, \omega_i=e\}}n=\alpha_e \, \, \text{and }\\
\tau_j((j_m)_{m \geqslant 1})=\ & \lim_{n\to\infty}\frac{\#\{1 \leqslant i \leqslant n \, : \, j_i=j\}}n=\alpha_j.
\end{split}\end{align*}Therefore,
\begin{equation*}\underline{d}_{m_{y,{\boldsymbol{\alpha}}}}(x)\geq\frac{\sum_{e\in\mathcal E}\alpha_e\log\alpha_e-\sum_{0 \leqslant j \lt J}\alpha_j\log\alpha_j}{\sum_{e\in\mathcal E}\alpha_e\log q_e}=\frac{h_{\mu_{\boldsymbol{\alpha}}}(\sigma)-h_{\mathbb P_{\boldsymbol{\alpha}}}(\sigma_J)}{-\sum_{e \in \mathcal E} \alpha_e \log q_e}.\end{equation*} Since this holds for 
$m_{y,{\boldsymbol{\alpha}}}$-a.e. 
$x \in F_y(\boldsymbol{\alpha})$ and 
$m_{y,{\boldsymbol{\alpha}}}(F_y({\boldsymbol{\alpha}}))=1$, it follows from lemma 4.2 and lemma 2.1(ii) that
\begin{align*} \dim_{\mathrm{H}}(F_y({\boldsymbol{\alpha}})) \geq\dim_{\mathrm{H}}(m_{y,{\boldsymbol{\alpha}}}) \geq\frac{h_{\mu_{\boldsymbol{\alpha}}}(\sigma)-h_{\mathbb P_{\boldsymbol{\alpha}}}(\sigma_J)}{-\sum_{e \in \mathcal E}\alpha_e \log q_e}. \end{align*}Remark 4.1. The additional condition on the frequency vector α that for each 
$0 \leqslant j \lt J$ there are 
$k,\ell \in B_j$ with 
$\alpha_{(j,k)} \gt 0$ and 
$\alpha_{(j,\ell)} \gt 0$ is used only to remove the infimum in (4.5). Another assumption that would allow us to remove the infimum is to assume that 
$\# \pi^{-1}\{(x,y)\} =1$ for all 
$(x,y) \in F(\boldsymbol{\alpha})$. This holds for example in the following cases.
(i) If
$B_j \geqslant 3$ for some $0 \leqslant j \lt J$ and there is a $1 \leqslant k\leqslant B_j-2$ with $\alpha_{(j,k)} \gt 0$, then for any (x, y) with $\#\pi^{-1}\{(x,y)\} \gt 1$ and any $\omega \in \#\pi^{-1}\{(x,y)\}$, we obtain $\tau_{(j,1)}(\omega)=0 \neq \alpha_{(j,k)}$ and thus $(x,y) \not \in F(\boldsymbol{\alpha})$.(ii) If
$\alpha_e \gt 0$ for each $e \in \mathcal E$ and $h_{(j,0)}(0)=0$ and $h_{(j,B_j-1)}(1)=1$ for all $0 \leqslant j \lt J$, then for any (x, y) with $\#\pi^{-1}\{(x,y)\} \gt 1$ any $\omega \in \#\pi^{-1}\{(x,y)\}$ will either end in an infinite string of digits from the set $\{(j,0) \, : \, 0 \leqslant j \lt J\}$ or in an infinite string of digits from the set $\{(j,B_j-1) \, : \, 0 \leqslant j \lt J\}$ and again there is at least one $e \in \mathcal E$ for which $\tau_e (\omega) = 0 \neq \alpha_e$.
5. Examples
Example 5.1. For a concrete example, let J = 2, 
$B_0=2$, 
$B_1=3$ so that
\begin{equation*}\mathcal E = \{(0,0), (0,1), (1,0), (1,1), (1,2)\}.\end{equation*}Let
\begin{align*} \begin{split}
h_{(0,k)}(x) =\ & \frac{x+k}{2}, \quad k=0,1,\\
h_{(1,k)}(x) =\ & \frac{x+k}{3}, \quad k=0,1,2,
\end{split}\end{align*} so 
$\varepsilon_e = 0$ for all 
$e \in \mathcal E$ and 
$r_{(0,1)}= \frac12$, 
$r_{(1,1)}=\frac13$ and 
$r_{(1,2)}= \frac23$. Take 
$p \in (0,1)$ arbitrary and let 
$p_0=p$, so 
$f_0(y) = py$ and 
$f_1(y) = (1-p)y+p$. This gives
\begin{align*} A_{(0,k)} = \left[\begin{array}{cc}
1/2 & 0 \\
0 & p
\end{array}\right], \, A_{(1,k)} = \left[\begin{array}{cc}
1/3 & 0 \\
0 & 1-p
\end{array}\right]\end{align*}and
\begin{align*} \boldsymbol{v}_{(0,0)} = \boldsymbol{v}_{(1,0)} = \left[\begin{array}{c}
0 \\
0\end{array} \right], \, \boldsymbol{v}_{(0,1)} = \left[\begin{array}{c}
1/2 \\
0\end{array} \right], \, \boldsymbol{v}_{(1,1)} = \left[\begin{array}{c}
1/3 \\
p\end{array} \right] , \, \boldsymbol{v}_{(1,2)} = \left[\begin{array}{c}
2/3 \\
p\end{array} \right].\end{align*} See figure 1(b) for an illustration of how this GLS IFS 
$\{A_e + \boldsymbol{v}_e\}_{e \in \mathcal E}$ acts on 
$[0,1]^2$. For the number expansions, if 
$(e_m)_{m \geqslant 1} \in \mathcal E^{\mathbb N}$, then for each 
$m \geqslant 1$ we get 
$s_m=0$, 
$K_m=2$ if 
$j_m=0$ and 
$K_m=3$ if 
$j_m=1$ and 
$t_m \in \{0,\frac12, \frac13, \frac23\}$ for all 
$m \geqslant 1$. So, in fact, for each 
$(e_m)_{m \geqslant 1} \in \mathcal E^{\mathbb N}$ if we set 
$\kappa(n) = \#\{1 \leqslant m \leqslant n \, : \, j_m = 0\}$, then (1.9) becomes
\begin{align*} \lim_{m \to \infty} h_{e_1} \circ \cdots \circ h_{e_m} (0) = \sum_{m \geqslant 1} \frac{t_m}{2^{\kappa(m)}3^{m-\kappa(m)}}.\end{align*} Hence, this GLS IFS produces for each 
$x \in [0,1]$ number expansions in mixed base 2 and 3. Note that for this IFS, 
$h_e \neq h_{e'}$ if 
$e \neq e'$. So from proposition 4.1, it follows that each 
$x \in [0,1]$ has uncountably many different expansions with mixed bases 2 and 3. If α satisfies the assumption of theorem 1.2, then we can apply theorem 1.2. We have 
$\alpha_0 = \alpha_{(0,0)}+ \alpha_{(0,1)}$ and 
$\alpha_1 = \alpha_{(1,0)}+ \alpha_{(1,1)}+ \alpha_{(1,2)}$ and obtain for 
$\nu_{\boldsymbol{\alpha}}$-a.e. 
$y \in W(\boldsymbol{\alpha})$ that
\begin{align*} \dim_{\mathrm H}(F_y(\boldsymbol{\alpha})) = \frac{\sum_{e \in \mathcal E} \alpha_e \log \alpha_e - \alpha_0\log \alpha_0 - \alpha_1 \log \alpha_1}{-\alpha_0\log 2 - \alpha_1 \log 3}. \end{align*} Note that if we consider the IFS 
$\{A_{(0,0)} + \boldsymbol{v}_{(0,0)}, A_{(1,2)} + \boldsymbol{v}_{(1,2)} \}$ with 
$A_{(j,k)}$ and 
$\boldsymbol{v}_{(j,k)}$ as in the example and we take 
$\frac12 \lt p \lt \frac23$, then we obtain a diagonally affine IFS that satisfies condition (Dʹ)(a). Hence, we can apply theorem 1.1 to this IFS to obtain for each 
$\alpha \in (0,1)$ an expression for the Hausdorff dimension of the set of points 
$(x,y) \in [0,1]^2$ that have a GLS expansion containing only the digits 
$(0,0)$ and 
$(1,2)$ and in which 
$(0,0)$ occurs with frequency α (and thus 
$(1,2)$ with frequency 
$1-\alpha$). Of course, here we can take any other combination of a digit from 
$\{(0,0), (0,1)\}$ and a digit from 
$\{(1,0), (1,1), (1,2)\}$ to obtain a similar result.
We can extend this example in the following sense. Fix some 
$J \in \mathbb N_{\geqslant 2}$ and different integers 
$M_0,M_1, \ldots, M_{J-1} \gt J$. Also fix some probability vector 
$(p_j)_{0 \leqslant j \lt J}$. So,
\begin{align*} \mathcal E = \{(j,k) \, : \, 0 \leqslant j \lt J , \, \, 0 \leqslant k \lt M_j\}. \end{align*} For 
$(j,k) \in \mathcal E$, set
\begin{align*} A_{(j,k)} = \left[\begin{array}{cc}
1/M_j & 0 \\
0 & p_j
\end{array}\right], \quad \boldsymbol{v}_{(j,k)} = \left[\begin{array}{c}
k/M_j \\
\sum_{i=0}^{j-1}p_i\end{array} \right].\end{align*} For each 
$x \in [0,1]$ and sequence 
$(j_m)_{m \geqslant 1} \in \{0,1, \ldots, J-1\}^{\mathbb N}$, the GLS expansion produced by this system has the form
\begin{equation}
x = \sum_{m \geqslant 1} \frac{d_m}{M_0^{c_{0,m}} M_1^{c_{1,m}} \cdots M_{J -1}^{c_{J-1,m}}},
\end{equation} with 
$d_m \in \{0, \ldots, M_{j_m}-1\}$ and 
$c_{j,m} = \# \{1 \leqslant i \leqslant m \, : \, j_i=j\}$. In other words, the system produces for each 
$x\in [0,1]$ uncountably many different mixed base expansions with bases 
$M_0, \ldots, M_{J-1}$. Here, we need to remark that we consider two GLS expansions produced by the system different if the two corresponding sequences in
\begin{align*} \mathcal A = \bigcup_{(j,k) \in \mathcal E}\left\{\left(0,M_j,\frac{k}{M_j}\right) \right\} \end{align*}are different. For the point 0, for example, this means that the GLS expansions generated by the system are all of the expansions of the form
\begin{align*} 0=\sum_{m \geqslant 1} \frac{0}{M_0^{c_{0,m}} \cdots M_{J -1}^{c_{J-1,m}}},\end{align*} with 
$(c_{0,m}, \ldots, c_{J-1,m}) \in \mathbb N^{J}$ satisfying 
$\sum_{\ell=0}^{J-1} c_{\ell,m} =m$.
Example 5.2. Fix an 
$N \in \mathbb N_{\geqslant 3}$ and a 
$0 \lt p \lt 1$ and let 
$\mathcal E = \{(j,k)\, : \, j=0,1, \, 0 \leqslant k \lt N\}$. For 
$0 \leqslant k \lt N$, set
\begin{align*} A_{(0,k)} = \left[\begin{array}{cc}
1/N & 0 \\
0 & p
\end{array}\right], \, A_{(1,k)} = \left[\begin{array}{cc}
-1/N & 0 \\
0 & 1-p
\end{array}\right]\end{align*}and
\begin{align*} \boldsymbol{v}_{(0,k)} = \left[\begin{array}{c}
k/N \\
0\end{array} \right], \, \boldsymbol{v}_{(1,k)} = \left[\begin{array}{c}
(k+1)/N \\
p\end{array} \right].\end{align*} Then for any 
$x \in [0,1]$ and any 
$(j_m)_{m \geqslant 1} \in \{0,1\}^{\mathbb N}$, the number expansion of x produced by this system has the form
\begin{align*} x = \sum_{m \geqslant 1} (-1)^{j_m} \frac{d_m}{N^m},\end{align*} for some 
$d_m \in \{0, \ldots, N-1\}$, 
$m \geqslant 1$. So, the system produces for each x a signed base N-expansion in which the signs of the terms correspond to a preset sequence of signs 
$(j_m)_{m \geqslant 1}$.
Also this system satisfies 
$h_e \neq h_{e'}$ whenever 
$e \neq e'$ and together with any frequency vector 
$\boldsymbol{\alpha} \in (0,1)^{2N}$ for which the conditions of theorem 1.2 are satisfied, the Hausdorff dimension of the Besicovitch–Eggleston set 
$F_y(\boldsymbol{\alpha})$ for 
$\nu_{\boldsymbol{\alpha}}$-a.e. 
$y \in W(\boldsymbol{\alpha})$ is given by theorem 1.2. For 
$\frac1N \lt p \lt \frac{N-1}{N}$ and any 
$0 \leqslant k,\ell \lt N$, the system 
$\{A_{(0,k)} + \boldsymbol{v}_{(0,k)}, A_{(1,\ell)} + \boldsymbol{v}_{(1,\ell)} \}$ satisfies (Dʹ)(a), so then also theorem 1.1 applies.
Acknowledgements
The first two authors acknowledge the hospitality of Uppsala University and the support of the Knut and Alice Wallenberg Foundation for their research visits. The third author acknowledges the hospitality of Leiden University and the support of the Knut and Alice Wallenberg Foundation and the Swedish Research Council under grant no. 2021-06594 while the author was in residence at Institut Mittag-Leffler in Djursholm, Sweden during the program ‘Two Dimensional Maps’. We also thank the anonymous referee for their valuable corrections and suggestions.
Open access
