Changeset 871f01 in git
- Timestamp:
- Mar 17, 2011, 12:08:31 PM (12 years ago)
- Branches:
- (u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', 'a800fe4b3e9d37a38c5a10cc0ae9dfa0c15a4ee6')
- Children:
- fbb23ba606a10b85b9967f56b46d70c6ed807282
- Parents:
- 3df116abe842435d60d8e8dc559ed7309d4ef38f
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
kernel/gfan.cc
r3df116 r871f01 924 924 } 925 925 926 LObject *P = new sLObject(); TODO What's the difference between sLObject and LObject?926 LObject *P = new sLObject();//TODO What's the difference between sLObject and LObject? 927 927 memset(P,0,sizeof(LObject)); 928 928 … … 938 938 if(P->p!=NULL) spoly non zero=? 939 939 { 940 poly p; NOTE Don't use pInit here. Evil memleak will follow940 poly p;//NOTE Don't use pInit here. Evil memleak will follow 941 941 poly q; 942 942 poly pDel,qDel; … … 1206 1206 dd_FreeMatrix(posRestr); 1207 1207 } 1208 Clean up but don't delete the return value!1208 //Clean up but don't delete the return value! 1209 1209 dd_FreeMatrix(ddineq); 1210 1210 set_free(ddredrows);check … … 1674 1674 TODO Not sufficient -> vol2 II/125&127 1675 1675 NOTE Sufficient according to cddlibs doc. These ARE rays 1676 What the hell... let's just take interior points1676 //What the hell... let's just take interior points 1677 1677 if(gcone::hasHomInput==FALSE) 1678 1678 { … … 1837 1837 ideal srcRing_HH; 1838 1838 srcRing_H=idrCopyR(H,tmpRing); 1839 H is needed further below, so don't idDelete here1839 //H is needed further below, so don't idDelete here 1840 1840 srcRing_HH=ffG(srcRing_H,this->gcBasis); 1841 1841 idDelete(&srcRing_H); … … 1859 1859 for (int ii=0;ii<IDELEMS(srcRing_HH);ii++) 1860 1860 { 1861 poly aktpoly=(poly)srcRing_HH->m[ii]; This is a pointer, so don't pDelete1861 poly aktpoly=(poly)srcRing_HH->m[ii];//This is a pointer, so don't pDelete 1862 1862 iPMatrixRows = iPMatrixRows+pLength(aktpoly); 1863 1863 } … … 1871 1871 { 1872 1872 markingsAreCorrect=FALSE; crucial to initialise here 1873 poly aktpoly=srcRing_HH->m[ii]; Only a pointer, so don't pDelete1873 poly aktpoly=srcRing_HH->m[ii]; //Only a pointer, so don't pDelete 1874 1874 /*Comparison of leading monomials is done via exponent vectors*/ 1875 1875 for (int jj=0;jj<IDELEMS(H);jj++) … … 1965 1965 dd_set_si(intPointMatrix->matrix[aktrow][jj],1); 1966 1966 } 1967 Let's make sure we compute interior points from the positive orthant1967 //Let's make sure we compute interior points from the positive orthant 1968 1968 dd_MatrixPtr posRestr=dd_CreateMatrix(this->numVars,this->numVars+1); 1969 1969 … … 2209 2209 ideal srcRing_HH; 2210 2210 srcRing_H=idrCopyR(H,tmpRing); 2211 H is needed further below, so don't idDelete here2211 //H is needed further below, so don't idDelete here 2212 2212 srcRing_HH=ffG(srcRing_H,this->gcBasis); 2213 2213 idDelete(&srcRing_H); … … 2287 2287 { 2288 2288 poly initialFormElement; 2289 poly aktpoly = (poly)gb->m[ii]; Ptr, so don't pDelete(aktpoly)2289 poly aktpoly = (poly)gb->m[ii];//Ptr, so don't pDelete(aktpoly) 2290 2290 int *leadExpV=(int *)omAlloc((this->numVars+1)*sizeof(int)); 2291 2291 pGetExpV(aktpoly,leadExpV); find the leading exponent in leadExpV[1],...,leadExpV[n], use pNext(p) … … 2942 2942 for(int k=1;k<n;k++) 2943 2943 { 2944 Let's find a j s.t. m[j][k]!=0 && c[j]=02944 //Let's find a j s.t. m[j][k]!=0 && c[j]=0 2945 2945 int condCtr=0;Check each row for zeroness 2946 2946 for(int j=1;j<m;j++)
Note: See TracChangeset
for help on using the changeset viewer.