Changeset 6d8734 in git


Ignore:
Timestamp:
Jul 26, 2021, 2:34:40 PM (3 years ago)
Author:
Hans Schoenemann <hannes@…>
Branches:
(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
Children:
62de185de0060170d8506f950869e3cb1d29e6ab
Parents:
fe2306067e5322b21216ff8f673910a0dd8dd110
git-author:
Hans Schoenemann <hannes@mathematik.uni-kl.de>2021-07-26 14:34:40+02:00
git-committer:
Hans Schoenemann <hannes@mathematik.uni-kl.de>2021-07-26 14:35:37+02:00
Message:
ringlist -> ring_list in primdec.lib
File:
1 edited

Legend:

Unmodified
Added
Removed
  • Singular/LIB/primdec.lib

    rfe23060 r6d8734  
    11////////////////////////////////////////////////////////////////////////////
    2 version="version primdec.lib 4.2.0.0 Mar_2021 "; // $Id$
     2version="version primdec.lib 4.2.1.1 Jul_2021 "; // $Id$
    33category="Commutative Algebra";
    44info="
     
    18821882   i=qr[1];
    18831883
    1884    def gnir=ring(ringlist(@P));
     1884   def gnir=ring(ring_list(@P));
    18851885   setring gnir;
    18861886
     
    20642064  int j;          // Counter
    20652065  def P0 = basering;
    2066   list Pl=ringlist(P0);
     2066  list Pl=ring_list(P0);
    20672067  intvec dp_w;
    20682068  for(j=nvars(P0);j>0;j--) {dp_w[j]=1;}
     
    24172417                                &&(find(ordstr(basering),"s")==0))
    24182418  {
    2419      def gnir=ring(ringlist(basering));
     2419     def gnir=ring(ring_list(basering));
    24202420     setring gnir;
    24212421     ideal i=imap(P,i);
     
    25222522                                &&(find(ordstr(basering),"s")==0))
    25232523  {
    2524      def gnir=ring(ringlist(basering));
     2524     def gnir=ring(ring_list(basering));
    25252525     setring gnir;
    25262526     ideal i=imap(P,i);
     
    25602560
    25612561  indep=maxIndependSet(j);
    2562   ring gnir1 = create_ring(ringlist(basering)[1],"("+indep[1][1]+")","("+indep[1][2]+")","no_minpoly");
     2562  ring gnir1 = create_ring(ring_list(basering)[1],"("+indep[1][1]+")","("+indep[1][2]+")","no_minpoly");
    25632563  if(homo==1)
    25642564  {
     
    31063106      newva[size(newva)]=")";
    31073107      newma[size(newma)]=";";
    3108       ring @deirf = create_ring(ringlist(gnir)[1],"("+newva,"lp","no_minpoly");
     3108      ring @deirf = create_ring(ring_list(gnir)[1],"("+newva,"lp","no_minpoly");
    31093109      execute("map @kappa=gnir,"+newma);
    31103110      ideal @j= @kappa(@j);
     
    32703270
    32713271  ideal jkeep=@j;
    3272   if((ordstr(@P)[1]=="w")&&(size(ringlist(@P)[3])==2)) // weighted ordering
    3273   {
    3274     list gnir_l=ringlist(gnir);
    3275     list @P_l=ringlist(@P);
     3272  if((ordstr(@P)[1]=="w")&&(size(ring_list(@P)[3])==2)) // weighted ordering
     3273  {
     3274    list gnir_l=ring_list(gnir);
     3275    list @P_l=ring_list(@P);
    32763276    gnir_l[3]=@P_l[3]; // ord
    32773277    def @Phelp=ring(gnir_l);
     
    32883288  {
    32893289    if(((ordstr(@P)[3]=="d")||(ordstr(@P)[1]=="d")||(ordstr(@P)[1]=="w")
    3290        ||(ordstr(@P)[3]=="w"))&&(size(ringlist(@P)[3])==2))
     3290       ||(ordstr(@P)[3]=="w"))&&(size(ring_list(@P)[3])==2))
    32913291    {
    32923292      ideal jwork=imap(@P,tras);
     
    33183318    //change the ring
    33193319    {
    3320       def gnir1=ring(ringlist(basering));
     3320      def gnir1=ring(ring_list(basering));
    33213321      setring gnir1;
    33223322      ideal @j=fetch(gnir,@j);
     
    33313331        @j=fetch(@P,i);
    33323332      }
    3333       ring gnir1 = create_ring(ringlist(basering)[1],"("+indep[@m][1]+")","("+indep[@m][2]+")","no_minpoly");
     3333      ring gnir1 = create_ring(ring_list(basering)[1],"("+indep[@m][1]+")","("+indep[@m][2]+")","no_minpoly");
    33343334      execute("map phi=gnir,"+@va+";");
    33353335      op=option(get);
     
    36313631           //change the ring
    36323632        {
    3633           def gnir1=ring(ringlist(basering));
     3633          def gnir1=ring(ring_list(basering));
    36343634          setring gnir1;
    36353635          ideal @j=fetch(gnir,jkeep);
     
    36393639        {
    36403640          @va=string(maxideal(1));
    3641           ring gnir1 = create_ring(ringlist(basering)[1],"("+restindep[@m][1]+")","("+restindep[@m][2]+")","no_minpoly");
     3641          ring gnir1 = create_ring(ring_list(basering)[1],"("+restindep[@m][1]+")","("+restindep[@m][2]+")","no_minpoly");
    36423642          execute("map phi=gnir,"+@va+";");
    36433643          op=option(get);
     
    46094609  def oldring=basering;
    46104610  string n=system("neworder",PS);
    4611   ring r = create_ring(ringlist(oldring)[1],"("+n+")","dp","no_minpoly");
     4611  ring r = create_ring(ring_list(oldring)[1],"("+n+")","dp","no_minpoly");
    46124612  ideal PS=imap(oldring,PS);
    46134613  matrix m=char_series(PS);  // We compute an irreducible
     
    46244624    for(i=1;i<nvars(oldring);i++) { n=n+","+string(var(i)); }
    46254625    kill r;
    4626     ring r = create_ring(ringlist(oldring)[1],"("+n+")","dp","no_minpoly");
     4626    ring r = create_ring(ring_list(oldring)[1],"("+n+")","dp","no_minpoly");
    46274627    ideal PS=imap(oldring,PS);
    46284628    matrix m=char_series(PS);
     
    52675267
    52685268      U[size(U)]=")";           // we compute the extractor of I (w.r.t. U)
    5269       ring RAU = create_ring(ringlist(basering)[1],"("+A+U,"(dp("+string(a)+"),dp)","no_minpoly");
     5269      ring RAU = create_ring(ring_list(basering)[1],"("+A+U,"(dp("+string(a)+"),dp)","no_minpoly");
    52705270      ideal I=imap(R,SI);
    52715271      //I=std(I,hv);            // the standard basis in (R[U])[A]
     
    58375837  int k=find(newvar,",",find(newvar,",")+1);
    58385838  newvar=newvar[k+1..size(newvar)];
    5839   list lR=ringlist(R);
     5839  list lR=ring_list(R);
    58405840  int i,de,ii;
    58415841  intvec vv=1:n;
     
    63746374  int j;
    63756375  def P0 = basering;
    6376   list Pl=ringlist(P0);
     6376  list Pl=ring_list(P0);
    63776377  intvec dp_w;
    63786378  for(j=nvars(P0);j>0;j--) {dp_w[j]=1;}
     
    67496749    //change the ring
    67506750    {
    6751       def gnir1=ring(ringlist(basering));
     6751      def gnir1=ring(ring_list(basering));
    67526752      setring gnir1;
    67536753      ideal @j = fetch(@P, I);
     
    67586758      @va = string(maxideal(1));
    67596759
    6760       ring gnir1 = create_ring(ringlist(basering)[1],"(" + indep[@m][1] + ")","(" + indep[@m][2] + ")","no_minpoly");
     6760      ring gnir1 = create_ring(ring_list(basering)[1],"(" + indep[@m][1] + ")","(" + indep[@m][2] + ")","no_minpoly");
    67616761      execute("map phi = @P," + @va + ";");
    67626762      ideal @j = groebner(phi(I));
     
    73197319    string(timer-time));
    73207320//---------------- compute characteristic polynomial of matrix --------------
    7321   ring P1 = create_ring(ringlist(R)[1],"T","dp","no_minpoly");
     7321  ring P1 = create_ring(ring_list(R)[1],"T","dp","no_minpoly");
    73227322  matrix n=imap(R,n);
    73237323  time = timer;
     
    73497349//Hence it is better to simultaneously reduce with I. For this we need a new
    73507350//ring.
    7351     ring P = create_ring(ringlist(R)[1],"(T,"+varstr(R)+")","(dp(1),dp)","no_minpoly");
     7351    ring P = create_ring(ring_list(R)[1],"(T,"+varstr(R)+")","(dp(1),dp)","no_minpoly");
    73527352    list rfac=imap(P1,fac);
    73537353    intvec ov=option(get);;
     
    76467646      newva[size(newva)]=")";
    76477647      newma[size(newma)]=";";
    7648       ring @deirf = create_ring(ringlist(gnir)[1],"("+newva,"lp","no_minpoly");
     7648      ring @deirf = create_ring(ring_list(gnir)[1],"("+newva,"lp","no_minpoly");
    76497649      execute("map @kappa=gnir,"+newma);
    76507650      ideal @j= @kappa(@j);
     
    78467846
    78477847  ideal jkeep=@j;
    7848   if((ordstr(@P)[1]=="w")&&(size(ringlist(@P)[3])==2)) // weighted ordering
    7849   {
    7850     def @Phelp=ring(ringlist(gnir));
     7848  if((ordstr(@P)[1]=="w")&&(size(ring_list(@P)[3])==2)) // weighted ordering
     7849  {
     7850    def @Phelp=ring(ring_list(gnir));
    78517851    setring @Phelp;
    78527852  }
     
    80298029          //change the ring
    80308030          {
    8031              def gnir1=ring(ringlist(basering));
     8031             def gnir1=ring(ring_list(basering));
    80328032             setring gnir1;
    80338033             ideal @j=fetch(gnir,jkeep);
     
    80378037          {
    80388038             @va=string(maxideal(1));
    8039              ring gnir1 = create_ring(ringlist(basering)[1],"(" +restindep[@m][1]+")","(" +restindep[@m][2]+")","no_minpoly");
     8039             ring gnir1 = create_ring(ring_list(basering)[1],"(" +restindep[@m][1]+")","(" +restindep[@m][2]+")","no_minpoly");
    80408040             execute("map phi=gnir,"+@va+";");
    80418041             op=option(get);
     
    83658365   //change the ring
    83668366   {
    8367      def gnir1=ring(ringlist(basering));
     8367     def gnir1=ring(ring_list(basering));
    83688368     setring gnir1;
    83698369     ideal @j = fetch(gnir, @j);
     
    83808380     //}
    83818381//=======================
    8382      ring gnir1 = create_ring(ringlist(basering)[1],"("+indepInfo[1]+")","(" +indepInfo[2]+")","no_minpoly");
     8382     ring gnir1 = create_ring(ring_list(basering)[1],"("+indepInfo[1]+")","(" +indepInfo[2]+")","no_minpoly");
    83838383     execute("map phi=gnir,"+@va+";");
    83848384     op=option(get);
Note: See TracChangeset for help on using the changeset viewer.