Polytope volume in Normaliz

We survey the computation of polytope volumes by the algorithms of Normaliz to which the Lawrence algorithm has recently been added. It has enabled us to master volume computations for polytopes from social choice in dimension $119$. This challenge required a sophisticated implementation of the Lawrence algorithm.


INTRODUCTION
About 20 years ago Amelia Taylor asked the author whether Normaliz [10] could compute polytope volumes.It was easy to include this computation goal into the triangulation based "primal algorithm".Since then, polytope volumes have played an important role in the development of Normaliz, and in recent years specific algorithms have been added.
Polytope volumes can be interpreted as degrees of projective toric varieties and multiplicities of monomial algebras.In 2011, Bogdan Ichim pointed out their applications in social choice.Since then they have been a driving challenge for the volume algorithms in Normaliz whose history we sketch briefly.Before going on, let us emphasize that Normaliz computes lattice normalized volumes that for rational polytopes are rational numbers.Euclidean volumes, if asked for by the user, are derived from them.
In social choice, polytope volumes are interpreted as probabilities of certain paradoxa and quality measures of voting schemes.See the books by Gehrlein and Lepelley [16], [17].These applications become rapidly very difficult since the relevant polytopes explode in dimension: for n candidates they are cross-sections of cones of dimension n!.In their paper [19, p. 382] of 2008 Lepelley, Louichi and Smaoui state: Consequently, it is not possible to analyze four candidate elections, where the total number of variables (possible preference rankings) is 24.We hope that further developments of these algorithms will enable the overcoming of this difficulty.With the efficient parallelization of Normaliz in 2012 and the addition of Schürmann's symmetrization method [23], Normaliz could compute a wide variety of social choice polytopes for 4 candidates; see Bruns, Ichim and Söger [9].
When the limitations of the available algorithms became visible in 2017, the author implemented an algorithm for polytope volumes by descent in the face lattice.It is implicitly based on a reverse lexicographic triangulation of the polytope, but does not compute the triangulation explicitly.It brought a significant improvement in computation times for polytopes defined by inequalities, and made more computations for elections with 4 candidates possible, as shown in Bruns and Ichim [6].
But the case n = 5, n! = 120, remained elusive.The breakthrough came with the Normaliz implementation of Lawrence's algorithm [18].It is based on a duality between "generic" triangulations of the dual cone and signed decompositions of the "primal" polytope into simplices.In principle, signed decompositions are as good as ordinary ones for volume computation, but they present hard numerical problems.The rational arithmetic of Normaliz can cope with them, but it must pay by computation time.The applications to 5 candidates elections have been documented by Bruns and Ichim [7].They would have been unreachable without the sophisticated implementation that we explain in Section 4.
Simultaneously with the Lawrence algorithm, we introduced a refinement of the descent algorithm that identifies isomorphic faces in the descent.Isomorphism classes, as explained in [4], are computed by McKay and Piperno's package nauty [21].Even for nauty, isomorphism classes are expensive, but they help in the volume computation of some classical polytopes.
Section 2 gives a very brief introduction to the terminology of this note.It explains basic results that are used in the computation of lattice normalized volume.Section 3 contains an overview of the volume algorithms and explains them, in particular Lawrence's algorithm.The final section 5 lists computation times, with emphasis on the new algorithms, signed decomposition and descent with the exploitation of isomorphism types.They are not only applied to polytopes from social choice, but also to classical polytopes, for example cubes, Birkhoff polytopes and linear ordering polytopes.
The package vinci [12] contains algorithms for polytope volumes.See [6] for a comparative study of computation times and memory usages.Because of its floating point arithmetic, vinci is often faster than Normaliz, but its results come without an error bound.Because of the extreme numerical difficulty, its implementation of the Lawrence algorithm fails reliable results already for polytopes coming from 4 candidates elections; see Remark 6.This note is dedicated to my friends Wolmer Vasconcelos and Rafael Villarreal.Their constant support has been very encouraging in the 25 years of the Normaliz project.One of the first third party publications citing Normaliz is their paper [14] with Delfino, Taylor and Weininger.The example collection of Normaliz still contains input files supplied by Rafael a quarter of a century ago, and his book [24] documents numerous applications.
Acknowledgement.The author was partially supported by the DFG grant BR 688/26-1.He thanks Ulrich von der Ohe for fruitful discussions.

