playground:extensions:max_mediated_sets

This is an old revision of the document!


Maximal Mediated Sets

max_mediated_sets is a polymake extension allowing creation and inspection of the maximal mediated set associated to a polynomial. It accompanies this paper. The parallel enumeration algorithm described there will be made available soon, too.

You'll need a running polymake installation with the bundled extension libnormaliz for computing lattice points. libnormaliz comes bundled with polymake, so if you did not do anything special like deliberately installing the minimal version, you should already have libnormaliz and do not have to do any extra work. See here for instructions in how to install polymake.

Use git to obtain your copy:

git clone https://git.polymake.org/extensions/max_mediated_sets

Assuming the directory you pulled the git to has the path ~/max_mediated_sets. To use the extension, fire up your polymake interactive shell and enter the following command:

>  import_extension "~/max_mediated_sets";

After this, you should be able to use the functionality provided by this extension.

This extension allows to compute maximal mediated sets from PointConfigurations, that is, all lattice points that are the midpoint of two all-even points.

Newton Configurations

This is interesting for Newton polytopes of polynomials, so there is a function producing the point configuration whose convex hull is the Newton polytope of a polynomial p. For example, the following creates said point configuration for the polynomial $1+x_0^2+x_1^6$:

> $p=new Polynomial('1+x_0^2+x_1^6');
> $n = newton_configuration($p);
> print $n->CONVEX_HULL->VERTICES;
1 0 0
1 0 1
1 2 0

Note that the vertices are handled in homogeneous coordinates. The polynomial gets stored in the attachment named POLYNOMIAL:

> print $n->get_attachment("POLYNOMIAL");
x_0^2 + x_1^6 + 1

The coefficients of the polynomial get stored in the property TRIANGULATION.WEIGHTS.

> print $n->TRIANGULATION->WEIGHTS;
1 1 1

Maximal Mediated Sets

You can compute the maximal mediated set of any point configuration whose convex hull has all-even vertices, like the one we created in the previous section.

> print $n->MAXIMAL_MEDIATED_SET;
1 0 0
1 0 2
1 0 4
1 0 6
1 2 0
1 1 3
1 1 2
1 0 3
1 1 0
1 0 1
1 1 1
1 0 5

You can look at the vanishing ideal, which in this case has two generators:

> print $n->VANISHING_IDEAL_GENERATORS->[1];
x_1^3 - x_3
> print $n->VANISHING_IDEAL_GENERATORS->[2];
x_0 - 1

You might also be interested in the dimension of the toric variety:

> print $n->TORIC_VARIETY_DIMENSION;
9
  • playground/extensions/max_mediated_sets.1561715303.txt.gz
  • Last modified: 2019/06/28 09:48
  • by oroehrig