Changeset c96539c in git


Ignore:
Timestamp:
Jan 28, 2014, 4:01:12 PM (9 years ago)
Author:
Adi Popescu <adi_popescum@…>
Branches:
(u'spielwiese', 'f6c3dc58b0df4bd712574325fe76d0626174ad97')
Children:
c385523e4771455b4fb306117540e4b11bbc011d
Parents:
4e4e24bd5220c690dc16529b1fb431b5ef67482c
Message:
Bug #540 was fixed by deactivating OPT_SB_1 (test(29)) when the coefficient ring a ring is
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • Singular/extra.cc

    r4e4e24b rc96539c  
    24942494        }
    24952495        else
     2496  /*==================== setsyzcomp ==================================*/
     2497      if(strcmp(sys_cmd,"setsyzcomp")==0)
     2498      {
     2499     
     2500      if ((h!=NULL) && (h->Typ()==INT_CMD))
     2501         {
     2502           int k = (int)(long)h->Data();
     2503           if ( currRing->order[0] == ringorder_s )
     2504           {
     2505                rSetSyzComp(k, currRing);
     2506           }
     2507          }
     2508     
     2509      }
    24962510  /*==================== ring debug ==================================*/
    24972511        if(strcmp(sys_cmd,"r")==0)

  • kernel/ideals.cc

    r4e4e24b rc96539c  
    12141214  int k2 = id_RankFreeModule(h2,currRing);
    12151215  tHomog   hom=isNotHomog;
    1216 
    12171216  k=si_max(k1,k2);
    12181217  if (k==0)
    12191218    k = 1;
    12201219  if ((k2==0) && (k>1)) *addOnlyOne = FALSE;
    1221 
     1220*addOnlyOne = FALSE;
    12221221  intvec * weights;
    12231222  hom = (tHomog)idHomModule(h1,currQuotient,&weights);
    … …  
    13041303    }
    13051304    h4->m[IDELEMS(h4)-1] = p;
     1305    #ifdef HAVE_RINGS
     1306    if(!rField_is_Ring)
     1307    #endif
    13061308    si_opt_1 |= Sy_bit(OPT_SB_1);
    13071309  }
     1310  idPrint(h4);
    13081311  idDelete(&temph1);
    13091312  return h4;

  • kernel/kstd1.cc

    r4e4e24b rc96539c  
    21342134    strat->syzComp = syzComp;
    21352135  if (TEST_OPT_SB_1)
     2136    #ifdef HAVE_RINGS
     2137    if(!rField_is_Ring)
     2138    #endif
    21362139    strat->newIdeal = newIdeal;
    21372140  if (rField_has_simple_inverse(currRing))
    … …  
    22782281    strat->syzComp = syzComp;
    22792282  if (TEST_OPT_SB_1)
     2283    #ifdef HAVE_RINGS
     2284    if(!rField_is_Ring)
     2285    #endif
    22802286    strat->newIdeal = newIdeal;
    22812287  if (rField_has_simple_inverse(currRing))
    … …  
    24022408    strat->syzComp = syzComp;
    24032409  if (TEST_OPT_SB_1)
     2410    #ifdef HAVE_RINGS
     2411    if(!rField_is_Ring)
     2412    #endif
    24042413    strat->newIdeal = newIdeal;
    24052414  if (rField_has_simple_inverse(currRing))

  • kernel/kstd2.cc

    r4e4e24b rc96539c  
    344344    j = kFindDivisibleByInT(strat->T, strat->sevT, strat->tl, h);
    345345    if (j < 0) return 1;
    346 
    347346    ksReducePoly(h, &(strat->T[j]), NULL, NULL, strat); // with debug output
     347    //pWrite(strat->T[j].p);rWrite(currRing);getchar();
    348348
    349349    if (h->GetLmTailRing() == NULL)
    … …  
    16081608  if (TEST_OPT_SB_1)
    16091609  {
    1610     int k=1;
    1611     int j;
    1612     while(k<=strat->sl)
    1613     {
    1614       j=0;
    1615       loop
    1616       {
    1617         if (j>=k) break;
    1618         clearS(strat->S[j],strat->sevS[j],&k,&j,strat);
    1619         j++;
    1620       }
    1621       k++;
     1610    #ifdef HAVE_RINGS
     1611    if(!rField_is_Ring)
     1612    #endif
     1613    {
     1614        int k=1;
     1615        int j;
     1616        while(k<=strat->sl)
     1617        {
     1618          j=0;
     1619          loop
     1620          {
     1621            if (j>=k) break;
     1622            clearS(strat->S[j],strat->sevS[j],&k,&j,strat);
     1623            j++;
     1624          }
     1625          k++;
     1626        }
    16221627    }
    16231628  }
    … …  
    22552260  if (TEST_OPT_SB_1)
    22562261  {
    2257     int k=1;
    2258     int j;
    2259     while(k<=strat->sl)
    2260     {
    2261       j=0;
    2262       loop
    2263       {
    2264         if (j>=k) break;
    2265         clearS(strat->S[j],strat->sevS[j],&k,&j,strat);
    2266         j++;
    2267       }
    2268       k++;
     2262    #ifdef HAVE_RINGS
     2263    if(!rField_is_Ring)
     2264    #endif
     2265    {
     2266        int k=1;
     2267        int j;
     2268        while(k<=strat->sl)
     2269        {
     2270          j=0;
     2271          loop
     2272          {
     2273            if (j>=k) break;
     2274            clearS(strat->S[j],strat->sevS[j],&k,&j,strat);
     2275            j++;
     2276          }
     2277          k++;
     2278        }
    22692279    }
    22702280  }

  • kernel/kutil.cc

    r4e4e24b rc96539c  
    75907590  if(TEST_OPT_SB_1)
    75917591  {
    7592     int i;
    7593     ideal P=idInit(IDELEMS(F)-strat->newIdeal,F->rank);
    7594     for (i=strat->newIdeal;i<IDELEMS(F);i++)
    7595     {
    7596       P->m[i-strat->newIdeal] = F->m[i];
    7597       F->m[i] = NULL;
    7598     }
    7599     initSSpecial(F,Q,P,strat);
    7600     for (i=strat->newIdeal;i<IDELEMS(F);i++)
    7601     {
    7602       F->m[i] = P->m[i-strat->newIdeal];
    7603       P->m[i-strat->newIdeal] = NULL;
    7604     }
    7605     idDelete(&P);
     7592    #ifdef HAVE_RINGS
     7593    if(!rField_is_Ring)
     7594    #endif
     7595    {
     7596        int i;
     7597        ideal P=idInit(IDELEMS(F)-strat->newIdeal,F->rank);
     7598        for (i=strat->newIdeal;i<IDELEMS(F);i++)
     7599        {
     7600          P->m[i-strat->newIdeal] = F->m[i];
     7601          F->m[i] = NULL;
     7602        }
     7603        initSSpecial(F,Q,P,strat);
     7604        for (i=strat->newIdeal;i<IDELEMS(F);i++)
     7605        {
     7606          F->m[i] = P->m[i-strat->newIdeal];
     7607          P->m[i-strat->newIdeal] = NULL;
     7608        }
     7609        idDelete(&P);
     7610    }
    76067611  }
    76077612  else
    … …  
    76147619  if (!TEST_OPT_SB_1)
    76157620  {
     7621    #ifdef HAVE_RINGS
     7622    if(!rField_is_Ring)
     7623    #endif
    76167624    updateS(TRUE,strat);
    76177625  }
    … …  
    77717779  if(TEST_OPT_SB_1)
    77727780  {
    7773     int i;
    7774     ideal P=idInit(IDELEMS(F)-strat->newIdeal,F->rank);
    7775     for (i=strat->newIdeal;i<IDELEMS(F);i++)
    7776     {
    7777       P->m[i-strat->newIdeal] = F->m[i];
    7778       F->m[i] = NULL;
    7779     }
    7780     initSSpecialSba(F,Q,P,strat);
    7781     for (i=strat->newIdeal;i<IDELEMS(F);i++)
    7782     {
    7783       F->m[i] = P->m[i-strat->newIdeal];
    7784       P->m[i-strat->newIdeal] = NULL;
    7785     }
    7786     idDelete(&P);
     7781    #ifdef HAVE_RINGS
     7782    if(!rField_is_Ring)
     7783    #endif
     7784    {
     7785        int i;
     7786        ideal P=idInit(IDELEMS(F)-strat->newIdeal,F->rank);
     7787        for (i=strat->newIdeal;i<IDELEMS(F);i++)
     7788        {
     7789          P->m[i-strat->newIdeal] = F->m[i];
     7790          F->m[i] = NULL;
     7791        }
     7792        initSSpecialSba(F,Q,P,strat);
     7793        for (i=strat->newIdeal;i<IDELEMS(F);i++)
     7794        {
     7795          F->m[i] = P->m[i-strat->newIdeal];
     7796          P->m[i-strat->newIdeal] = NULL;
     7797        }
     7798        idDelete(&P);
     7799    }
    77877800  }
    77887801  else
    … …  
    77957808  if (!TEST_OPT_SB_1)
    77967809  {
     7810    #ifdef HAVE_RINGS
     7811    if(!rField_is_Ring)
     7812    #endif
    77977813    updateS(TRUE,strat);
    77987814  }
    … …  
    90249040  if(TEST_OPT_SB_1)
    90259041  {
    9026     int i;
    9027     ideal P=idInit(IDELEMS(F)-strat->newIdeal,F->rank);
    9028     for (i=strat->newIdeal;i<IDELEMS(F);i++)
    9029     {
    9030       P->m[i-strat->newIdeal] = F->m[i];
    9031       F->m[i] = NULL;
    9032     }
    9033     initSSpecial(F,Q,P,strat);
    9034     for (i=strat->newIdeal;i<IDELEMS(F);i++)
    9035     {
    9036       F->m[i] = P->m[i-strat->newIdeal];
    9037       P->m[i-strat->newIdeal] = NULL;
    9038     }
    9039     idDelete(&P);
     9042    #ifdef HAVE_RINGS
     9043    if(!rField_is_Ring)
     9044    #endif
     9045    {
     9046        int i;
     9047        ideal P=idInit(IDELEMS(F)-strat->newIdeal,F->rank);
     9048        for (i=strat->newIdeal;i<IDELEMS(F);i++)
     9049        {
     9050          P->m[i-strat->newIdeal] = F->m[i];
     9051          F->m[i] = NULL;
     9052        }
     9053        initSSpecial(F,Q,P,strat);
     9054        for (i=strat->newIdeal;i<IDELEMS(F);i++)
     9055        {
     9056          F->m[i] = P->m[i-strat->newIdeal];
     9057          P->m[i-strat->newIdeal] = NULL;
     9058        }
     9059        idDelete(&P);
     9060    }
    90409061  }
    90419062  else
    … …  
    90499070  {
    90509071    /* the only change: we do not fill the set T*/
     9072    #ifdef HAVE_RINGS
     9073    if(!rField_is_Ring)
     9074    #endif
    90519075    updateS(FALSE,strat);
    90529076  }
Note: See TracChangeset for help on using the changeset viewer.