PRELIMINARIES
We refer the reader to [5] for discrete convex geometry.Here we content ourselves to a very brief overview.

Cones and polytopes.
A cone C in the real space R n is the intersection of finitely many linear halfspaces: and for each i the halfspace H + i is the set {x ∈ R n : σ i (x) ≥ 0} for a linear form σ i in the dual space (R n ) * .By the theorem of Minkowski-Weyl, one can equivalently describe cones as the conical set generated by finitely many vectors v j ∈ R n , Since we want to deal only with polytopes and cones derived from them, we can restrict our cones to a subclass: They intersect C in its facets.More generally, a face of C is the intersection of C with a hyperplane that has C inside one of the two closed halfspaces it defines.A face of C is again a cone.
In Sections 3.5 and 4 the dual cone C * will play a central role.Its definition does not only depend on the intrinsic structure of C, but also on the ambient space.Therefore we will then assume that C is a full dimensional pointed cone: dimC = n.Consequently the halfspaces H + i in an irredundant representation (2.1) are uniquely determined, and there is a unique primitive choice for σ i .These linear forms σ i are called the support forms of C. In this note the hyperplanes H i are the support hyperplanes of C. The dual cone is full dimensional and pointed as well.Under the natural identification R n = (R n ) * * the bidual cone C * * is identified with C: the extreme rays of C * are the linear forms defining the facets of C, and vice versa.In the rational case the extreme integral generators of C * are the support forms of C, and vice versa.
A polytope P is the convex hull of finitely many points in a real space R n .Our polytopes will be rational: such polytopes have vertices in Q n .Computationally, polytopes are treated as compact intersections of pointed cones and hyperplanes.The hyperplane is defined by a linear form with integral coefficients, called degree, such that (2. 3) The intersection P is compact (and nonempty) if and only if C = 0 and deg x > 0 for x ∈ C, x = 0.This is not a restriction of generality: if P ⊂ R n is not given as in (2.3), then we can easily re-embed it suitably: we identify P with P ′ = P × {1} ⊂ R n+1 , and choose This definition is independent of the choice of w since Vol is translation invariant on A 0 .Note that the polytope containing a single point x ∈ Q n has lattice volume 1.If desired, the definition of lattice volume can be extended to arbitrary measurable subsets of A, and Normaliz does it for algebraic polytopes.
If P is a lattice polytope, i.e., a polytope with vertices in Z n , then Vol P is an integer.For an arbitrary rational polytope we have Vol P ∈ Q.As a consequence, Vol P can be computed precisely by rational arithmetic.This is not true for Euclidean volume in general: the diagonal of the unit square has length √ 2. A second invariant we need is the lattice height of a rational point x over a rational subspace H = / 0.More generally, one can consider points x such that aff(H, x) is again rational; for example, this is the case if Assume first that 0 ∈ H. Then H is cut out from A by an equation λ (y) = 0 with a primitive Z-linear form λ on L = A ∩ Z n .With this choice of λ , Ht H (x) = |λ (x)| is called the lattice height of x over H. (There are exactly two choices for λ , differing by the factor −1.) If 0 / ∈ H, then we choose an auxiliary point v ∈ H, replace H by H − v, A by A − v and x by x − v.In the algorithms we will only have to deal with the case 0 ∈ H.If P is a rational polytope and F is a facet or, more generally, a face of P, then we set Ht F (x) = Ht H (x) where H = aff F.
The following proposition relates lattice volume and lattice height.
Proposition 1.Let P be a rational polytope and v ∈ P a vertex of P such that there is a single facet F of P with v / ∈ F. Then This is part of [6, Prop.1], to which we refer for the proof.The next basic result tells us how to compute the volume of a simplex, which is a polytope of dimension d with d + 1 vertices.
This follows immediately from the transformation formula of Lebesgue measure.See [5, 2.C] for an algebraic proof.
As mentioned already, we present rational polytopes P in the form P = C ∩ H where C is a pointed cone and H is defined by the condition deg x = 1 with a Z-linear form deg.This brings a second polytope into play, namely P = conv(0, P) as in Figure 1.
This follows immediately from Proposition 2 if we set v 0 = 0, observing that v 1 /g 1 , . . ., v d /g d are the remaining vertices of S.

