Skip to content

Commit 58079dd

Browse files
authored
Remove _NmzPrintSomeConeProperties (#119)
1 parent f416e45 commit 58079dd

7 files changed

Lines changed: 320 additions & 266 deletions

File tree

tst/descent.tst

Lines changed: 30 additions & 25 deletions
Original file line numberDiff line numberDiff line change
@@ -137,31 +137,36 @@ gap> Display(NmzSublattice(cone));
137137
[ [ 1, 0, 0, 0, 0, 0, 0 ], [ 0, 1, 0, 0, 0, 0, 0 ], [ 0, 0, 1, 0, 0, 0, 0 ],
138138
[ 0, 0, 0, 1, 0, 0, 0 ], [ 0, 0, 0, 0, 1, 0, 0 ],
139139
[ 0, 0, 0, 0, 0, 1, 0 ], [ 0, 0, 0, 0, 0, 0, 1 ] ], 1 ]
140-
gap> _NmzPrintSomeConeProperties(cone, [
141-
> "Generators",
142-
> "ExtremeRays",
143-
> "SupportHyperplanes",
144-
> "HilbertBasis",
145-
> "Deg1Elements",
146-
> "Sublattice",
147-
> "NumberLatticePoints",
148-
> "OriginalMonoidGenerators",
149-
> ]);
150-
BasicTriangulation = fail
151-
EmbeddingDim = 7
152-
Grading = [ 1, 1, 1, 1, 1, 1, -2 ]
153-
GradingDenom = 1
154-
IsDeg1ExtremeRays = true
155-
IsDeg1HilbertBasis = false
156-
IsInhomogeneous = false
157-
IsPointed = true
158-
IsTriangulationNested = false
159-
IsTriangulationPartial = true
160-
MaximalSubspace = [ ]
161-
Multiplicity = 72
162-
Rank = 7
163-
TriangulationDetSum = 2
164-
TriangulationSize = 1
140+
gap> NmzConeProperty(cone, "BasicTriangulation");
141+
fail
142+
gap> NmzConeProperty(cone, "EmbeddingDim");
143+
7
144+
gap> NmzConeProperty(cone, "Grading");
145+
[ 1, 1, 1, 1, 1, 1, -2 ]
146+
gap> NmzConeProperty(cone, "GradingDenom");
147+
1
148+
gap> NmzConeProperty(cone, "IsDeg1ExtremeRays");
149+
true
150+
gap> NmzConeProperty(cone, "IsDeg1HilbertBasis");
151+
false
152+
gap> NmzConeProperty(cone, "IsInhomogeneous");
153+
false
154+
gap> NmzConeProperty(cone, "IsPointed");
155+
true
156+
gap> NmzConeProperty(cone, "IsTriangulationNested");
157+
false
158+
gap> NmzConeProperty(cone, "IsTriangulationPartial");
159+
true
160+
gap> NmzConeProperty(cone, "MaximalSubspace");
161+
[ ]
162+
gap> NmzConeProperty(cone, "Multiplicity");
163+
72
164+
gap> NmzConeProperty(cone, "Rank");
165+
7
166+
gap> NmzConeProperty(cone, "TriangulationDetSum");
167+
2
168+
gap> NmzConeProperty(cone, "TriangulationSize");
169+
1
165170

166171
#
167172
gap> NmzVolume(cone);

tst/dual.tst

