Changeset 2e4757c in git
- Timestamp:
- Apr 13, 2011, 3:39:30 PM (13 years ago)
- Branches:
- (u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', 'fc741b6502fd8a97288eaa3eba6e5220f3c3df87')
- Children:
- ee33a7ced725d6bd7bcfb344c937981897a8a06a
- Parents:
- 91a72fa06f22b35afc357574afa7c13204ff616e
- git-author:
- Hans Schoenemann <hannes@mathematik.uni-kl.de>2011-04-13 15:39:30+02:00
- git-committer:
- Mohamed Barakat <mohamed.barakat@rwth-aachen.de>2011-11-09 12:12:31+01:00
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
libpolys/polys/simpleideals.cc
r91a72f r2e4757c 430 430 * for idSort: compare a and b revlex inclusive module comp. 431 431 */ 432 static int p_Comp_RevLex(poly a, poly b,BOOLEAN nolex, const ring r)432 static int p_Comp_RevLex(poly a, poly b,BOOLEAN nolex, const ring R) 433 433 { 434 434 if (b==NULL) return 1; … … 437 437 if (nolex) 438 438 { 439 int r=p_LmCmp(a,b, r);439 int r=p_LmCmp(a,b,R); 440 440 if (r!=0) return r; 441 number h=n_Sub(pGetCoeff(a),pGetCoeff(b), r->cf);442 r = -1+n_IsZero(h, r->cf)+2*n_GreaterZero(h,r->cf); /* -1: <, 0:==, 1: > */443 n_Delete(&h, r->cf);441 number h=n_Sub(pGetCoeff(a),pGetCoeff(b),R->cf); 442 r = -1+n_IsZero(h,R->cf)+2*n_GreaterZero(h,R->cf); /* -1: <, 0:==, 1: > */ 443 n_Delete(&h, R->cf); 444 444 return r; 445 445 } 446 int l=rVar( r);447 while ((l>0) && (p_GetExp(a,l, r)==pGetExp(b,l,r))) l--;446 int l=rVar(R); 447 while ((l>0) && (p_GetExp(a,l,R)==p_GetExp(b,l,R))) l--; 448 448 if (l==0) 449 449 { 450 if (p_GetComp(a, r)==p_GetComp(b,r))451 { 452 number h=n_Sub(pGetCoeff(a),pGetCoeff(b), r->cf);453 int r = -1+n_IsZero(h, r->cf)+2*n_GreaterZero(h,r->cf); /* -1: <, 0:==, 1: > */454 n_Delete(&h, r->cf);450 if (p_GetComp(a,R)==p_GetComp(b,R)) 451 { 452 number h=n_Sub(pGetCoeff(a),pGetCoeff(b),R->cf); 453 int r = -1+n_IsZero(h,R->cf)+2*n_GreaterZero(h,R->cf); /* -1: <, 0:==, 1: > */ 454 n_Delete(&h,R->cf); 455 455 return r; 456 456 } 457 if (p_GetComp(a, r)>p_GetComp(b,r)) return 1;458 } 459 else if (p_GetExp(a,l, r)>p_GetExp(b,l,r))457 if (p_GetComp(a,R)>p_GetComp(b,R)) return 1; 458 } 459 else if (p_GetExp(a,l,R)>p_GetExp(b,l,R)) 460 460 return 1; 461 461 return -1;
Note: See TracChangeset
for help on using the changeset viewer.