next | previous | forward | backward | up | top | index | toc | Macaulay2 web site
FourTiTwo :: toricGraverDegrees

toricGraverDegrees -- displays the degrees of all Graver basis elements for the toric ideal I_A; equivalent to "output --degrees foo.gra" in 4ti2

Synopsis

Description

Very often the Graver basis consits of too many binomials, and one is only interested in their degrees. In this case, instead of looking at the Graver basis of IA, we may just want to look for the degrees of binomials which show up:

i1 : A = matrix "1,1,1,1; 1,2,3,4"

o1 = | 1 1 1 1 |
     | 1 2 3 4 |

              2        4
o1 : Matrix ZZ  <--- ZZ
i2 : toricGraver(A) -- the Graver basis

o2 = | 1 -2 1  0 |
     | 2 -3 0  1 |
     | 1 -1 -1 1 |
     | 0 1  -2 1 |
     | 1 0  -3 2 |

              5        4
o2 : Matrix ZZ  <--- ZZ
i3 : toricGraverDegrees(A) -- just the degrees
-------------------------------------------------
4ti2 version 1.6.7
Copyright 1998, 2002, 2006, 2015 4ti2 team.
4ti2 comes with ABSOLUTELY NO WARRANTY.
This is free software, and you are welcome
to redistribute it under certain conditions.
For details, see the file COPYING.
-------------------------------------------------
Norm = 4,   number of elements = 3
Norm = 6,   number of elements = 2

Note that these are all 1-norms of the vectors. Since IA is homogeneous, there are 3 binomials of degree 2 (norm 4) and 2 binomials of degree 3 (norm 6).

Here is a more complicated example, where one may not want to see the Graver basis elements explicitly. The toric ideal IM is the ideal of the rational normal scroll S(3,2,3):

i4 : M = matrix "1,1,1,1,1,1,1,1,1,1,1; 1,1,1,1,0,0,0,0,0,0,0; 0,0,0,0,1,1,1,0,0,0,0; 0,0,0,0,0,0,0,1,1,1,1; 1,2,3,4,1,2,3,1,2,3,4"

o4 = | 1 1 1 1 1 1 1 1 1 1 1 |
     | 1 1 1 1 0 0 0 0 0 0 0 |
     | 0 0 0 0 1 1 1 0 0 0 0 |
     | 0 0 0 0 0 0 0 1 1 1 1 |
     | 1 2 3 4 1 2 3 1 2 3 4 |

              5        11
o4 : Matrix ZZ  <--- ZZ
i5 : toricGraverDegrees(M)
-------------------------------------------------
4ti2 version 1.6.7
Copyright 1998, 2002, 2006, 2015 4ti2 team.
4ti2 comes with ABSOLUTELY NO WARRANTY.
This is free software, and you are welcome
to redistribute it under certain conditions.
For details, see the file COPYING.
-------------------------------------------------
Norm = 4,   number of elements = 37
Norm = 6,   number of elements = 115
Norm = 8,   number of elements = 76
Norm = 10,   number of elements = 15

Here is another example where with many Graver basis elements. The following matrix is a design matrix for a particular statistical model for a 4-node p1 network; see Fienberg-Petrovic-Rinaldo.

i6 : A = matrix "1,0,1,1,0,1,1,0,1,0,0,0,0,0,0,0,0,0;0,1,1,0,0,0,0,0,0,1,0,1,1,0,1,0,0,0;0,0,0,0,1,1,0,0,0,0,1,1,0,0,0,1,0,1;0,0,0,0,0,0,0,1,1,0,0,0,0,1,1,0,1,1;0,1,1,0,1,1,0,1,1,0,0,0,0,0,0,0,0,0;1,0,1,0,0,0,0,0,0,0,1,1,0,1,1,0,0,0;0,0,0,1,0,1,0,0,0,1,0,1,0,0,0,0,1,1;0,0,0,0,0,0,1,0,1,0,0,0,1,0,1,1,0,1;0,0,1,0,0,1,0,0,1,0,0,1,0,0,1,0,0,1";

              9        18
o6 : Matrix ZZ  <--- ZZ
i7 : toricGraverDegrees(A)
-------------------------------------------------
4ti2 version 1.6.7
Copyright 1998, 2002, 2006, 2015 4ti2 team.
4ti2 comes with ABSOLUTELY NO WARRANTY.
This is free software, and you are welcome
to redistribute it under certain conditions.
For details, see the file COPYING.
-------------------------------------------------
Norm = 4,   number of elements = 9
Norm = 6,   number of elements = 67
Norm = 8,   number of elements = 378
Norm = 10,   number of elements = 792
Norm = 12,   number of elements = 1800
Norm = 14,   number of elements = 1872
Norm = 16,   number of elements = 1224
Norm = 18,   number of elements = 792
Norm = 20,   number of elements = 432
Norm = 22,   number of elements = 96

Ways to use toricGraverDegrees :