VOLUME ALGORITHMS IN NORMALIZ
There are three basic algorithms: (1) the primal volume algorithm: Normaliz computes a lexicographic triangulation, and finds the volume as the sum of the volumes of the simplices in the triangulation; (2) volume by descent in the face lattice: there is a reverse lexicographic triangulation in the background, but it is not computed explicitly; (3) volume by signed decomposition, the Lawrence algorithm: Normaliz computes a triangulation of the dual cone and converts it into a signed decomposition of the polytope.Normaliz also computes the exact volume of full dimensional polytopes defined over real algebraic number fields.For them only ( 1) is implemented at present.One could extend (3) to them, whereas ( 2) is not suitable.The algorithms (1) and ( 3) are also used in the computations of integrals of rational polynomials over polytopes.
By rule of thumb one can say that the best choice is (1) if the polytope has few vertices, but potentially many facets; (2) if the number of vertices and the number of facets are of the same order of magnitude; (3) if there are very few facets and many vertices.This recommendation will be confirmed by the computational data in Section 5.There are variants: (a) exploitation of isomorphism types of faces in the descent algorithm; (b) symmetrization as explained below.Normaliz checks the default conditions of the algorithms in the order signed decomposition → descent → symmetrization.If the default conditions are not satisfied for any of them, the primal triangulation algorithm is used.These decisions must often be made on the basis of partial information.Therefore it can be useful to choose a certain variant explicitly or to exclude others.The exploitation of isomorphism types must always be asked for by the user.
Normaliz uses OpenMP for parallelization.Unless the user insists on computations with GMP integers, Normaliz tries 64 bit arithmetic first, and restarts the computation with GMP integers if it recognizes an overflow.
3.1.The primal volume algorithm.Mathematically there is not much to say: if a polytope P is decomposed into simplices with non-overlapping interiors, then its volume is the sum of the volumes of the simplices forming the decomposition.Since the computa-FIGURE 2. A triangulation tion of Hilbert bases and Hilbert series is based on (lexicographic) triangulations as well, Normaliz has a sophisticated algorithm for them, using pyramid decomposition; see [8].Normaliz tries to avoid determinant computations by the "exploitation of unimodularity"; see [8,Prop. 7].

3.2.
Volume by descent in the face lattice.The idea is to exploit the following proposition: Proposition 5. Let P ⊂ R n be a rational polytope, and v ∈ P. Then Proposition 5 follows immediately from Proposition 1 since the polytopes conv(v, F) constitute a polyhedral decomposition of P. Usually v is a vertex of the polytope P with as few opposite facets F i as possible, as illustrated by Figure 3.The recursive application Pyramid decomposition of a polytope results in building a descent system, i.e., a subset F of the face lattice so that for each face F ∈ F , to which (3.1) is applied, all facets of F that are opposite to the selected vertex are contained in F .However, if a face is simplicial, its multiplicity is computed by the standard determinant formula.The algorithm is implemented in such a way that all data are collected in the descent and no backtracking is necessary.The RAM usage is essentially determined by the two largest layers.For a detailed discussion we refer the reader to [6].

