7 const gfan::ZVector &interiorPoint,
const gfan::ZMatrix &normalVectors)
9 int k = normalVectors.getHeight();
10 std::vector<bool> needToFlip(k,
true);
12 int n = normalVectors.getWidth();
13 gfan::ZMatrix testVectors(k,n);
14 gfan::ZVector bigInteriorPoint = 1000*interiorPoint;
15 for (
int i=0;
i<
k;
i++)
16 testVectors[
i] = bigInteriorPoint+normalVectors[
i];
18 for (groebnerCones::iterator sigma = tropicalVariety.begin(); sigma!=tropicalVariety.end(); sigma++)
20 if (sigma->contains(interiorPoint))
22 for (
int i=0;
i<
k;
i++)
24 if (needToFlip[
i] && sigma->contains(testVectors[
i]))
26 needToFlip[
i] =
false;
33 for (groebnerCones::iterator sigma = workingList.begin(); sigma!=workingList.end(); sigma++)
35 if (sigma->contains(interiorPoint))
37 for (
int i=0;
i<
k;
i++)
39 if (needToFlip[
i] && sigma->contains(testVectors[
i]))
41 needToFlip[
i] =
false;
60 workingList.insert(startingCone);
62 std::set<gfan::ZVector> finishedInteriorPoints;
63 while(!workingList.empty())
71 for (
int i=0;
i<interiorPoints.getHeight();
i++)
76 gfan::ZVector interiorPoint = interiorPoints[
i];
89 for (
int j=0;
j<normalVectors.getHeight();
j++)
94 workingList.insert(neighbour);
98 finishedInteriorPoints.insert(interiorPoint);
102 workingList.erase(sigma);
103 tropicalVariety.insert(sigma);
105 Print(
"cones finished: %lu cones in working list: %lu\n",
106 (
unsigned long)tropicalVariety.size(), (
unsigned long)workingList.size());
113 const gfan::ZMatrix &interiorPoints)
115 int k = interiorPoints.getHeight();
116 std::vector<bool> needToFlip(k,
true);
118 for (groebnerCones::iterator sigma = groebnerFan.begin(); sigma!=groebnerFan.end(); sigma++)
120 for (
int i=0;
i<
k;
i++)
122 if (needToFlip[
i] && sigma->contains(interiorPoints[
i]))
123 needToFlip[i] =
false;
127 for (groebnerCones::iterator sigma = workingList.begin(); sigma!=workingList.end(); sigma++)
129 for (
int i=0;
i<
k;
i++)
131 if (needToFlip[
i] && sigma->contains(interiorPoints[
i]))
132 needToFlip[i] =
false;
146 workingList.insert(startingCone);
147 std::set<gfan::ZVector> finishedInteriorPoints;
150 while(!workingList.empty())
157 workingList.erase(workingList.begin());
160 gfan::ZMatrix interiorPoints = interiorPointsAndOuterFacetNormals.first;
161 gfan::ZMatrix outerFacetNormals = interiorPointsAndOuterFacetNormals.second;
164 for (
int i=0;
i<interiorPoints.getHeight();
i++)
166 gfan::ZVector interiorPoint = interiorPoints[
i];
168 if (needToFlip[
i]==
true)
171 workingList.insert(neighbour);
173 finishedInteriorPoints.insert(interiorPoints[
i]);
177 groebnerFan.insert(sigma);
179 Print(
"cones finished: %lu cones in working list: %lu\n",
180 (
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)
ideal gfanlib_satStd_wrapper(ideal I, ring r, tHomog h=testHomog)
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