6 const gfan::ZVector &interiorPoint,
const gfan::ZMatrix &normalVectors)
8 int k = normalVectors.getHeight();
9 std::vector<bool> needToFlip(k,
true);
11 int n = normalVectors.getWidth();
12 gfan::ZMatrix testVectors(k,n);
13 gfan::ZVector bigInteriorPoint = 1000*interiorPoint;
14 for (
int i=0;
i<
k;
i++)
15 testVectors[
i] = bigInteriorPoint+normalVectors[
i];
17 for (groebnerCones::iterator sigma = tropicalVariety.begin(); sigma!=tropicalVariety.end(); sigma++)
19 if (sigma->contains(interiorPoint))
21 for (
int i=0;
i<
k;
i++)
23 if (needToFlip[
i] && sigma->contains(testVectors[
i]))
24 needToFlip[i] =
false;
29 for (groebnerCones::iterator sigma = workingList.begin(); sigma!=workingList.end(); sigma++)
31 if (sigma->contains(interiorPoint))
33 for (
int i=0;
i<
k;
i++)
35 if (needToFlip[
i] && sigma->contains(testVectors[
i]))
36 needToFlip[i] =
false;
53 workingList.insert(startingCone);
55 std::set<gfan::ZVector> finishedInteriorPoints;
56 while(!workingList.empty())
64 for (
int i=0;
i<interiorPoints.getHeight();
i++)
69 gfan::ZVector interiorPoint = interiorPoints[
i];
80 for (
int j=0;
j<normalVectors.getHeight();
j++)
85 workingList.insert(neighbour);
89 finishedInteriorPoints.insert(interiorPoint);
93 workingList.erase(sigma);
94 tropicalVariety.insert(sigma);
96 Print(
"cones finished: %lu cones in working list: %lu\n",
97 (
unsigned long)tropicalVariety.size(), (
unsigned long)workingList.size());
104 const gfan::ZMatrix &interiorPoints)
106 int k = interiorPoints.getHeight();
107 std::vector<bool> needToFlip(k,
true);
109 for (groebnerCones::iterator sigma = groebnerFan.begin(); sigma!=groebnerFan.end(); sigma++)
111 for (
int i=0;
i<
k;
i++)
113 if (needToFlip[
i] && sigma->contains(interiorPoints[
i]))
114 needToFlip[i] =
false;
118 for (groebnerCones::iterator sigma = workingList.begin(); sigma!=workingList.end(); sigma++)
120 for (
int i=0;
i<
k;
i++)
122 if (needToFlip[
i] && sigma->contains(interiorPoints[
i]))
123 needToFlip[i] =
false;
137 workingList.insert(startingCone);
138 std::set<gfan::ZVector> finishedInteriorPoints;
141 while(!workingList.empty())
148 workingList.erase(workingList.begin());
151 gfan::ZMatrix interiorPoints = interiorPointsAndOuterFacetNormals.first;
152 gfan::ZMatrix outerFacetNormals = interiorPointsAndOuterFacetNormals.second;
155 for (
int i=0;
i<interiorPoints.getHeight();
i++)
157 gfan::ZVector interiorPoint = interiorPoints[
i];
159 if (needToFlip[
i]==
true)
162 workingList.insert(neighbour);
164 finishedInteriorPoints.insert(interiorPoints[
i]);
168 groebnerFan.insert(sigma);
170 Print(
"cones finished: %lu cones in working list: %lu\n",
171 (
unsigned long)groebnerFan.size(), (
unsigned long)workingList.size());
gfan::ZCone getPolyhedralCone() const
ring getPolynomialRing() const
gfan::ZFan * groebnerFan(const tropicalStrategy currentStrategy)
std::vector< bool > checkNecessaryTropicalFlips(const groebnerCones &tropicalVariety, const groebnerCones &workingList, const gfan::ZVector &interiorPoint, const gfan::ZMatrix &normalVectors)
void id_Delete(ideal *h, ring r)
deletes an ideal/module/matrix
std::pair< gfan::ZMatrix, gfan::ZMatrix > interiorPointsAndNormalsOfFacets(const gfan::ZCone zc, const std::set< gfan::ZVector > &exceptThesePoints, const bool onlyLowerHalfSpace)
groebnerCones groebnerTraversal(const groebnerCone startingCone)
std::set< groebnerCone, groebnerCone_compare > groebnerCones
const tropicalStrategy * getTropicalStrategy() const
void deletePolynomialData()
groebnerCones tropicalTraversalMinimizingFlips(const groebnerCone startingCone)
poly initial(const poly p, const ring r, const gfan::ZVector w)
Returns the initial form of p with respect to w.
bool isValuationTrivial() const
BOOLEAN tropicalVariety(leftv res, leftv args)
bool restrictToLowerHalfSpace() const
returns true, if valuation non-trivial, false otherwise
implementation of the class groebnerCone
std::vector< bool > checkNecessaryGroebnerFlips(const groebnerCones &groebnerFan, const groebnerCones &workingList, const gfan::ZMatrix &interiorPoints)
gfan::ZMatrix raysOfTropicalStar(ideal I, const ring r, const gfan::ZVector &u, const tropicalStrategy *currentStrategy)
groebnerCone flipCone(const gfan::ZVector &interiorPoint, const gfan::ZVector &facetNormal) const
Given an interior point on the facet and the outer normal factor on the facet, returns the adjacent g...
gfan::ZMatrix interiorPointsOfFacets(const gfan::ZCone &zc, const std::set< gfan::ZVector > &exceptThese)
ideal getPolynomialIdeal() const