Changeset fb87b82 in git for kernel/GBEngine/kutil.cc


Ignore:
Timestamp:
Oct 19, 2015, 7:00:58 PM (8 years ago)
Author:
Adi Popescu <adi_popescum@…>
Branches:
(u'spielwiese', 'a719bcf0b8dbc648b128303a49777a094b57592c')
Children:
f6db0d20f2a2a1e61303f9831ef6201e89faad2f
Parents:
a07db8714e0495d86b2f4758df4475f17535d842
Message:
partially fixed chainCrit over Rings
File:
1 edited

Legend:

Unmodified
Added
Removed
  • kernel/GBEngine/kutil.cc

    ra07db8 rfb87b82  
    13931393      }
    13941394      strat->initEcartPair(&Lp,strat->S[i],p,strat->ecartS[i],ecart);
    1395       l = strat->posInL(strat->L,strat->Ll,&Lp,strat);
    1396       enterL(&strat->L,&strat->Ll,&strat->Lmax,Lp,l);
     1395      l = strat->posInL(strat->B,strat->Bl,&Lp,strat);
     1396      enterL(&strat->B,&strat->Bl,&strat->Bmax,Lp,l);
    13971397    }
    13981398  #else
     
    14491449              pWrite(strat->B[j].p1);
    14501450              pWrite(strat->B[j].p2);
     1451              pWrite(strat->B[j].lcm);
    14511452              printf("\nh - neue Paar\n");
    14521453              pWrite(h.p);
     
    14551456              pWrite(h.lcm);
    14561457              #endif
    1457         if ((n_DivBy(strat->B[j].lcm->coef, h.lcm->coef, currRing->cf)) && ((strat->fromQ==NULL) || (isFromQ==0) || (strat->fromQ[i]==0)))
    1458         {
     1458        if ((n_DivBy(h.lcm->coef, strat->B[j].lcm->coef,  currRing->cf)) && ((strat->fromQ==NULL) || (isFromQ==0) || (strat->fromQ[i]==0)))
     1459        {
     1460          #if ADIDEBUG
     1461          printf("\nGelöscht neue h\n");
     1462          #endif
    14591463          strat->c3++;
    14601464          pLmDelete(h.lcm);
     
    14711475              pWrite(strat->B[j].p1);
    14721476              pWrite(strat->B[j].p2);
     1477              pWrite(strat->B[j].lcm);
    14731478              printf("\nh - neue Paar\n");
    14741479              pWrite(h.p);
     
    14771482              pWrite(h.lcm);
    14781483              #endif
    1479           if(n_DivBy( h.lcm->coef, strat->B[j].lcm->coef, currRing->cf))
    1480           {
     1484          //Ist schon im coeffCompare
     1485          //if(n_DivBy( h.lcm->coef, strat->B[j].lcm->coef, currRing->cf))
     1486          //{
    14811487          #if ADIDEBUG
    14821488                        printf("\nGelöscht: B[j]\n");
     
    14841490            deleteInL(strat->B,&strat->Bl,j,strat);
    14851491            strat->c3++;
    1486           }
     1492          //}
    14871493      }
    14881494    }
     
    14971503              pWrite(strat->B[j].p1);
    14981504              pWrite(strat->B[j].p2);
     1505              pWrite(strat->B[j].lcm);
    14991506              printf("\nh - neue Paar\n");
    15001507              pWrite(h.p);
     
    15301537              pWrite(h.lcm);
    15311538              #endif
    1532         if(n_DivBy( h.lcm->coef, strat->B[j].lcm->coef,currRing->cf))
    1533         {
     1539        //if(n_DivBy( h.lcm->coef, strat->B[j].lcm->coef,currRing->cf))
     1540        //{
    15341541          #if ADIDEBUG
    15351542          printf("\nGelöscht B[j]\n");
     
    15371544          deleteInL(strat->B,&strat->Bl,j,strat);
    15381545          strat->c3++;
    1539         }
     1546        //}
    15401547      }
    15411548    }
     
    15691576        pDelete(&h.lcm);
    15701577      h.Clear();
     1578      if (strat->pairtest==NULL) initPairtest(strat);
     1579      strat->pairtest[i] = TRUE;
     1580      strat->pairtest[strat->sl+1] = TRUE;
    15711581      return;
    15721582    }
     
    16001610  h.i_r = -1;
    16011611  if(h.p == NULL)
     1612  { 
     1613    if (strat->pairtest==NULL) initPairtest(strat);
     1614    strat->pairtest[i] = TRUE;
     1615    strat->pairtest[strat->sl+1] = TRUE;
    16021616    return;
     1617  }
    16031618  h.tailRing = strat->tailRing;
    16041619  int posx;
    1605   h.pCleardenom();
    1606   if(h.p == NULL)
    1607     return;
    1608   pSetm(h.p);
     1620  //h.pCleardenom();
     1621  //pSetm(h.p);
    16091622  #if ADIDEBUG
    16101623  printf("\nThis is afterwards:\n");
     
    17491762      h.i_r2 = -1;
    17501763    }
    1751     if (strat->Ll==-1)
     1764    if (strat->Bl==-1)
    17521765      posx =0;
    17531766    else
    1754       posx = strat->posInL(strat->L,strat->Ll,&h,strat);
    1755     enterL(&strat->L,&strat->Ll,&strat->Lmax,h,posx);
     1767      posx = strat->posInL(strat->B,strat->Bl,&h,strat);
     1768    enterL(&strat->B,&strat->Bl,&strat->Bmax,h,posx);
    17561769  }
    17571770  else
     
    34343447                wrp(strat->S[j]);
    34353448                PrintS("  strat->B[i].lcm:");
     3449                wrp(strat->B[i].lcm);PrintLn();
     3450                pWrite(strat->B[i].p);
     3451                pWrite(strat->B[i].p1);
     3452                pWrite(strat->B[i].p2);
    34363453                wrp(strat->B[i].lcm);
    34373454                PrintLn();
Note: See TracChangeset for help on using the changeset viewer.