Changeset a13cae in git


Ignore:
Timestamp:
Jul 27, 2008, 12:21:59 PM (16 years ago)
Author:
Hans Schönemann <hannes@…>
Branches:
(u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', '45e267b9942dec9429fe249ce3e5e44ab80a6a3a')
Children:
37d43ec78deb09e94d2b41b08283e6fe466687d7
Parents:
e4df0ba39c30b0d3f0f545e4cf34a78b4b74fbea
Message:
*hannes: OPT_DEBUG fuor ratGB


git-svn-id: file:///usr/local/Singular/svn/trunk@10924 2c84dea3-7e68-4137-9b89-c4e89433aadc
Location:
kernel
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • kernel/gr_kstd2.cc

    re4df0ba ra13cae  
    22*  Computer Algebra System SINGULAR     *
    33****************************************/
    4 /* $Id: gr_kstd2.cc,v 1.15 2008-07-26 18:11:42 Singular Exp $ */
     4/* $Id: gr_kstd2.cc,v 1.16 2008-07-27 10:21:59 Singular Exp $ */
    55/*
    66*  ABSTRACT -  Kernel: noncomm. alg. of Buchberger
     
    229229        }
    230230        //poly hh = nc_CreateSpoly(strat->S[c_j],(*h).p, currRing);
    231         poly hh=c_p; c_p=NULL;
    232         pDelete(&((*h).p));
    233         (*h).p=hh;
     231        poly hh=c_p; c_p=NULL;
     232        pDelete(&((*h).p));
     233        (*h).p=hh;
    234234        if (!TEST_OPT_INTSTRATEGY)
    235235        {
     
    254254        /*- try to reduce the s-polynomial again -*/
    255255        pass++;
    256         j=0;
     256        j=0;
    257257      }
    258258      else
     
    265265    {
    266266      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      }
    268272      if ((c_e==-1)||(c_e>a_e))
    269273      {
    270274        c_e=a_e; c_j=j;
    271         pDelete(&c_p);
     275        pDelete(&c_p);
    272276        c_p = nc_CreateSpoly(pCopy(strat->S[c_j]),pCopy((*h).p), currRing);
    273277      }
    274278      else if (c_e == a_e)
    275279      {
    276         poly cc_pp= nc_CreateSpoly(pCopy(strat->S[j]),pCopy((*h).p), currRing);
    277         if (((cc_pp==NULL)&&(c_p!=NULL)) || (pCmp(cc_pp,c_p)==-1))
     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))
    278282         {
    279            assume(pTotaldegree(cc_pp)<=pTotaldegree(c_p));
     283           assume(pTotaldegree(cc_pp)<=pTotaldegree(c_p));
    280284           c_e=a_e; c_j=j;
    281            pDelete(&c_p);
     285           pDelete(&c_p);
    282286           c_p = cc_pp;
    283         }
     287        }
    284288      }
    285289      /*computes the ecart*/
     
    295299    else
    296300    {
    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      }
    298306    }
    299307    j++;
  • kernel/kutil.cc

    re4df0ba ra13cae  
    22*  Computer Algebra System SINGULAR     *
    33****************************************/
    4 /* $Id: kutil.cc,v 1.104 2008-07-25 15:40:59 Singular Exp $ */
     4/* $Id: kutil.cc,v 1.105 2008-07-27 10:21:59 Singular Exp $ */
    55/*
    66* ABSTRACT: kernel: utils for kStd
     
    19801980               currRing->real_var_start,currRing->real_var_end))
    19811981            {
     1982              if(TEST_OPT_DEBUG)
     1983              {
     1984                 Print("chain-crit-part: S[%d]=",j);
     1985                 p_wrp(strat->S[j],currRing);
     1986                 Print(" divide B[%d].lcm=",i);
     1987                 p_wrp(strat->B[i].lcm,currRing);
     1988                 PrintLn();
     1989              }
    19821990              deleteInL(strat->B,&strat->Bl,i,strat);
    19831991              strat->c3++;
     
    20072015          if (strat->L[j].p == strat->tail)
    20082016          {
     2017              if(TEST_OPT_DEBUG)
     2018              {
     2019                 PrintS("chain-crit-part: pCompareChainPart p=");
     2020                 p_wrp(p,currRing);
     2021                 Print(" delete L[%d]",j);
     2022                 p_wrp(strat->L[j].lcm,currRing);
     2023                 PrintLn();
     2024              }
    20092025              deleteInL(strat->L,&strat->Ll,j,strat);
    20102026              strat->c3++;
     
    20302046            if (sugarDivisibleBy(strat->B[j].ecart,strat->B[i].ecart))
    20312047            {
     2048              if(TEST_OPT_DEBUG)
     2049              {
     2050                 Print("chain-crit-part: sugar B[%d].lcm=",j);
     2051                 p_wrp(strat->B[j].lcm,currRing);
     2052                 Print(" delete B[%d]",i);
     2053                 p_wrp(strat->B[i].lcm,currRing);
     2054                 PrintLn();
     2055              }
    20322056              deleteInL(strat->B,&strat->Bl,i,strat);
    20332057              j--;
     
    20352059            else
    20362060            {
     2061              if(TEST_OPT_DEBUG)
     2062              {
     2063                 Print("chain-crit-part: sugar B[%d].lcm=",i);
     2064                 p_wrp(strat->B[i].lcm,currRing);
     2065                 Print(" delete B[%d]",j);
     2066                 p_wrp(strat->B[j].lcm,currRing);
     2067                 PrintLn();
     2068              }
    20372069              deleteInL(strat->B,&strat->Bl,j,strat);
    20382070              break;
     
    20572089          if ((pNext(strat->L[j].p) == strat->tail)||(pOrdSgn==1))
    20582090          {
     2091              if(TEST_OPT_DEBUG)
     2092              {
     2093                 PrintS("chain-crit-part: sugar:pCompareChainPart p=");
     2094                 p_wrp(p,currRing);
     2095                 Print(" delete L[%d]",j);
     2096                 p_wrp(strat->L[j].lcm,currRing);
     2097                 PrintLn();
     2098              }
    20592099            deleteInL(strat->L,&strat->Ll,j,strat);
    20602100            strat->c3++;
     
    20752115          if (pLmEqual(strat->B[j].lcm,strat->B[i].lcm))
    20762116          {
     2117              if(TEST_OPT_DEBUG)
     2118              {
     2119                 Print("chain-crit-part: equal lcm B[%d].lcm=",j);
     2120                 p_wrp(strat->B[j].lcm,currRing);
     2121                 Print(" delete B[%d]\n",i);
     2122              }
    20772123            strat->c3++;
    20782124            deleteInL(strat->B,&strat->Bl,i,strat);
     
    21052151        if ((pNext(strat->L[j].p) == strat->tail)||(pOrdSgn==1))
    21062152        {
     2153              if(TEST_OPT_DEBUG)
     2154              {
     2155                 PrintS("chain-crit-part: pCompareChainPart p=");
     2156                 p_wrp(p,currRing);
     2157                 Print(" delete L[%d]",j);
     2158                 p_wrp(strat->L[j].lcm,currRing);
     2159                 PrintLn();
     2160              }
    21072161          deleteInL(strat->L,&strat->Ll,j,strat);
    21082162          strat->c3++;
     
    21672221              *so we mark it with "tail"
    21682222              */
     2223              if(TEST_OPT_DEBUG)
     2224              {
     2225                 PrintS("chain-crit-part: divisible_by p=");
     2226                 p_wrp(p,currRing);
     2227                 Print(" delete L[%d]",l);
     2228                 p_wrp(strat->L[l].lcm,currRing);
     2229                 PrintLn();
     2230              }
    21692231              deleteInL(strat->L,&strat->Ll,l,strat);
    21702232              i--;
     
    21722234            else
    21732235            {
     2236              if(TEST_OPT_DEBUG)
     2237              {
     2238                 PrintS("chain-crit-part: divisible_by(2) p=");
     2239                 p_wrp(p,currRing);
     2240                 Print(" delete L[%d]",i);
     2241                 p_wrp(strat->L[i].lcm,currRing);
     2242                 PrintLn();
     2243              }
    21742244              deleteInL(strat->L,&strat->Ll,i,strat);
    21752245            }
Note: See TracChangeset for help on using the changeset viewer.