--- src/halfopencone.cpp.orig 2017-06-20 14:47:37 UTC +++ src/halfopencone.cpp @@ -602,7 +602,7 @@ void tropicalHyperSurfaceIntersectionWithProcessor(int intersectionOld.print(&P); */ } - log2 fprintf(Stderr,"Halfopen intersection completed.\n"); + gfan_log2 fprintf(Stderr,"Halfopen intersection completed.\n"); } PolyhedralFan tropicalHyperSurfaceIntersectionClosed(int dimension, PolynomialSet const &g, PolyhedralCone *restrictingCone, bool expand, bool saveResult, int intervalLow, int intervalHigh) @@ -1296,7 +1296,7 @@ struct RecursionData (public) if(index == fans.size()) { - log2 fprintf(Stderr,"ADDING CONE\n"); + gfan_log2 fprintf(Stderr,"ADDING CONE\n"); //ret.push_back(current); processor.process(current,chosenFans,chosen); numberOfUsefulCalls++; @@ -1339,7 +1339,7 @@ struct RecursionData (public) static int iterationNumber; if(!(iterationNumber++ & 31)) - log2 + gfan_log2 { fprintf(Stderr,"Iteration level:%i, Chosen fan:%i, Number of candidates:%i, Iteration Number:%i, Useful (%i/%i)=%f\n",index,bestIndex,bestNumberOfCandidates,iterationNumber,numberOfUsefulCalls,totalNumberOfCalls,float(numberOfUsefulCalls)/totalNumberOfCalls); fprintf(Stderr,"Chosen fans vector: "); @@ -1490,7 +1490,7 @@ struct RecursionData (public) } bool closure() { - log2 cerr<<"computing closure"<<endl; + gfan_log2 cerr<<"computing closure"<<endl; bool ret=false; int a=0; for(int f1=0;f1<fans.size();f1++) @@ -1532,8 +1532,8 @@ struct RecursionData (public) } } } - log2 fprintf(Stderr,"%i FOR FREE\n",a); - log2 cerr<<"done computing closure"<<endl; + gfan_log2 fprintf(Stderr,"%i FOR FREE\n",a); + gfan_log2 cerr<<"done computing closure"<<endl; return ret; } @@ -1607,8 +1607,8 @@ struct RecursionData (public) freeLPs++; } // table.print(); -// log2 fprintf(Stderr,"LPs solved:%i for relation table\n",data2.table.numberOfSolvedLPs); - log2 fprintf(Stderr,"Number of infeassible LPs discovered from higherdimensional cones: %i\n",freeLPs); +// gfan_log2 fprintf(Stderr,"LPs solved:%i for relation table\n",data2.table.numberOfSolvedLPs); + gfan_log2 fprintf(Stderr,"Number of infeassible LPs discovered from higherdimensional cones: %i\n",freeLPs); } }; @@ -1630,9 +1630,9 @@ void tropicalHyperSurfaceIntersection(int dimension, P HalfOpenConeList l,lf; buildFanFromPolynomial(*i, &lf, &l, &(parents[I]),restrictingCone); - log2 AsciiPrinter(Stderr).printVector(parents[I]); + gfan_log2 AsciiPrinter(Stderr).printVector(parents[I]); - log2 fprintf(Stderr,"\n"); + gfan_log2 fprintf(Stderr,"\n"); vector<HalfOpenCone> L; for(HalfOpenConeList::const_iterator i=l.begin();i!=l.end();i++) @@ -1658,7 +1658,7 @@ void tropicalHyperSurfaceIntersection(int dimension, P // data.table.print();//HERE IntegerVectorList empty; data.rek(0, HalfOpenCone(dimension,empty,empty,empty)); - log2 fprintf(Stderr,"LPs solved:%i for relation table\n",data.table.numberOfSolvedLPs); + gfan_log2 fprintf(Stderr,"LPs solved:%i for relation table\n",data.table.numberOfSolvedLPs); } @@ -1861,7 +1861,7 @@ void tropicalHyperSurfaceIntersectionInSubspace(int di are gotten as the pull-back of the lower dimensional equations by this normal form map. */ - log2 cerr<<"Projecting Newton polytopes modulo the homogeneity space."; + gfan_log2 cerr<<"Projecting Newton polytopes modulo the homogeneity space."; int N=G.getRing().getNumberOfVariables(); IntegerVectorList w=wallInequalities(G); @@ -1887,7 +1887,7 @@ void tropicalHyperSurfaceIntersectionInSubspace(int di restrictingCone=&restrictedConeNew; } - log2 cerr<<"Done projecting Newton polytopes modulo the homogeneity space."; + gfan_log2 cerr<<"Done projecting Newton polytopes modulo the homogeneity space."; /* Now do the computation with the new set of polynomials. */ @@ -1895,7 +1895,7 @@ void tropicalHyperSurfaceIntersectionInSubspace(int di /* Here follows restriction to subspace cut out by the binomials. */ - log2 cerr<<"Restricting to subspace determined by binomials and computing tropical hypersurfaces."; + gfan_log2 cerr<<"Restricting to subspace determined by binomials and computing tropical hypersurfaces."; int n=g.getRing().getNumberOfVariables(); IntegerVectorList equations; @@ -1937,12 +1937,12 @@ void tropicalHyperSurfaceIntersectionInSubspace(int di for(PolynomialSet::const_iterator i=g.begin();i!=g.end();i++) if(i->numberOfTerms()!=2) { - log2 cerr << I; + gfan_log2 cerr << I; HalfOpenConeList l,lf; - log2 cerr<< "Building fan"<<endl; + gfan_log2 cerr<< "Building fan"<<endl; buildFanFromPolynomial(*i, &lf, &l, &(parents[I]),restrictingCone); - log2 cerr<< "Number of cones:" << lf.size()<<","<< l.size()<<endl; - log2 cerr<< "rewriting"<<endl; + gfan_log2 cerr<< "Number of cones:" << lf.size()<<","<< l.size()<<endl; + gfan_log2 cerr<< "rewriting"<<endl; vector<HalfOpenCone> L; for(HalfOpenConeList::const_iterator i=l.begin();i!=l.end();i++) @@ -1953,10 +1953,10 @@ void tropicalHyperSurfaceIntersectionInSubspace(int di fullDimFanList.push_back(F); coDimOneFanList.push_back(L); - log2 cerr<< "Done Building fan"<<endl; + gfan_log2 cerr<< "Done Building fan"<<endl; I++; } - log2 cerr<<"Done computing tropical hypersurfaces."; + gfan_log2 cerr<<"Done computing tropical hypersurfaces."; /* Now we must create a new HalfOpenConeProcessor and tell it how to expand a cone. We must insert the pivot columns from A and add in the equtions gotten from the binomials. @@ -1974,12 +1974,12 @@ void tropicalHyperSurfaceIntersectionInSubspace(int di IntegerVectorList empty; - log2 cerr<<"Doing intersection."; + gfan_log2 cerr<<"Doing intersection."; if(intervalHigh!=-1)data.setInterval(intervalLow,intervalHigh); data.rek(0, HalfOpenCone(nonPivots.size(),empty,empty,empty)); - log2 cerr<<"Done doing intersection."; + gfan_log2 cerr<<"Done doing intersection."; // data.rek(0, HalfOpenCone(dimension,empty,empty,empty)); - // log2 fprintf(Stderr,"LPs solved:%i for relation table\n",data.table.numberOfSolvedLPs); + // gfan_log2 fprintf(Stderr,"LPs solved:%i for relation table\n",data.table.numberOfSolvedLPs); /* HalfOpenConeList ret; @@ -2151,7 +2151,7 @@ bool hasMixedCellOfDimension(PolynomialSet const &g, i IntegerVectorList empty; data.rek(0, HalfOpenCone(g.getRing().getNumberOfVariables(),empty,empty,empty)); - log2 fprintf(Stderr,"LPs solved:%i for relation table\n",data.table.numberOfSolvedLPs); + gfan_log2 fprintf(Stderr,"LPs solved:%i for relation table\n",data.table.numberOfSolvedLPs); if(maximalSeen)*maximalSeen=p.getMaximalDimensionFound(); @@ -2185,7 +2185,7 @@ bool hasMixedCellOfDimension(list<list<IntegerVector> IntegerVectorList empty; data.rek(0, HalfOpenCone(/*g.getRing().getNumberOfVariables()*/d,empty,empty,empty)); - log2 fprintf(Stderr,"LPs solved:%i for relation table\n",data.table.numberOfSolvedLPs); + gfan_log2 fprintf(Stderr,"LPs solved:%i for relation table\n",data.table.numberOfSolvedLPs); if(maximalSeen)*maximalSeen=p.getMaximalDimensionFound();