Changeset 9441be in git for kernel/kutil.cc


Ignore:
Timestamp:
Feb 27, 2009, 12:21:24 PM (15 years ago)
Author:
Hans Schönemann <hannes@…>
Branches:
(u'spielwiese', '8e0ad00ce244dfd0756200662572aef8402f13d5')
Children:
bc772d8c1eb752fbedb6e31bf0a94ac364666dc1
Parents:
e371bf4ee8fd6db22ae0856bc12a6c78fc50df7e
Message:
*hannes: avoid equal LM after NF (with redSB)


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

Legend:

Unmodified
Added
Removed

  • kernel/kutil.cc

    re371bf r9441be  
    22*  Computer Algebra System SINGULAR     *
    33****************************************/
    4 /* $Id: kutil.cc,v 1.130 2009-02-27 11:05:15 Singular Exp $ */
     4/* $Id: kutil.cc,v 1.131 2009-02-27 11:21:24 Singular Exp $ */
    55/*
    66* ABSTRACT: kernel: utils for kStd
    … …  
    19901990              if(TEST_OPT_DEBUG)
    19911991              {
    1992                 Print("chain-crit-part: S[%d]=",j);
    1993                 p_wrp(strat->S[j],currRing);
    1994                 Print(" divide B[%d].lcm=",i);
    1995                 p_wrp(strat->B[i].lcm,currRing);
    1996                 PrintLn();
     1992                Print("chain-crit-part: S[%d]=",j);
     1993                p_wrp(strat->S[j],currRing);
     1994                Print(" divide B[%d].lcm=",i);
     1995                p_wrp(strat->B[i].lcm,currRing);
     1996                PrintLn();
    19971997              }
    19981998              deleteInL(strat->B,&strat->Bl,i,strat);
    … …  
    20252025              if(TEST_OPT_DEBUG)
    20262026              {
    2027                 PrintS("chain-crit-part: pCompareChainPart p=");
    2028                 p_wrp(p,currRing);
    2029                 Print(" delete L[%d]",j);
    2030                 p_wrp(strat->L[j].lcm,currRing);
    2031                 PrintLn();
     2027                PrintS("chain-crit-part: pCompareChainPart p=");
     2028                p_wrp(p,currRing);
     2029                Print(" delete L[%d]",j);
     2030                p_wrp(strat->L[j].lcm,currRing);
     2031                PrintLn();
    20322032              }
    20332033              deleteInL(strat->L,&strat->Ll,j,strat);
    … …  
    20562056              if(TEST_OPT_DEBUG)
    20572057              {
    2058                 Print("chain-crit-part: sugar B[%d].lcm=",j);
    2059                 p_wrp(strat->B[j].lcm,currRing);
    2060                 Print(" delete B[%d]",i);
    2061                 p_wrp(strat->B[i].lcm,currRing);
    2062                 PrintLn();
     2058                Print("chain-crit-part: sugar B[%d].lcm=",j);
     2059                p_wrp(strat->B[j].lcm,currRing);
     2060                Print(" delete B[%d]",i);
     2061                p_wrp(strat->B[i].lcm,currRing);
     2062                PrintLn();
    20632063              }
    20642064              deleteInL(strat->B,&strat->Bl,i,strat);
    … …  
    20692069              if(TEST_OPT_DEBUG)
    20702070              {
    2071                 Print("chain-crit-part: sugar B[%d].lcm=",i);
    2072                 p_wrp(strat->B[i].lcm,currRing);
    2073                 Print(" delete B[%d]",j);
    2074                 p_wrp(strat->B[j].lcm,currRing);
    2075                 PrintLn();
     2071                Print("chain-crit-part: sugar B[%d].lcm=",i);
     2072                p_wrp(strat->B[i].lcm,currRing);
     2073                Print(" delete B[%d]",j);
     2074                p_wrp(strat->B[j].lcm,currRing);
     2075                PrintLn();
    20762076              }
    20772077              deleteInL(strat->B,&strat->Bl,j,strat);
    … …  
    20992099              if(TEST_OPT_DEBUG)
    21002100              {
    2101                 PrintS("chain-crit-part: sugar:pCompareChainPart p=");
    2102                 p_wrp(p,currRing);
    2103                 Print(" delete L[%d]",j);
    2104                 p_wrp(strat->L[j].lcm,currRing);
    2105                 PrintLn();
     2101                PrintS("chain-crit-part: sugar:pCompareChainPart p=");
     2102                p_wrp(p,currRing);
     2103                Print(" delete L[%d]",j);
     2104                p_wrp(strat->L[j].lcm,currRing);
     2105                PrintLn();
    21062106              }
    21072107            deleteInL(strat->L,&strat->Ll,j,strat);
    … …  
    21252125              if(TEST_OPT_DEBUG)
    21262126              {
    2127                 Print("chain-crit-part: equal lcm B[%d].lcm=",j);
    2128                 p_wrp(strat->B[j].lcm,currRing);
    2129                 Print(" delete B[%d]\n",i);
     2127                Print("chain-crit-part: equal lcm B[%d].lcm=",j);
     2128                p_wrp(strat->B[j].lcm,currRing);
     2129                Print(" delete B[%d]\n",i);
    21302130              }
    21312131            strat->c3++;
    … …  
    21522152              if(TEST_OPT_DEBUG)
    21532153              {
    2154                 PrintS("chain-crit-part: pCompareChainPart p=");
    2155                 p_wrp(p,currRing);
    2156                 Print(" delete L[%d]",j);
    2157                 p_wrp(strat->L[j].lcm,currRing);
    2158                 PrintLn();
     2154                PrintS("chain-crit-part: pCompareChainPart p=");
     2155                p_wrp(p,currRing);
     2156                Print(" delete L[%d]",j);
     2157                p_wrp(strat->L[j].lcm,currRing);
     2158                PrintLn();
    21592159              }
    21602160          deleteInL(strat->L,&strat->Ll,j,strat);
    … …  
    22162216              if(TEST_OPT_DEBUG)
    22172217              {
    2218                 PrintS("chain-crit-part: divisible_by p=");
    2219                 p_wrp(p,currRing);
    2220                 Print(" delete L[%d]",l);
    2221                 p_wrp(strat->L[l].lcm,currRing);
    2222                 PrintLn();
     2218                PrintS("chain-crit-part: divisible_by p=");
     2219                p_wrp(p,currRing);
     2220                Print(" delete L[%d]",l);
     2221                p_wrp(strat->L[l].lcm,currRing);
     2222                PrintLn();
    22232223              }
    22242224              deleteInL(strat->L,&strat->Ll,l,strat);
    … …  
    22292229              if(TEST_OPT_DEBUG)
    22302230              {
    2231                 PrintS("chain-crit-part: divisible_by(2) p=");
    2232                 p_wrp(p,currRing);
    2233                 Print(" delete L[%d]",i);
    2234                 p_wrp(strat->L[i].lcm,currRing);
    2235                 PrintLn();
     2231                PrintS("chain-crit-part: divisible_by(2) p=");
     2232                p_wrp(p,currRing);
     2233                Print(" delete L[%d]",i);
     2234                p_wrp(strat->L[i].lcm,currRing);
     2235                PrintLn();
    22362236              }
    22372237              deleteInL(strat->L,&strat->Ll,i,strat);
    … …  
    58515851    int q;
    58525852    poly p;
     5853    BOOLEAN reduction_found=FALSE;
    58535854    for (l=IDELEMS(r)-1;l>=0;l--)
    58545855    {
    … …  
    58655866              r->m[l]=kNF(Q,NULL,p);
    58665867              pDelete(&p);
     5868              reduction_found=TRUE;
    58675869            }
    58685870            else
    … …  
    58715873            }
    58725874            break;
     5875          }
     5876        }
     5877      }
     5878    }
     5879    if (/*TEST_OPT_REDSB &&*/ reduction_found)
     5880    {
     5881      for (l=IDELEMS(r)-1;l>=0;l--)
     5882      {
     5883        if (r->m[l]!=NULL)
     5884        {
     5885          for(q=IDELEMS(r)-1;q>=0;q--)
     5886          {
     5887            if ((l!=q)
     5888            && (r->m[q]!=NULL)
     5889            &&(pLmDivisibleBy(r->m[l],r->m[q])))
     5890            {
     5891              pDelete(&r->m[q]);
     5892            }
    58735893          }
    58745894        }
Note: See TracChangeset for help on using the changeset viewer.