Exploitation of isomorphism classes.
If the integral automorphism group of the cone C over the polytope P is large enough, one can expect that each face in the descent system F has many isomorphic copies in F .These can be detected and identified so that only one representative of every isomorphism class must be kept in F .This reduces F in size and can significantly speed up the volume computation.It must be used with care since the computation of isomorphism classes is rather expensive.See [4] for a discussion of their computation.
If the polytope is specified by generators and the number of facets is large, then the first step in the descent system is built differently.Normaliz computes the automorphism group of the polytope and selects one representative in each orbit of facets.The vertex v above is replaced by a fix point of the automorphism group, and the first step in the volume computation is the formula where F i represents one of the c orbits and O(F i ) is the number of facets in the orbit of F i .Then F 1 , . .., F c form the first layer in the descent system.This allows the application of descent in cases where the number of facets is too large for a successful computation without exploitation of isomorphism classes.If the user does not prohibit it, Normaliz encodes isomorphism classes by their SHA256 checksums.
3.4.Symmetrization.To understand the computation of volumes through symmetrization one must take a detour through Ehrhart series.As usual, assume that our polytope P is given as the intersection P = C ∩ H where C ⊂ R d is a pointed rational cone and H = {x ∈ R d : deg x = 1} is the hyperplane of degree 1 points.For symmetrization we assume that deg is primitive.
Under certain conditions one can count lattice points of degree k, k ∈ N, in C by mapping C to a cone C ′ of lower dimension and then counting each degree k lattice point y in C ′ with the number of its lattice preimages.This approach works well if the number of preimages is given by a polynomial in the coordinates of y.Since C ′ has lower dimension, one can hope that its combinatorial structure is much simpler than that of C. One must of course pay a price: instead of counting each lattice point with the weight 1, one must count it with a polynomial weight.
The availability of this approach depends on symmetries in the coordinates of C, and therefore we call it symmetrization.Normaliz tries symmetrization under the following condition: C and the relevant lattice are given by constraints (inequalities, equations, congruences) and the inequalities contain the sign conditions x i ≥ 0 for all coordinates x i of C. Then Normaliz groups coordinates that appear in all constraints and the grading (!) with the same coefficients, and, roughly speaking, replaces them by their sum.The number of preimages that one must count for the vector y of sums is then a product of binomial coefficients-a polynomial as desired.More precisely, if y j , j = 1, . . ., m, is the sum of u j variables x i then is the number of preimages of (y 1 , . . ., y m ).
Since the Lebesgue measure can be approximated by scaled counting measures, one obtains where h is the highest homogeneous component of f with respect to total degree, and λ is the suitably scaled Lebesgue measure.We learnt this approach from Schürmann [23].The Normaliz algorithm for integrals is described in [11].This note contains a complete elementary treatment and several references to advanced aspects.
Symmetrization can have stunning effects.Nevertheless we do not include it in the computations of Section 5 since it does not help for any of them, at least not in the present implementation.Plenty of examples are contained in [9], where it is often very useful in the computation of Hilbert series.

3.5.
Volume by signed decomposition.This algorithm uses that a generic triangulation of the dual cone induces a signed decomposition of the primal polytope, as we will now explain.
Let C ⊂ R d be a pointed cone of dimension d (it is important that C is full dimensional).The polytope P is the intersection of C with the hyperplane H defined by a grading deg: H = {x : deg(x) = 1}.The grading is an interior element of the dual cone C * = {λ ∈ (R d ) * : λ (x) ≥ 0 for all x ∈ C}.In order to visualize the situation we take an auxiliary (irrelevant) cross-section Q of the dual cone as in Figure 4.
A square and a cross-section of the dual cone Now suppose that we have a generic triangulation ∆ of the dual cone where genericity is defined as follows: deg is not contained in any hyperplane through a facet of any δ ∈ ∆.Let δ ∈ ∆ be given, and denote the linear forms on (R d ) * defining its facets by ℓ 1 , . . .ℓ d ∈ (R d ) * * = R d . (ℓ 1 , . . .ℓ d are the extreme rays of the dual of δ .)The hyperplanes defined by the vanishing of ℓ 1 , . . .ℓ d decompose (R d ) * into "orthants" that can be labeled by a sign vector σ = (s 1 , . .., s d ) ∈ {±1} d : By the assumption on deg, there is exactly one sign vector σ such that deg lies in the interior of D(δ , σ ).Consequently the hyperplane H intersects the dual D(δ , σ ) * in a polytope R δ .Set e(δ ) = |{i : s i = −1}|.
Let ι X denote the indicator function of a subset X ⊂ R d .Then for all x ∈ R d outside a union of finitely many hyperplanes.Since volume (lattice normalized or Euclidean) is additive on indicator functions, this formula can be used for the computation Vol P, and more generally for the computation of integrals over P.
In order to find a generic triangulation, Normaliz first computes a triangulation ∆ of C * and saves the induced hollow triangulation Γ that ∆ induces on the boundary of C * .Then it finds a generic element ω ∈ C * such that the star triangulation Σ of C * , in which every simplicial cone is generated by the center ω and a facet of the hollow triangulation, is generic.Figure 5 illustrates the signed decomposition of a square into 4 simplices.
The algorithm has been developed by Lawrence [18] in the language of linear programming.We have learnt it from Filliman's paper [15], which contains a proof of equation . Generic triangulation of the dual and signed decomposition (3.2).See Beck, Haase and Sottile [1] for an elementary approach and the relationship to other duality theorems.

