Changeset a13cae in git for kernel/gr_kstd2.cc
- Timestamp:
- Jul 27, 2008, 12:21:59 PM (15 years ago)
- Branches:
- (u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', 'a800fe4b3e9d37a38c5a10cc0ae9dfa0c15a4ee6')
- Children:
- 37d43ec78deb09e94d2b41b08283e6fe466687d7
- Parents:
- e4df0ba39c30b0d3f0f545e4cf34a78b4b74fbea
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
kernel/gr_kstd2.cc
re4df0ba ra13cae 2 2 * Computer Algebra System SINGULAR * 3 3 ****************************************/ 4 /* $Id: gr_kstd2.cc,v 1.1 5 2008-07-26 18:11:42Singular Exp $ */4 /* $Id: gr_kstd2.cc,v 1.16 2008-07-27 10:21:59 Singular Exp $ */ 5 5 /* 6 6 * ABSTRACT - Kernel: noncomm. alg. of Buchberger … … 229 229 } 230 230 //poly hh = nc_CreateSpoly(strat->S[c_j],(*h).p, currRing); 231 232 233 231 poly hh=c_p; c_p=NULL; 232 pDelete(&((*h).p)); 233 (*h).p=hh; 234 234 if (!TEST_OPT_INTSTRATEGY) 235 235 { … … 254 254 /*- try to reduce the s-polynomial again -*/ 255 255 pass++; 256 256 j=0; 257 257 } 258 258 else … … 265 265 { 266 266 int a_e=(pTotaldegree(strat->S[j],currRing)-pFDeg(strat->S[j],currRing)); 267 if (TEST_OPT_DEBUG) { Print("j=%d, e=%d\n",j,a_e); } 267 if(TEST_OPT_DEBUG) 268 { 269 p_wrp(h->p,currRing); Print(" divisibly by S[%d]=",j); 270 p_wrp(strat->S[j],currRing); Print(" e0%d\n",a_e); 271 } 268 272 if ((c_e==-1)||(c_e>a_e)) 269 273 { 270 274 c_e=a_e; c_j=j; 271 275 pDelete(&c_p); 272 276 c_p = nc_CreateSpoly(pCopy(strat->S[c_j]),pCopy((*h).p), currRing); 273 277 } 274 278 else if (c_e == a_e) 275 279 { 276 277 280 poly cc_pp= nc_CreateSpoly(pCopy(strat->S[j]),pCopy((*h).p), currRing); 281 if (((cc_pp==NULL)&&(c_p!=NULL)) || (pCmp(cc_pp,c_p)==-1)) 278 282 { 279 283 assume(pTotaldegree(cc_pp)<=pTotaldegree(c_p)); 280 284 c_e=a_e; c_j=j; 281 285 pDelete(&c_p); 282 286 c_p = cc_pp; 283 287 } 284 288 } 285 289 /*computes the ecart*/ … … 295 299 else 296 300 { 297 if (TEST_OPT_DEBUG) PrintS("-\n"); 301 if(TEST_OPT_DEBUG) 302 { 303 p_wrp(h->p,currRing); Print(" not divisibly by S[%d]=",j); 304 p_wrp(strat->S[j],currRing); PrintLn(); 305 } 298 306 } 299 307 j++;
Note: See TracChangeset
for help on using the changeset viewer.