Changeset d4f1b95 in git
- Timestamp:
- Nov 24, 2009, 10:19:36 AM (14 years ago)
- Branches:
- (u'spielwiese', '828514cf6e480e4bafc26df99217bf2a1ed1ef45')
- Children:
- ceaa049b15f01877c8cd64ac61771aff342c1da7
- Parents:
- d945f30c526b2f3c4c937f5acb82e282284898d9
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
kernel/gfan.cc
rd945f3 rd4f1b95 126 126 facet::facet(const facet& f) 127 127 { 128 // this->setFacetNormal(f.getFacetNormal()); 129 // this->setUCN(f.getUCN()); 130 // this->isFlippable=f.isFlippable; 131 // facet* f2Copy; 132 // f2Copy=f.codim2Ptr; 133 // facet* f2Act; 134 // f2Act=this->codim2Ptr; 135 // while(f2Copy!=NULL) 136 // { 137 // if(f2Act==NULL) 138 // { 139 // f2Act=new facet(2); 140 // } 141 // else 142 // { 143 // f2Act->next = new facet(2); 144 // } 145 // f2Act->setFacetNormal(f2Copy->getFacetNormal()); 146 // f2Copy = f2Copy->next; 147 // } 128 148 } 129 149 … … 220 240 } 221 241 222 /** Return the flipped GB*/ 242 /** Return the flipped GB 243 Seems not be used 244 Anyhow idCopy would make sense here. 245 */ 223 246 ideal facet::getFlipGB() 224 247 { … … 245 268 int facet::getUCN() 246 269 { 247 if(this!=NULL && ( this!=(facet *)0xfbfbfbfb || this!=(facet *)0xfbfbfbfbfbfbfbfb) )270 if(this!=NULL && this!=(facet *)0xfbfbfbfbfbfbfbfb)// || this!=(facet *)0xfbfbfbfb) ) 248 271 // if(this!=NULL && ( this->fNormal!=(intvec *)0xfbfbfbfb || this->fNormal!=(intvec *)0xfbfbfbfbfbfbfbfb) ) 249 272 return this->UCN; … … 409 432 intvec *gcone::getIntPoint() 410 433 { 411 return this->ivIntPt;434 return ivCopy(this->ivIntPt); 412 435 } 413 436 … … 1724 1747 gcTmp->normalize(); 1725 1748 #ifdef gfan_DEBUG 1726 gcTmp->showFacets(1);1749 // gcTmp->showFacets(1); 1727 1750 #endif 1728 1751 /*add facets to SLA here*/ … … 1735 1758 for(int ii=0;ii<IDELEMS(gcTmp->gcBasis);ii++) 1736 1759 { 1737 gcTmp->gcBasis->m[ii]=(poly)NULL; 1760 // gcTmp->gcBasis->m[ii]=(poly)NULL; 1761 pDelete(&gcTmp->gcBasis->m[ii]); 1738 1762 } 1739 1763 } 1740 1764 1741 1765 #if gfan_DEBUG 1742 if(SearchListRoot!=NULL)1743 gcTmp->showSLA(*SearchListRoot);1766 // if(SearchListRoot!=NULL) 1767 // gcTmp->showSLA(*SearchListRoot); 1744 1768 #endif 1745 1769 rChangeCurrRing(gcAct->baseRing); … … 1988 2012 slEndStatic = slEnd; 1989 2013 /*1st step: compare facetNormals*/ 1990 intvec *fNormal=NULL; //= new intvec(this->numVars);1991 intvec *f2Normal=NULL; //= new intvec(this->numVars);1992 intvec *slNormal=NULL; //= new intvec(this->numVars);1993 intvec *sl2Normal=NULL; //= new intvec(this->numVars);2014 intvec *fNormal=NULL; 2015 intvec *f2Normal=NULL; 2016 intvec *slNormal=NULL; 2017 intvec *sl2Normal=NULL; 1994 2018 1995 2019 while(fAct!=NULL) … … 2055 2079 //while(slAct!=slEndStatic->next) 2056 2080 { 2057 //if(deleteMarker!=NULL)2058 //{2081 if(deleteMarker!=NULL) 2082 { 2059 2083 // delete deleteMarker; 2060 2084 // deleteMarker=NULL; 2061 //}2085 } 2062 2086 removalOccured=FALSE; 2063 2087 slNormal = slAct->getFacetNormal(); … … 2072 2096 { 2073 2097 notParallelCtr++; 2074 // slAct = slAct->next;2075 2098 } 2076 2099 else //fN and slN are parallel … … 2108 2131 if(slAct==slHead) //We want to delete the first element of SearchList 2109 2132 { 2110 deleteMarker = sl Head;2133 deleteMarker = slAct; 2111 2134 slHead = slAct->next; 2112 2135 if(slHead!=NULL) 2113 2136 slHead->prev = NULL; 2114 /*delete deleteMarker*/; 2115 2137 // delete deleteMarker; deleteMarker=NULL; 2138 //set a bool flag to mark slAct as to be deleted 2116 2139 }//NOTE find a way to delete without affecting slAct = slAct->next 2117 2140 else if(slAct==slEndStatic) … … 2129 2152 slEndStatic->next->prev = slEndStatic->prev; 2130 2153 slEndStatic = slEndStatic->prev; 2131 2154 //slEnd = slEndStatic; 2132 2155 } 2133 2156 } … … 2137 2160 slAct->prev->next = slAct->next; 2138 2161 slAct->next->prev = slAct->prev; 2139 } 2140 2162 } 2141 2163 //update lengthOfSearchList 2142 lengthOfSearchList--; 2143 //delete slAct; 2144 //slAct=NULL; 2145 //slAct = slAct->next; //not needed, since facets are equal 2146 // delete deleteMarker; 2147 //deleteMarker=NULL; 2148 //fAct = fAct->next; 2164 lengthOfSearchList--; 2149 2165 break; 2150 2166 }//if(ctr==fAct->numCodim2Facets) … … 2152 2168 //is no other facet later in SLA that might be equal. 2153 2169 { 2154 maybe=TRUE; 2155 // if(slAct->next==NULL && maybe==TRUE) 2156 // { 2157 // doNotAdd=FALSE; 2158 // slAct = slAct->next; 2159 // break; 2160 // } 2161 // else 2162 // slAct=slAct->next; 2163 } 2164 //slAct = slAct->next; 2165 //delete deleteMarker; 2170 maybe=TRUE;// 2171 } 2166 2172 }//if(!isParallel(fNormal, slNormal)) 2167 2173 if( (slAct->next==NULL) && (maybe==TRUE) ) … … 2176 2182 be a way it gets called twice thus ommiting one facet: 2177 2183 slAct = slAct->next;*/ 2178 //delete deleteMarker;2179 //deleteMarker=NULL;2184 // delete deleteMarker; 2185 // deleteMarker=NULL; 2180 2186 //if slAct was marked as to be deleted, delete it here! 2181 2187 }//while(slAct!=NULL) … … 2206 2212 slEnd->prev = marker; 2207 2213 //Copy codim2-facets 2208 intvec *f2Normal; // = new intvec(this->numVars);2214 intvec *f2Normal; 2209 2215 while(f2Act!=NULL) 2210 2216 { … … 2235 2241 }//while(fAct!=NULL) 2236 2242 return slHead; 2237 // delete sl2Normal;2238 // delete slNormal;2239 // delete f2Normal;2240 // delete fNormal;2241 2243 }//addC2N 2242 2244 … … 2574 2576 ring gcone::getBaseRing() 2575 2577 { 2576 return this->baseRing;2578 return rCopy(this->baseRing); 2577 2579 } 2578 2580 /** \brief Gather the output … … 2740 2742 gcDel = gcRoot; 2741 2743 gcAct = gcRoot; 2742 do2744 while(gcAct!=NULL) 2743 2745 { 2744 2746 gcDel = gcAct; 2745 2747 gcAct = gcAct->next; 2746 2748 delete gcDel; 2747 } while(gcAct!=NULL);2749 } 2748 2750 } 2749 2751 dd_free_global_constants();
Note: See TracChangeset
for help on using the changeset viewer.