Pascal Chossat

Grégory Faye

Olivier Faugeras

Bifurcation of hyperbolic planforms

Keywords: Equivariant bifurcation analysis, neural fields, Poincaré disc, periodic lattices, hyperbolic planforms, irreductible representations, Laplace-Beltrami operator

Motivated by a model for the perception of textures by the visual cortex in primates, we analyse the bifurcation of periodic patterns for nonlinear equations describing the state of a system defined on the space of structure tensors, when these equations are further invariant with respect to the isometries of this space. We show that the problem reduces to a bifurcation problem in the hyperbolic plane D (Poincaré disc). We make use of the concept of periodic lattice in D to further reduce the problem to one on a compact Riemann surface D/Γ, where Γ is a cocompact, torsion-free Fuchsian group. The knowledge of the symmetry group of this surface allows to carry out the machinery of equivariant bifurcation theory. Solutions which generically bifurcate are called "Hplanforms", by analogy with the "planforms" introduced for pattern formation in Euclidean space. This concept is applied to the case of an octagonal periodic pattern, where we are able to classify all possible H-planforms satisfying the hypotheses of the Equivariant Branching Lemma. These patterns are however not straightforward to compute, even numerically, and in the last section we describe a method for computation illustrated with a selection of images of octagonal H-planforms.

Introduction

In a recent paper [START_REF] Chossat | Hyperbolic planforms in relation to visual edges and textures perception[END_REF] a model for the visual perception of textures by the cortex was proposed, which assumes that populations of neurons in each so-called "hypercolumn" of the visual cortex layer V1 are sensitive to informations carried by the structure tensor of the image (in fact, of a small part of it since each hypercolumn is dedicated to a specific part of the visual field). The structure tensor is a 2 × 2 positive definite, symmetric matrix, the eigenvalues of which characterize to some extent the image geometric properties like the presence and position of an edge, the contrast, etc. This led to write the average membrane potential V of the populations of neurons in a hypercolumn as a function of the structure tensor (and time). The equations that govern the evolution of the average membrane potential are ∂V ∂τ (T ,τ)=-αV (T ,τ)+ H w(T , T )S(V (T ,τ)) dT + I(T ,τ)

where T , T ∈ SDP(2) (the space of structure tensors), S is a smooth function R → R of sigmoid type (S(x) →± 1 as x →± ∞ ), I corresponds to some input signal coming from different brain areas such as the thalamus and w(T , T ) is a function which expresses the interaction between the populations of neurons of types T and T in the hypercolumn. We can introduce a distance on the set of structure tensors by noting that it has the structure of a homogeneous space for the group of invertible matrices GL(2, R) acting by coordinate change on quadratic forms associated with the structure tensors [START_REF] Weil | Groupes des formes quadratiques indéfinies et des formes bilinéaires alternées[END_REF]: for any T∈SDP [START_REF] Aurich | Periodic-orbit sum rules for the hadamard-gutzwiller model[END_REF] and G ∈ GL(2, R), we define G •T = t G T G. This induces a Riemannian metric on the tangent space: g T (A, B)=tr(T -1 AT -1 B) and the corresponding distance in SDP( 2) is given by (see [START_REF] Moakher | A differential geometric approach to the geometric mean of symmetric positive-definite matrices[END_REF])

d 0 (T 1 , T 2 )= log 2 λ 1 + log 2 λ 2 (2)
where λ 1 , λ 2 are the eigenvalues of T -1 1 T 2 . Because of its invariance this distance is biologically plausible since the neurons have no obvious way of knowing in which coordinate system are expressed the components of the structure tensors and this leads to the assumption that the function w is insensitive to coordinate changes in SDP [START_REF] Aurich | Periodic-orbit sum rules for the hadamard-gutzwiller model[END_REF]. This implies that w is in fact a function of the distance d 0 : w(T , T )=f (d 0 (T , T ))

(3)

As a consequence, the equation ( 1) is invariant under the action of the isometry group GL(2, R).

Assume now that there is no input: I =0 , and that V =0is a solution of (1) (we can assume w.l.o.g. that S(0) = 0). If the slope µ = S (0) is small enough, small perturbations of this basic state are damped to 0 and the basic state is stable. We can however expect that if µ exceeds a critical value an instability will grow and lead to a new state which need not be invariant under the action of GL(2, R). In other word we expect pattern formation through spontaneous symmetry breaking, a phenomenon which has been widely studied in other contexts in the last 40 years (see [START_REF] Hoyle | Pattern formation: an introduction to methods[END_REF] for a review). From the mathematical point of view the theoretical framework is the well-established bifurcation theory with symmetry [START_REF] Golubitsky | Singularities and Groups in Bifurcation Theory[END_REF], [START_REF] Chossat | Methods in Equivariant Bifurcations and Dynamical Systems[END_REF]. There is however a basic difference between the problem at hand and previous studies in pattern formation and symmetry-breaking bifurcation: the symmetries (isometries) are not Euclidean. Moreover the group GL(2, R) being non compact, the spectrum of the equivariant linear operators is continuous and eigenvalues have infinite multiplicity in general. This is also true when the symmetry group is the group of displacements and reflections in Euclidean space. In this case a rational approach is to resctrict the analysis to the bifurcation of patterns which are periodic (in space) with a period equal to 2π/k where k is the wave number of the most unstable Fourier modes. This hypothesis is supported by observations (to some extent) and it provides a framework which is suitable for the bifurcation analysis. Namely, the spectrum of operators restricted to functions which are invariant on a periodic lattice is discrete and consists of eigenvalues with finite multiplicity. This allows to fully exploit the power of equivariant bifurcation theory and to describe the different types of bifurcated solutions by their (remaining) symmetries.

We would like to proceed in a similar way for equation [START_REF] Allaire | Analyse numérique et optimisation[END_REF], that is to consider solutions which are periodic in the space of structure tensors (one has to give a precise meaning to this statement) and then apply equivariant bifurcation methods to describe the set of solutions. We are not addressing here the question of the relevance of looking for such solutions, neither from the neurophysiological point of view, nor even from the point of view of their stability under perturbations which do not have the same periodicity (or no periodicity at all). Let us mention the fact that in [START_REF] Chossat | Hyperbolic planforms in relation to visual edges and textures perception[END_REF] we identified families of subgroups of the group of isometries of the set of structure tensors that naturally arose from the analysis of the retinal input to the hypercolumns in visual area V1. These subgroups, which we called neuronal Fuchsian groups, have the property that the membrane potential functions are invariant under their action. This is one example of the neurophysiological relevance of looking for solutions of (1) that are periodic in this sense but it certainly does not give the final answer: These difficult questions will be addressed subsequently. Let us point out that, on the other hand, a classification of possible periodic patterns should be largely independent of the model equations as long as these equations share some basic properties (like GL(2, R) invariance). In this paper we shall therefore not focus on equation [START_REF] Allaire | Analyse numérique et optimisation[END_REF] and our results would apply equally to other models such as reaction-diffusion or Swift-Hohenberg equations with Laplace-Beltrami operator in SDP(2)).

The structure of the paper is as follows:

1. In section 2 we introduce the necessary basic material from hyperbolic geometry. Periodic lattices and functions are defined in SDP(2) and it will be shown that the problem can be decoupled and reduced to looking for periodic patterns in the hyperbolic plane (we shall work with its Poincaré disc representation D).

2. In section 3 we set the bifurcation problem for periodic patterns in D. Applying the equivariant branching lemma [START_REF] Golubitsky | Singularities and Groups in Bifurcation Theory[END_REF] (or isolated stratum principle in the variational case), this leads us to defining hyperbolic planforms, which are the hyperbolic counterparts of the planforms defined in the Euclidean case (see [START_REF] Dionne | Planforms in two and three dimensions[END_REF]).

3. In section 4 the methods of section 3 are applied to the case of a regular octagonal pattern for which we can describe all planforms which result from the applicaiton of the equivariant branching lemma.

4.

In section 5 we compute eigenvalues and eigenfunctions of the Laplace-Beltrami operator in the hyperbolic octagon in order to exhibit hyperbolic planforms which satisfy certain isotropy conditions.

2 Periodic lattices and functions in the space of structure tensors 2.1 The space SDP(2)

The following decomposition will prove to be very convenient. Let ξ 2 be the determinant of T∈SDP [START_REF] Aurich | Periodic-orbit sum rules for the hadamard-gutzwiller model[END_REF].W e set ξ>0 by convention. Writing T = ξT , we see that

SDP(2) = R + * × SSDP(2)
where SSDP(2) denotes the subspace of tensors with determinant 1. In the open cone SDP(2) the surface SSDP( 2) is a hyperboloid sheet and it can be shown that it carries a metric induced by the metric g which is just the usual metric of the hyperbolic plane. The isometry group of displacements in SSDP(2) is the special linear group SL(2, R), and indeed we may also write the group of orientation preserving isometries

GL + (2, R) = SL(2, R) × R + * acting on a tensor T =(T ,ξ) by (Γ,α) • (T ,ξ)=( t ΓT Γ,α 2 ξ). (4)
We now identify SDP(2) with the "half" open cylinder D×R + * , where D is the Poincaré disc, through the following change of variables.

Let T = x 1 x 3 x 3 x 2 with x 1 > 0 and x 1 x 2 -x 2 3 > 0.
We set

x i = z 3 xi ,i=1, 2, 3,z 3 > 0, x1 x2 -x2 3 =1.
Now the hyperbolic plane SSDP( 2) is further identified with D through the change of coordinates

       x1 = (1+z1) 2 +z 2 2 1-z 2 1 -z 2 2 x2 = (1-z1) 2 +z 2 2 1-z 2 1 -z 2 2 x3 = 2z2 1-z 2 1 -z 2 2 (5)
These formulas define a diffeomorphism