THE IMPLEMENTATION OF THE LAWRENCE ALGORITHM
The complete course of the volume computation consists of 4 steps that can be clearly delimited from each other: (1) computation of a triangulation ∆ of C * ; (2) computation of the induced hollow triangulation Γ; (3) choice of the generic element ω; (4) evaluation of the star triangulation Σ.Each of the 4 steps is highly critical if one wants to reach the applications in social choice that were our driving challenge.For (1) we could essentially rely on the standard triangulation algorithm of Normaliz.Step (2) and the arithmetic for (3) and (4) are described in the following.Both (3) and ( 4) are iterations over the hollow triangulation and star triangulations derived from it.
4.1.The hollow triangulation.Suppose the triangulation ∆ of C * has been computed.For each simplicial cone δ ∈ ∆ we must now find the facets of δ that lie in the boundary of C * .There are various solutions for this task.The first that comes to mind is to compute the facets of C * and match the facets of δ with it.But C * can have an enormous number of facets that one does better not compute since they can easily exhaust RAM.The facets of C * are of course extreme rays of the cone C over P, but for signed decomposition Normaliz only computes them if asked for by the user.A second approach that is much better in terms of RAM is to compute the facets of δ and select those that have all extreme rays of C * on the same side as δ .However, this requires an enormous number of scalar products that in high dimension are expensive.
Instead of these geometric approaches, Normaliz uses only the combinatorics of triangulations of cones.It forms the set of all pairs (F, δ ) where F is a facet of δ , δ ∈ ∆.Both components are encoded as 0-1-vectors that indicate the extreme rays of C * spanning F and δ , respectively.From this set one must discard all pairs (F, δ ) for which there exists a pair (F, δ ′ ), δ ′ ∈ ∆, δ ′ = δ .In principle one could eliminate all F that appear a second time without remembering the "mother" δ , but δ helps in several ways.The first is that one can store the hollow triangulation as a set of pairs (δ , φ (δ )) where φ (δ ) is a second 0-1-vector indicating those extreme rays of δ whose omission yields a facet of the hollow triangulation.
In order not to blow up memory for large ∆, the pairs (δ , φ (δ )) are computed in small portions controlled by "patterns".Each pattern is an increasing sequence (p 1 , . .., p r ) of indices, and it is required that the facet F satisfies the following condition: if q 1 , . .., q d−1 , d = dimC * , are the indices of the extreme rays of F in ascending order, then q i = p i for i = 1, . . ., r.

