Changeset c0e2c36 in git for kernel/gr_kstd2.cc


Ignore:
Timestamp:
Feb 26, 2009, 4:57:01 PM (14 years ago)
Author:
Hans Schönemann <hannes@…>
Branches:
(u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', 'a800fe4b3e9d37a38c5a10cc0ae9dfa0c15a4ee6')
Children:
3ec38a567b6025adb57223e0a734c42f135941f9
Parents:
84d8b09321d52723c68edd2e14a97c84d72a8022
Message:
*hannes: pCleardenom changes


git-svn-id: file:///usr/local/Singular/svn/trunk@11482 2c84dea3-7e68-4137-9b89-c4e89433aadc
File:
1 edited

Legend:

Unmodified
Added
Removed

  • kernel/gr_kstd2.cc

    r84d8b09 rc0e2c36  
    22*  Computer Algebra System SINGULAR     *
    33****************************************/
    4 /* $Id: gr_kstd2.cc,v 1.31 2009-02-26 12:17:38 levandov Exp $ */
     4/* $Id: gr_kstd2.cc,v 1.32 2009-02-26 15:57:01 Singular Exp $ */
    55/*
    66*  ABSTRACT -  Kernel: noncomm. alg. of Buchberger
    … …  
    103103      {
    104104        if (rField_is_Zp_a()) pContent(h->p);
    105         else pCleardenom(h->p);// also does a pContent
     105        else h->pCleardenom();// also does a pContent
    106106      }
    107107      /*computes the ecart*/
    … …  
    231231  {
    232232    if (rField_is_Zp_a()) pContent(h->p);
    233     else pCleardenom(h->p);// also does a pContentRat
     233    else h->pCleardenom();// also does a pContentRat
    234234  }
    235235  loop
    … …  
    264264        {
    265265          if (rField_is_Zp_a()) pContent(h->p);
    266           else pCleardenom(h->p);// also does a pContent
     266          else h->pCleardenom();// also does a pContent
    267267        }
    268268
    … …  
    463463        {
    464464          if (rField_is_Zp_a()) pContent(h->p);
    465           else pCleardenom(h->p);// also does a pContent
     465          else h->pCleardenom();// also does a pContent
    466466        }
    467467        if (strat->syzComp!=0)
    … …  
    487487        {
    488488          if (rField_is_Zp_a()) pContent(h->p);
    489           else pCleardenom(h->p);// also does a pContent
     489          else h->pCleardenom();// also does a pContent
    490490        }
    491491/*
    … …  
    607607        {
    608608          if (rField_is_Zp_a()) pContent(h->p);
    609           else pCleardenom(h->p);// also does a pContent
     609          else h->pCleardenom();// also does a pContent
    610610        }
    611611        enterT((*h),strat);
    … …  
    724724      {
    725725        //pContent(h->p);
    726         pCleardenom(h->p);// also does a pContent
     726        h->pCleardenom();// also does a pContent
    727727      }
    728728      /* compute the ecart */
    … …  
    790790        {
    791791          //pContent(h->p);
    792           pCleardenom(h->p);// also does a pContent
     792          h->pCleardenom();// also does a pContent
    793793        }
    794794        enterT((*h),strat);
    … …  
    945945        {
    946946          //pContent(h->p);
    947           pCleardenom(h->p);// also does a pContent
     947          h->pCleardenom();// also does a pContent
    948948        }
    949949        enterT((*h),strat);
    … …  
    10591059#endif
    10601060
    1061 
    1062 
    10631061  assume(pOrdSgn != -1); // no mora!!! it terminates only for global ordering!!! (?)
    10641062
    … …  
    11771175      }
    11781176#endif
    1179 
    11801177    }
    11811178    if (strat->P.p != NULL)
    11821179    {
    1183           /* statistic */
    1184           if (TEST_OPT_PROT)
     1180      if (TEST_OPT_PROT)
     1181      {
     1182        PrintS("s\n");
     1183      }
     1184      /* enter P.p into s and L */
     1185      {
     1186        strat->P.sev=0;
     1187        int pos=posInS(strat,strat->sl,strat->P.p, strat->P.ecart);
     1188        {
     1189          if (TEST_OPT_INTSTRATEGY)
    11851190          {
    1186             PrintS("s\n");
     1191            if ((strat->syzComp==0)||(!strat->homog))
     1192            {
     1193              #ifdef HAVE_RATGRING
     1194              if(!rIsRatGRing(currRing))
     1195              #endif
     1196                strat->P.p = redtailBba(strat->P.p,pos-1,strat);
     1197            }
     1198
     1199            strat->P.p=pCleardenom(strat->P.p);
    11871200          }
    1188           /* enter P.p into s and L */
     1201          else
    11891202          {
    1190             strat->P.sev=0;
    1191             int pos=posInS(strat,strat->sl,strat->P.p, strat->P.ecart);
     1203            pNorm(strat->P.p);
     1204            if ((strat->syzComp==0)||(!strat->homog))
    11921205            {
    1193               if (TEST_OPT_INTSTRATEGY)
    1194               {
    1195                 if ((strat->syzComp==0)||(!strat->homog))
    1196                 {
    1197                   strat->P.p = redtailBba(strat->P.p,pos-1,strat);
    1198                 }
    1199                 pCleardenom(strat->P.p);
    1200               }
    1201               else
    1202               {
    1203                 pNorm(strat->P.p);
    1204                 if ((strat->syzComp==0)||(!strat->homog))
    1205                 {
    1206                   strat->P.p = redtailBba(strat->P.p,pos-1,strat);
    1207                 }
    1208               }
    1209               // PLURAL debug
    1210               /* should be used only internally!!! */
    1211 
    1212               //pWrite(strat->P.p);
    1213 
    1214               if (TEST_OPT_DEBUG)
    1215               {
    1216                 PrintS("new s:"); wrp(strat->P.p);
    1217                 PrintLn();
     1206              strat->P.p = redtailBba(strat->P.p,pos-1,strat);
     1207            }
     1208          }
     1209          if (TEST_OPT_DEBUG)
     1210          {
     1211            PrintS("new s:"); wrp(strat->P.p);
     1212            PrintLn();
    12181213#if MYTEST
    1219                 Print("s: "); pWrite(strat->P.p);
    1220 #endif
    1221 
    1222               }
    1223               // kTest(strat);
    1224               //
    1225               enterpairs(strat->P.p,strat->sl,strat->P.ecart,pos,strat);
    1226 
    1227               if (strat->sl==-1) pos=0;
    1228               else pos=posInS(strat,strat->sl,strat->P.p,strat->P.ecart);
    1229 
    1230               strat->enterS(strat->P,pos,strat,-1);
    1231             }
    1232 //            if (hilb!=NULL) khCheck(Q,w,hilb,hilbeledeg,hilbcount,strat);
     1214            Print("s: "); pWrite(strat->P.p);
     1215#endif
     1216
    12331217          }
    1234           if (strat->P.lcm!=NULL) pLmFree(strat->P.lcm);
     1218          // kTest(strat);
     1219          //
     1220          enterpairs(strat->P.p,strat->sl,strat->P.ecart,pos,strat);
     1221
     1222          if (strat->sl==-1) pos=0;
     1223          else pos=posInS(strat,strat->sl,strat->P.p,strat->P.ecart);
     1224
     1225          strat->enterS(strat->P,pos,strat,-1);
     1226        }
     1227//      if (hilb!=NULL) khCheck(Q,w,hilb,hilbeledeg,hilbcount,strat);
     1228      }
     1229      if (strat->P.lcm!=NULL) pLmFree(strat->P.lcm);
    12351230      if (strat->sl>srmax) srmax = strat->sl;
    12361231    }
Note: See TracChangeset for help on using the changeset viewer.