Lines changed: 38 additions & 30 deletions
Original file line numberDiff line numberDiff line change
@@ -496,38 +496,46 @@ gap> Display(NmzSublattice(cone));
496496
[ [ 1, 0, 0, 0, 0, 0, 0 ], [ 0, 1, 0, 0, 0, 0, 0 ], [ 0, 0, 1, 0, 0, 0, 0 ],
497497
[ 0, 0, 0, 1, 0, 0, 0 ], [ 0, 0, 0, 0, 1, 0, 0 ],
498498
[ 0, 0, 0, 0, 0, 1, 0 ], [ 0, 0, 0, 0, 0, 0, 1 ] ], 1 ]
499-
gap> _NmzPrintSomeConeProperties(cone, [
500-
> "Generators",
501-
> "ExtremeRays",
502-
> "SupportHyperplanes",
503-
> "HilbertBasis",
504-
> "Deg1Elements",
505-
> "Sublattice",
506-
> "NumberLatticePoints",
507-
> "OriginalMonoidGenerators",
508-
> ]);
509-
BasicTriangulation = fail
510-
ClassGroup = [ 17 ]
511-
EhrhartQuasiPolynomial = [ [ 60, 194, 284, 245, 130, 41, 6 ], 60 ]
512-
EmbeddingDim = 7
513-
Grading = [ 1, 1, 1, 1, 1, 1, -2 ]
514-
GradingDenom = 1
515-
HilbertQuasiPolynomial =
499+
gap> NmzConeProperty(cone, "BasicTriangulation");
500+
fail
501+
gap> NmzConeProperty(cone, "ClassGroup");
502+
[ 17 ]
503+
gap> NmzConeProperty(cone, "EhrhartQuasiPolynomial");
504+
[ [ 60, 194, 284, 245, 130, 41, 6 ], 60 ]
505+
gap> NmzConeProperty(cone, "EmbeddingDim");
506+
7
507+
gap> NmzConeProperty(cone, "Grading");
508+
[ 1, 1, 1, 1, 1, 1, -2 ]
509+
gap> NmzConeProperty(cone, "GradingDenom");
510+
1
511+
gap> NmzConeProperty(cone, "HilbertQuasiPolynomial");
516512
[ 1/10*t^6+41/60*t^5+13/6*t^4+49/12*t^3+71/15*t^2+97/30*t+1 ]
517-
HilbertQuasiPolynomial =
513+
gap> NmzConeProperty(cone, "HilbertQuasiPolynomial");
518514
[ 1/10*t^6+41/60*t^5+13/6*t^4+49/12*t^3+71/15*t^2+97/30*t+1 ]
519-
HilbertSeries = [ 6*t^4+25*t^3+31*t^2+9*t+1, [ [ 1, 7 ] ] ]
520-
IsDeg1ExtremeRays = true
521-
IsDeg1HilbertBasis = false
522-
IsInhomogeneous = false
523-
IsPointed = true
524-
IsTriangulationNested = false
525-
IsTriangulationPartial = false
526-
MaximalSubspace = [ ]
527-
Multiplicity = 72
528-
Rank = 7
529-
TriangulationDetSum = 72
530-
TriangulationSize = 69
515+
gap> NmzConeProperty(cone, "HilbertSeries");
516+
[ 6*t^4+25*t^3+31*t^2+9*t+1, [ [ 1, 7 ] ] ]
517+
gap> NmzConeProperty(cone, "IsDeg1ExtremeRays");
518+
true
519+
gap> NmzConeProperty(cone, "IsDeg1HilbertBasis");
520+
false
521+
gap> NmzConeProperty(cone, "IsInhomogeneous");
522+
false
523+
gap> NmzConeProperty(cone, "IsPointed");
524+
true
525+
gap> NmzConeProperty(cone, "IsTriangulationNested");
526+
false
527+
gap> NmzConeProperty(cone, "IsTriangulationPartial");
528+
false
529+
gap> NmzConeProperty(cone, "MaximalSubspace");
530+
[ ]
531+
gap> NmzConeProperty(cone, "Multiplicity");
532+
72
533+
gap> NmzConeProperty(cone, "Rank");
534+
7
535+
gap> NmzConeProperty(cone, "TriangulationDetSum");
536+
72
537+
gap> NmzConeProperty(cone, "TriangulationSize");
538+
69
531539

532540
#
533541
gap> STOP_TEST("dual.tst", 0);

tst/project.tst

