293 | | #27637 |
294 | | |
| 293 | Sage 9.6 allows for the calculation of the H*-series (also known as the φ-series) of a rational polytope which is invariant under the linear action of a finite group, as introduced by Alan Stapledon as part of equivariant Ehrhart theory #27637. |
| 294 | The computation is made with the method `Hstar_function` of a rational polytope. |
| 295 | The fixed subset of a polytope under the action of a group element may also be computed via the methods `fixed_subpolytope` or `fixed_subpolytopes`. |
| 296 | |
| 297 | As an example, consider the action of the symmetric group S3 on the 2-dimensional permutahedron in 3-dimensional space, given by permuting the three basis vectors. |
| 298 | As shown by Ardila, Supina, and Vindas Meléndez, the corresponding H*-series is polynomial and effective: |
| 299 | {{{ |
| 300 | sage: p3 = polytopes.permutahedron(3, backend = 'normaliz') |
| 301 | sage: G = p3.restricted_automorphism_group(output='permutation') |
| 302 | sage: reflection12 = G([(0,2),(1,4),(3,5)]) |
| 303 | sage: reflection23 = G([(0,1),(2,3),(4,5)]) |
| 304 | sage: S3 = G.subgroup(gens=[reflection12, reflection23]) |
| 305 | sage: S3.is_isomorphic(SymmetricGroup(3)) |
| 306 | True |
| 307 | sage: p3.Hstar_function(S3, output='complete') |
| 308 | {'Hstar': chi_0*t^2 + (chi_0 + chi_1 + chi_2)*t + chi_0, |
| 309 | 'Hstar_as_lin_comb': (t^2 + t + 1, t, t), |
| 310 | 'conjugacy_class_reps': [(), (0,1)(2,3)(4,5), (0,3,4)(1,5,2)], |
| 311 | 'character_table': [ 1 1 1] |
| 312 | [ 1 -1 1] |
| 313 | [ 2 0 -1], |
| 314 | 'is_effective': True} |
| 315 | sage: p3.fixed_subpolytope(reflection12) |
| 316 | A 1-dimensional polyhedron in QQ^3 defined as the convex hull of 2 vertices |
| 317 | sage: p3.fixed_subpolytope(reflection12).vertices() |
| 318 | (A vertex at (3/2, 3/2, 3), A vertex at (5/2, 5/2, 1)) |
| 319 | }}} |