Changeset 67d4b5 in git
- Timestamp:
- Jul 20, 2011, 1:34:49 PM (12 years ago)
- Branches:
- (u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', 'a800fe4b3e9d37a38c5a10cc0ae9dfa0c15a4ee6')
- Children:
- 9d4b8cd59b57658f59577acd8fb7501642b875a8
- Parents:
- 64606013a8ec045f8331f1c67fea0078d7e5148b
- git-author:
- Burcin Erocal <burcin@erocal.org>2011-07-20 13:34:49+02:00
- git-committer:
- Mohamed Barakat <mohamed.barakat@rwth-aachen.de>2011-11-09 13:01:22+01:00
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
kernel/kpolys.cc
r646060 r67d4b5 8 8 * LE(p2, j) != LE(lcm, j) ==> LCM(p2, p) != lcm 9 9 */ 10 BOOLEAN pCompareChain (poly p,poly p1,poly p2,poly lcm )10 BOOLEAN pCompareChain (poly p,poly p1,poly p2,poly lcm, const ring R = currRing) 11 11 { 12 12 int k, j; … … 14 14 if (lcm==NULL) return FALSE; 15 15 16 for (j=( currRing->N); j; j--)17 if ( p GetExp(p,j) > pGetExp(lcm,j)) return FALSE;16 for (j=(R->N); j; j--) 17 if ( p_GetExp(p,j, R) > p_GetExp(lcm,j, R)) return FALSE; 18 18 if ( pGetComp(p) != pGetComp(lcm)) return FALSE; 19 for (j=( currRing->N); j; j--)19 for (j=(R->N); j; j--) 20 20 { 21 if (p GetExp(p1,j)!=pGetExp(lcm,j))21 if (p_GetExp(p1,j, R)!=p_GetExp(lcm,j, R)) 22 22 { 23 if (p GetExp(p,j)!=pGetExp(lcm,j))23 if (p_GetExp(p,j, R)!=p_GetExp(lcm,j, R)) 24 24 { 25 for (k=( currRing->N); k>j; k--)25 for (k=(R->N); k>j; k--) 26 26 { 27 if ((p GetExp(p,k)!=pGetExp(lcm,k))28 && (p GetExp(p2,k)!=pGetExp(lcm,k)))27 if ((p_GetExp(p,k, R)!=p_GetExp(lcm,k, R)) 28 && (p_GetExp(p2,k, R)!=p_GetExp(lcm,k, R))) 29 29 return TRUE; 30 30 } 31 31 for (k=j-1; k; k--) 32 32 { 33 if ((p GetExp(p,k)!=pGetExp(lcm,k))34 && (p GetExp(p2,k)!=pGetExp(lcm,k)))33 if ((p_GetExp(p,k, R)!=p_GetExp(lcm,k, R)) 34 && (p_GetExp(p2,k, R)!=p_GetExp(lcm,k, R))) 35 35 return TRUE; 36 36 } … … 38 38 } 39 39 } 40 else if (p GetExp(p2,j)!=pGetExp(lcm,j))40 else if (p_GetExp(p2,j, R)!=p_GetExp(lcm,j, R)) 41 41 { 42 if (p GetExp(p,j)!=pGetExp(lcm,j))42 if (p_GetExp(p,j, R)!=p_GetExp(lcm,j, R)) 43 43 { 44 for (k=( currRing->N); k>j; k--)44 for (k=(R->N); k>j; k--) 45 45 { 46 if ((p GetExp(p,k)!=pGetExp(lcm,k))47 && (p GetExp(p1,k)!=pGetExp(lcm,k)))46 if ((p_GetExp(p,k, R)!=p_GetExp(lcm,k, R)) 47 && (p_GetExp(p1,k, R)!=p_GetExp(lcm,k, R))) 48 48 return TRUE; 49 49 } 50 50 for (k=j-1; k!=0 ; k--) 51 51 { 52 if ((p GetExp(p,k)!=pGetExp(lcm,k))53 && (p GetExp(p1,k)!=pGetExp(lcm,k)))52 if ((p_GetExp(p,k, R)!=p_GetExp(lcm,k, R)) 53 && (p_GetExp(p1,k, R)!=p_GetExp(lcm,k, R))) 54 54 return TRUE; 55 55 } … … 61 61 } 62 62 #ifdef HAVE_RATGRING 63 BOOLEAN pCompareChainPart (poly p,poly p1,poly p2,poly lcm) 63 BOOLEAN pCompareChainPart (poly p,poly p1,poly p2,poly lcm, 64 const ring R = currRing) 64 65 { 65 66 int k, j; … … 67 68 if (lcm==NULL) return FALSE; 68 69 69 for (j= currRing->real_var_end; j>=currRing->real_var_start; j--)70 if ( p GetExp(p,j) > pGetExp(lcm,j)) return FALSE;70 for (j=R->real_var_end; j>=R->real_var_start; j--) 71 if ( p_GetExp(p,j, R) > p_GetExp(lcm,j, R)) return FALSE; 71 72 if ( pGetComp(p) != pGetComp(lcm)) return FALSE; 72 for (j= currRing->real_var_end; j>=currRing->real_var_start; j--)73 for (j=R->real_var_end; j>=R->real_var_start; j--) 73 74 { 74 if (p GetExp(p1,j)!=pGetExp(lcm,j))75 if (p_GetExp(p1,j, R)!=p_GetExp(lcm,j, R)) 75 76 { 76 if (p GetExp(p,j)!=pGetExp(lcm,j))77 if (p_GetExp(p,j, R)!=p_GetExp(lcm,j, R)) 77 78 { 78 for (k=( currRing->N); k>j; k--)79 for (k= currRing->real_var_end; k>j; k--)79 for (k=(R->N); k>j; k--) 80 for (k=R->real_var_end; k>j; k--) 80 81 { 81 if ((p GetExp(p,k)!=pGetExp(lcm,k))82 && (p GetExp(p2,k)!=pGetExp(lcm,k)))82 if ((p_GetExp(p,k, R)!=p_GetExp(lcm,k, R)) 83 && (p_GetExp(p2,k, R)!=p_GetExp(lcm,k, R))) 83 84 return TRUE; 84 85 } 85 for (k=j-1; k>= currRing->real_var_start; k--)86 for (k=j-1; k>=R->real_var_start; k--) 86 87 { 87 if ((p GetExp(p,k)!=pGetExp(lcm,k))88 && (p GetExp(p2,k)!=pGetExp(lcm,k)))88 if ((p_GetExp(p,k, R)!=p_GetExp(lcm,k, R)) 89 && (p_GetExp(p2,k, R)!=p_GetExp(lcm,k, R))) 89 90 return TRUE; 90 91 } … … 92 93 } 93 94 } 94 else if (p GetExp(p2,j)!=pGetExp(lcm,j))95 else if (p_GetExp(p2,j, R)!=p_GetExp(lcm,j, R)) 95 96 { 96 if (p GetExp(p,j)!=pGetExp(lcm,j))97 if (p_GetExp(p,j, R)!=p_GetExp(lcm,j, R)) 97 98 { 98 for (k= currRing->real_var_end; k>j; k--)99 for (k=R->real_var_end; k>j; k--) 99 100 { 100 if ((p GetExp(p,k)!=pGetExp(lcm,k))101 && (p GetExp(p1,k)!=pGetExp(lcm,k)))101 if ((p_GetExp(p,k, R)!=p_GetExp(lcm,k, R)) 102 && (p_GetExp(p1,k, R)!=p_GetExp(lcm,k, R))) 102 103 return TRUE; 103 104 } 104 for (k=j-1; k>= currRing->real_var_start; k--)105 for (k=j-1; k>=R->real_var_start; k--) 105 106 { 106 if ((p GetExp(p,k)!=pGetExp(lcm,k))107 && (p GetExp(p1,k)!=pGetExp(lcm,k)))107 if ((p_GetExp(p,k, R)!=p_GetExp(lcm,k, R)) 108 && (p_GetExp(p1,k, R)!=p_GetExp(lcm,k, R))) 108 109 return TRUE; 109 110 }
Note: See TracChangeset
for help on using the changeset viewer.