Θ: (x 1 ,x 2 ,x 3 ) ∈ SDP(2) → (z 1 ,z 2 ,z 3 ) ∈D×R + * (6)
The group of direct isometries (displacements) in D×R

g 1 (z 1 ,z 2 ,z 3 )= 4(dz 2 1 + dz 2 2 ) (1 -z 2 1 -z 2 2 ) 2 + dz 2 3 z 2 3 ( 7
)
Proof. The fact that the decomposition is diagonal follows directly from the direct product decomposition and group action [START_REF] Balazs | Chaos on the pseudosphere[END_REF]. In order to compute the precise expression for the metric, let us first express the form g in SDP(2) in suitable coordinates of the tangent space (which is the space of symmetric 2 × 2 matrices). Such a basis is given, at any tensor T , by

∂ ∂x 1 = 10 00 ∂ ∂x 2 = 00 01 ∂ ∂x 3 = 01 10
Then a straightforward calculation shows that

g(x 1 ,x 2 ,x 3 )=z -4 3 [x 2 2 dx 2 1 +2x 2 3 dx 1 dx 2 -4x 2 x 3 dx 1 dx 3 +x 2 1 dx 2 2 -4x 1 x 3 dx 2 dx 3 + 2(x 1 x 2 + x 2 3 )dx 2 3 ].
The proposition follows by applying the pull-back Θ * to this form.

Corollary 1. The volume element in the z j coordinates is

dV = 4 dz 1 dz 2 dz 3 (1 -z 2 1 -z 2 2 ) 2 z 3 (8)
Corollary 2. The Laplace -Beltrami operator in D×R + * in z j coordinates is

= (1 -z 2 1 -z 2 2 ) 2 4 ∂ 2 ∂z 2 1 + ∂ 2 ∂z 2 2 + z 3 ∂ ∂z 3 + z 2 3 ∂ 2 ∂z 2 3 (9)
We note D the first term on the r.h.s.

Let us now compute the distance d in D×R + * . We set z =(z 1 ,z 2 ,z 3 ) and z =(z 1 ,z 2 ,z 3 ). Then Proposition 2.

