Changeset 110e721 in git
- Timestamp:
- Dec 4, 2009, 1:17:19 PM (14 years ago)
- Branches:
- (u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
- Children:
- 1a68d1d2c26c94af9b8ee4cf98aad93763403437
- Parents:
- 60eb1ddb931cc15fe0ad0a8376dc9a7b83de0508
- Location:
- kernel
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
kernel/gfan.cc
r60eb1d r110e721 17 17 #include "ideals.h" 18 18 #include "kmatrix.h" 19 #include "fast_maps.h" //Mapping of ideals19 //#include "fast_maps.h" //Mapping of ideals 20 20 #include "maps.h" 21 21 #include "ring.h" … … 149 149 f2Act->prev = marker; 150 150 } 151 f2Act->setFacetNormal(f2Copy->getFacetNormal()); 151 intvec *f2Normal; 152 f2Normal = f2Copy->getFacetNormal(); 153 // f2Act->setFacetNormal(f2Copy->getFacetNormal()); 154 f2Act->setFacetNormal(f2Normal); 155 delete f2Normal; 152 156 f2Act->setUCN(f2Copy->getUCN()); 153 157 f2Copy = f2Copy->next; … … 181 185 idDelete((ideal *)&this->flipGB); 182 186 if(this->flipRing!=NULL && this->flipRing->idroot!=(idhdl)0xfbfbfbfbfbfbfbfb) 183 rDelete(this->flipRing);187 // rDelete(this->flipRing); 184 188 // this->flipRing=NULL; 185 189 this->prev=NULL; … … 288 292 int facet::getUCN() 289 293 { 290 if(this!=NULL && this!=(facet *)0xfbfbfbfbfbfbfbfb)// || this!=(facet *)0xfbfbfbfb) )294 if(this!=NULL)// && this!=(facet *)0xfbfbfbfbfbfbfbfb)// || this!=(facet *)0xfbfbfbfb) ) 291 295 // if(this!=NULL && ( this->fNormal!=(intvec *)0xfbfbfbfb || this->fNormal!=(intvec *)0xfbfbfbfbfbfbfbfb) ) 292 296 return this->UCN; … … 584 588 int gcone::getUCN() 585 589 { 586 if( this!=NULL && ( this!=(gcone * const)0xfbfbfbfbfbfbfbfb && this!=(gcone * const)0xfbfbfbfb ) )590 if( this!=NULL)// && ( this!=(gcone * const)0xfbfbfbfbfbfbfbfb && this!=(gcone * const)0xfbfbfbfb ) ) 587 591 return this->UCN; 588 592 else … … 1372 1376 { 1373 1377 //NOTE: Can't this be done without new? 1374 intvec *res = new intvec(iv->length());1378 intvec *res;// = new intvec(iv->length()); 1375 1379 res=ivCopy(iv); 1376 1380 *res *= (int)-1; -
kernel/gfan.h
r60eb1d r110e721 166 166 void showFacets(short codim=1); 167 167 volatile void showSLA(facet &f); 168 void idDebugPrint( ideal const&I);169 void invPrint( ideal const&I);170 bool isMonomial( ideal const&I);168 void idDebugPrint(const ideal &I); 169 void invPrint(const ideal &I); 170 bool isMonomial(const ideal &I); 171 171 intvec *ivNeg(const intvec *iv); 172 172 int dotProduct(const intvec &iva, const intvec &ivb); 173 bool isParallel( intvec const &a, intvec const&b);174 bool areEqual( intvec const &a, intvec const&b);173 bool isParallel(const intvec &a, const intvec &b); 174 bool areEqual(const intvec &a, const intvec &b); 175 175 bool areEqual(facet *f, facet *g); 176 176 int intgcd(int a, int b); 177 void writeConeToFile( gcone const&gc, bool usingIntPoints=FALSE);177 void writeConeToFile(const gcone &gc, bool usingIntPoints=FALSE); 178 178 void readConeFromFile(int gcNum, gcone *gc); 179 179 intvec f2M(gcone *gc, facet *f, int n=1); 180 180 181 181 //The real stuff 182 void getConeNormals( ideal const&I, bool compIntPoint=FALSE);183 void getCodim2Normals( gcone const&gc);182 void getConeNormals(const ideal &I, bool compIntPoint=FALSE); 183 void getCodim2Normals(const gcone &gc); 184 184 void flip(ideal gb, facet *f); 185 185 void computeInv(ideal &gb, ideal &inv, intvec &f); 186 186 // poly restOfDiv(poly const &f, ideal const &I); removed with r12286 187 ideal ffG( ideal const &H, ideal const&G);187 ideal ffG(const ideal &H, const ideal &G); 188 188 void getGB(ideal const &inputIdeal); 189 void interiorPoint( dd_MatrixPtr const&M, intvec &iv);190 ring rCopyAndAddWeight( ring const &r, intvec const*ivw);191 ring rCopyAndChangeWeight( ring const&r, intvec *ivw);189 void interiorPoint(const dd_MatrixPtr &M, intvec &iv); 190 ring rCopyAndAddWeight(const ring &r, const intvec *ivw); 191 ring rCopyAndChangeWeight(const ring &r, intvec *ivw); 192 192 // void reverseSearch(gcone *gcAct); //NOTE both removed from r12286 193 193 // bool isSearchFacet(gcone &gcTmp, facet *testfacet); 194 194 void noRevS(gcone &gcRoot, bool usingIntPoint=FALSE); 195 void makeInt( dd_MatrixPtr const &M, int const line, intvec &n);195 void makeInt(const dd_MatrixPtr &M, const int line, intvec &n); 196 196 void normalize(); 197 197 facet * enqueueNewFacets(facet *f); 198 dd_MatrixPtr facets2Matrix( gcone const&gc);198 dd_MatrixPtr facets2Matrix(const gcone &gc); 199 199 // static void gcone::idPrint(ideal &I); 200 200 friend class facet;
Note: See TracChangeset
for help on using the changeset viewer.