Piggyback simplices.
After the purely combinatorial computation of the hollow triangulation, arithmetic must be used in steps ( 3) and ( 4) above, namely in finding a generic linear form ω ∈ C * and then in the volume computation.Both tasks are accelerated significantly if one takes advantage of the fact that simplices G and G ′ of the star triangulation are in "piggyback" position to each other, if the facets F and F ′ of the hollow triangulation that define them belong to the same simplex δ ∈ ∆.By "piggyback" position we mean that the simplices share a facet and lie on different sides of it, as indicated in Figure 6.
Since the computation of the k i and −ℓ d from λ 1 , . . ., λ d−1 , µ alone amounts to the inversion of a matrix, it is clear that the use of (4.1) is a significant advantage, even if the computation of the values ℓ i (µ) needs d 2 multiplications.There is actually no need to compute the k i completely.We will only need their values on elements in the dual space, for which the values of the ℓ i are known, for example the degree.
In the primal space, ℓ 1 , . . ., ℓ d and k 1 , . . ., k d−1 , −ℓ d are extreme rays of the simplices dual to G and G ′ .For the volume of the corresponding simplices we need the determinants.By standard rules So the piggyback relation between G and G ′ pays off a second time.
In dealing with the simplices of the star triangulation that belong to the same simplicial cone δ of ∆, we pick one of them, say G 1 , and take all others piggyback.For G 1 we must indeed invert the matrix M with rows λ 1 , . .., λ d over Q, using D = | det M| as the denominator: M −1 = (1/D)N with a matrix N ∈ Z d×d .After extraction of their greatest common divisors, the columns ℓ 1 , . . ., ℓ d of N are the support forms of G 1 , equivalently, the extreme integral generators of the dual cone of G 1 .For the volume computation we need the determinant of the matrix N ′ with columns ℓ 1 , . . ., ℓ d .There is no need to compute it directly: since MN ′ is the diagonal matrix with entries ℓ i (λ i ), i = 1, . .., d, one has and det M has already been computed.

4.3.
Finding a generic element.The generic element ω of C * must satisfy the following condition: the grading deg does not lie on any hyperplane through a facet of a simplicial cone in the hollow triangulation and ω.But this condition is symmetric in ω and deg!It is much better to first take the star triangulation with center deg and check that ω is not on any of the critical hyperplanes.The main difference between deg and ω is the size of the coordinates: those of deg are usually very small and those of ω very large.So, in working with deg as the center, there is a very good chance to get away with 64 bit arithmetic.The computation based on ω which is necessary for the volume, must very often be done with GMP integers.Instead of choosing one vector ω at random and verifying that it is generic, Normaliz takes two vectors ω 1 and ω 2 and checks that not both of them lie together on a critical hyperplane.If this condition is satisfied, then a suitable linear combination ω = a 1 ω 1 + a 2 ω 2 , a 1 , a 2 ∈ Z, a 1 , a 2 ≥ 0, is generic.For the check we apply the piggyback trick of Section 4.2, and especially (4.1).If ω 1 and ω 2 do not work, then their coordinates are increased.

The addition of fractions.
There is one more extremely critical aspect, namely the addition of fractions.If one accumulates the volume as an alternating sum of simplex volumes by successive addition of fractions, one can easily spend 99% of the computation time on this addition, or, in extreme situations, not finish at all, creating fractions whose numerators and denominators fill gigabytes.The reason is that the extreme integral generators of the duals of cones involving the generic element can be very large: we must divide by them.
In extreme cases there is no other choice but to work with fixed precision.If the user asks for it, then the simplex volumes are still computed exactly, but for addition they are truncated to a fixed number of decimal digits.The default choice is truncation to an integral multiple of 10 −100 .With this choice the final volume is computed up to an error ≤ |Γ| • 10 −100 .The user can set a higher or lower precision.
For computations with full precision, Normaliz uses an addition pyramid, an extensible vector (q 0 , . . ., q n ) of rational numbers.It serves as an accumulator.A new summand s is added to q 0 , provided the number of summands that have already been accumulated in q 0 is smaller than the capacity.If the capacity has been reached, then q 0 is added to q 1 and set to s.The addition of q 0 to q 1 is handled in the same way as that of s to q 0 etc.At present the capacity is 8.This scheme has proved to be very efficient.Of course, at the very end, all entries of the pyramid must be added.Remark 6. vinci [12] contains the Lawrence algorithm in a floating point implementation.As the authors state in [13], it is numerically unstable, and our preceding discussion should also indicate this problem.
If one has a look at the volumes whose alternating sum must be formed, then they easily reach absolute values of 10 100 .For any precision of their alternating sum, which may be of order 10 −6 , one therefore needs a very high number of significant digits that a standard floating point format does not offer.
For example, this becomes visible already in the comparison of 4 voting schemes for 4 candidates [6, Sect.6.1].Let P be the corresponding polytope.With its algorithm HOT (based on the same principle as descent) vinci correctly computes the Euclidean volume of 1.260510232743 • 10 −25 , for P, whereas the vinci Lawrence algorithm yields 9.287423132835 • 10 −8 .(We are grateful to Bogdan Ichim for these computations.)For this reason we are not comparing the Lawrence algorithms in Normaliz and vinci.

