1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
|
#include <stdlib.h>
#include <vector>
#include <fstream>
#include <string>
#ifdef _OPENMP
#include <omp.h>
#endif
using namespace std;
#include "libnormaliz/libnormaliz.h"
#include "libnormaliz/cone.h"
#include "libnormaliz/vector_operations.h"
#include "libnormaliz/cone_property.h"
#include "libnormaliz/integer.h"
#include "libnormaliz/matrix.h"
using namespace libnormaliz;
typedef long long Integer;
int main(int argc, char* argv[]){
Matrix<Integer> Gens=readMatrix<Integer>(string("small_gens.mat"));
/*Cone<Integer> C(Type::cone,Gens);
C. setVerbose(true);
C.compute(ConeProperty::DefaultMode);
C.getVerticesFloatMatrix().pretty_print(cout);
C.setExpansionDegree(25);
cout << C.getHilbertSeries().getExpansion();
cout << "========================" << endl;
C.setExpansionDegree(50);
cout << C.getHilbertSeries().getExpansion();
cout << "========================" << endl;
C.compute(ConeProperty::HSOP);
cout << "HSOP " << C.getHilbertSeries().getHSOPNum();
C.setNrCoeffQuasiPol(2);
cout << "========================" << endl;
Matrix<mpz_class> Q(C.getHilbertSeries().getHilbertQuasiPolynomial());
Q.pretty_print(cout);
cout << "========================" << endl;
C.setNrCoeffQuasiPol(5);
Q=C.getHilbertSeries().getHilbertQuasiPolynomial();
Q.pretty_print(cout);*/
vector<Cone<Integer> > ParCones(16);
#pragma omp parallel for
for(size_t i=0;i<ParCones.size();++i){
ParCones[i]=Cone<Integer>(Type::cone,Gens);
// ParCones[i].setVerbose(true);
switch(i%8){
case 0: ParCones[i].compute(ConeProperty::DefaultMode);
break;
case 1: ParCones[i].compute(ConeProperty::DualMode, ConeProperty::Deg1Elements);
break;
case 2: ParCones[i].compute(ConeProperty::Projection, ConeProperty::Deg1Elements);
break;
case 3: ParCones[i].compute(ConeProperty::ProjectionFloat, ConeProperty::Deg1Elements);
break;
case 4: ParCones[i].compute(ConeProperty::Approximate, ConeProperty::Deg1Elements, ConeProperty::IsGorenstein);
break;
case 5: ParCones[i].compute(ConeProperty::SupportHyperplanes);
break;
case 6: ParCones[i].compute(ConeProperty::IntegerHull);
break;
case 7: ParCones[i].compute(ConeProperty::IsIntegrallyClosed);
break;
default: break;
}
}
} //end main
|