d(z, z )= (d 2 ((z 1 ,z 2 ), (z 1 ,z 2 )) 2 + log 2 ( z 3 z 3 ) ( 10
)
where d 2 denotes the distance in the Poincaré disc.

Proof. We start from expression (2) which gives the distance between two tensors T 1 and T 2 . Let T = z 3 T , T = z

3

T and α = z 3 /z 3 . We note λ, 1/λ the eigenvalues of T -1 T . Then

d 0 (T , T ) 2 = log 2 αλ + log 2 α λ = 2 log 2 α + 2 log 2 λ = 2 log 2 α + d 0 ( T , T ) 2 . But d 0 ( T , T )=d 2 ((z 1 ,z 2 ), (z 1 ,z 2
)), hence the result.

Note that the distance in D is given by the following formula where z, z are complex numbers in the unit disc.

d 2 (z, z ) = 2 arctanh |z -z | |1 -zz | , (11)

Periodic lattices in D×R + *

Let us recall first some basic properties of the group of isometries in the hyperbolic plane (which we shall always identify with the Poincaré disc D in the sequel). We refer to textbooks in hyperbolic geometry for details and proofs. The direct isometries (preserving the orientation) in D are the elements of the special unitary group, noted SU(1, 1), of 2 × 2 Hermitian matrices with determinant equal to 1. Given

γ = αβ β α such that |α| 2 -|β| 2 =1,
the corresponding isometry γ is defined by

γ • z = αz + β βz + α ,z ∈D (12)
Orientation reversing isometries of D are obtained by composing any transformation [START_REF] Chossat | Methods in Equivariant Bifurcations and Dynamical Systems[END_REF] with the reflection κ : z → z. The full symmetry group of the Poincaré disc is therefore

U(1, 1) = SU(1, 1) ∪ κ • SU(1, 1).
These isometries preserve angles, however they do not transform straight lines into straight lines. Given two points z = z in D, there is a unique geodesic passing through them: the portion in D of the circle containing z and z and intersecting the unit circle at right angles. This circle degenerates to a straight line when the two points lie on the same diameter. Any geodesic uniquely defines the reflection through it. Reflections are orientation reversing, one representative being the complex conjugation: κ • z = z.

We distinguish three different kinds of direct isometries in D, according to which conjugacy class of the following one parameter subgroups it does belong:

Definition 1.                        K = {r ϕ = e iϕ/2 0 0 e -iϕ/2 ,ϕ ∈ S 1 } A = {a t = cosh t/2s i n h t/2 sinh t/2 cosh t/2 ,t ∈ R} N = {n s = 1+is -is is 1 -is ,s ∈ R}
Note that r ϕ • z = e iϕ z for z ∈Dand also, a t • 0 = tanh(t). The elements of A are sometimes called "boosts" [START_REF] Balazs | Chaos on the pseudosphere[END_REF]. The following theorem is called the Iwasawa decomposition (see [START_REF] Iwaniec | Spectral methods of automorphic forms[END_REF]).

Theorem 1. SU(1, 1) = K • A • N
The orbits of A converge to the same limit points b ±1 = ±1 on the unit circle when t →± ∞ . In particular the diameter (b -1 ,b 1 ) is the orbit {x = tanh(t),t∈ R}. The orbits of N are the circles tangent to the unit circle at b 1 . These circles are called horocycles with limit point b 1 (N is called the horocyclic group). The orbits of K are circles inside D (they coincide with Euclidean circles only when they are centered at the origin). Isometries of the types K, A, N are respectively called elliptic, hyperbolic and parabolic. Elliptic isometries have one fixed point in D while hyperbolic isometries have two fixed points on the boundary ∂D and parabolic isometries have one fixed point on ∂D (hence "at infinity"). Hyperbolic isometries play the same role as translations in Euclidean space while parabolic isometries have no counterpart. There is however one important difference with Euclidean translations: two hyperbolic translations do not commute in general. The fundamental difference from this point of view between D and the Euclidean plane R 2 is that the latter is itself an Abelian group while DSU(1, 1)/SO(2, R) is not a group. This makes its analysis, especially its Fourier analysis, harder and less intuitive.

We can now define a (periodic) lattice in D and in D×R * . Let Γ be a discrete subgroup of SU(1, 1) such that the orbits of points in D under the action of Γ have no accumulation point in D. This is a Fuchsian group. To any Fuchsian group we can associate a fundamental domain which is the closure, noted F Γ , of an open set o F Γ ⊂D with the following properties [START_REF] Katok | Fuchsian Groups[END_REF]:

(i) if γ = Id ∈ Γ, then γF Γ ∩ o F Γ = ∅; (ii) γ∈Γ γF Γ = D.
Hence F Γ generates a periodic tiling (or tessellation) of D. A fundamental domain need not be a compact subset of D (it may have vertices on the circle at infinity ∂D). When it does, Γ is called a cocompact Fuchsian group. In this case Γ contains no parabolic element, its area is finite and a fundamental domain can always be built as a polygon (the Dirichlet region, see [START_REF] Katok | Fuchsian Groups[END_REF]). As mentioned in the introduction we identified in [START_REF] Chossat | Hyperbolic planforms in relation to visual edges and textures perception[END_REF] a family of Fuchsian groups, some of them compact, that naturally arose from the analysis of the symmetries in the spatial distribution of the photoreceptors in the retina.

The following definition is just a translation to the hyperbolic plane of the definition of an Euclidean lattice.

Definition 2. A lattice group of D is a cocompact Fuchsian group which contains no elliptic element.

The action of a lattice group has no fixed point, therefore the quotient surface D/Γ is a (compact) manifold and it is in fact a Riemann surface. A remarkable theorem states that any compact Riemann surface is isomorphic to a lattice fundamental domain of D if and only if it has genus g ≥ 2 [START_REF] Katok | Fuchsian Groups[END_REF]. The case g =1corresponds to lattices in the Euclidean plane (in this case there are three kinds of fundamental domains: rectangles, squares and hexagons). The simplest lattice in D, with genus 2, is generated by an octagon and will be studied in detail in Section 4.

Given a lattice, we may ask what is the symmetry group of the fundamental domain F Γ , identified with the quotient surface D/Γ. Indeed, this information will play a fundamental role in the subsequent bifurcation analysis. In the case of Euclidean lattice, the quotient R 2 /Γ is a torus T (genus one surface), and the group of automorphisms is H T where H is the holohedry of the lattice: H = D 2 ,D 4 or D 6 for the rectangle, square and hexagonal lattices respectively. In the hyperbolic case the group of automorphisms of the surface is finite. In order to build this group we need first to introduce some additional definitions.

Tilings of the hyperbolic plane can be generated by reflections through the edges of a triangle τ with vertices P , Q, R and angles π/, π/m and π/n respectively, where , m, n are integers such that 1/ +1/m +1/n < 1 [START_REF] Katok | Fuchsian Groups[END_REF].

Remember that reflections are orientation-reversing isometries. We note κ, κ and κ the reflections thro ugh the edges PQ, QR and RP respectively (Figure 1). The group generated by these reflections contains an index 2 Fuchsian subgroup Λ called a triangle group, which always contains elliptic elements because the product of the reflections through two adjacent edges of a polygon is elliptic with fixed point at the corresponding vertex. One easily shows that Λ is generated by the rotations of angles 2π/l, 2π/m and 2π/n around the vertices P , Q, R re spectively. A fundamental domain of Λ is the "quadrangle" F Λ = τ ∪ κτ [START_REF] Katok | Fuchsian Groups[END_REF]. Note that F Λ D /Λ is a sphere (genus 0 surface) obtained by identifying the three edges of τ . The subgroup of hyperbolic translations in Λ is a lattice group Γ, normal in Λ, whose fundamental domain is filled with copies of the basic tile τ . The group of orientation-preserving automorphisms of F Γ D/Γ is therefore G =Λ/Γ. From the algebraic point of view, G is generated by three elements a, b, c satisfy therelations a = b m = c n =1and a • b • c =1 . We say that G is an (l, m, n) group. Taking account of orientation-reversing isometries, the full symmetry group of

F Γ is G * = G ∪ κG = G Z 2 (κ)
. This is also a tiling group of F Γ with tile τ : the orbit G * τ fills F Γ and its elements can only intersect at their edges.

Given a lattice, how to determine the groups G and G * ? The following theorem gives conditions for this, see [START_REF] Broughton | Triangular surface tiling groups for low genus[END_REF]. Theorem 2. An (l, m, n) group G is the tiling rotation group of a compact Riemann surface of genus g if and only if its order satisfies the Riemann-Hurwitz relation

|G| = 2g -2 1 -( 1 + 1 m + 1 n )
.

Tables of triangle groups for surfaces of genus up to 13 can be found in [START_REF] Broughton | Triangular surface tiling groups for low genus[END_REF]. These definitions extend naturally to SDP(2) D×R * as follows.

Definition 3. A lattice group of D×R * is a subgroup of the form Γ × Ξ where Γ is a lattice group acting in D and Ξ is a non trivial discrete subgroup of R * .

Any discrete subgroup of R * is generated by a positive number a and can be further identified with

Z: Ξ= {a n ,n∈ Z}. A fundamental domain for Γ × Ξ is a "box" F Γ × [1,a].

Plane waves in SDP(2)

Let us first recall the Euclidean setting for functions defined in R 3 (we could take any other R n , n>0). In this case every function of the form e λk•r where k ∈ R 3 is a unit vector, is an eigenfunction of the Laplace operator in R 3 :

∆e λk•r = -λ 2 e λk•r , r ∈ R 3 .
The fact that the eigenvalues do not depend upon the direction of the wave vector k reflects the rotational invariance of the Laplace operator. Moreover if we take λ = iα, α ∈ R, then e λk•r is invariant under translations in R 2 by any vector e satisfying the condition k • e =2 nπ where n ∈ Z (it clearly does not depend upon the coordinate along the axis orthogonal to k). The functions e iαk•r are elementary spatial waves in R 3 . Now, given α>0 and a basis of unit vectors {k 1 , k 2 , k 3 } of R 3 we can define the translation group L spanned by e i , i =1 , 2, 3, such that k i • e j =2 π/αδ ij . Hence L is a lattice group of R 3 . It defines a periodic tiling, the fundamental domain of which is a compact cell which we may identify with the quotient space R 3 /L and which we can identify with a 3-torus. Any smooth enough function in R 3 which is invariant under the action of L can be expanded in a Fourier series of elementary spatial waves e iα(mk1+nk2+pk3)•r , m, n, p ∈ Z. The Laplace operator in the space of square-integrable functions in R 3 /L is self-adjoint and its spectrum consists of real isolated eigenvalues with finite multiplicities. The multiplicity depends upon the holohedry of the lattice, which we defined in the previous section (the largest subgroup of O(3) leaving invariant the lattice). There are finitely many holohedries (see [START_REF] Miller | Symmetry groups and their applications[END_REF] for details). It follows from the above considerations that by restricting the analysis to classes of functions which are invariant under the action of a lattice group, one can apply standard techniques of equivariant bifurcation theory to assert the generic existence of branches of solutions of Euclidean invariant bifurcation problems, which are spatially periodic with respect to lattice groups and whose properties are largely determined by the holohedry of the lattice [START_REF] Golubitsky | Singularities and Groups in Bifurcation Theory[END_REF], [START_REF] Dionne | Planforms in two and three dimensions[END_REF]. Note also this was the approach of [START_REF] Bressloff | What Geometric Visual Hallucinations Tell Us about the Visual Cortex[END_REF] for the analysis of the occurence of visual hallucinations in the cortex. We wish to apply the same idea to bifurcation problems defined in SDP [START_REF] Aurich | Periodic-orbit sum rules for the hadamard-gutzwiller model[END_REF]. For this we need to define elementary eigenfunctions of the Laplace-Beltrami operator such that spatially periodic functions (in a sense to be defined later) can be expanded in series of these elementary "waves" in SDP(2). In the sequel will denote the Laplace-Beltrami operator in SDP(2) or, equivalently, in D×R + * . Let b be a point on the circle ∂D, which we may take equal to b 1 =1after a suitable rotation. For z ∈D , we define the "inner product" z, b as the algebraic distance to the origin of the (unique) horocycle based at b and passing through z. This distance is defined as the hyperbolic signed length of the segment Oξ where ξ is the intersection point of the horocycle and the line (geodesic) Ob. Note that z, b does not depend on the position of z on the horocycle. In other words, z, b is invariant under the action of the one-parameter group N (see definition above). The "hyperbolic plane waves"

e ρ,b (z)=e (iρ+ 1 2 )z,b ,ρ∈ C, satisfy -D e ρ,b =(ρ 2 + 1 4 )e ρ,b .
where D is defined in Corollary 2. These are the elementary eigenfunctions with which Helgason built a Fourier transform theory for the Poincaré disc, see [START_REF] Helgason | Groups and geometric analysis[END_REF]. It follows from Helgason's theory that any eigenfunction of D can be expressed as an integral over the boundary elements:

Theorem 3. [24]
Any eigenfunction of the operator -D admits a decomposition of the form

∂D e ρ,b (z)dT ρ (b)
where T ρ is a distribution defined on the circle ∂D and the eigenvalue is

ρ 2 + 1 4 . Real eigenvalues -(ρ 2 + 1 4
) of D correspond to taking ρ real or ρ ∈ iR. The latter case is irrelevant for our study as it corresponds to exponentially diverging eigenfunctions. Therefore the real spectrum of is continuous and is bounded from above by -1/4. By using ( 9) we extend Theorem 3 to D×R + * :

Corollary 3. (i) Let us note z =(z 1 ,z 2 ,z 3 ) ∈D×R + * and z = z 1 + iz 2 . The function ψ ρ,b,β (z)=e ρ,b (z) e i log β log z3 ( 13
)
satisfies the relation ψ ρ,b,β = -(ρ 2 + 1 4 + log 2 β)ψ ρ,b,β . (

ii) Any eigenfunction of admits a decomposition of the form

e i log β log z3 ∂D e ρ,b (z) dT ρ (b). (14)
From there we can extend the definition (and properties) of the Fourier transform in the hyperbolic plane given by Helgason [START_REF] Helgason | Groups and geometric analysis[END_REF] to the space D×R + * :

Definition 4. Given a function f on D×R + * , its Fourier transform is defined by

f (ρ, b, β)= D×R + * f (z)e -ρ,b (z)e -i log β log z3 dz (15)
In the following we will look for solutions of bifurcation problems in D×R + * , which are invariant under the action of a lattice group:

(γ, ξ) • u(z, z 3 )=u(γ -1 z, ξ -1 z 3 )=u(z, z 3 ) for γ ∈ Γ, ξ ∈ Ξ.
This boils down to looking for the problem restricted to a fundamental domain with suitable boundary conditions imposed by the Γperiodicity, or, equivalently, to looking for the solutions of the problem projected onto the orbit space D/Γ×R + * /Ξ (which inherits a Riemannian structure from D×R + * ). Because the fundamental domain is compact, it follows from general spectral theory thatis self-adjoint, non negative and has compact resolvent in L 2 (D/Γ × R + * /Ξ) [START_REF] Buser | Geometry and spectra of compact Riemann surfaces[END_REF]. Hence its spectrum consists of real positive and isolated eigenvalues of finite multiplicity.

Coming back to Theorem 3, we observe that those eigenvalues λ of -D which correspond to Γ-invariant eigenfunctions, must have ρ ∈ R or ρ ∈ iR. The case ρ real corresponds to the Euclidean situation of planar waves with a given wave number, the role of which is played by ρ in D. In this case the eigenvalues of -D satisfy 1/4 <λ. On the other hand there is no Euclidean equivalent of the case ρ ∈ iR, for which the eigenvalues 0 <λ≤ 1/4 are in finite number. It turns out that such "exceptional" eigenvalues do not occur for "simple" groups such as the octagonal group to be considered in more details in the Section 4. This follows from formulas which give lower bounds for these eigenvalues. Let us give two examples of such estimates (derived by Buser [START_REF] Buser | Geometry and spectra of compact Riemann surfaces[END_REF], see also [START_REF] Iwaniec | Spectral methods of automorphic forms[END_REF]): (i) if g is the genus of the surface D/Γ, there are at most 3g -2 exceptional eigenvalues; (ii) if d is the diameter of the fundamental domain, then the smallest (non zero) eigenvalue is bounded from below by 4π sinh d 2 -2 .

Suppose now that the eigenfunction in Theorem 3 is Γ-periodic. Then the distribution T ρ satisfies the following equivariance relation [START_REF] Pollicott | Distributions at infinity for riemann surfaces[END_REF]. Let γ(θ) denote the image of θ ∈ ∂D under the action of γ ∈ Γ. Then

T ρ (γ • θ)=|γ (θ)| 1 2 +iρ T ρ (θ).
As observed by [START_REF] Series | Some geometrical models of chaotic dynamics[END_REF], this condition is not compatible with T ρ being a "nice" function. In fact, not only does there not exist any explicit formula for these eigenfunctions, but their approximate computation is itself an uneasy task. We shall come back to this point in subsequent sections.

3 Bifurcation of patterns in SDP [START_REF] Aurich | Periodic-orbit sum rules for the hadamard-gutzwiller model[END_REF] We now consider again equation ( 1), which we set in D×R + * by the change of coordinates [START_REF] Bressloff | Geometric visual hallucinations, euclidean symmetry and the functional architecture of striate cortex[END_REF]. Assuming I =0 (no external input) and the invariance hypothesis (3) for the connectivity function w, and after a choice of time scale such that α =1, the equation reads

∂V ∂τ = -V + µw * V + R(V ) (16)
where:

• µ = S (0),
• w * V denotes the convolution product

D×R + * w(z, z )V (z ) dz (with w(z, z )=f (d(z, z ))),
• R(V ) stands for the remainder terms in the integral part of (3). This implies R (0) = 0.

It is further assumed that f is a "Mexican hat" function, typically of the form

f (x)= 1 2πσ 2 1 e -x 2 2σ 2 1 -θ 1 2πσ 2 2 e -x 2 2σ 2 2 ,
where σ 1 <σ 2 and θ ≤ 1.

Let us look at the linear stability of the trivial solution of ( 16) against perturbations in the form of hyperbolic waves [START_REF]Handbook of Numerical Analysis[END_REF] with ρ ∈ R. This comes back to looking for σ's such that σ = -1+µ ŵ where ŵ is the hyperbolic Fourier transform of w as defined in Definition 15. The numerical calculation shows that for each value of ρ and β, there exists a value µ(ρ, β) such that if µ<µ (ρ, β) then all σ's are negative, while σ =0at µ = µ(ρ, β). The "neutral stability surface" defined by µ(ρ, β) is typically convex and reaches a minimum µ c at some values ρ c ,β c . Therefore when µ<µ c the trivial state V =0is stable against such perturbations while it becomes marginally stable when µ = µ c with critical modes ψ ρc,b,βc , for any b ∈ S 1 (rotational invariance). Therefore a bifurcation takes place at this critical value.

The situation is absolutely similar if instead of equation 16 we consider systems of PDEs in D×R + * with pattern selection behavior and with U (1, 1) × R + * invariance. A paradigm for such systems is the "Laplace-Beltrami" version of Swift-Hohenberg equation

∂u ∂t = µu -( + α) 2 u + u 2 ,α ∈ R + * .
with as in [START_REF] Broughton | Triangular surface tiling groups for low genus[END_REF]. It is unconceivable to solve the bifurcation problem at this level of generality, because the fact that the spectrum is continuous plus that each eigenvalue σ has an infinite multiplicity (indifference to b) makes impossible the use of the classical tools of bifurcation theory. As in the Euclidean case of pattern formation, we therefore want to look for solutions in the restricted class of patterns which are spatially periodic. In the present framework, this means looking for bifurcating patterns which are invariant under the action of a lattice group

Γ × Z + * in SU(1, 1) × R + * .
There is however an immediate big difference with the Euclidean case. While in the latter any critical wave number α c can be associated with a periodic lattice (of period 2π/α c ), in the hyperbolic case not every value of ρ c can be associated with a lattice in D. More precisely it is not generic for the Laplace-Beltrami operator in D to possess eigenfunctions, as defined in Theorem 3, which are invariant under a lattice group. We can therefore look for the bifurcation of spatially periodic solutions associated with a given lattice, but these patterns will not in general correspond to the most unstable perturbations unless the parameters in the equations are tuned so that it happens this way. The question of the observability of such patterns is therefore completely open.

We henceforth look for patterns in D×R + * which are invariant under a lattice Γ in D and which are periodic, with period 2π/β c , in R + * . This boils down to looking for solutions in the space

L 2 (D/Γ × R + * /β c Z + * ). Note that R + * /β c Z + * S 1 .
With a suitable inner product this space admits an orthonormal Hilbert basis which is made of functions of the form Ψ(z)e ni log(βc)log(z3) ,z ∈D,n ∈ N where Ψ are the eigenfunctions of in L 2 (D/Γ). As we mentionned in the previous section, these eigenfunctions are not known explicitely. By restricting the "neutral stability surface" µ(ρ, β) to those values which correspond to eigenfunctions with Γ × β c Z + * periodicity, we obtain a discrete set of points on this surface with one minimum µ 0 associated with a value ρ 0 of ρ and β 0 of β. In general this minimum is unique. Moreover the multiplicity of the 0 eigenvalue is now finite and this eigenvalue is semi-simple. Let us call X the eigenspace associated with the 0 eigenvalue (therefore X is the kernel of the critical linear operator).

The full symmetry group of

D/Γ × S 1 is equal to G * × O(2)
where G * is the (finite) group of automorphisms in U (1, 1) of the Riemann surface D/Γ (notation of Section 2.2) and O(2) is the symmetry group of the circle (generated by S 1 and by reflection across a diameter). The equation restricted to this class of Γ × β c Z + * -periodic patterns is invariant under the action of G * × O(2). We can therefore apply an equivariant Lyapunov-Schmidt reduction to this bifurcation problem [START_REF] Chossat | Methods in Equivariant Bifurcations and Dynamical Systems[END_REF], leading to a bifurcation equation in

X f (x, µ)=0,x ∈ X ( 17
)
where

f : X × R → X is smooth, f (0, 0) = 0, ∂ x f (0, 0) is not invertible and f (•,µ) commutes with the action of G * × O(2) in X.
Now the methods of equivariant bifurcation theory can be applied to [START_REF] Dionne | Planforms in two and three dimensions[END_REF]. In particular we can apply the Equivariant Branching lemma (see [START_REF] Golubitsky | Singularities and Groups in Bifurcation Theory[END_REF] for a detailed exposition): 2) is absolutely irreducible in X (i.e. real equivariant linear maps in X are scalar multiple of the identity). Let H be an isotropy subgroup of G * × O(2) such that the subspace X H = {x ∈ X | H • x = x} is one dimensional. Then generically a branch of solutions of [START_REF] Dionne | Planforms in two and three dimensions[END_REF] bifurcates in X H . The conjugacy class of H (or isotropy type) is called "symmetry breaking".