COMPUTATIONAL DATA
All computations have been done on the Dell R640 server of the Institute of Mathematics at Osnabrück.It is equipped with two Intel Xeon Gold 6152 cards (a total of 44 cores) and 1 TB of RAM.The computations use 32 parallel threads (of the maximum of 88).The listed times are "wall clock" times.In order to avoid overloading the tables, information about RAM usage has been inserted into the text.
That we allow 32 threads for a computation does of course not mean that they can be used.The percentage of CPU that the computations got varies from ∼ 400% to almost 3200%.

5.1.
Polytopes defined by vertices.The computation times for several polytopes defined by vertices are listed in Table 1.In the tables, dim is the dimension of the cone C over the polytope, #ext the number of extreme rays of C and sup the number of its support hyperplanes.The size of the triangulation computed by the primal triangulation is to be found in the column #tri.The number of determinants computed is usually smaller.[20].
The maximum n reachable is 7.For n = 8 not even the number of facets is known.It is however > 800 • 10 6 .The computation of the volumes is surprisingly fast if one exploits the isomorphism classes of faces.Note that for n = 7 the computation of the support hyperplanes alone takes > 20 h so that the computation of the volume needs ∼ 1 h.The maximum RAM usage of lo-6 is 2.1 GB for descent, the other two algorithms need < 1.5 GB.The computation for lo-7 takes 14.7 GB. 5.1.2.Cross polytopes.cr-<n> is the unit cross polytope of dimension n.We have computed their (known) volumes for n = 20, 24, 28.They have only 2n vertices, but 2 n facets.But all facets are simplices, and therefore the descent algorithm and its variant exploiting isomorphism classes are applicable.However, the primal algorithm behaves better for two reasons: (i) it avoids the administrative overhead of the descent algorithm, and (ii) the formation of the single orbit of facets takes rather long-it cannot be parallelized.So the saving in the computations of determinants is overcompensated.Since all facets of a cross polytope are simplicial, there is only one descent step, namely from the full polytope to the facets opposite to the chosen vertex.
For n = 20 the primal algorithm gets away with 713 MB, whereas the two descent algorithms need about 1 GB.For n = 24 the numbers are 4.2 GB and 21 GB.For n = 28 they rise to 194 GB and ∼ 300 GB.It takes a lot of space to accommodate the 2 28 extreme rays.5.1.3.Other polytopes.A543 and A553 are taken from the Ohsugi-Hibi classification [22] of polytopes related to contingency tables.A553 shows that descent with isomorphism types can be favorable if the automorphism group of the polytope is sufficiently large.For A543 this effect is already visible, but still small.This applies to cy-60 as well, the cyclotomic polytope of order 60 defined by Beck and Hos ¸ten [2].That the pure descent algorithm is not suitable for this type of polytope is shown by A543.
The RAM usage of A543 is about 1.6 GB for the primal algorithm and 1 GB for descent with isomorphism classes.For A553 the corresponding numbers are 4.3 GB and 101 GB.For cy-60 they are 316 MB and 1.4 GB.

