Changeset 861529 in git for kernel/ideals.cc


Ignore:
Timestamp:
Jul 26, 2011, 1:21:51 PM (12 years ago)
Author:
Hans Schoenemann <hannes@…>
Branches:
(u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', 'f875bbaccd0831e36aaed09ff6adeb3eb45aeb94')
Children:
92058150caf53dc60a028c63bafb64f39bfa73a8
Parents:
b7cfaff1d9a4ec8364a64a54d2b7e009f91418e2
git-author:
Hans Schoenemann <hannes@mathematik.uni-kl.de>2011-07-26 13:21:51+02:00
git-committer:
Mohamed Barakat <mohamed.barakat@rwth-aachen.de>2011-11-09 12:53:34+01:00
Message:
fix pShift -> p_Shift
File:
1 edited

Legend:

Unmodified
Added
Removed

  • kernel/ideals.cc

    rb7cfaf r861529  
    274274  rChangeCurrRing(r);
    275275  // create 1-t, t
    276   poly omt=pOne();
    277   pSetExp(omt,r->N,1);
    278   poly t=pCopy(omt);
    279   pSetm(omt);
    280   omt=pNeg(omt);
    281   omt=pAdd(omt,pOne());
     276  poly omt=p_One(currRing);
     277  p_SetExp(omt,r->N,1,currRing);
     278  poly t=p_Copy(omt,currRing);
     279  p_Setm(omt,currRing);
     280  omt=p_Neg(omt,currRing);
     281  omt=p_Add_q(omt,pOne(),currRing);
    282282  // compute (1-t)*h1
    283   h1=(ideal)mpMultP((matrix)h1,omt);
     283  h1=(ideal)mp_MultP((matrix)h1,omt,currRing);
    284284  // compute t*h2
    285   h2=(ideal)mpMultP((matrix)h2,pCopy(t));
     285  h2=(ideal)mp_MultP((matrix)h2,pCopy(t),currRing);
    286286  // (1-t)h1 + t*h2
    287287  h=idInit(IDELEMS(h1)+IDELEMS(h2),1);
    … …  
    347347
    348348  ring orig_ring=currRing;
    349   ring syz_ring=rCurrRingAssure_SyzComp();
     349  ring syz_ring=rAssure_SyzComp(orig_ring,TRUE);
     350  rChangeCurrRing(syz_ring);
    350351  rSetSyzComp(length, syz_ring);
    351352
    … …  
    360361        temp->m[k] = pCopy(first->m[i]);
    361362      else
    362         temp->m[k] = prCopyR(first->m[i], orig_ring);
     363        temp->m[k] = prCopyR(first->m[i], orig_ring, syz_ring);
    363364      q = pOne();
    364365      pSetComp(q,i+1+length);
    365366      pSetmComp(q);
    366       if (flength==0) pShift(&(temp->m[k]),1);
     367      if (flength==0) p_Shift(&(temp->m[k]),1,currRing);
    367368      p = temp->m[k];
    368369      while (pNext(p)!=NULL) pIter(p);
    … …  
    378379        temp->m[k] = pCopy(second->m[i]);
    379380      else
    380         temp->m[k] = prCopyR(second->m[i], orig_ring);
    381       if (slength==0) pShift(&(temp->m[k]),1);
     381        temp->m[k] = prCopyR(second->m[i], orig_ring,currRing);
     382      if (slength==0) p_Shift(&(temp->m[k]),1,currRing);
    382383      k++;
    383384    }
    … …  
    486487
    487488  ring orig_ring=currRing;
    488   ring syz_ring=rCurrRingAssure_SyzComp();
     489  ring syz_ring=rAssure_SyzComp(orig_ring,TRUE);
     490  rChangeCurrRing(syz_ring);
    489491  rSetSyzComp(syzComp, syz_ring);
    490492
    … …  
    515517            bigmat->m[i] = pCopy(arg[j]->m[l]);
    516518          else
    517             bigmat->m[i] = prCopyR(arg[j]->m[l], orig_ring);
    518           pShift(&(bigmat->m[i]),k*maxrk+isIdeal);
     519            bigmat->m[i] = prCopyR(arg[j]->m[l], orig_ring,currRing);
     520          p_Shift(&(bigmat->m[i]),k*maxrk+isIdeal,currRing);
    519521          i++;
    520522        }
    … …  
    542544      else
    543545        p = prCopyR(tempstd->m[j], syz_ring);
    544       pShift(&p,-syzComp-isIdeal);
     546      p_Shift(&p,-syzComp-isIdeal,currRing);
    545547      result->m[k] = p;
    546548      k++;
    … …  
    578580  if (k == 0)
    579581  {
    580     for (j=0; j<=i; j++) pShift(&(h2->m[j]),1);
     582    for (j=0; j<=i; j++) p_Shift(&(h2->m[j]),1,currRing);
    581583    k = 1;
    582584  }
    … …  
    678680  {
    679681    ideal result=idFreeModule(idElemens_h1/*IDELEMS(h1)*/);
    680     int curr_syz_limit=rGetCurrSyzLimit();
     682    int curr_syz_limit=rGetCurrSyzLimit(currRing);
    681683    if (curr_syz_limit>0)
    682684    for (ii=0;ii<idElemens_h1/*IDELEMS(h1)*/;ii++)
    683685    {
    684686      if (h1->m[ii]!=NULL)
    685         pShift(&h1->m[ii],curr_syz_limit);
     687        p_Shift(&h1->m[ii],curr_syz_limit,currRing);
    686688    }
    687689    return result;
    … …  
    692694  assume(currRing != NULL);
    693695  ring orig_ring=currRing;
    694   ring syz_ring=rCurrRingAssure_SyzComp();
     696  ring syz_ring=rAssure_SyzComp(orig_ring.TRUE);
     697  rChangeCurrRing(syz_ring);
    695698
    696699  if (setSyzComp)
    … …  
    723726      {
    724727        if (p_MinComp(s_h3->m[j],syz_ring) > k)
    725           pShift(&s_h3->m[j], -k);
     728          p_Shift(&s_h3->m[j], -k,currRing);
    726729        else
    727730          pDelete(&s_h3->m[j]);
    … …  
    814817  assume(currRing != NULL);
    815818  ring orig_ring=currRing;
    816   ring syz_ring=rCurrRingAssure_SyzComp();
     819  ring syz_ring=rAssure_SyzComp(orig_ring,TRUE);
     820  rChangeCurrRing(syz_ring);
    817821
    818822  rSetSyzComp(k,syz_ring);
    … …  
    869873    {
    870874      *syz=idFreeModule(IDELEMS(h1));
    871       int curr_syz_limit=rGetCurrSyzLimit();
     875      int curr_syz_limit=rGetCurrSyzLimit(currRing);
    872876      if (curr_syz_limit>0)
    873877      for (int ii=0;ii<IDELEMS(h1);ii++)
    874878      {
    875879        if (h1->m[ii]!=NULL)
    876           pShift(&h1->m[ii],curr_syz_limit);
     880          p_Shift(&h1->m[ii],curr_syz_limit,currRing);
    877881      }
    878882    }
    … …  
    887891
    888892  ring orig_ring = currRing;
    889   ring syz_ring = rCurrRingAssure_SyzComp();
     893  ring syz_ring = rAssure_SyzComp(orig_ring,TRUE);
     894  rChangeCurrRing(syz_ring);
    890895  rSetSyzComp(k,syz_ring);
    891896
    … …  
    930935          }
    931936        }
    932         if (!inputIsIdeal) pShift(&(s_h3->m[j]), -1);
     937        if (!inputIsIdeal) p_Shift(&(s_h3->m[j]), -1,currRing);
    933938      }
    934939      else
    … …  
    937942        if (lift3)
    938943        {
    939           pShift(&s_h3->m[j], -k);
     944          p_Shift(&s_h3->m[j], -k,currRing);
    940945          (*syz)->m[j]=s_h3->m[j];
    941946          s_h3->m[j]=NULL;
    … …  
    10941099
    10951100  ring orig_ring=currRing;
    1096   ring syz_ring=rCurrRingAssure_SyzComp();
     1101  ring syz_ring=rAssure_SyzComp(orig_ring,TRUE);
     1102  rChangeCurrRing(syz_ring);
    10971103  rSetSyzComp(k,syz_ring);
    10981104
    … …  
    11321138    {
    11331139      if (s_temp->m[j-1]!=NULL)
    1134         pShift(&(s_temp->m[j-1]),1);
     1140        p_Shift(&(s_temp->m[j-1]),1,currRing);
    11351141    }
    11361142  }
    … …  
    11851191        }
    11861192      }
    1187       pShift(&(s_result->m[j]),-k);
     1193      p_Shift(&(s_result->m[j]),-k,currRing);
    11881194      pNeg(s_result->m[j]);
    11891195    }
    … …  
    11951201      if (s_rest->m[j-1]!=NULL)
    11961202      {
    1197         pShift(&(s_rest->m[j-1]),-1);
     1203        p_Shift(&(s_rest->m[j-1]),-1,currRing);
    11981204        s_rest->m[j-1] = s_rest->m[j-1];
    11991205      }
    … …  
    12451251        }
    12461252      }
    1247       pShift(&s_result->m[i],-comps_to_add);
     1253      p_Shift(&s_result->m[i],-comps_to_add,currRing);
    12481254    }
    12491255  }
    … …  
    13511357      p = pCopy(h2->m[i]);
    13521358      if (k2 == 0)
    1353         pShift(&p,j*k+1);
     1359        p_Shift(&p,j*k+1,currRing);
    13541360      else
    1355         pShift(&p,j*k);
     1361        p_Shift(&p,j*k,currRing);
    13561362      q = pAdd(q,p);
    13571363      j++;
    … …  
    13811387      {
    13821388        p = pCopy_noCheck(h4->m[i-1]);
    1383         pShift(&p,1);
     1389        p_Shift(&p,1,currRing);
    13841390        h4->m[i] = p;
    13851391      }
    … …  
    13971403        p = pCopy(temph1->m[l]);
    13981404        if (k1 == 0)
    1399           pShift(&p,ll*k+1);
     1405          p_Shift(&p,ll*k+1,currRing);
    14001406        else
    1401           pShift(&p,ll*k);
     1407          p_Shift(&p,ll*k,currRing);
    14021408        if (kkk >= IDELEMS(h4))
    14031409        {
    … …  
    14541460
    14551461  ring orig_ring=currRing;
    1456   ring syz_ring=rCurrRingAssure_SyzComp();
     1462  ring syz_ring=rAssure_SyzComp(orig_ring.TRUE);
     1463  rChangeCurrRing(syz_ring);
    14571464  rSetSyzComp(kmax-1,syz_ring);
    14581465  if (orig_ring!=syz_ring)
    … …  
    14951502    {
    14961503      if (resultIsIdeal)
    1497         pShift(&s_h3->m[i],-kmax);
     1504        p_Shift(&s_h3->m[i],-kmax,currRing);
    14981505      else
    1499         pShift(&s_h3->m[i],-kmax+1);
     1506        p_Shift(&s_h3->m[i],-kmax+1,currRing);
    15001507    }
    15011508    else
    … …  
    21382145    else
    21392146    {
    2140       pShift(&(temp1->m[i]),-length);
     2147      p_Shift(&(temp1->m[i]),-length,currRing);
    21412148    }
    21422149  }
    … …  
    21992206    if(temp->m[i]!=NULL)
    22002207    {
    2201       if (slength==0) pShift(&(temp->m[i]),1);
     2208      if (slength==0) p_Shift(&(temp->m[i]),1,currRing);
    22022209      p = temp->m[i];
    22032210      while (pNext(p)!=NULL) pIter(p);
    … …  
    22172224      {
    22182225        temp->m[k] = pCopy(h1->m[i]);
    2219         if (flength==0) pShift(&(temp->m[k]),1);
     2226        if (flength==0) p_Shift(&(temp->m[k]),1,currRing);
    22202227        k++;
    22212228      }
    … …  
    22242231
    22252232  ring orig_ring=currRing;
    2226   ring syz_ring=rCurrRingAssure_SyzComp();
     2233  ring syz_ring=rAssure_SyzComp(orig_ring.TRUE);
     2234  rChangeCurrRing(syz_ring);
    22272235  rSetSyzComp(length, syz_ring);
    22282236  ideal s_temp;
    … …  
    22592267    else
    22602268    {
    2261       pShift(&(s_temp1->m[i]),-length);
     2269      p_Shift(&(s_temp1->m[i]),-length,currRing);
    22622270    }
    22632271  }
Note: See TracChangeset for help on using the changeset viewer.