Theorem 4. Suppose the action of G * × O(
Let us briefly recall the meaning of this theorem. By equivariance of f , any subspace of X defined as X H , H a (closed) subgroup of G * × O(2), is invariant under f . By the irreducibility assumption, if

H = G * × O(2), then {x ∈ X | H • x = x} = {0}.
Therefore f (0,µ)=0for all µ. Now the assumption of absolute irreducibility implies that ∂ x f (0,µ)=a(µ)Id X where a is a smooth real function such that a(0) = 0 and generically a (0) =0. It follows that if now H is a subgroup such that dim X H =1, then equation [START_REF] Dionne | Planforms in two and three dimensions[END_REF] restricted to this subspace reduces to a scalar equation 0=a (0)µx + xh(x, µ) (with h(0, 0) = 0), which has a branch of non trivial solutions by the implicit function theorem. Remark 1: the word "generically" can be interpreted as follows: the result will fail only if additional degeneracies are introduced in the equations. See [START_REF] Chossat | Methods in Equivariant Bifurcations and Dynamical Systems[END_REF] and [START_REF] Golubitsky | Singularities and Groups in Bifurcation Theory[END_REF] for a rigorous definition and proof. Remark 2: the assumption of absolute irreducibility is itself generic (in the above stated sense) for one-parameter steady-state bifurcation problems. A given irreducible representation of a compact or finite group need not be absolutely irreducible, this fact has to be proven. Remark 3: Theorem 4 does not necessarily give an account of all possible branches of solutions of [START_REF] Dionne | Planforms in two and three dimensions[END_REF], see [START_REF] Chossat | Methods in Equivariant Bifurcations and Dynamical Systems[END_REF]. It gives nevertheless a large set of generic ones. To go further it is necessary to compute the equivariant structure of f , or at least of its Taylor expansion to a sufficient order. The same is true if one wants to determine the stability of the bifurcated solutions, within the class of Γ × S 1 periodic solutions of the initial evolution equation.

This theorem, together with the knowledge of the lattices and the (absolutely) irreducible representations of the groups G * , gives us a mean to classify the periodic patterns which can occur in D×R + * . By analogy with the Euclidean case (bifurcation of spatially periodic solutions in the Euclidean space), we call H-planforms the solutions of a U (1, 1) (resp. GL(2R)) invariant bifurcation problem in D (resp. D×R + * ), which are invariant by a lattice group Γ (resp. Γ × S 1 ).

Being interested here in the classification of solutions rather than in their actual computation for a specific equation, all remains to do is to determine the absolutely real irreducible representations of the group G * × O(2) and the computation of the dimensions of the subspaces X H ). For this purpose we can get rid of the S 1 component of the domain of periodicity. Indeed let H = H 1 × H 2 be an isotropy subgroup for the representation R of G * × O(2) acting in X. Then X = V ⊗ W and R = S ⊗ T , where S is an irreducible representation G * in V and T is an irreducible representation of O(2) in W [START_REF] Serre | Représentations linéaires des groupes finis[END_REF], and therefore H 1 acts in V and H 2 acts in W . Now we have the following lemma, the proof of which is straightforward: Lemma 1. dim(X H )=1if and only if dim(V H1 )=dim(W H2 )=1. Now, the irreducible real representations of O(2) are well-known: they are either one dimensional (in which case every point is rotationally invariant) or two-dimensional, and in the latter case the only possible one dimensional subspaces W H2 are the reflection symmetry axes in R 2 (which are all equivalent under rotations in O(2)). It follows that the classification is essentially obtained from the classification of the isotropy subgroups H 1 of the irreducible representations of G * .

Once the irreducible representations are known, this can be achieved by applying the "trace formula" [START_REF] Golubitsky | Singularities and Groups in Bifurcation Theory[END_REF][START_REF] Chossat | Methods in Equivariant Bifurcations and Dynamical Systems[END_REF]:

Proposition 3. Let H be a subgroup of G * acting in a space V by a representation ρ : G * → Aut(V ), then dim(V H )= 1 |H| h∈H tr(ρ(h)) (18)
Note that tr(ρ) is the character of the representation ρ (a homomorphism G * → C). What is really needed to apply the Equivaraint Branching Lemma is therefore the character table of the representations. In the next section we investigate this classification in the case when the lattice is the regular octogonal group.