5.2.
Polytopes defined by inequalities.We now turn to polytopes defined by inequalities and equations.Among them we have chosen Birkhoff polytopes, cubes and polytopes from social choice-as said already, the latter were the driving challenge for our implementation of the Lawrence algorithm.In the tables, #tri dual is the size of the triangulation of the dual cone, and #hollow that of the associated hollow triangulation.5.2.1.Birkhoff polytopes.The Birkhoff polytope of order n is the set of doubly stochastic n × n matrices.Its vertices are the n × n permutation matrices, and their number n! is rapidly growing.Their volumes have been computed for n ≤ 10 by Beck and Pixton [3] with residue methods that are not (yet) available in Normaliz.
For n ≤ 5 any of the Normaliz algorithm does the job very quickly, but for n = 6 the primal algorithm must already give up since the lexicographic triangulation becomes too large.In the table we start with this case.The Lawrence algorithm reaches n = 8.As one can see, even the triangulations of the dual cone grow too quickly for the next step.The bulk of the computation time for n = 8 is taken by the computation of the hollow triangulation, namely ∼ 13 h.The coordinates of the generic element ω are small enough to allow 64 bit arithmetic for the volume computation (< 3 h).RAM usage for bi-6 is 44 MB and 473 MB.The computations for bi-7 take 7.9 GB and 12.5 GB.That for bi-8 needs 216 GB. 5.2.2.Cubes.cu-<n> is the unit cube of dimension n.It is a good test object since its volume is known.Since faces of the same codimension are isomorphic, descent with isomorphism types is expected to be fast, and it is indeed.However, it must use the huge number of vertices explicitly, and for this reason signed decomposition is even faster.The Normaliz binary in the distribution does never reach any of these algorithms since Normaliz recognizes parallelotopes P, computes the volume of a "corner simplex" and multiplies it by n!, n = dim P.This takes ∼ 0.01 s, even for n = 28.It would certainly be possible to go to n = 32 with the Lawrence algorithm or descent with isomorphism types.
In [6] the reader can find performance data for the descent algorithm applied to cu-20 and cu-24.In addition, more general parallelotopes of the same dimensions are computed there.The computation times show that the arithmetic is secondary and the times are dominated by the combinatorial complexity.
RAM usage for cu-20 is ∼ 1.1 GB for both algorithms, and for cu-24 we need 23 GB and 2 GB.cu-28 takes 428 GB and 127 GB.That signed decomposition is so much better, is due to the fact that it does not store the extreme rays.5.2.3.Polytopes from social choice.Computational data for polytopes from social choice are contained in [9] for the primal algorithm and symmetrization and in [6] for descent in the face lattice.Voting schemes with 5 candidates are essentially inaccessible to them.The Lawrence algorithm has now reached them, and [7] contains data for them.Tables 3 4 are imported from there.The names of the polytopes are explained in [7].
For two polytopes we have included the number of extreme rays to show the order of magnitude.For the Lawrence algorithm they are not needed explicitly, and in particular they need not be stored.The stages (1)-( 3) of all computations could be done by 64 bit arithmetic, and this holds even for the volume computations of the first and third polytope.The volume computations of the two largest had to be done with fixed precision.
For the two largest examples it was necessary to use distributed computation on a high performance cluster (indicated by HPC).For this reason we have split the computation times.For CondEffPlur on the HPC of the University of Osnabrück the time was < 9 h.We refer the reader to [7]

FIGURE 6 .
FIGURE 6. Piggyback simplices Suppose that G has extreme rays λ 1 , . . ., λ d ∈ (R d ) * and its facets are given by linear forms ℓ 1 , . . ., ℓ d ∈ R d = (R d ) * * where ℓ d defines the common facet, λ 1 , . . ., λ d−1 are the extreme rays shared by G and G ′ and µ is the remaining extreme ray of G ′ .Then the facets of G ′ are determined by −ℓ d and [6,.Lattice normalized volume.Normaliz computes lattice normalized volume.We review this notion with emphasis of its computation.The reader can find more details in[6, Sect.3].Let P ⊂ R n be a rational polytope.The affine hull A = aff P is a rational affine subspace of R n .First assume that 0 ∈ A. Then L = (aff P) ∩ Z n is a subgroup of Z n of rank d = dim P (and Z n /L is torsionfree).Choose a Z-basis v 1 , . .., v d of L. The lattice (normalized) volume Vol on A is the Lebesgue measure on A scaled in such a way that the simplex conv(0, v 1 , . . ., v d ) has measure 1.The definition is independent of the choice of v 1 , . . ., v d since all invertible d × d matrices over Z have determinant ±1.If 0 / ∈ A, then we replace A by a translate A 0

TABLE 1 .
Polytopes defined by vertices 5.1.1.Linear ordering polytopes.Lo<n> is the linear ordering polytope for a set of n elements.These polytopes have been investigated in combinatorial optimization; see

TABLE 4 .
for more information.Memory usage and times for parallelized volume computations