Lines changed: 62 additions & 67 deletions
Original file line numberDiff line numberDiff line change
@@ -12,17 +12,6 @@ gap> proj:=NmzProjectCone(cone);
1212
# check what was computed for the input cone
1313
gap> tmp := NmzKnownConeProperties(cone);;
1414
gap> RemoveSet(tmp, "NumberLatticePoints");
15-
gap> Perform(tmp, Display);
16-
EmbeddingDim
17-
ExtremeRays
18-
Generators
19-
IsInhomogeneous
20-
IsPointed
21-
MaximalSubspace
22-
ProjectCone
23-
Rank
24-
Sublattice
25-
SupportHyperplanes
2615
gap> Display(NmzTriangulation(cone));
2716
[ [ rec(
2817
Excluded := [ ],
@@ -41,33 +30,38 @@ gap> Display(NmzDeg1Elements(cone));
4130
gap> Display(NmzSublattice(cone));
4231
[ [ [ 1, 0, 0, 0 ], [ 0, 1, 0, 0 ] ],
4332
[ [ 1, 0 ], [ 0, 1 ], [ 0, 0 ], [ 0, 0 ] ], 1 ]
44-
gap> _NmzPrintSomeConeProperties(cone, [
45-
> "Generators",
46-
> "ExtremeRays",
47-
> "SupportHyperplanes",
48-
> "HilbertBasis",
49-
> "Deg1Elements",
50-
> "Sublattice",
51-
> "NumberLatticePoints",
52-
> "OriginalMonoidGenerators",
53-
> ]);
54-
BasicTriangulation = fail
55-
EmbeddingDim = 4
56-
Grading = [ 0, 1, 0, 0 ]
57-
GradingDenom = 1
58-
IsDeg1ExtremeRays = true
59-
IsDeg1HilbertBasis = true
60-
IsInhomogeneous = false
61-
IsPointed = false
62-
IsTriangulationNested = false
63-
IsTriangulationPartial = true
64-
MaximalSubspace =
65-
[ [ 1, 0, 0, 0 ] ]
66-
Multiplicity = 1
67-
ProjectCone = <object>
68-
Rank = 2
69-
TriangulationDetSum = 0
70-
TriangulationSize = 0
33+
gap> NmzConeProperty(cone, "BasicTriangulation");
34+
fail
35+
gap> NmzConeProperty(cone, "EmbeddingDim");
36+
4
37+
gap> NmzConeProperty(cone, "Grading");
38+
[ 0, 1, 0, 0 ]
39+
gap> NmzConeProperty(cone, "GradingDenom");
40+
1
41+
gap> NmzConeProperty(cone, "IsDeg1ExtremeRays");
42+
true
43+
gap> NmzConeProperty(cone, "IsDeg1HilbertBasis");
44+
true
45+
gap> NmzConeProperty(cone, "IsInhomogeneous");
46+
false
47+
gap> NmzConeProperty(cone, "IsPointed");
48+
false
49+
gap> NmzConeProperty(cone, "IsTriangulationNested");
50+
false
51+
gap> NmzConeProperty(cone, "IsTriangulationPartial");
52+
true
53+
gap> NmzConeProperty(cone, "MaximalSubspace");
54+
[ [ 1, 0, 0, 0 ] ]
55+
gap> NmzConeProperty(cone, "Multiplicity");
56+
1
57+
gap> NmzConeProperty(cone, "ProjectCone");
58+
<a Normaliz cone>
59+
gap> NmzConeProperty(cone, "Rank");
60+
2
61+
gap> NmzConeProperty(cone, "TriangulationDetSum");
62+
0
63+
gap> NmzConeProperty(cone, "TriangulationSize");
64+
0
7165

7266
# check what was computed for the projected cone
7367
gap> Perform(NmzKnownConeProperties(proj), Display);
@@ -102,35 +96,36 @@ gap> Display(NmzDeg1Elements(proj));
10296
[ [ 0, 1, 0 ] ]
10397
gap> Display(NmzSublattice(proj));
10498
[ [ [ 1, 0, 0 ], [ 0, 1, 0 ] ], [ [ 1, 0 ], [ 0, 1 ], [ 0, 0 ] ], 1 ]
105-
gap> _NmzPrintSomeConeProperties(proj, [
106-
> "Generators",
107-
> "ExtremeRays",
108-
> "SupportHyperplanes",
109-
> "HilbertBasis",
110-
> "Deg1Elements",
111-
> "Sublattice",
112-
> "NumberLatticePoints",
113-
> "OriginalMonoidGenerators",
114-
> ]);
115-
BasicTriangulation = fail
116-
EmbeddingDim = 3
117-
Grading = [ 0, 1, 0 ]
118-
GradingDenom = 1
119-
InternalIndex = 1
120-
IsDeg1ExtremeRays = true
121-
IsDeg1HilbertBasis = true
122-
IsInhomogeneous = false
123-
IsIntegrallyClosed = true
124-
IsPointed = false
125-
IsTriangulationNested = false
126-
IsTriangulationPartial = true
127-
MaximalSubspace =
128-
[ [ 1, 0, 0 ] ]
129-
Multiplicity = 1
130-
Rank = 2
131-
TriangulationDetSum = 0
132-
TriangulationSize = 0
133-
UnitGroupIndex = 1
99+
gap> NmzConeProperty(cone, "BasicTriangulation");
100+
fail
101+
gap> NmzConeProperty(cone, "EmbeddingDim");
102+
4
103+
gap> NmzConeProperty(cone, "Grading");
104+
[ 0, 1, 0, 0 ]
105+
gap> NmzConeProperty(cone, "GradingDenom");
106+
1
107+
gap> NmzConeProperty(cone, "IsDeg1ExtremeRays");
108+
true
109+
gap> NmzConeProperty(cone, "IsDeg1HilbertBasis");
110+
true
111+
gap> NmzConeProperty(cone, "IsInhomogeneous");
112+
false
113+
gap> NmzConeProperty(cone, "IsPointed");
114+
false
115+
gap> NmzConeProperty(cone, "IsTriangulationNested");
116+
false
117+
gap> NmzConeProperty(cone, "IsTriangulationPartial");
118+
true
119+
gap> NmzConeProperty(cone, "MaximalSubspace");
120+
[ [ 1, 0, 0, 0 ] ]
121+
gap> NmzConeProperty(cone, "Multiplicity");
122+
1
123+
gap> NmzConeProperty(cone, "Rank");
124+
2
125+
gap> NmzConeProperty(cone, "TriangulationDetSum");
126+
0
127+
gap> NmzConeProperty(cone, "TriangulationSize");
128+
0
134129

135130
#
136131
gap> STOP_TEST("project.tst", 0);

tst/rational.tst

Lines changed: 54 additions & 44 deletions
Original file line numberDiff line numberDiff line change
@@ -73,51 +73,61 @@ gap> Display(NmzOriginalMonoidGenerators(cone));
7373
[ [ 1, 1, 2 ],
7474
[ -1, -1, 3 ],
7575
[ 1, -2, 4 ] ]
76-
gap> _NmzPrintSomeConeProperties(cone, [
77-
> "Generators",
78-
> "ExtremeRays",
79-
> "SupportHyperplanes",
80-
> "HilbertBasis",
81-
> "Deg1Elements",
82-
> "Sublattice",
83-
> "NumberLatticePoints",
84-
> "OriginalMonoidGenerators",
85-
> ]);
86-
BasicTriangulation = fail
87-
ClassGroup = [ 0, 3, 15 ]
88-
EhrhartQuasiPolynomial = [ [ 48, 28, 15 ], [ 11, 22, 15 ], [ -20, 28, 15 ],
89-
[ 39, 22, 15 ], [ 32, 28, 15 ], [ -5, 22, 15 ], [ 12, 28, 15 ],
90-
[ 23, 22, 15 ], [ 16, 28, 15 ], [ 27, 22, 15 ], [ -4, 28, 15 ],
91-
[ 7, 22, 15 ], 48 ]
92-
EmbeddingDim = 3
93-
Grading = [ 0, 0, 1 ]
94-
GradingDenom = 1
95-
HilbertQuasiPolynomial = [ 5/16*t^2+7/12*t+1, 5/16*t^2+11/24*t+11/48,
96-
5/16*t^2+7/12*t-5/12, 5/16*t^2+11/24*t+13/16, 5/16*t^2+7/12*t+2/3,
97-
5/16*t^2+11/24*t-5/48, 5/16*t^2+7/12*t+1/4, 5/16*t^2+11/24*t+23/48,
98-
5/16*t^2+7/12*t+1/3, 5/16*t^2+11/24*t+9/16, 5/16*t^2+7/12*t-1/12,
99-
5/16*t^2+11/24*t+7/48 ]
100-
HilbertQuasiPolynomial = [ 5/16*t^2+7/12*t+1, 5/16*t^2+11/24*t+11/48,
101-
5/16*t^2+7/12*t-5/12, 5/16*t^2+11/24*t+13/16, 5/16*t^2+7/12*t+2/3,
102-
5/16*t^2+11/24*t-5/48, 5/16*t^2+7/12*t+1/4, 5/16*t^2+11/24*t+23/48,
103-
5/16*t^2+7/12*t+1/3, 5/16*t^2+11/24*t+9/16, 5/16*t^2+7/12*t-1/12,
104-
5/16*t^2+11/24*t+7/48 ]
105-
HilbertSeries = [ 2*t^12+t^11+t^10+t^9+t^8+2*t^7+2*t^6-t^5+2*t^4+3*t^3+1,
76+
gap> NmzConeProperty(cone, "BasicTriangulation");
77+
fail
78+
gap> NmzConeProperty(cone, "ClassGroup");
79+
[ 0, 3, 15 ]
80+
gap> NmzConeProperty(cone, "EhrhartQuasiPolynomial");
81+
[ [ 48, 28, 15 ], [ 11, 22, 15 ], [ -20, 28, 15 ], [ 39, 22, 15 ],
82+
[ 32, 28, 15 ], [ -5, 22, 15 ], [ 12, 28, 15 ], [ 23, 22, 15 ],
83+
[ 16, 28, 15 ], [ 27, 22, 15 ], [ -4, 28, 15 ], [ 7, 22, 15 ], 48 ]
84+
gap> NmzConeProperty(cone, "EmbeddingDim");
85+
3
86+
gap> NmzConeProperty(cone, "Grading");
87+
[ 0, 0, 1 ]
88+
gap> NmzConeProperty(cone, "GradingDenom");
89+
1
90+
gap> NmzConeProperty(cone, "HilbertQuasiPolynomial");
91+
[ 5/16*t^2+7/12*t+1, 5/16*t^2+11/24*t+11/48, 5/16*t^2+7/12*t-5/12,
92+
5/16*t^2+11/24*t+13/16, 5/16*t^2+7/12*t+2/3, 5/16*t^2+11/24*t-5/48,
93+
5/16*t^2+7/12*t+1/4, 5/16*t^2+11/24*t+23/48, 5/16*t^2+7/12*t+1/3,
94+
5/16*t^2+11/24*t+9/16, 5/16*t^2+7/12*t-1/12, 5/16*t^2+11/24*t+7/48 ]
95+
gap> NmzConeProperty(cone, "HilbertQuasiPolynomial");
96+
[ 5/16*t^2+7/12*t+1, 5/16*t^2+11/24*t+11/48, 5/16*t^2+7/12*t-5/12,
97+
5/16*t^2+11/24*t+13/16, 5/16*t^2+7/12*t+2/3, 5/16*t^2+11/24*t-5/48,
98+
5/16*t^2+7/12*t+1/4, 5/16*t^2+11/24*t+23/48, 5/16*t^2+7/12*t+1/3,
99+
5/16*t^2+11/24*t+9/16, 5/16*t^2+7/12*t-1/12, 5/16*t^2+11/24*t+7/48 ]
100+
gap> NmzConeProperty(cone, "HilbertSeries");
101+
[ 2*t^12+t^11+t^10+t^9+t^8+2*t^7+2*t^6-t^5+2*t^4+3*t^3+1,
106102
[ [ 1, 1 ], [ 2, 1 ], [ 12, 1 ] ] ]
107-
InternalIndex = 15
108-
IsDeg1ExtremeRays = false
109-
IsDeg1HilbertBasis = false
110-
IsInhomogeneous = false
111-
IsIntegrallyClosed = false
112-
IsPointed = true
113-
IsTriangulationNested = false
114-
IsTriangulationPartial = false
115-
MaximalSubspace = [ ]
116-
Multiplicity = 5/8
117-
Rank = 3
118-
TriangulationDetSum = 15
119-
TriangulationSize = 1
120-
UnitGroupIndex = 1
103+
gap> NmzConeProperty(cone, "InternalIndex");
104+
15
105+
gap> NmzConeProperty(cone, "IsDeg1ExtremeRays");
106+
false
107+
gap> NmzConeProperty(cone, "IsDeg1HilbertBasis");
108+
false
109+
gap> NmzConeProperty(cone, "IsInhomogeneous");
110+
false
111+
gap> NmzConeProperty(cone, "IsIntegrallyClosed");
112+
false
113+
gap> NmzConeProperty(cone, "IsPointed");
114+
true
115+
gap> NmzConeProperty(cone, "IsTriangulationNested");
116+
false
117+
gap> NmzConeProperty(cone, "IsTriangulationPartial");
118+
false
119+
gap> NmzConeProperty(cone, "MaximalSubspace");
120+
[ ]
121+
gap> NmzConeProperty(cone, "Multiplicity");
122+
5/8
123+
gap> NmzConeProperty(cone, "Rank");
124+
3
125+
gap> NmzConeProperty(cone, "TriangulationDetSum");
126+
15
127+
gap> NmzConeProperty(cone, "TriangulationSize");
128+
1
129+
gap> NmzConeProperty(cone, "UnitGroupIndex");
130+
1
121131
gap> Display(NmzConeDecomposition(cone));
122132
[ [ rec(
123133
Excluded := [ false, false, false ],

0 commit comments

Comments
 (0)