A case study: the octagonal lattice

According to the comment following Lemma 1, in all of this section we only consider the classification of isotropy subgroups satisfying the conditions of the equivaraint branching lemma in the Poincaré disc D.

The octagonal lattice and its symmetries

Among all lattices in the hyperbolic plane, the octagonal lattice is the simplest one. As before we use the Poincaré disc representation of the hyperbolic plane. Then the octagonal lattice group Γ is generated by the following four hyperbolic translations (boosts), see [START_REF] Balazs | Chaos on the pseudosphere[END_REF]:

g 0 = 1+ √ 2 2+2 √ 2 2+2 √ 21 + √ 2 (19)
and g j = r jπ/4 g 0 r -jπ/4 , j =1 , 2, 3, where r ϕ indicates the rotation of angle ϕ around the origin in D. The fundamental domain of the lattice is a regular octagon O as shown in Figure . The opposite sides of the octagon are identified by periodicity, so that the corresponding quotient surface D/Γ is isomorphic to a "double doughnut" (genus two surface) [START_REF] Balazs | Chaos on the pseudosphere[END_REF]. Note that the same octagon is also the fundamental domain of another group, not isomorphic to Γ, obtained by identifying not the opposite sides but pairs of sides as indicated in Figure. This is called the Gutzwiller octagon. A procedure of classification of the lattices using graphs is presented in [START_REF] Sausset | Periodic boundary conditions on the pseudosphere[END_REF]. For us however there is no difference between the two kinds of octagons because we are really interested in the full symmetry group of the pattern generated by Γ, which includes the rotations r jπ/4 , j =1 , ••• , 8, and therefore the boosts r π/2 g -1 0 , g 1 r -π/2 and their conjugates by the rotation r π , which are precisely the generators of the Gutzwiller lattice group.

We now determine what is the full symmetry group G * of the octagonal lattice, or equivalently, of the surface D/Γ. Clearly the symmetry group of the octagon itself is part of it. This is the dihedral group D 8 generated by the rotation r π/4 and by the reflection κ through the real axis, but there is more. We have seen in Section 2.2 that the group G * =Λ/Γ, Λ being the triangle group generated by reflections through the edges of a triangle τ which tiles (by the action of Λ/Γ) the surface D/Γ. The smallest triangle (up to symmetry) with these properties is the one shown in Figure 1. It has angles π/8, π/2 and π/3 at vertices P = O (the center of D), Q, R respectively, and its area is, by Gauss-Bonnet formula, equal to π/24. There are exactly 96 copies of τ filling the octagon, hence |G * | = 96. The index two subgroup G of orientation-preserving transformations in G * has therefore 48 elements.

In [START_REF] Broughton | Classifying finite group actions on surfaces of low genus[END_REF] it has been found that G GL(2, 3), the group of invertible 2 × 2 matrices over the 3 elements field Z 3 . In summary: The isomorphism between GL(2, 3) and G can be built as follows. We use the notation Z 3 = {0, 1, 2} and we call ρ the rotation by π/4 centered at P (mod Γ), σ the rotation by π centered at Q (mod Γ) and the rotation by 2π/3 centered at R (mod Γ). In the notations of Section 2.2, a = σ, b = , c = ρ, and ρσ =1. Then we can take

ρ = 02 22 ,σ= 20 01 ,= 21 20
since these matrices satisfy the conditions ρ 8 = σ 2 = 3 = Id and ρσ = Id. Note that ρ 4 = -Id where Id is the identity matrix. We shall subsequently use this notation. The group GL(2, 3), therefore the group G, is made of 8 conjugacy classes which we list in Table 1, indicating one representative, the number of elements in each class and their order. This result is classical and can be found, e.g., in [START_REF] Lang | Algebra. Addison-Wesley[END_REF].

representative Id ρ ρ We now turn to the full symmetry group G * which is generated by G and κ, the reflection through the real axis in D and which maps the octagon O to itself. We write κ = ρκ the reflection through the side PR of the triangle τ . Note that (i) κ preserves also O, (ii) κ = κ = σκ is the reflection through the third side QR.

In what follows we rely on the group algebra software GAP [START_REF] Gap | Groups, Algorithms and programming[END_REF]. For this we have first identified a presentation for G GL(2, 3) considered as an abstract group, then a presentation for G * . The presentation for GL(2, 3) can be obtained with the command "P := PresentationViaCosetTable(GL(2,3))" and the relations are shown with the command "TzPrintRelators(P)": Applying the above lemma we find with GAP that the 96 elements group G * has 13 conjugacy classes which are listed in table 2 for direct isometries and in table 3 for isometries which reverse orientation. GAP gives representatives of the conjugacy classes in the abstract presentation, which in general have complicated expressions. In some cases we have chosen other representatives, using in particular the 8-fold generator. To simplify some expressions in the tables 2 to 5 we also use the notations

σ = σ -1 , σ = ρ 2 σρ -2 ,
where σ is the rotation by π centered at S (mod Γ) and σ is the rotation by π centered at S (mod Γ), see figure 2. We shall also need in Section 4.3 the list of subgroups of G * together with their decomposition in conjugacy classes (in G * ) in order to apply the trace formula [START_REF] Ermentrout | Neural networks as spatio-temporal pattern-forming systems[END_REF]. Here again we rely on GAP to obtain the necessary informations.

Then representatives of each class are determined by inspection. These datas are listed in tables 4 (subgroups of G) and 5 (subgroups containing orientation reversing elements). The subgroups are listed up to conjugacy in G * , the subgroups of order two are not listed. The rationale for the notations is as follows:

class
} C 4 4 <ρ 2 > {Id (1), -Id (1), ρ 2 (2)} D 2 4 < -Id,σ > {Id (1), -Id (1), σ (2)} C 3 3 <> {Id (1), (2)} C 2 2 < -Id > {Id (1), -Id (1)} C 2 2 <σ> {Id (1), σ (1)} Table 4: Subgroups of G ⊂ G * (up to conjugacy). (3)
The last column provides datas about their conjugacy subclasses (in G * ).

• G 0 is an index 2 subgroup of G. Seen as a subgroup of GL(2, 3) it is SL(2, 3), the subgroup of determinant 1 matrices. It contains no order 2 elements except -Id and no order 8 elements.

• C n , C n , C n , denote order n cyclic groups. The notation C n is standard for the n-fold rotation group centered at the origin.

• D n denotes a group isomorphic to the dihedral group of order 2n, generated by an n-fold rotation and a reflection. Hence D 8 is the symmetry group of the octagon. The notation D n is used for a 2n element group which has an n-fold rotation and a 2-fold rotation as generators. For example D 8 =< ρ, σ> , and one can verify that σρ σ -1 = ρ 3 , which makes D 8 a quasidihedral group (see [START_REF] Gorenstein | Finite groups[END_REF]).

• Q 8 is a usual notation for the 8 elements quaternionic group.

• The notation H nκ indicates a group generated by the group H n and κ. Same thing if replacing κ by κ . For example C 3κ is the 6 elements group generated by C 3 and κ .

The irreducible representations of G *

There are 13 conjugacy classes and therefore we know there are 13 complex irreducible representations of G * , the characters of which will be denoted χ j , j =1,...,13. The character table, as computed by GAP, is shown in table 6.

The character of the identity is equal to the dimension of the corresponding representation. It follows from table 6 that there are 4 irreducible representations of dimension 1, 2 of dimension 2, 4 of dimension 3 and 3 of dimension 4. In the following we shall denote the irreducible representations by their character: χ j is the representation with this character.

Lemma 3. All irreduclible representations of G

} D 4κ 16 < D 4 ,κ > D 4 ∪{κ (2), κ (4), ρ σκ (2)} C 12 12 <κ> C 6 ∪{κ (2), -κ (2), ρ σκ (2)} C 6κ 12 < C 6 ,κ > C 6 ∪{κ (6)} C 8 8 < σκ > C 4 ∪{ σκ (4)} C 4κ (= D 4 ) 8 <C 4 ,κ > C 4 ∪{κ (4)} C 4κ 8 <C 4 ,κ > C 4 ∪{κ (4)} D 2κ 8 < D 2 ,κ > D 2 ∪{κ (2), κ (2)} C 4κ 8 <C 4 ,κ > C 4 ∪{ρ 2 (2), κ (2)} C 4κ 8 <C 4 ,κ > C 4 ∪{σ (2), κ (2)} C 3κ 6 < C 3 ,κ > C 3 ∪{κ (3)} C 4 4 <ρ σκ > {Id (1), -Id (1), ρ σκ (2)} C 2κ 4 < -Id,κ > {Id (1), -Id (1), κ (2)} C 2k 4 < -Id,κ > {Id (1), -Id (1), κ (2)} C 2k 4 <σ,κ> {Id (1), σ (1), κ (1), κ (1)} C 2κ 4 < σ, κ > {Id (1), σ (1), κ (1), κ (1)} C 1κ 2 <κ> {Id (1), κ (1)} C 1κ 2 <κ > {Id (1), κ (1)} (4)
-Id σ -κ κ σκ ρ σκ κ -κ χ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 χ 2 1 -1 1 1 -1 1 1 1 -1 -1 1 1 1 χ 3 1 -1 1 1 -1 1 1 -1 1 1 -1 -1 -1 χ 4 1 1 1 1 1 1 1 -1 -1 -1 -1 -1 -1 χ 5 2 0 2 2 0 -1 -1 -2 0 0 -2 1 1 χ 6 2 0 2 2 0 -1 -1 2 0 0 2 -1 -1 χ 7 3 1 -1 3 -1 0 0 -1 -1 1 3 0 0 χ 8 3 1 -1 3 -1 0 0 1 1 -1 -3 0 0 χ 9 3 -1 -1 3 1 0 0 1 -1 1 -3 0 0 χ 10 3 -1 -1 3 1 0 0 -1 1 -1 3 0 0 χ 11 4 0 0 -4 0 -2 2 0 0 0 0 0 0 χ 12 4 0 0 -4 0 1 -1 0 0 0 0 √ 3 - √ 3 χ 13 4 0 0 -4 0 1 -1 0 0 0 0 - √ 3 √ 3 Table 6: Irreducible characters of G *
Proof. This is clear for the one dimensional representations whose characters are real.

