- Timestamp:
- Jul 17, 2009, 11:51:42 AM (15 years ago)
- Branches:
- (u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', 'b4f17ed1d25f93d46dbe29e4b499baecc2fd51bb')
- Children:
- 3c7278bb58cb6b216292adf91871509a8d8455da
- Parents:
- 9446f84181f539db0f2689bb27e177084073a62f
- Location:
- kernel
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
kernel/kstd1.cc
r9446f84 rd95cee 2 2 * Computer Algebra System SINGULAR * 3 3 ****************************************/ 4 /* $Id: kstd1.cc,v 1.5 3 2009-07-10 15:13:56Singular Exp $ */4 /* $Id: kstd1.cc,v 1.54 2009-07-17 09:51:42 Singular Exp $ */ 5 5 /* 6 6 * ABSTRACT: … … 323 323 if (h->pTotalDeg()+h->ecart >= strat->tailRing->bitmask) 324 324 { 325 325 strat->overflow=TRUE; 326 326 //Print("OVERFLOW in redEcart d=%ld, max=%ld",d,strat->tailRing->bitmask); 327 328 329 330 327 h->GetP(); 328 at = strat->posInL(strat->L,strat->Ll,h,strat); 329 enterL(&strat->L,&strat->Ll,&strat->Lmax,*h,at); 330 h->Clear(); 331 331 return -1; 332 332 } … … 454 454 if (h->pTotalDeg()+h->ecart >= strat->tailRing->bitmask) 455 455 { 456 456 strat->overflow=TRUE; 457 457 //Print("OVERFLOW in redFirst d=%ld, max=%ld",d,strat->tailRing->bitmask); 458 459 460 461 458 h->GetP(); 459 at = strat->posInL(strat->L,strat->Ll,h,strat); 460 enterL(&strat->L,&strat->Ll,&strat->Lmax,*h,at); 461 h->Clear(); 462 462 return -1; 463 463 } … … 1665 1665 if (i<=kModW->length()) 1666 1666 return o+(*kModW)[i-1]; 1667 return o; 1667 return o; 1668 1668 } 1669 1669 long kHomModDeg(poly p, ring r) … … 2273 2273 if (strat->Ll > lrmax) lrmax =strat->Ll;/*stat.*/ 2274 2274 #ifdef KDEBUG 2275 #ifdef HAVE_RINGS2276 if (TEST_OPT_DEBUG) PrintS("--- next step ---\n");2277 #endif2278 2275 if (TEST_OPT_DEBUG) messageSets(strat); 2279 2276 #endif … … 2354 2351 if (strat->sl>srmax) srmax = strat->sl; 2355 2352 if (pos<strat->sl) 2356 need_retry++; 2353 { 2354 // move all "larger" elements fromS to L 2355 // remove them from T 2356 int ii=pos+1; 2357 for(;ii<=strat->sl;ii++) 2358 { 2359 LObject h; 2360 memset(&h,0,sizeof(h)); 2361 h.tailRing=strat->tailRing; 2362 h.p=strat->S[ii]; strat->S[ii]=NULL; 2363 strat->initEcart(&h); 2364 h.sev=strat->sevS[ii]; 2365 int jj=strat->tl; 2366 while (jj>=0) 2367 { 2368 if (strat->T[jj].p==h.p) 2369 { 2370 strat->T[jj].p=NULL; 2371 if (jj<strat->tl) 2372 { 2373 memmove(&(strat->T[jj]),&(strat->T[jj+1]), 2374 (strat->tl-jj)*sizeof(strat->T[jj])); 2375 memmove(&(strat->sevT[jj]),&(strat->sevT[jj+1]), 2376 (strat->tl-jj)*sizeof(strat->sevT[jj])); 2377 } 2378 strat->tl--; 2379 break; 2380 } 2381 jj--; 2382 } 2383 int lpos=strat->posInL(strat->L,strat->Ll,&h,strat); 2384 enterL(&strat->L,&strat->Ll,&strat->Lmax,h,lpos); 2385 } 2386 strat->sl=pos; 2387 } 2357 2388 } 2358 2389 … … 2360 2391 memset(&(strat->P), 0, sizeof(strat->P)); 2361 2392 #endif 2362 kTest_TS(strat);2393 //kTest_TS(strat);: i_r out of sync in kInterRedBba, but not used! 2363 2394 } 2364 2395 #ifdef KDEBUG … … 2413 2444 return kInterRedOld(F,Q); 2414 2445 2446 //return kInterRedOld(F,Q); 2447 2415 2448 BITSET save=test; 2416 2449 //test|=Sy_bit(OPT_NOT_SUGAR); -
kernel/kstd2.cc
r9446f84 rd95cee 2 2 * Computer Algebra System SINGULAR * 3 3 ****************************************/ 4 /* $Id: kstd2.cc,v 1.9 7 2009-07-13 16:36:35Singular Exp $ */4 /* $Id: kstd2.cc,v 1.98 2009-07-17 09:51:42 Singular Exp $ */ 5 5 /* 6 6 * ABSTRACT - Kernel: alg. of Buchberger … … 1000 1000 #ifdef KDEBUG 1001 1001 loop_count++; 1002 #ifdef HAVE_RINGS1003 if (TEST_OPT_DEBUG) PrintS("--- next step ---\n");1004 #endif1005 1002 if (TEST_OPT_DEBUG) messageSets(strat); 1006 1003 #endif
Note: See TracChangeset
for help on using the changeset viewer.