Changeset 19addd1 in git
- Timestamp:
- Oct 25, 2013, 6:22:01 PM (9 years ago)
- Branches:
- (u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', 'f875bbaccd0831e36aaed09ff6adeb3eb45aeb94')
- Children:
- 69b2c1754cbeecd7c1b844e71e2101e1b3794062
- Parents:
- b7b2c516034df879b01c7c8666bbdf0928da91dc
- git-author:
- Yue Ren <ren@mathematik.uni-kl.de>2013-10-25 18:22:01+02:00
- git-committer:
- Yue Ren <ren@mathematik.uni-kl.de>2013-10-25 18:25:28+02:00
- Location:
- gfanlib
- Files:
-
- 10 edited
Legend:
- Unmodified
- Added
- Removed
-
gfanlib/gfanlib_matrix.h
rb7b2c5 r19addd1 69 69 void appendRow(Vector<typ> const &v) 70 70 { 71 assert( v.size()==width);71 assert((int)v.size()==width); 72 72 rows.push_back(v); 73 73 height++; … … 410 410 { 411 411 assert(typ::isField()); 412 assert( v.size()==getWidth());412 assert((int)v.size()==getWidth()); 413 413 414 414 int pivotI=-1; -
gfanlib/gfanlib_polyhedralfan.cpp
rb7b2c5 r19addd1 352 352 353 353 bool notAll=false; 354 for( intj=0;j<theCone.indices.size();j++)354 for(unsigned j=0;j<theCone.indices.size();j++) 355 355 if(dot(rays[theCone.indices[j]],facetCandidates[i]).sign()==0) 356 356 indices.insert(theCone.indices[j]); … … 395 395 void addFacesToSymmetricComplex(SymmetricComplex &c, ZCone const &cone, ZMatrix const &facetCandidates, ZMatrix const &generatorsOfLinealitySpace) 396 396 { 397 ZMatrix const &rays=c.getVertices();397 // ZMatrix const &rays=c.getVertices(); 398 398 std::set<int> indices; 399 399 … … 502 502 } 503 503 504 std::string PolyhedralFan::toString(int flags)const504 std::string PolyhedralFan::toString(int /*flags*/)const 505 505 //void PolyhedralFan::printWithIndices(class Printer *p, bool printMultiplicities, SymmetryGroup *sym, bool group, bool ignoreCones, bool xml, bool tPlaneSort, vector<string> const *comments)const 506 506 { -
gfanlib/gfanlib_polymakefile.cpp
rb7b2c5 r19addd1 53 53 list<int> ret; 54 54 int c=s.peek(); 55 while(( c>='0') && (c<='9')|| (c==' '))55 while(((c>='0') && (c<='9')) || (c==' ')) 56 56 { 57 57 // fprintf(Stderr,"?\n"); … … 66 66 namespace gfan{ 67 67 PolymakeProperty::PolymakeProperty(const std::string &name_, const std::string &value_): 68 name(name_),69 value(value_)68 value(value_), 69 name(name_) 70 70 { 71 71 } … … 271 271 272 272 273 bool PolymakeFile::readBooleanProperty(const char * p)273 bool PolymakeFile::readBooleanProperty(const char */*p*/) 274 274 { 275 275 return false; … … 277 277 278 278 279 void PolymakeFile::writeBooleanProperty(const char * p, bool n)279 void PolymakeFile::writeBooleanProperty(const char */*p*/, bool /*n*/) 280 280 { 281 281 } … … 322 322 stringstream t; 323 323 324 if(comments)assert( comments->size()>=m.getHeight());324 if(comments)assert((int)comments->size()>=m.getHeight()); 325 325 if(isXml) 326 326 { … … 389 389 { 390 390 t<<"<incidence_matrix>"; 391 for( inti=0;i<m.size();i++)391 for(unsigned i=0;i<m.size();i++) 392 392 { 393 393 t<<"<set>"; … … 405 405 else 406 406 { 407 for( inti=0;i<m.size();i++)407 for(unsigned i=0;i<m.size();i++) 408 408 { 409 409 t<<'{'; … … 445 445 { 446 446 t<<"<vector>"; 447 for( inti=0;i<v.size();i++)447 for(unsigned i=0;i<v.size();i++) 448 448 { 449 449 if(i!=0)t<<" "; … … 454 454 else 455 455 { 456 for( inti=0;i<v.size();i++)456 for(unsigned i=0;i<v.size();i++) 457 457 { 458 458 if(i!=0)t<<" "; -
gfanlib/gfanlib_q.h
rb7b2c5 r19addd1 156 156 return mpq_sgn(value); 157 157 } 158 static Rational gcd(Rational const &a, Rational const &b, Rational &s, Rational &t)158 static Rational gcd(Rational const &a, Rational const /*&b*/, Rational /*&s*/, Rational /*t*/) 159 159 { 160 160 /* mpz_t r; -
gfanlib/gfanlib_symmetriccomplex.cpp
rb7b2c5 r19addd1 15 15 16 16 SymmetricComplex::Cone::Cone(std::set<int> const &indices_, int dimension_, Integer multiplicity_, bool sortWithSymmetry, SymmetricComplex const &complex): 17 isKnownToBeNonMaximalFlag(false), 17 18 dimension(dimension_), 18 19 multiplicity(multiplicity_), 19 isKnownToBeNonMaximalFlag(false),20 20 sortKeyPermutation(complex.n) 21 21 { … … 27 27 ZMatrix const &vertices=complex.getVertices(); 28 28 ZVector sum(vertices.getWidth()); 29 for( inti=0;i<indices.size();i++)29 for(unsigned i=0;i<indices.size();i++) 30 30 sum+=vertices[indices[i]]; 31 31 … … 56 56 ZMatrix const &vertices=complex.getVertices(); 57 57 ZVector sum(vertices.getWidth()); 58 for( inti=0;i<indices.size();i++)58 for(unsigned i=0;i<indices.size();i++) 59 59 sum+=vertices[indices[i]]; 60 60 … … 62 62 Permutation const &bestPermutation=sortKeyPermutation; 63 63 64 assert( bestPermutation.size()==n);64 assert((int)bestPermutation.size()==n); 65 65 66 66 IntVector indicesNew(indices.size()); 67 67 int I=0; 68 for( inti=0;i<indices.size();i++,I++)68 for(unsigned i=0;i<indices.size();i++,I++) 69 69 { 70 70 ZVector ny=bestPermutation.apply(complex.vertices[indices[i]]); … … 80 80 { 81 81 std::set<int> ret; 82 for( inti=0;i<indices.size();i++)82 for(unsigned i=0;i<indices.size();i++) 83 83 ret.insert(indices[i]); 84 84 … … 89 89 { 90 90 int next=0; 91 for( inti=0;i<indices.size();i++)91 for(unsigned i=0;i<indices.size();i++) 92 92 { 93 93 while(1) 94 94 { 95 if(next>= c.indices.size())return false;95 if(next>=(int)c.indices.size())return false; 96 96 if(indices[i]==c.indices[next])break; 97 97 next++; … … 105 105 { 106 106 std::set<int> r; 107 for( inti=0;i<indices.size();i++)107 for(unsigned i=0;i<indices.size();i++) 108 108 { 109 109 ZVector ny=permutation.apply(complex.vertices[indices[i]]); … … 139 139 { 140 140 ZMatrix l; 141 for( inti=0;i<indices.size();i++)141 for(unsigned i=0;i<indices.size();i++) 142 142 l.appendRow(complex.vertices[indices[i]]); 143 143 … … 150 150 SymmetricComplex::SymmetricComplex(ZMatrix const &rays, ZMatrix const &linealitySpace_, SymmetryGroup const &sym_): 151 151 n(rays.getWidth()), 152 linealitySpace(canonicalizeSubspace(linealitySpace_)), 152 153 sym(sym_), 153 dimension(-1), 154 linealitySpace(canonicalizeSubspace(linealitySpace_)) 154 dimension(-1) 155 155 { 156 156 assert(rays.getWidth()==linealitySpace.getWidth()); … … 222 222 } 223 223 224 /* 224 #if 0 225 225 IntVector SymmetricComplex::dimensionsAtInfinity()const 226 226 { … … 228 228 dimension of the intersection of each cone in the complex with 229 229 the plane x_0=0 */ 230 /*231 230 IntVector ret(cones.size()); 232 231 … … 247 246 return ret; 248 247 } 249 */ 248 #endif 250 249 251 250 void SymmetricComplex::buildConeLists(bool onlyMaximal, bool compressed, std::vector<std::vector<IntVector > >*conelist/*, ZMatrix *multiplicities*/)const … … 259 258 int numberOfOrbitsOutput=0; 260 259 int numberOfOrbitsOfThisDimension=0; 261 bool newDimension=true;260 // bool newDimension=true; 262 261 { 263 262 int I=0; 264 263 for(ConeContainer::const_iterator i=cones.begin();i!=cones.end();i++,I++) 265 266 267 264 if(i->dimension==d) 265 { 266 numberOfOrbitsOfThisDimension++; 268 267 if(!onlyMaximal || isMaximal(*i)) 269 268 { 270 269 numberOfOrbitsOutput++; 271 bool isMax=isMaximal(*i);272 bool newOrbit=true;270 // bool isMax=isMaximal(*i); 271 // bool newOrbit=true; 273 272 std::set<std::set<int> > temp; 274 275 273 for(SymmetryGroup::ElementContainer::const_iterator k=sym.elements.begin();k!=sym.elements.end();k++) 274 { 276 275 Cone temp1=i->permuted(*k,*this,false); 277 276 temp.insert(temp1.indexSet()); … … 291 290 *multiplicities << std::endl; 292 291 }*/ 293 newOrbit=false;294 newDimension=false;292 // newOrbit=false; 293 // newDimension=false; 295 294 } 296 }297 298 } 299 } 300 301 } 302 303 std::string SymmetricComplex::toStringJustCones(int dimLow, int dimHigh, bool onlyMaximal, bool group, std::ostream *multiplicities, bool compressed, bool tPlaneSort)const295 } 296 } 297 } 298 } 299 300 } 301 302 std::string SymmetricComplex::toStringJustCones(int dimLow, int dimHigh, bool onlyMaximal, bool group, std::ostream *multiplicities, bool compressed, bool /*tPlaneSort*/)const 304 303 { 305 304 std::stringstream ret; … … 392 391 { 393 392 bool isBounded=true; 394 for( intj=0;j<i->indices.size();j++)393 for(unsigned j=0;j<i->indices.size();j++) 395 394 if(vertices[i->indices[j]][0].sign()==0)isBounded=false; 396 395 doAdd=isBounded; … … 692 691 { 693 692 ZMatrix generators(indices.size(),getAmbientDimension()); 694 for( inti=0;i<indices.size();i++)693 for(unsigned i=0;i<indices.size();i++) 695 694 generators[i]=vertices[indices[i]]; 696 695 return ZCone::givenByRays(generators,linealitySpace); -
gfanlib/gfanlib_symmetry.cpp
rb7b2c5 r19addd1 25 25 TrieNode(IntVector const &v, int i) 26 26 { 27 if(i<v.size())27 if(i<(int)v.size()) 28 28 m[v[i]]=TrieNode(v,i+1); 29 29 } … … 31 31 { 32 32 int ret=0; 33 if(i== v.size())return 1;33 if(i==(int)v.size())return 1; 34 34 for(Map::const_iterator j=m.begin();j!=m.end();j++) 35 35 { … … 41 41 void search(ZVector const &v, ZVector &building, Permutation &tempPerm, Permutation &ret, ZVector &optimal, int i, bool &isImproving)const 42 42 { 43 if(i==v.size()){ret=tempPerm;optimal=building;isImproving=false;return;}43 if(i==(int)v.size()){ret=tempPerm;optimal=building;isImproving=false;return;} 44 44 if(isImproving) 45 45 building[i]=-0x7fffffff; … … 61 61 void searchStabalizer(ZVector const &v, ZVector &building, Permutation &tempPerm, Permutation &ret, ZVector &optimal, int i, bool &isImproving, ZVector const &toBeFixed)const 62 62 { 63 if(i==v.size())63 if(i==(int)v.size()) 64 64 if(!(tempPerm.apply(v)<optimal)) 65 65 { … … 102 102 void insert(Permutation const &v, int i) 103 103 { 104 if(i==v.size())return;104 if(i==(int)v.size())return; 105 105 if(m.count(v[i])) 106 106 m[v[i]].insert(v,i+1); … … 132 132 int n; 133 133 Trie(int n_): 134 n(n_),135 theTree(Permutation(n_),0)134 theTree(Permutation(n_),0), 135 n(n_) 136 136 { 137 137 } … … 276 276 IntVector ret(size()); 277 277 assert(size()==b.size()); 278 for( inti=0;i<size();i++)ret[i]=b[(*this)[i]];278 for(unsigned i=0;i<size();i++)ret[i]=b[(*this)[i]]; 279 279 return Permutation(ret); 280 280 } … … 284 284 IntVector ret(size()); 285 285 assert(size()==b.size()); 286 for( inti=0;i<size();i++)ret[(*this)[i]]=b[i];286 for(unsigned i=0;i<size();i++)ret[(*this)[i]]=b[i]; 287 287 return Permutation(ret); 288 288 } … … 292 292 IntVector ret(size()); 293 293 assert(size()==v.size()); 294 for( inti=0;i<size();i++)ret[i]=v[(*this)[i]];294 for(unsigned i=0;i<size();i++)ret[i]=v[(*this)[i]]; 295 295 return ret; 296 296 } … … 300 300 ZVector ret(size()); 301 301 assert(size()==v.size()); 302 for( inti=0;i<size();i++)ret[i]=v[(*this)[i]];302 for(unsigned i=0;i<size();i++)ret[i]=v[(*this)[i]]; 303 303 return ret; 304 304 } … … 309 309 ZVector ret(size()); 310 310 assert(size()==v.size()); 311 for( inti=0;i<size();i++)ret[(*this)[i]]=v[i];311 for(unsigned i=0;i<size();i++)ret[(*this)[i]]=v[i]; 312 312 return ret; 313 313 } … … 431 431 ZVector Permutation::fundamentalDomainInequality()const 432 432 { 433 for( inti=0;i<size();i++)434 if((*this)[i]!= i)433 for(unsigned i=0;i<size();i++) 434 if((*this)[i]!=(int)i) 435 435 return ZVector::standardVector(size(),i)-ZVector::standardVector(size(),(*this)[i]); 436 436 return ZVector(size()); -
gfanlib/gfanlib_vector.h
rb7b2c5 r19addd1 101 101 if(size()<b.size())return true; 102 102 if(size()>b.size())return false; 103 for( inti=0;i<size();i++)103 for(unsigned i=0;i<size();i++) 104 104 { 105 105 if(v[i]<b[i])return true; … … 156 156 return pq*pq==pp*qq; 157 157 */ 158 intn=p.size();158 unsigned n=p.size(); 159 159 assert(n==q.size()); 160 inti;160 unsigned i; 161 161 for(i=0;i<n;i++) 162 162 { … … 167 167 typ a=p.v[i]; 168 168 typ b=q.v[i]; 169 for( intj=0;j<n;j++)169 for(unsigned j=0;j<n;j++) 170 170 if(a*q.v[j]!=b*p.v[j])return false; 171 171 return true; … … 175 175 // Arithmetic slow 176 176 //----------------- 177 inline friend Vector operator*(typ s, const Vector& q){Vector p=q;for( inti=0;i<q.size();i++)p[i]*=s;return p;}177 inline friend Vector operator*(typ s, const Vector& q){Vector p=q;for(unsigned i=0;i<q.size();i++)p[i]*=s;return p;} 178 178 // inline friend Vektor operator/(const Vektor& q, typ s){Vektor p=q;for(int i=0;i<q.size();i++)p[i]/=s;return p;} 179 179 /* inline friend Vector operator*(const Vektor& p, const Vektor& q){assert(p.size()==q.size());Vektor p1=p;for(int i=0;i<p.size();i++)p1.v[i]*=q.v[i];return p1;} … … 181 181 */ 182 182 inline friend Vector operator/(const Vector& p, typ const &s){Vector ret(p.size());for(unsigned i=0;i<p.size();i++)ret[i]=p[i]/s;return ret;} 183 inline friend Vector operator+(const Vector& p, const Vector& q){assert(p.size()==q.size());Vector p1=p;for( inti=0;i<p.size();i++)p1[i]+=q[i];return p1;}184 inline friend Vector operator-(const Vector& p, const Vector& q){assert(p.size()==q.size());Vector p1=p;for( inti=0;i<p.size();i++)p1[i]-=q[i];return p1;}185 inline friend Vector max(const Vector& p, const Vector& q){assert(p.size()==q.size());Vector p1=p;for( inti=0;i<p.size();i++)if(p1[i]<q[i])p1[i]=q[i];return p1;}186 inline friend Vector min(const Vector& p, const Vector& q){assert(p.size()==q.size());Vector p1=p;for( inti=0;i<p.size();i++)if(p1[i]>q[i])p1[i]=q[i];return p1;}183 inline friend Vector operator+(const Vector& p, const Vector& q){assert(p.size()==q.size());Vector p1=p;for(unsigned i=0;i<p.size();i++)p1[i]+=q[i];return p1;} 184 inline friend Vector operator-(const Vector& p, const Vector& q){assert(p.size()==q.size());Vector p1=p;for(unsigned i=0;i<p.size();i++)p1[i]-=q[i];return p1;} 185 inline friend Vector max(const Vector& p, const Vector& q){assert(p.size()==q.size());Vector p1=p;for(unsigned i=0;i<p.size();i++)if(p1[i]<q[i])p1[i]=q[i];return p1;} 186 inline friend Vector min(const Vector& p, const Vector& q){assert(p.size()==q.size());Vector p1=p;for(unsigned i=0;i<p.size();i++)if(p1[i]>q[i])p1[i]=q[i];return p1;} 187 187 188 188 friend Vector operator-(const Vector &b) 189 189 { 190 190 Vector ret(b.size()); 191 for( inti=0;i<b.size();i++)ret[i]=-b[i];191 for(unsigned i=0;i<b.size();i++)ret[i]=-b[i]; 192 192 return ret; 193 193 } … … 226 226 { 227 227 assert(begin>=0); 228 assert(end<= size());228 assert(end<=(int)size()); 229 229 assert(end>=begin); 230 230 Vector ret(end-begin); … … 293 293 typ temp1,temp2; 294 294 typ ret(1); 295 for( inti=0;i<size();i++)295 for(unsigned i=0;i<size();i++) 296 296 ret=typ::gcd(ret,v[i],temp1,temp2); 297 297 return ret; -
gfanlib/gfanlib_zcone.cpp
rb7b2c5 r19addd1 44 44 dd_colrange d_input,j; 45 45 dd_RepresentationType rep=dd_Inequality; 46 dd_boolean found=dd_FALSE, newformat=dd_FALSE, successful=dd_FALSE;47 char command[dd_linelenmax], comsave[dd_linelenmax];46 // dd_boolean found=dd_FALSE, newformat=dd_FALSE, successful=dd_FALSE; 47 // char command[dd_linelenmax], comsave[dd_linelenmax]; 48 48 dd_NumberType NT; 49 49 50 50 (*Error)=dd_NoError; 51 51 52 rep=dd_Inequality; newformat=dd_TRUE;52 rep=dd_Inequality; // newformat=dd_TRUE; 53 53 54 54 m_input=g.getHeight(); … … 70 70 } 71 71 72 successful=dd_TRUE;72 // successful=dd_TRUE; 73 73 74 74 return M; … … 108 108 { 109 109 bool ret; 110 dd_MatrixPtr M=NULL ,M2=NULL,M3=NULL;111 dd_colrange d;110 dd_MatrixPtr M=NULL/*,M2=NULL,M3=NULL*/; 111 // dd_colrange d; 112 112 dd_ErrorType err=dd_NoError; 113 dd_rowset redrows,linrows,ignoredrows, basisrows;114 dd_colset ignoredcols, basiscols;115 dd_DataFileType inputfile;116 FILE *reading=NULL;113 // dd_rowset redrows,linrows,ignoredrows, basisrows; 114 // dd_colset ignoredcols, basiscols; 115 // dd_DataFileType inputfile; 116 // FILE *reading=NULL; 117 117 118 118 cddinitGmp(); … … 121 121 if (err!=dd_NoError) goto _L99; 122 122 123 d=M->colsize;123 // d=M->colsize; 124 124 125 125 static dd_Arow temp; … … 339 339 if(numberOfRows==0)return;//the full space, so description is already irredundant 340 340 341 dd_rowset r=NULL;341 // dd_rowset r=NULL; 342 342 ZMatrix g=inequalities; 343 343 g.append(equations); 344 344 345 dd_LPSolverType solver=dd_DualSimplex;345 // dd_LPSolverType solver=dd_DualSimplex; 346 346 dd_MatrixPtr A=NULL; 347 347 dd_ErrorType err=dd_NoError; … … 402 402 int numberOfRows=numberOfEqualities+numberOfInequalities; 403 403 404 dd_rowset r=NULL;404 // dd_rowset r=NULL; 405 405 ZMatrix g=inequalities; 406 406 g.append(equations); … … 448 448 void dual(ZMatrix const &inequalities, ZMatrix const &equations, ZMatrix &dualInequalities, ZMatrix &dualEquations) 449 449 { 450 int result;450 // int result; 451 451 452 452 dd_MatrixPtr A=NULL; … … 472 472 473 473 return; 474 _L99:475 assert(0);474 // _L99: 475 // assert(0); 476 476 } 477 477 // this procedure is take from cddio.c. … … 556 556 int dim2=inequalities.getHeight(); 557 557 if(dim2==0)return std::vector<std::vector<int> >(); 558 int dimension=inequalities.getWidth();558 // int dimension=inequalities.getWidth(); 559 559 560 560 dd_MatrixPtr A=NULL; … … 603 603 604 604 return ret; 605 _L99:606 assert(0);607 return std::vector<std::vector<int> >();605 // _L99: 606 // assert(0); 607 // return std::vector<std::vector<int> >(); 608 608 } 609 609 … … 707 707 inequalities=LpSolver::fastNormals(inequalities2); 708 708 goto noFallBack; 709 fallBack://alternativ (disabled)710 lpSolver.removeRedundantRows(inequalities,equations,true);709 // fallBack://alternativ (disabled) 710 // lpSolver.removeRedundantRows(inequalities,equations,true); 711 711 noFallBack:; 712 712 } … … 730 730 } 731 731 732 std::ostream &operator<<(std::ostream &f, ZCone const &c)732 void operator<<(std::ostream &f, ZCone const &c) 733 733 { 734 734 f<<"Ambient dimension:"<<c.n<<std::endl; … … 741 741 742 742 ZCone::ZCone(int ambientDimension): 743 preassumptions(PCP_impliedEquationsKnown|PCP_facetsKnown), 744 state(1), 745 n(ambientDimension), 746 multiplicity(1), 747 linearForms(ZMatrix(0,ambientDimension)), 743 748 inequalities(ZMatrix(0,ambientDimension)), 744 749 equations(ZMatrix(0,ambientDimension)), 745 n(ambientDimension), 746 state(1), 747 preassumptions(PCP_impliedEquationsKnown|PCP_facetsKnown), 750 haveExtremeRaysBeenCached(false) 751 { 752 } 753 754 755 ZCone::ZCone(ZMatrix const &inequalities_, ZMatrix const &equations_, int preassumptions_): 756 preassumptions(preassumptions_), 757 state(0), 758 n(inequalities_.getWidth()), 748 759 multiplicity(1), 749 haveExtremeRaysBeenCached(false), 750 linearForms(ZMatrix(0,ambientDimension)) 751 { 752 } 753 754 755 ZCone::ZCone(ZMatrix const &inequalities_, ZMatrix const &equations_, int preassumptions_): 760 linearForms(ZMatrix(0,inequalities_.getWidth())), 756 761 inequalities(inequalities_), 757 762 equations(equations_), 758 state(0), 759 preassumptions(preassumptions_), 760 multiplicity(1), 761 haveExtremeRaysBeenCached(false), 762 n(inequalities_.getWidth()), 763 linearForms(ZMatrix(0,inequalities_.getWidth())) 763 haveExtremeRaysBeenCached(false) 764 764 { 765 765 assert(preassumptions_<4);//OTHERWISE WE ARE DOING SOMETHING STUPID LIKE SPECIFYING AMBIENT DIMENSION … … 1042 1042 std::vector<std::vector<int> > indices=lpSolver.extremeRaysInequalityIndices(inequalities); 1043 1043 1044 for( inti=0;i<indices.size();i++)1044 for(unsigned i=0;i<indices.size();i++) 1045 1045 { 1046 1046 /* At this point we know lineality space, implied equations and … … 1082 1082 1083 1083 std::vector<int> asVector(inequalities.getHeight()); 1084 for( intj=0;j<indices[i].size();j++){asVector[indices[i][j]]=1;}1084 for(unsigned j=0;j<indices[i].size();j++){asVector[indices[i][j]]=1;} 1085 1085 ZMatrix equations=this->equations; 1086 1086 ZVector theInequality; 1087 1087 1088 for( intj=0;j<asVector.size();j++)1088 for(unsigned j=0;j<asVector.size();j++) 1089 1089 if(asVector[j]) 1090 1090 equations.appendRow(inequalities[j]); … … 1229 1229 ZCone ZCone::faceContaining(ZVector const &v)const 1230 1230 { 1231 assert(n== v.size());1231 assert(n==(int)v.size()); 1232 1232 assert(contains(v)); 1233 1233 ZMatrix newEquations=equations; -
gfanlib/gfanlib_zcone.h
rb7b2c5 r19addd1 352 352 */ 353 353 // PolyhedralCone projection(int newn)const; 354 friend std::ostream &operator<<(std::ostream &f, ZCone const &c);354 friend void operator<<(std::ostream &f, ZCone const &c); 355 355 }; 356 356 -
gfanlib/gfanlib_zfan.cpp
rb7b2c5 r19addd1 16 16 { 17 17 assert(dimension>=0); 18 if(dimension>= T.size())return 0;18 if(dimension>=(int)T.size())return 0; 19 19 return T[dimension].size(); 20 20 } … … 162 162 163 163 // log2 cerr<< "Number of orbits to expand "<<cones.size()<<endl; 164 for( inti=0;i<cones.size();i++)164 for(unsigned i=0;i<cones.size();i++) 165 165 // if(coneIndices==0 || coneIndices->count(i)) 166 166 { … … 202 202 } 203 203 ZFan::ZFan(ZFan const& f): 204 coneCollection(0), 204 205 complex(0), 205 coneCollection(0),206 206 cones(f.table(0,0)), 207 207 maximalCones(f.table(0,1)),
Note: See TracChangeset
for help on using the changeset viewer.