For the two dimensional representations, let us consider the dihedral subgroup D 3 generated by the 3-fold symmetry and the reflection κ . The representation of D 3 in either representations planes of χ 5 and χ 6 have characters χ j ()=-1 and χ j (κ )=0(j =5or 6). These are the characters of the 2D irreducible representation of D 3 , which is absolutely irreducible, and therefore the representations χ 5 and χ 6 of G * are also absolutely irreducible. Indeed if any real linear map which commutes with the elements of a subgroup is a scalar multiple of the identity, this is a fortiori true for the maps which commute with the full group.

For the three dimensional representations χ 7 to χ 10 , let us first remark that if we write C 2 = {Id,-Id}, then G/C 2 O, the octahedral group. Its subgroup T (tetrahedral group) can easily be identified with the 12 elements group generated by the "pairs" {Id,-Id}, {, -} and {ρ 2 , -ρ 2 }. Now we consider the representation of G defined by the action of χ j restricted to G (for each 3D χ j ). One can check easily from the character table that it projects onto a representation of G * /C 2 , the character of which is given by the value of χ j on the corresponding conjugacy classes, and in particular the character for the representation of the group T is given, for any j =7to 10, by χ j ({Id,-Id})=3, χ j ({, -})=0and χ j ({ρ 2 , -ρ 2 })=-1. But this is the character of the irreducible representation of T [START_REF] Miller | Symmetry groups and their applications[END_REF], which is absolutely irreducible (natural action of T in R 3 ). Hence the three dimensional representations of G * are absolutely irreducible by the same argument as above.

It remains to prove the result for the four dimensional representations χ 11 , χ 12 and χ 13 . For this we consider the action of the group D 8 generated by ρ and κ, as defined by either one of these 4D irreducible representations of G * . We observe from the character table that in all cases, the character of this action is χ(ρ)=0 , χ(ρ 2 )=0 , χ(-Id)=-4, χ(ρ 3 )=0(ρ and ρ 3 are conjugate in G * ), and χ(κ)=χ(κ )=0. We can determine the isotypic decomposition for this action of D 8 from these character values. The character tables of the four one dimensional and three two dimensional irreducible representations of D 8 can be computed easily either by hand (see [START_REF] Miller | Symmetry groups and their applications[END_REF] for the method) or using a computer group algebra software like GAP. For all one dimensional characters the value at -Id is 1, while for all two dimensional characters, the value at -Id is -2. Since χ(-Id)=-4, it is therefore not possible to have one dimensional representations in this isotypic decomposition. It must therefore be the sum of two representations of dimension 2. Moreover, since χ(ρ)=χ(ρ 2 )=χ(ρ 3 )=0 , it can't be twice the same representation. In fact it must be the sum of the representations whose character values at ρ are √ 2 and -√ 2 respectively. Now, these representations are absolutely irreducible (well-know fact which is straightforward to check), hence any D 8 -equivariant matrix which commutes with this action decomposes into a direct sum of two scalar 2×2 matrices λI 2 and µI 2 where λ and µ are real. But the representation of G * is irreducible, hence λ = µ, which proves that it is also absolutely irreducible.

The octagonal H-planforms

We can now apply Lemma 1 in order to determine the H-planforms for the octagonal lattice.

Theorem 5. The irreducible representations of G * admit H-planforms with the following isotropy types:

• χ 1 : G * ; • χ 2 : G 0κ ; • χ 3 : G 0κ ; • χ 4 : G GL(2, 3); • χ 5 : D 8 , Q 8κ ; • χ 6 : D 8κ ; • χ 7 : C 8κ , C 12 , C 4κ ; • χ 8 : C 8κ , C 6κ , D 2κ ; • χ 9 : D 6 , D 4κ ; • χ 10 : D 6κ , D 4κ ; • χ 11 : C 2κ , C 2κ ; • χ 12 : D 3 , C 3κ , C 2κ , C 2κ ; • χ 13 : D 3 , C 3κ , C 2κ , C 2κ ;
Proof. For the one dimensional representations of G * this is straightforward: each element whose character image is +1 belongs to the isotropy group. The result follows therefore directly from the character table and list of subgroups of G * . For the higher dimensional irreducible representations we need to find those isotropy subgroups H such that (see [START_REF] Ermentrout | Neural networks as spatio-temporal pattern-forming systems[END_REF]):

1=dim(V H )= 1 |H| h∈H χ j (h)
where χ j denotes the character of the j-th irreducible representation. This can be done in a systematic way by using the character table 6 and applying the datas on subgroups and their conjugacy classes listed in 4 and 5. The calculations are cumbersome but can be slightly simplified by noting that if H ⊂ H and dim(V H )=0(a case which occurs many times), then not only H is not symmetry breaking but also H , since

H ⊂ H ⇒ V H ⊂ V H .
The following lemma is also useful as it eliminates most candidates in the case of 4D representations:

Lemma 4. If a subgroup H contains -Id, then for j = 11, 12 or 13, one has h∈H χ j (h)=0. Proof of the lemma. In all three cases the result follows from the relations: (i) χ j (-Id)=-χ j (Id), (ii) χ(-)= -χ() and χ(-κ)=-χ(κ), (iii) χ j (s)=0for all s which is not conjugate to one in (i) or (ii).

For the lower dimensional representations it is possible to reduce the problem to known situations and to provide bifurcation diagrams without any further calculations. The next theorem provides these informations. Stability of the solutions has to be understood here with respect to perturbations with the same octagonal periodicity in D and under the condition that, in L 2 (D/Γ), the corresponding representation corresponds to the "most unstable" modes ("neutral modes" at bifurcation). Theorem 6. For the one and two dimensional representations, the generic bifurcation diagrams have the following properties:

• χ 1 : transcritical branch, exchange of stability principle holds.

• χ 2 , χ 3 and χ 4 : pictchfork bifurcation, exchange of stability principle holds.

• χ 5 : same as bifurcation with hexagonal symmetry in the plane, see Figure 3.

• χ 6 : same as bifurcation with triangular symmetry in the plane. In particular H-planforms are always unstable on both sides of the bifurcation point unless the subcritical branch bends back sufficiently near the bifurcation point (see Figure 4).

Proof. (i) For the one dimensional representations, this follows from classical bifurcation theory: in χ 1 there is no symmetry breaking, hence generically the bifurcation is of transcritical type and the trivial and bifurcated solutions exchange stability at the bifurcation point. In the three other cases, a symmetry exists which acts by reversing direction on the axis as can be seen from Table 6. For example in χ 2 this can be taken as σ (but also κ does the same thing). Hence the bifurcation is of pitchfork type and exchange of stability holds.

(ii) For χ 5 , note that the subgroup Q 8 acts trivially on any point of this plane (dim(V Q8 )=2 ). In fact Q 8 is the isotropy group of the principal stratum in this group action. Now,

G GL(2, 3) = Q 8 D 3 , hence G * /Q 8 D 3 Z 2 D 6
, the symmetry group of an hexagon. This group action is isomorphic to the natural action of D 6 in the plane. It follows that the problem reduces in this case to a bifurcation problem with the action of D 6 in the plane, see [START_REF] Golubitsky | Singularities and Groups in Bifurcation Theory[END_REF] for details.

(iii) In the case of χ 6 , the maximal subgroup which keeps every point in the plane fixed is the 16 element group Q 8κ . It follows that the problem reduces to a bifurcation problem in the plane with symmetry G * /Q 8κ D 3 . Details on this bifurcation can be found in [START_REF] Golubitsky | Singularities and Groups in Bifurcation Theory[END_REF].

Remark 1: A similar reduction can be made with the 3 dimensional representations. Indeed it can be seen that the principal isotropy type (which keeps all points in the three dimensional representation space fixed) is C 4 for χ 7 and χ 10 , and C 2 for χ 8 and χ 9 . In the first case, this leads to reducing the problem to one with G * /C

4 O symmetry, where O is the 24 elements group of direct symmetries (rotations) of a cube. However there are two irreducible representations of dimension 3 of O [START_REF] Miller | Symmetry groups and their applications[END_REF], and it turns out that χ 7 corresponds to one of them (the "natural" action of O in R 3 ) while χ 10 corresponds to the other representation. This explains why there are 3 types of H-planforms for χ 7 and only 2 for χ 10 . Similarly, the principal isotropy type for χ 8 and χ 9 is the two element group C 2 , and G * /C 2 OZ 2 . Then the same remark holds for these cases as for the previous ones.

Remark 2: In the 4 dimensional cases, the principal isotropy type is the trivial group, hence no reduction can be made. Bifurcation in this case (and in the 3 dimensional cases as well) will be the subject of a forthcoming paper.

Computing the H-planforms

It follows from the definition that H-planforms are eigenfunctions of the Laplace-Beltrami operator in D which satisfy certain isotropy conditions: (i) being invariant under a lattice group Γ and (ii) being invariant under the action of an isotropy subgroup of the symmetry group of the fundamental domain D/Γ (mod Γ). Therefore in order to exhibit H-planforms, we need first to compute eigenvalues and eigenfunctions of in D, and second to find those eigenfunctions which satisfy the desired isotropy conditions. In this section we tackle this question in the case where the lattice has the regular octagon as a fundamental domain.

