Changeset 3bae69 in git


Ignore:
Timestamp:
Dec 2, 2014, 6:41:48 PM (10 years ago)
Author:
Hans Schoenemann <hannes@…>
Branches:
(u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', '4188d308699580d975efd0f6cca8dcb41c396f70')
Children:
8fb1468d9b21997e9955d5e5f6903a1337853915
Parents:
a62aaf77249f98fad26ba201fd0f646f392e5fcc
Message:
revisiting kutil/for ringcf: cleanup
Location:
kernel/GBEngine
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • kernel/GBEngine/kInline.h

    ra62aaf r3bae69  
    533533  {
    534534    int l = GetpLength();
    535     if (use_bucket && l > 1)
     535    if (use_bucket && (l > 1))
    536536    {
    537537      poly tp = GetLmTailRing();
  • kernel/GBEngine/kutil.cc

    ra62aaf r3bae69  
    335335
    336336#ifdef HAVE_RINGS
    337     if (rField_is_Ring(currRing) && (currRing->OrdSgn == -1))
    338                   lc = p_GetCoeff(p,r);
     337    if (rField_is_Ring(r) /*&& (!rHasGlobalOrdering(r))*/)
     338      lc = p_GetCoeff(p,r);
    339339#endif
    340340
     
    359359      if (!inNF)
    360360      {
    361              number eins;
    362               #ifdef HAVE_RINGS
    363               if (rField_is_Ring(currRing) && (currRing->OrdSgn == -1))
    364                               eins = nCopy(lc);
    365                     else
    366                      #endif
    367               eins=nInit(1);
    368               if (L->p != NULL)  pSetCoeff(L->p,eins);
    369         else if (L->t_p != NULL) nDelete(&pGetCoeff(L->t_p));
    370         if (L->t_p != NULL) pSetCoeff0(L->t_p,eins);
     361        number eins;
     362#ifdef HAVE_RINGS
     363        if (rField_is_Ring(r) /*&& (!rHasGlobalOrdering(r))*/)
     364          eins = nCopy(lc);
     365        else
     366#endif
     367          eins=nInit(1);
     368        if (L->p != NULL)
     369        {
     370          pSetCoeff(L->p,eins);
     371          if (L->t_p != NULL)
     372            pSetCoeff0(L->t_p,eins);
     373        }
     374        else
     375          pSetCoeff(L->t_p,eins);
     376        /* p and t_p share the same coeff, if both are !=NULL */
     377        /* p==NULL==t_p cannot happen here */
    371378      }
    372379      L->ecart = 0;
     
    377384
    378385      if (L->t_p != NULL && pNext(L->t_p) != NULL)
    379         pNext(L->t_p) = NULL;
     386        p_Delete(&pNext(L->t_p),r);
    380387      if (L->p != NULL && pNext(L->p) != NULL)
    381388        pNext(L->p) = NULL;
     
    391398      // Note: As long as qring j forbidden if j contains integer (i.e. ground rings are
    392399      //       domains), no zerodivisor test needed  CAUTION
    393       if (rField_is_Ring(currRing) && currRing->OrdSgn == -1)
    394               if(n_DivBy(p_GetCoeff(h,r->cf),lc,r->cf) == 0)
    395                       return;
     400      if (rField_is_Ring(r) /*&&(!rHasGlobalOrdering(r)) */)
     401        if(n_DivBy(p_GetCoeff(h,r->cf),lc,r->cf) == 0)
     402          return;
    396403      #endif
    397404      if (i == r->N) break; // does divide, try next monom
     
    18291836//#if 1
    18301837#if DEBUGF5
    1831   printf("IN PAIR GENERATION - COMPARING SIGS: %d\n",sigCmp);
     1838  Print("IN PAIR GENERATION - COMPARING SIGS: %d\n",sigCmp);
    18321839  pWrite(pSigMult);
    18331840  pWrite(sSigMult);
     
    19972004    pDelete (&m2);
    19982005#if DEBUGF5
    1999     printf("SIGNATURE OF PAIR:  ");
     2006    PrintS("SIGNATURE OF PAIR:  ");
    20002007    pWrite(Lp.sig);
    20012008#endif
     
    35843591void superenterpairs (poly h,int k,int ecart,int pos,kStrategy strat, int atR)
    35853592{
    3586   #if ADIDEBUG
    3587   PrintLn();
    3588   PrintS("Enter superenterpairs");
    3589   PrintLn();
     3593#if ADIDEBUG
     3594  PrintS("\nEnter superenterpairs\n");
    35903595  int iii = strat->Ll;
    3591   #endif
     3596#endif
    35923597  assume (rField_is_Ring(currRing));
    35933598  // enter also zero divisor * poly, if this is non zero and of smaller degree
    35943599  if (!(rField_is_Domain(currRing))) enterExtendedSpoly(h, strat);
    3595    #if ADIDEBUG
     3600#if ADIDEBUG
    35963601  if(iii==strat->Ll)
    35973602  {
     3603    PrintS("\n                enterExtendedSpoly has not changed the list L.\n");
     3604  }
     3605  else
     3606  {
    35983607    PrintLn();
    3599     PrintS("                enterExtendedSpoly has not changed the list L.");
    3600     PrintLn();
    3601   }
    3602   else
    3603   {
    3604     PrintLn();
    3605     PrintS("                enterExtendedSpoly changed the list L: ");
    3606     PrintLn();
     3608    PrintS("\n                enterExtendedSpoly changed the list L:\n");
    36073609    for(iii=0;iii<=strat->Ll;iii++)
    36083610    {
    3609       PrintLn();
    3610       PrintS("                L[");printf("%d",iii);PrintS("]:");PrintLn();
     3611      Print("\n                L[%d]:\n",iii);
    36113612      PrintS("                     ");p_Write(strat->L[iii].p1,strat->tailRing);
    36123613      PrintS("                     ");p_Write(strat->L[iii].p2,strat->tailRing);
     
    36153616  }
    36163617  iii = strat->Ll;
    3617   #endif
     3618#endif
    36183619  initenterpairs(h, k, ecart, 0, strat, atR);
    3619   #if ADIDEBUG
     3620#if ADIDEBUG
    36203621  if(iii==strat->Ll)
    36213622  {
    3622     PrintLn();
    3623     PrintS("                initenterpairs has not changed the list L.");
    3624     PrintLn();
     3623    PrintS("\n                initenterpairs has not changed the list L.\n");
    36253624  }
    36263625  else
    36273626  {
    3628     PrintLn();
    3629     PrintS("                initenterpairs changed the list L: ");
    3630     PrintLn();
     3627    PrintS("\n                initenterpairs changed the list L:\n");
    36313628    for(iii=0;iii<=strat->Ll;iii++)
    36323629    {
    3633       PrintLn();
    3634       PrintS("                L[");printf("%d",iii);PrintS("]:");PrintLn();
     3630      Print("\n                L[%d]:\n",iii);
    36353631      PrintS("                     ");p_Write(strat->L[iii].p1,strat->tailRing);
    36363632      PrintS("                     ");p_Write(strat->L[iii].p2,strat->tailRing);
     
    36393635  }
    36403636  iii = strat->Ll;
    3641   #endif
     3637#endif
    36423638  initenterstrongPairs(h, k, ecart, 0, strat, atR);
    3643   #if ADIDEBUG
     3639#if ADIDEBUG
    36443640  if(iii==strat->Ll)
    36453641  {
    3646     PrintLn();
    3647     PrintS("                initenterstrongPairs has not changed the list L.");
    3648     PrintLn();
     3642    PrintS("\n                initenterstrongPairs has not changed the list L.\n");
    36493643  }
    36503644  else
    36513645  {
    3652     PrintLn();
    3653     PrintS("                initenterstrongPairs changed the list L: ");
    3654     PrintLn();
     3646    PrintS("\n                initenterstrongPairs changed the list L:\n");
    36553647    for(iii=0;iii<=strat->Ll;iii++)
    36563648    {
    3657       PrintLn();
    3658       PrintS("                L[");printf("%d",iii);PrintS("]:");PrintLn();
     3649      Print("\n                L[%d]:\n",iii);
    36593650      PrintS("                     ");p_Write(strat->L[iii].p1,strat->tailRing);
    36603651      PrintS("                     ");p_Write(strat->L[iii].p2,strat->tailRing);
     
    36623653    }
    36633654  }
    3664   PrintLn();
    3665   PrintS("End of superenterpairs");
    3666   PrintLn();
    3667   #endif
     3655  PrintS("\nEnd of superenterpairs\n");
     3656#endif
    36683657  clearSbatch(h, k, pos, strat);
    36693658}
     
    50935082//#if 1
    50945083#ifdef DEBUGF5
    5095       printf("DELETE!\n");
     5084      PrintS("DELETE!\n");
    50965085#endif
    50975086      //printf("- T -\n\n");
     
    51345123    {
    51355124#ifdef F5DEBUG
    5136       printf("COMP %d/%d - MIN %d - MAX %d - SYZL %ld\n",comp,strat->currIdx,min,max,strat->syzl);
     5125      Print("COMP %d/%d - MIN %d - MAX %d - SYZL %ld\n",comp,strat->currIdx,min,max,strat->syzl);
    51375126      Print("checking with: %d --  ",k);
    51385127      pWrite(pHead(strat->syz[k]));
     
    51535142//#if 1
    51545143#ifdef DEBUGF5
    5155   printf("rewritten criterion checks:  ");
     5144  PrintS("rewritten criterion checks:  ");
    51565145  pWrite(sig);
    51575146#endif
     
    51685157//#if 1
    51695158#ifdef DEBUGF5
    5170       printf("DELETE!\n");
     5159      PrintS("DELETE!\n");
    51715160#endif
    51725161      return TRUE;
     
    55735562        else
    55745563        {
    5575           if ((Ln.t_p!=NULL)&&(Ln.p==NULL))
    5576             Ln.GetP();
     5564          if ((Ln.t_p!=NULL)&&(Ln.p==NULL))
     5565            Ln.GetP();
    55775566          if (Ln.p!=NULL)
    5578           {
    5579             Ln.p=pAdd(Ln.p,mm);
    5580             if (Ln.t_p!=NULL)
    5581             {
    5582               pNext(Ln.t_p)=NULL;
    5583               p_LmDelete(Ln.t_p,strat->tailRing);
    5584             }
    5585           }
     5567          {
     5568            Ln.p=pAdd(Ln.p,mm);
     5569            if (Ln.t_p!=NULL)
     5570            {
     5571              pNext(Ln.t_p)=NULL;
     5572              p_LmDelete(Ln.t_p,strat->tailRing);
     5573            }
     5574          }
    55865575        }
    55875576      }
     
    61216110    strat->syzidxmax  = comp;
    61226111#if defined(DEBUGF5) || defined(DEBUGF51)
    6123     printf("------------- GENERATING SYZ RULES NEW ---------------\n");
     6112    PrintS("------------- GENERATING SYZ RULES NEW ---------------\n");
    61246113#endif
    61256114    i = 1;
     
    62156204//#if 1
    62166205#ifdef DEBUGF5
    6217     Print("Principal syzygies:\n");
    6218     printf("syzl   %d\n",strat->syzl);
    6219     printf("syzmax %d\n",strat->syzmax);
    6220     printf("ps     %d\n",ps);
    6221     Print("--------------------------------\n");
     6206    PrintS("Principal syzygies:\n");
     6207    Print("syzl   %d\n",strat->syzl);
     6208    Print("syzmax %d\n",strat->syzmax);
     6209    Print("ps     %d\n",ps);
     6210    PrintS("--------------------------------\n");
    62226211    for(i=0;i<=strat->syzl-1;i++)
    62236212    {
    6224       printf("%d - ",i);
     6213      Print("%d - ",i);
    62256214      pWrite(strat->syz[i]);
    62266215    }
    62276216    for(i=0;i<strat->currIdx;i++)
    62286217    {
    6229       printf("%d - %d\n",i,strat->syzIdx[i]);
    6230     }
    6231     Print("--------------------------------\n");
    6232 #endif
    6233 
     6218      Print("%d - %d\n",i,strat->syzIdx[i]);
     6219    }
     6220    PrintS("--------------------------------\n");
     6221#endif
    62346222  }
    62356223}
     
    72957283//#if 1
    72967284#ifdef DEBUGF5
    7297     Print("--- Syzygies ---\n");
    7298     printf("syzl   %d\n",strat->syzl);
    7299     printf("syzmax %d\n",strat->syzmax);
    7300     Print("--------------------------------\n");
     7285    PrintS("--- Syzygies ---\n");
     7286    Print("syzl   %d\n",strat->syzl);
     7287    Print("syzmax %d\n",strat->syzmax);
     7288    PrintS("--------------------------------\n");
    73017289    for(i=0;i<=strat->syzl-1;i++)
    73027290    {
    7303       printf("%d - ",i);
     7291      Print("%d - ",i);
    73047292      pWrite(strat->syz[i]);
    73057293    }
    7306     Print("--------------------------------\n");
     7294    PrintS("--------------------------------\n");
    73077295#endif
    73087296}
Note: See TracChangeset for help on using the changeset viewer.