Changeset 6f4f9c in git


Ignore:
Timestamp:
Aug 15, 2019, 5:02:21 PM (5 years ago)
Author:
Sachin <sachinkm308@…>
Branches:
(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
Children:
d08c429701793a3ea40aacddaa1284f49fbf6167
Parents:
a910024f22f240ce40b16d4375550fdd82c28a29
git-author:
Sachin <sachinkm308@gmail.com>2019-08-15 17:02:21+02:00
git-committer:
Sachin <sachinkm308@gmail.com>2019-08-27 15:42:13+02:00
Message:
replacing execute with create_ring (6)- using for loops
Location:
Singular/LIB
Files:
7 edited

Legend:

Unmodified
Added
Removed

  • Singular/LIB/chern.lib

    ra91002 r6f4f9c  
    190190  def br@=basering; // remember the base ring
    191191  // add additional variables to the base ring
    192   execute("ring r@= ("+ charstr(basering) +"),("+varstr(basering)+",c@(1..nV),A@(1..nV)), dp;" );
     192  int ii;
     193  list l2 = ringlist(basering)[2];
     194  for (ii = 1; ii <= nV; ii++)
     195  {
     196   l2[size(l2)+1] = "c@("+string(ii)+")";
     197  }
     198    for (ii = 1; ii <= nV; ii++)
     199  {
     200   l2[size(l2)+1] = "A@("+string(ii)+")";
     201  }
     202  ring r@ = create_ring(ringlist(basering)[1], l2, "dp", "no_minpoly");
    193203  execute( "map F= br@,"+varstr(br@)+";" ); // define the corresponding inclusion of rings
    194204  ideal V=F(V);
    … …  
    688698    def br@=basering; // remember the base ring
    689699    // add additional variables c@, a@(1..n) to the base ring
    690     execute("ring r@= (" + charstr(basering) + "),(c@,"+varstr(basering)+", a@(1..n)), lp;" );
     700    list l3 = "c@";
     701    l3 = l3+ ringlist(basering)[2];
     702    for (int ii = 1; ii <= n; ii++)
     703    {
     704     l3[size(l3)+1] = "a@("+string(ii)+")";
     705    }
     706    ring r@ = create_ring(ringlist(basering)[1], l3, "lp", "no_minpoly");
    691707    execute( "map F= br@,"+varstr(br@)+";" ); // define the corresponding inclusion of rings
    692708    list c=F(c); // embed c in the bigger ring
    … …  
    10181034  def br@=basering; // remember the base ring
    10191035  // add additional variables a@(1..r), b@(1..R), x@ to the base ring
    1020   execute("ring r@=("+ charstr(basering) +"),(x@,"+varstr(basering)+",a@(1..r),b@(1..R)), lp;");
     1036  int ii;
     1037  list l4 = "x@";
     1038  l4 = l4+ ringlist(basering)[2];
     1039  for (ii = 1; ii <= r; ii++)
     1040  {
     1041   l4[size(l4)+1] = "a@("+string(ii)+")";
     1042  }
     1043  for (ii = 1; ii <= R; ii++)
     1044  {
     1045   l4[size(l4)+1] = "b@("+string(ii)+")";
     1046  }
     1047  ring r@ = create_ring(ringlist(basering)[1], l4, "lp", "no_minpoly");
    10211048  execute( "map F= br@,"+varstr(br@)+";" ); // define the corresponding inclusion of rings
    10221049  list c=F(c); // embed c in the bigger ring
    … …  
    11011128  def br@=basering; // remember the base ring
    11021129  // add  additional variables to the base ring
    1103   execute("ring r@=("+charstr(basering)+"), ("+varstr(basering)+", t@, c@(1..r), C@(1..R)), dp;");
     1130  int ii;
     1131  list l5 = ringlist(basering)[2]; 
     1132  l5[size(l5)+1] = "t@";
     1133  for (ii = 1; ii <= r; ii++)
     1134  {
     1135   l5[size(l5)+1] = "c@("+string(ii)+")";
     1136  }
     1137  for (ii = 1; ii <= R; ii++)
     1138  {
     1139   l5[size(l5)+1] = "C@("+string(ii)+")";
     1140  }
     1141  ring r@ = create_ring(ringlist(basering)[1], l5, "dp", "no_minpoly");
    11041142  execute( "map F= br@,"+varstr(br@)+";" ); // define the corresponding inclusion of rings
    11051143  list c, C;
    … …  
    12351273  def br@=basering; // remember the base ring
    12361274  // add  additional variables to the base ring
    1237   execute("ring r@=("+charstr(basering)+"), ("+varstr(basering)+", t@, c@(1..r), C@(1..R)), dp;");
     1275  int ii;
     1276  list l6 = ringlist(basering)[2]; 
     1277  l6[size(l6)+1] = "t@";
     1278  for (ii = 1; ii <= r; ii++)
     1279  {
     1280   l6[size(l6)+1] = "c@("+string(ii)+")";
     1281  }
     1282  for (ii = 1; ii <= R; ii++)
     1283  {
     1284   l6[size(l6)+1] = "C@("+string(ii)+")";
     1285  }
     1286  ring r@ = create_ring(ringlist(basering)[1], l6, "dp", "no_minpoly");
    12381287  execute( "map F= br@,"+varstr(br@)+";" ); // define the corresponding inclusion of rings
    12391288  list c, C;
    … …  
    15711620  def br@=basering; // remember the base ring
    15721621  // add additional variables to the base ring
    1573   execute("ring r@=(" + charstr(basering) +  "),(x@,"+varstr(basering)+", a@(1..r)), lp;" );
     1622  int ii;
     1623  list l7 = "x@";
     1624  l7 = l7+ ringlist(basering)[2];
     1625  for (ii = 1; ii <= r; ii++)
     1626  {
     1627   l7[size(l7)+1] = "a@("+string(ii)+")";
     1628  }
     1629  ring r@ = create_ring(ringlist(basering)[1], l7, "lp", "no_minpoly");
    15741630  execute( "map F= br@,"+varstr(br@)+";" ); // define the corresponding inclusion of rings
    15751631  list c=F(c); // embed c into the bigger ring
    … …  
    16971753  def br@=basering; // remember the base ring
    16981754  // add  additional variables to the base ring
    1699   execute("ring r@=(" + charstr(basering) +  "), ("+varstr(basering)+", t@, c@(1..r)), dp;" );
     1755  list l8 = ringlist(basering)[2]; 
     1756  l8[size(l8)+1] = "t@";
     1757  for (int ii = 1; ii <= r; ii++)
     1758  {
     1759   l8[size(l8)+1] = "c@("+string(ii)+")";
     1760  }
     1761  ring r@ = create_ring(ringlist(basering)[1], l8, "dp", "no_minpoly");
    17001762  execute( "map F= br@,"+varstr(br@)+";" ); // define the corresponding inclusion of rings
    17011763  list c;
    … …  
    18631925  def br@=basering; // remember the base ring
    18641926  // add additional variables a@(1..r), x@ to the base ring
    1865   execute("ring r@= (" + charstr(basering) + "), (x@,"+varstr(basering)+", a@(1..r)), lp;" );
     1927  list l9 = "x@";
     1928  l9 = l9+ ringlist(basering)[2];
     1929  for (int ii = 1; ii <= r; ii++)
     1930  {
     1931   l9[size(l9)+1] = "a@("+string(ii)+")";
     1932  }
     1933  ring r@ = create_ring(ringlist(basering)[1], l9, "lp", "no_minpoly");
    18661934  execute( "map F= br@,"+varstr(br@)+";" ); // define the corresponding inclusion of rings
    18671935  list c = F(c); // embed c into the bigger ring
    … …  
    19912059  def br@=basering; // remember the base ring
    19922060  // add  additional variables to the base ring
    1993   execute("ring r@=(" + charstr(basering) +  "), ("+varstr(basering)+", t@, c@(1..r)), dp;" );
     2061  list l10 = ringlist(basering)[2]; 
     2062  l10[size(l10)+1] = "t@";
     2063  for (int ii = 1; ii <= r; ii++)
     2064  {
     2065   l10[size(l10)+1] = "c@("+string(ii)+")";
     2066  }
     2067  ring r@ = create_ring(ringlist(basering)[1], l10, "dp", "no_minpoly");
    19942068  execute( "map F= br@,"+varstr(br@)+";" ); // define the corresponding inclusion of rings
    19952069  list c;
    … …  
    21132187    def br@=basering; // remember the base ring
    21142188    // add  additional variables to the base ring
    2115     execute("ring r@=(" + charstr(basering) +  "), ("+varstr(basering)+", a@, c@(1..n)), dp;" );
     2189    list l11 = ringlist(basering)[2]; 
     2190    l11[size(l11)+1] = "a@";
     2191    for (int ii = 1; ii <= n; ii++)
     2192    {
     2193     l11[size(l11)+1] = "c@("+string(ii)+")";
     2194    }
     2195    ring r@ = create_ring(ringlist(basering)[1], l11, "dp", "no_minpoly");
    21162196    execute( "map F= br@,"+varstr(br@)+";" ); // define the corresponding inclusion of rings
    21172197    list c=F(c); // embed c into the bigger ring
    … …  
    21952275    def br@=basering; // remember the base ring
    21962276    // add additional variables a@(1..n), x@ to the base ring
    2197     execute("ring r@=(" + charstr(basering) + "), (x@,"+varstr(basering)+", a@(1..n)), lp;" );
     2277    list l12 = "x@";
     2278    l12 = l12+ ringlist(basering)[2];
     2279    for (int ii = 1; ii <= n; ii++)
     2280    {
     2281     l12[size(l12)+1] = "a@("+string(ii)+")";
     2282    }
     2283    ring r@ = create_ring(ringlist(basering)[1], l12, "lp", "no_minpoly");
    21982284    execute( "map F= br@,"+varstr(br@)+";" ); // define the corresponding inclusion of rings
    21992285    list c=F(c); // embed c into the bigger ring
    … …  
    30763162  def br@=basering; // remember the base ring
    30773163  // add additional variables t@(1), ... , t@(sz) and u@ to the base ring
    3078   execute("ring r@=("+ charstr(basering) +"),("+varstr(basering)+",t@(1..sz),u@), dp;");
     3164  list l13 = ringlist(basering)[2];
     3165  for (int ii = 1; ii <= sz; ii++)
     3166  {
     3167   l13[size(l13)+1] = "t@("+string(ii)+")";
     3168  }
     3169  l13[size(l13)+1] = "u@";
     3170  ring r@ = create_ring(ringlist(basering)[1], l13, "dp", "no_minpoly");
    30793171  execute( "map F= br@,"+varstr(br@)+";" ); // define the corresponding inclusion of rings
    30803172  ideal I=F(I); // the ideal generated by I in the new ring
    … …  
    31703262  def br@=basering; // remember the base ring
    31713263  // add additional variables t@(1), ... , t@(sz) and u@ to the base ring
    3172   execute("ring r@=("+ charstr(basering) +"),("+varstr(basering)+",t@(1..sz),u@), dp;");
     3264  list l14 = ringlist(basering)[2];
     3265  for (int ii = 1; ii <= sz; ii++)
     3266  {
     3267   l14[size(l14)+1] = "t@("+string(ii)+")";
     3268  }
     3269  l14[size(l14)+1] = "u@";
     3270  ring r@ = create_ring(ringlist(basering)[1], l14, "dp", "no_minpoly");
    31733271  execute( "map F= br@,"+varstr(br@)+";" ); // define the corresponding inclusion of rings
    31743272  ideal I=F(I); // the ideal generated by I in the new ring
    … …  
    38833981  int sz=ncols(I);
    38843982  def br@=basering; // remember the base ring
    3885   execute("ring r@=("+ charstr(basering) +"),("+varstr(basering)+",t@(1..sz),u@), dp;");
     3983  list l15 = ringlist(basering)[2];
     3984  for (int ii = 1; ii <= sz; ii++)
     3985  {
     3986   l15[size(l15)+1] = "t@("+string(ii)+")";
     3987  }
     3988  l15[size(l15)+1] = "u@";
     3989  ring r@ = create_ring(ringlist(basering)[1], l15, "dp", "no_minpoly");
    38863990  execute( "map F= br@,"+varstr(br@)+";" ); // define the corresponding inclusion of rings
    38873991  ideal I=F(I);

  • Singular/LIB/finvar.lib

    ra91002 r6f4f9c  
    71927192    int m=ncols(F);
    71937193    string mp=string(minpoly);
    7194     execute("ring R=("+charstr(br)+"),("+varstr(br)+",y(1..m)),dp;");
     7194    list l2 = ringlist(br)[2];
     7195    for (int ii = 1; ii <= m; ii++)
     7196    {
     7197     l2[size(l2)+1] = "y("+string(ii)+")";
     7198    }
     7199    ring R = create_ring(ringlist(br)[1], l2, "dp", "no_minpoly");
    71957200    if (mp!="0")
    71967201    { execute("minpoly=number("+mp+");"); }
    … …  
    74147419    int m=ncols(F);
    74157420    string mp=string(minpoly);
    7416     execute("ring R=("+charstr(br)+"),("+varstr(br)+",y(1..m)),lp;");
     7421    list l2 = ringlist(br)[2];
     7422    for (int ii = 1; ii <= m; ii++)
     7423    {
     7424     l2[size(l2)+1] = "y("+string(ii)+")";
     7425    }
     7426    ring R = create_ring(ringlist(br)[1], l2, "lp", "no_minpoly");
    74177427    if (mp!="0")
    74187428    { execute("minpoly=number("+mp+");"); }

  • Singular/LIB/graal.lib

    ra91002 r6f4f9c  
    404404  string ostring = ordstr(basering);
    405405  ostring = ostring[1..size(ostring)-2];
    406   execute("ring Q0y = ("+charstr(basering)+"),(Y(1..s),"+varstr(basering)+"),(ds(s),c,"+ostring+")");
     406  list l2;
     407  for (int ii = 1; ii <= s; ii++)
     408  {
     409   l2[ii] = "Y("+string(ii)+")";
     410  }
     411  l2[s+1] = ringlist(basering)[2];
     412  ring Q0y = create_ring(ringlist(basering)[1], l2, "(ds(s),c,"+ostring+")", "no_minpoly");
    407413  setring Q0y;
    408414  ideal H = imap(Gr.A,H);

  • Singular/LIB/hnoether.lib

    ra91002 r6f4f9c  
    27362736   string strmip=string(minpoly);
    27372737   string strf=string(f);
    2738    execute("ring tempr=("+string(char(basering))+","+parstr(basering)+"),("
    2739            +varstr(basering)+"),dp;");
     2738   ring tempr = create_ring("("+string(char(basering))+","+parstr(basering)+")", "("+varstr(basering)+")", "dp");
    27402739   execute("minpoly="+strmip+";");
    27412740   execute("poly f="+strf+";");

  • Singular/LIB/ncalg.lib

    ra91002 r6f4f9c  
    595595  int @d = nrows(A);
    596596  int n = ncols(A);
    597   execute("ring r1=0,(d(1..n)),"+sord+";");
     597  list l2;
     598  for (int ii = 1; ii <= n; ii++)
     599  {
     600   l2[ii] = "d("+string(ii)+")";
     601  }
     602  ring r1 = create_ring(0, l2, sord);
    598603  ideal I0;
    599604  if (size(#)==0)
    … …  
    615620  }
    616621  string sord2 = "(a(0:" + string(n) + ")," + sord + ")";
    617   execute("ring GR = (0,b(1..@d)),(x(1..n),d(1..n)),"+sord2+";");
     622  string l3 = "(0";
     623  int ii;
     624  for (ii = 1; ii <= @d; ii++)
     625  {
     626   l3 = l3+",b("+string(ii)+")";
     627  }
     628  l3 = l3+")";
     629  list l4;
     630  for (ii = 1; ii <= n; ii++)
     631  {
     632   l4[ii] = "x("+string(ii)+")";
     633   l4[n+ii] = "d("+string(ii)+")";
     634  }
     635  ring GR = create_ring(l3, l4, sord2);
    618636  def W=Weyl(); setring W;
    619637  kill GR;

  • Singular/LIB/poly.lib

    ra91002 r6f4f9c  
    740740  else
    741741  {
    742     execute("ring R=("+charstr(br)+"),(@z,@x(1..n)),(dp(1),"+ordstr(br)+");");
     742    list l2 = "@z";
     743    for (int ii = 1; ii <= n; ii++)
     744    {
     745     l2[ii+1] = "@x("+string(ii)+")";
     746    }
     747    ring R = create_ring(ringlist(br)[1], l2, "(dp(1),"+ordstr(br)+")", "no_minpoly");
    743748  }
    744749  if (mp!="0")

  • Singular/LIB/resolve.lib

    ra91002 r6f4f9c  
    666666//---------------------------------------------------------------------------
    667667//--- control the names of the variables
    668    execute("ring R=("+charstr(R0)+"),(x(1..n)),dp;");
     668   list l2;
     669   for (int ii = 1; ii <= n; ii++)
     670   {
     671    l2[ii] = "x("+string(ii)+")";
     672   }
     673   ring R = create_ring(ringlist(R0)[1], l2, "dp", "no_minpoly");
    669674   list BO=fetch(R1,BO);
    670675   ideal C=fetch(R1,C);
    … …  
    706711         //to have the list depending on R in case BO[4] is empty
    707712
    708    execute("ring S=("+charstr(R)+"),("+varstr(R)+",y(0..s-1)),dp;");
     713   list l2 = ringlist(R)[2];
     714   for (int ii = 0; ii <= s-1; ii++)
     715   {
     716    l2[size(l2)+1] = "y("+string(ii)+")";
     717   }
     718   ring S = create_ring(ringlist(R)[1], l2, "dp", "no_minpoly");
    709719   list resu;
    710720   list B;
    711    execute("ring T=("+charstr(R)+"),("+varstr(R)+",y(0..s-1),t),dp;");
     721   list l2 = ringlist(R)[2];
     722   for (int ii = 0; ii <= s-1; ii++)
     723   {
     724    l2[size(l2)+1] = "y("+string(ii)+")";
     725   }
     726   l2[size(l2)+1] = "t";
     727   ring S = create_ring(ringlist(R)[1], l2, "dp", "no_minpoly");
    712728   ideal C=imap(R,C);
    713729   ideal W=imap(R,W);
Note: See TracChangeset for help on using the changeset viewer.