Over the past decades, computing the eigenmodes of the Laplace-Beltrami operator on compact manifolds has received much interest from physicists. The main applications are certainly in quantum chaos [START_REF] Balazs | Chaos on the pseudosphere[END_REF][START_REF] Aurich | Periodic-orbit sum rules for the hadamard-gutzwiller model[END_REF][START_REF] Aurich | Statistical properties of highly excited quantum eigenstates of a strongly chaotic system[END_REF][START_REF] Schmit | Quantum and classical properties of some billiards on the hyperbolic plane[END_REF][START_REF] Cornish | Ringing the eigenmodes from compact manifolds[END_REF] and in cosmology [START_REF] Inoue | Computation of eigenmodes on a compact hyperbolic 3-space[END_REF][START_REF] Cornish | On the eigenmodes of compact hyperbolic 3-manifolds[END_REF][START_REF] Lehoucq | Eigenmodes of 3-dimensional spherical spaces and their application to cosmology[END_REF].

To our knowledge, the interest in such computation was sparkled by the study of classical and quantum mechanics on surfaces of constant negative curvature, and the connections between them (for an overview on the subject see [START_REF] Balazs | Chaos on the pseudosphere[END_REF]). To be more precise, quantum chaology can be defined as the study of the semiclassical behaviour characteristic of systems whose classical motion exhibits chaos, for example the classical free motion of a mass point on a compact surface of constant negative curvature (as it is the most chaotic possible). In [START_REF] Balazs | Chaos on the pseudosphere[END_REF][START_REF] Aurich | Periodic-orbit sum rules for the hadamard-gutzwiller model[END_REF][START_REF] Cornish | Ringing the eigenmodes from compact manifolds[END_REF], authors studied the time-independent Schroedinger equation on the compact Riemannian surface of constant curvature -1 and genus 2, which is topologically equivalent to the regular octagon with four periodic boundary conditions. This is the same as solving the eigenvalue problem for Γ invariant eigenmodes in D. The first computations have been performed using the finite element method on "desymmetrised" domains of the hyperbolic octagon with a mixture of Dirichlet and Neumann boundary conditions [START_REF] Balazs | Chaos on the pseudosphere[END_REF][START_REF] Schmit | Quantum and classical properties of some billiards on the hyperbolic plane[END_REF]. We explain the procedure of desymmetrisation in the next subsection. Aurich and Steiner in [START_REF] Aurich | Periodic-orbit sum rules for the hadamard-gutzwiller model[END_REF] were the first to compute the eigenmodes on the whole octagon with periodic boundary conditions. They began with the finite element method of type P 2 and were able to exhibit the first 100 eigenvalues. In [START_REF] Aurich | Statistical properties of highly excited quantum eigenstates of a strongly chaotic system[END_REF], the same authors used the direct boundary-element method on an asymmetric octagon to reach the 20 000th eigenvalue.

There is also a strong interest of cosmologists for ringing the eigenmodes of the Laplace-Beltrami operator on compact surfaces. Indeed this is necessary in order to evaluate the cosmic microwave background anisotropy in multiply-connected compact cosmological models. For some models, this computation is performed on a compact hyperbolic 3-space called the Thurston manifold, and Inoue computed the first eigenmodes of Thurston space such that each corresponding eigenvalue λ satisfies λ ≤ 10 with the direct boundary-element method [START_REF] Inoue | Computation of eigenmodes on a compact hyperbolic 3-space[END_REF]. For 3-dimensional spherical spaces, several methods have been proposed: the "ghosts method" [START_REF] Cornish | On the eigenmodes of compact hyperbolic 3-manifolds[END_REF], the averaging method and the projection method. All these methods are explained and summarized in [START_REF] Lehoucq | Eigenmodes of 3-dimensional spherical spaces and their application to cosmology[END_REF].

Our aim is different in that we do not want to compute all the eigenvalues of the Laplace-Beltrami operator, but instead to calculate the H-planforms with the isotropy types listed in theorem 5. The methods of numerical computation are however similar, and one question is to choose the method best suited to our goal. For the H-planforms associated to irreductible representations of dimension 1 (i.e. for (χ i ) i=1•••4 ), we use a desymmetrization of the octagon with a reformulation of the boundary conditions. For H-planforms associated with irreductible representations of dimension ≥ 2, the desymmetrization of the octagon is also possible but much more complicated as noticed by Balazs and Voros [4]. This will therefore be the subject of another paper. Here we only identify some H-planforms of specific isotropy types. In order to find these H-planforms, we use the finite-element method with periodic boundary conditions. This choice is dictated by the fact that this method will allow us to compute all the first n eigenmodes and among all these we will indentify those which correspond to a given isotropy group. As explained before, if we have used the direct boundary-element method, we would have reached any eigenmode but the search for H-planforms would also have became more random. Indeed, each iteration of this method gives only one eigenmode while the finite-element method provides n eigenmodes depending on the precision of the discretization. This is why we prefer to use this last method in order to find some H-planforms associated with irreductible representations of dimension ≥ 2, although it is quite more complicated to implement because of the periodic boundary conditions.

Desymmetrization of the octagon

We have already seen that the fundamental domain T (2, 3, 8) of the group G * generates a tiling of the octagon. Desymmetrization consists in separating the individual solutions according to the symmetry classes of G * . This entails solving the eigenvalue problem in certain irreducible subregions of the fundamental domain, such as T (2, 3, 8), using special boundary conditions for these subregions. In effect the periodicity conditions in the original domain (octagon) may produce Dirichlet or Neumann conditions on the boundaries of these subregions. The symmetry group G * has a smaller fundamental domain, the triangle T (2, 3, 8) see figure 1, which is 1 96 th of the original octagon. The method of desymmetrization can be applied to many other techniques than finite element methods (see [START_REF] Fässler | Group theoretical methods and their applications[END_REF]).

We now focus on the four one-dimensional irreductible representations (χ i ) i=1•••4 acting upon the generators as indicated in table 6. One can find in the book of Fässler and Stiefel [19,Chapter 3] the principle of desymmetrization in the context of dihedral symmetry. We follow their method in the case of the symmetry group G * . The first step is to attribute one number (value) to each of the 96 triangles that tesselate the octagon under the action of G * (see 2), according to the character values obtained from table 6, i.e. ±1 depending on the conjugacy class (remember we restrict ourselves to the one dimensional representations χ 1 to χ 4 ). Let us take the example of the first irreductible representation χ 1 and explain how we obtain the domain and the boundary conditions depicted in figure 5. Table 6 shows that all 96 triangles end up with the same value, 1. This means that the eigenfunction we are looking for is even under all the 96 elements in G * and it follows that it must satisfy Neumann boundary conditions on all the edges of the tesselation of the hyperbolic octagon [START_REF] Balazs | Chaos on the pseudosphere[END_REF][START_REF] Aurich | Periodic-orbit sum rules for the hadamard-gutzwiller model[END_REF]. Finally, it is sufficient to solve the eigenproblem on the reduced domain T (2, 3, 8) with Neumann boundary conditions on its three edges. For the four one-dimensional representations one has to choose the correct combination of Neumann and Dirichlet boundary conditions as shown in figure 5.

The representations of dimension ≥ 2 require the same number of values as their dimension. For example there are two basis vectors determining the function values in the case of an irreductible representation of dimension 2. The table 6 is then no longer sufficient to set the values of the function on each triangles and one has to explicitly write the matrices of the irreductible representation in order to obtain the suitable conditions. This is why we have restricted ourselves to the four one-dimensional representations as described in the previous paragraphs.

Numerical experiments

As there exists an extensive literature on the finite element methods (see for an overview [START_REF]Handbook of Numerical Analysis[END_REF][START_REF] Allaire | Analyse numérique et optimisation[END_REF]) and as numerical analysis is not the main goal of this article, we do not detail the method itself but rather focus on the way to actually compute the eigenmodes of the Laplace-Beltrami operator.

Desymmetrized problem : For the four problems depicted in figure 5, we use the mesh generator Mesh2D from Matlab to tesselate the triangle T (2, 3, 8) with 2995 nodes and we implement the finite element method of order 1.

Our results are presented in figure 6 and are in a good agreement with those obtained by Balazs-Voros in [START_REF] Balazs | Chaos on the pseudosphere[END_REF] and Aurich-Steiner in [START_REF] Aurich | Periodic-orbit sum rules for the hadamard-gutzwiller model[END_REF]. Once we have computed the eigenfunction in T (2, 3, 8), we extend it to the whole octagon by applying the generators of G * . We superimpose in figure 6(a) the tesselation of the octagon by the 96 triangles in order to allow the reader to see the symmetry class of G * .

Non desymmetrized problem : As discussed previously, we also present some H-planforms of higher dimension. We mesh the full octagon with 3641 nodes in such a way that the resulting mesh enjoys a D 8 -symmetry, see figure 7. We implement, in the finite element method of order 1, the periodic boundary conditions of the eigenproblem and obtain the first 100 eigenvalues of the octagon. Our results are in agreement with those of Aurich and Steiner reported in [START_REF] Aurich | Periodic-orbit sum rules for the hadamard-gutzwiller model[END_REF]. Instead of giving a table of all the eigenvalues, we prefer to plot the staircase function N (λ)={λ n |λ n ≤ λ} for comparison with Weyl's law. Weyl's law is, in its simplest version, a statement on the asymptotic growth of the eigenvalues of the Laplace-Beltrami operator on bounded domains. If Ω is a given bounded domain of R 2 , then the staircase function has the following asymptotic behaviour: N (λ)= |Ω| 4π λ + o(λ) as λ →∞. We recall that in the case of the hyperbolic octagon, one has |Ω| =4π and hence N (λ) ∼ λ as λ →∞.

As can be seen in figure 8 the asymptotic law describes the staircase well down to the smallest eigenvalues, which confirms the validity of our numerical results.

We show in figure 9 two H-planforms, with D 8 and D 8κ isotropy respectively. These two H-planforms belong to irreductible representations of dimension 2: χ 5 for 9(a) and χ 6 for 9(b) and 10(a).

