Changeset e2af3af in git for Singular/LIB/primdec.lib


Ignore:
Timestamp:
Jul 1, 2019, 4:07:30 PM (5 years ago)
Author:
Hans Schoenemann <hannes@…>
Branches:
(u'spielwiese', '2a584933abf2a2d3082034c7586d38bb6de1a30a')
Children:
2050058328d464a197ff417b127e85fe08b3a3cd
Parents:
77006ca97061dca7e90bdacf1fbaff76644280a24142c906ad243d17ff54542115db04d9e8264de7
git-author:
Hans Schoenemann <hannes@mathematik.uni-kl.de>2019-07-01 16:07:30+02:00
git-committer:
GitHub <noreply@github.com>2019-07-01 16:07:30+02:00
Message:
Merge pull request #936 from steenpass/execute

replace some occurrences of execute()
File:
1 edited

Legend:

Unmodified
Added
Removed
  • Singular/LIB/primdec.lib

    r77006c re2af3af  
    25642564
    25652565  indep=maxIndependSet(j);
    2566 
    2567   execute("ring gnir1 = ("+charstr(basering)+"),("+indep[1][1]+"),("
    2568                               +indep[1][2]+");");
     2566  ring gnir1 = create_ring(ringlist(basering)[1],"("+indep[1][1]+")","("+indep[1][2]+")","no_minpoly");
    25692567  if(homo==1)
    25702568  {
     
    31123110      newva[size(newva)]=")";
    31133111      newma[size(newma)]=";";
    3114       execute("ring @deirf=("+charstr(gnir)+"),("+newva+",lp;");
     3112      ring @deirf = create_ring(ringlist(gnir)[1],"("+newva,"lp","no_minpoly");
    31153113      execute("map @kappa=gnir,"+newma);
    31163114      ideal @j= @kappa(@j);
     
    33373335        @j=fetch(@P,i);
    33383336      }
    3339       execute("ring gnir1 = ("+charstr(basering)+"),("+indep[@m][1]+"),("
    3340                               +indep[@m][2]+");");
     3337      ring gnir1 = create_ring(ringlist(basering)[1],"("+indep[@m][1]+")","("+indep[@m][2]+")","no_minpoly");
    33413338      execute("map phi=gnir,"+@va+";");
    33423339      op=option(get);
     
    36463643        {
    36473644          @va=string(maxideal(1));
    3648           execute("ring gnir1 = ("+charstr(basering)+"),("+
    3649                       restindep[@m][1]+"),(" +restindep[@m][2]+");");
     3645          ring gnir1 = create_ring(ringlist(basering)[1],"("+restindep[@m][1]+")","("+restindep[@m][2]+")","no_minpoly");
    36503646          execute("map phi=gnir,"+@va+";");
    36513647          op=option(get);
     
    46144610  def oldring=basering;
    46154611  string n=system("neworder",PS);
    4616   execute("ring r=("+charstr(oldring)+"),("+n+"),dp;");
     4612  ring r = create_ring(ringlist(oldring)[1],"("+n+")","dp","no_minpoly");
    46174613  ideal PS=imap(oldring,PS);
    46184614  matrix m=char_series(PS);  // We compute an irreducible
     
    46294625    for(i=1;i<nvars(oldring);i++) { n=n+","+string(var(i)); }
    46304626    kill r;
    4631     execute("ring r=("+charstr(oldring)+"),("+n+"),dp;");
     4627    ring r = create_ring(ringlist(oldring)[1],"("+n+")","dp","no_minpoly");
    46324628    ideal PS=imap(oldring,PS);
    46334629    matrix m=char_series(PS);
     
    52725268
    52735269      U[size(U)]=")";           // we compute the extractor of I (w.r.t. U)
    5274       execute("ring RAU=("+charstr(basering)+"),("+A+U+",(dp("+string(a)+"),dp);");
     5270      ring RAU = create_ring(ringlist(basering)[1],"("+A+U,"(dp("+string(a)+"),dp)","no_minpoly");
    52755271      ideal I=imap(R,SI);
    52765272      //I=std(I,hv);            // the standard basis in (R[U])[A]
     
    67636759      @va = string(maxideal(1));
    67646760
    6765       execute("ring gnir1 = (" + charstr(basering) + "), (" + indep[@m][1] + "),("
    6766                               + indep[@m][2] + ");");
     6761      ring gnir1 = create_ring(ringlist(basering)[1],"(" + indep[@m][1] + ")","(" + indep[@m][2] + ")","no_minpoly");
    67676762      execute("map phi = @P," + @va + ";");
    67686763      ideal @j = groebner(phi(I));
     
    73257320    string(timer-time));
    73267321//---------------- compute characteristic polynomial of matrix --------------
    7327   execute("ring P1=("+charstr(R)+"),T,dp;");
     7322  ring P1 = create_ring(ringlist(R)[1],"T","dp","no_minpoly");
    73287323  matrix n=imap(R,n);
    73297324  time = timer;
     
    73557350//Hence it is better to simultaneously reduce with I. For this we need a new
    73567351//ring.
    7357     execute("ring P=("+charstr(R)+"),(T,"+varstr(R)+"),(dp(1),dp);");
     7352    ring P = create_ring(ringlist(R)[1],"(T,"+varstr(R)+")","(dp(1),dp)","no_minpoly");
    73587353    list rfac=imap(P1,fac);
    73597354    intvec ov=option(get);;
     
    76527647      newva[size(newva)]=")";
    76537648      newma[size(newma)]=";";
    7654       execute("ring @deirf=("+charstr(gnir)+"),("+newva+",lp;");
     7649      ring @deirf = create_ring(ringlist(gnir)[1],"("+newva,"lp","no_minpoly");
    76557650      execute("map @kappa=gnir,"+newma);
    76567651      ideal @j= @kappa(@j);
     
    80438038          {
    80448039             @va=string(maxideal(1));
    8045              execute("ring gnir1 = ("+charstr(basering)+"),("+
    8046                      restindep[@m][1]+"),(" +restindep[@m][2]+");");
     8040             ring gnir1 = create_ring(ringlist(basering)[1],"(" +restindep[@m][1]+")","(" +restindep[@m][2]+")","no_minpoly");
    80478041             execute("map phi=gnir,"+@va+";");
    80488042             op=option(get);
     
    83878381     //}
    83888382//=======================
    8389      execute("ring gnir1 = ("+charstr(basering)+"),("+indepInfo[1]+"),("
    8390                               +indepInfo[2]+");");
     8383     ring gnir1 = create_ring(ringlist(basering)[1],"("+indepInfo[1]+")","(" +indepInfo[2]+")","no_minpoly");
    83918384     execute("map phi=gnir,"+@va+";");
    83928385     op=option(get);
Note: See TracChangeset for help on using the changeset viewer.