We finally present in figure 10 In figures 6,9 and 10, we have plotted, for convenience, the corresponding H-planforms in the octagon only. Nevertheless, H-planforms are periodic in the Poincaré disc, as stated before, and in figure 11, we plot the Hplanform with G 0κ isotropy type of figure 6(c). We recall that the octagonal lattice group Γ is generated by the four boosts g j of subsection 4.1. Then, once the H-planform is calculated, we report it periodically in the whole Poincaré disc by the actions of there four boosts and obtain figure 11. Note that it is arduous to tesselate the entire disc and this is why there remains some untesselated areas in the figure .   Remark : The finite element method relies on a variational principle and thus yields upper bounds for the eigenvalues. Note that it is not well-suited to compute high eigenvalues. Thus, if one wants to compute, for example, the first 125 eigenvalues, then the finite element method provides a very good approximation to the true eigenvalues. If however one wants to reach the 2000th eigenvalue by the finite element method we need to use matrices of size at least 10 6 × 10 6 which is untractable for desktop or laptop computers and requires going to more powerful architectures. This is why the direct boundary element method should be preferred for such computations (see [START_REF] Aurich | Statistical properties of highly excited quantum eigenstates of a strongly chaotic system[END_REF]).

Conclusion

In this paper, we have analysed the bifurcation of periodic patterns for neural field equations describing the state of a system defined on the space of structure tensors, when these equations are further invariant under the isometries in this space. We have made use of the concept of periodic lattice in D (Poincaré disc) to further reduce the problem to one on a compact Riemann surface D/Γ, where Γ is a cocompact, torsion-free Fuchsian group. Successfully, we have applied the machinery of equivariant bifurcation theory in the case of an octagonal periodic pattern, where we have been able to classify all possible H-planforms satisfying the hypotheses of the Equivariant Branching Lemma. In the last section, we have described a method to compute these patterns and illustrated it with a selection of images of octagonal H-planforms.

There are several questions which are raised by this analysis.

(i) The first one is that of the interpretation of these H-planforms for the modelling of the visual cortex. At this moment we believe that they could be involved in the process of defining texture tuning curves in a way that would resemble the definition of orientation tuning curves in the related ring model of orientation selectivity [START_REF] Hansel | Modeling feature selectivity in local cortical circuits[END_REF][START_REF] Shriki | Rate models for conductance-based cortical neuronal networks[END_REF][START_REF] Ermentrout | Neural networks as spatio-temporal pattern-forming systems[END_REF][START_REF] Dayan | Theoretical Neuroscience : Computational and Mathematical Modeling of Neural Systems[END_REF][START_REF] Bressloff | Dynamical mechanism for sharp orientation tuning in an integrate-and-fire model of a cortical hypercolumn[END_REF][START_REF] Bressloff | Geometric visual hallucinations, euclidean symmetry and the functional architecture of striate cortex[END_REF][START_REF] Veltz | Local/global analysis of the stationary solutions of some neural field equations[END_REF]. Another fascinating possibility is that they could be related to neural illusions caused by the existence of several stable stationary solutions to (1) when, e.g., the slope of the sigmoid at the origin becomes larger than that at which several branches of solutions bifurcate from the trivial one. These neural illusions would be functions of membrane potential values that would not correspond to the actual thalamic input and could be expressed as combinations of these H-planforms. This is of course still very speculative but very much worth investigating.

(ii) The second one is about the observability of such patterns in a natural system or under direct simulation of the evolution equations. Indeed, not only there is a high degeneracy of the bifurcation problem if one removes the assumption that all perturbations respect the periodicity of the pattern, which is also the case for patterns in Euclidean space, but in addition the fact that such patterns would be neutral modes for the bifurcation problem posed in full generality in D is non generic. We may however imagine mechanisms such as "spatial frequency locking" by which periodic patterns could become stable, hence observable.

(iii) The third one is that of a more effective computation of H-planforms for a given isotropy type. As we have seen in the last section, a desymmetrization of the domain allowed us to calculate all the H-planforms with isotropy types associated to irreductible representations of dimension 1. For irreductible representations of dimension ≥ 2 the computation becomes more intricate as the desymmetrization method is no longer straightforward and remains misunderstood. Naturally, one solution would be to elaborate a general algorithm which, for a given isotropy type, computes systematically the associated H-planform. We think that such an algorithm would be of interests for quantum physisists and cosmologists. These questions will be the subject of further studies.

Figure 1 :Figure 2 :

 12 Figure 1: The triangle τ , also noted T (2, 3, 8). The values of l, m, and n are l =8, m =2and n =3.

Proposition 4 .

 4 The full symmetry group G * of D/Γ is G ∪ κG where G GL(2, 3) has 48 elements.

Lemma 2 .

 2 (i) As an abstract group, G is presented with two generators a and b and three relations a 2 =1 , b 3 =1and (abab -1 ab -1 ) 2 =1 . (ii) As an abstract group, G * is presented with three generators a, b and c and six relations: the three relations for G plus the three relations c 2 =1, (ca) 2 =1and (cb) 2 =1. (iii) These abstract elements can be identified with automorphisms of D/Γ as follows: a = σ, b = and c = κ .

Figure 3 :

 3 Figure 3: Bifurcation diagram for the case χ 5 . Dotted lines: unstable branches, .

Figure 4 :

 4 Figure 4: Bifurcation diagram for the case χ 6 . Dotted lines: unstable branches.

(a) χ 1 :

 1 G * . (b) χ 2 : G 0κ . (c) χ 3 : G 0κ . (d) χ 4 : G.

Figure 5 :

 5 Figure 5: Boundary conditions for the one-dimensional irreducible representations. Top left: boundary conditions for χ 1 , corresponding to the isotropy group G * . Top right: boundary conditions for χ 2 corresponding to the isotropy group G Oκ . Bottom left: boundary conditions for χ 3 corresponding to the isotropy group G 0κ . Bottom right: boundary conditions for χ 4 corresponding to the isotropy group G.

(a) χ 1 :

 1 G * , the corrseponding eigenvalue is λ = 23.0790. (b) χ 2 : G 0κ , the corrseponding eigenvalue is λ = 91.4865. (c) χ 3 : G 0κ , the corrseponding eigenvalue is λ = 32.6757. (d) χ 4 : G, the corresponding eigenvalue is λ = 222.5434.

Figure 6 :

 6 Figure 6: The four H-planforms with their corresponding eigenvalue associated with the four irreductible representations of dimension 1, see text.

Figure 7 :

 7 Figure 7: A mesh of the octagon with 3641 nodes which is used for the method of finite elements. It leads to matrices of dimension 3641 × 3641.

  three H-planforms, with C 8κ , D 4κ and D 4κ isotropy. These three H-Planforms belong to irreductible representations of dimension 3: χ 8 for 10(a), χ 9 for 10(b) and χ 10 for 10(c).

Figure 8 :

 8 Figure 8: The staircase function N (λ), in dark, is shown in comparison with Weyl's law N (λ) ∼ λ as λ →∞.

(a) χ 5 :

 5 D 8 , the corrseponding eigenvalue is λ = 73.7323. (b) χ 6 : D 8κ , the corrseponding eigenvalue is λ = 8.2501.

Figure 9 :

 9 Figure 9: Two H-planforms with their corresponding eigenvalue associated to the two irreductible representations of dimension 2, see text.

(a) χ 8 :

 8 C 8κ , the corrseponding eigenvalue is λ = 3.8432. (b) χ 9 : D 4κ , the corresponding eigenvalue is λ = 28.0888. (c) χ 10 : D 4κ , the corrseponding eigenvalue is λ = 15.0518.

Figure 10 :

 10 Figure 10: Three H-planforms with their corresponding eigenvalue associated to three irreductible representations of dimension 3, see text.

Figure 11 :

 11 Figure 11: Extension of the G 0κ H-planform on the Poincaré disk.

Table 1 :

 1 Conjugacy classes of G GL[START_REF] Aurich | Periodic-orbit sum rules for the hadamard-gutzwiller model[END_REF][START_REF] Aurich | Statistical properties of highly excited quantum eigenstates of a strongly chaotic system[END_REF]

	order	2 -Id ρ 5 σ -1 8 4 2 8 2 3 6
	# elements	1 6 6	1	6 12 8	8

Table 2 :

 2 Conjugacy classes of G * , orientation preserving transformations

	number	1	2	3	4	5 6	7
	representative Id ρ ρ 2 -Id σ -order 1 8 4 2 2 3 6
	# elements	1 12 6	1	12 8	8

Table 3 :

 3 Conjugacy classes of G * , orientation reversing transformations

	Subgroup	Order	Generators	Subclasses: representatives (# elements)
	G 0 SL(2, 3) D 8 D 6 C 8 Q 8 D 4 C 6 D 3	24 16 12 8 8 8 6 6	<ρ 2 ,> < ρ, σ> < -, σ> <ρ> <ρ 2 ,σ ρ 2 σ> <ρ 2 , σ> < -> <, σ>	{Id (1), -Id (1), ρ 2 (6), (8), -(8)} σ (4)} {Id (1), -Id (1), ρ (4), ρ 2 (6), {Id (1), -Id (1), σ (6), (2), -(2)} {Id (1), -Id (1), ρ (4), ρ 2 (2)} {Id (1), -Id (1), ρ 2 (6)} σ (4)} {Id (1), -Id (1), ρ 2 (2), {Id (1), -Id (1), (2), -(2)} {Id (1), (2), σ

Table 5 :

 5 Subgroups of G

	Class #	1	2	3	4	5	6	7	8	9 10	11	12	13
	Representative Id ρ ρ 2

*

, not in G (up to conjugacy). The last column provides datas about their conjugacy subclasses (in G * ).

C 4 is a subgroup conjugate to C 4 with generator (ρ 2 σ) 2 .

Acknowledgments

This work was partially funded by the ERC advanced grant NerVi number 227747.