Changeset 85c8c60 in git


Ignore:
Timestamp:
Jun 3, 2020, 10:54:58 PM (4 years ago)
Author:
Sachin <sachinkm308@…>
Branches:
(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
Children:
e71151310cf4ea51a0b407491e6e91880cf38900
Parents:
6645dd011a1c3421072bb043781cd8f48ac790f2
git-author:
Sachin <sachinkm308@gmail.com>2020-06-03 22:54:58+02:00
git-committer:
Sachin <sachinkm308@gmail.com>2020-06-06 10:23:12+02:00
Message:
replacing execute with create_ring(16)
Location:
Singular/LIB
Files:
8 edited

Legend:

Unmodified
Added
Removed
  • Singular/LIB/atkins.lib

    r6645dd r85c8c60  
    567567          s3=string(t);
    568568
    569           execute("a1="+s1+";");
    570           execute("b1="+s2+";");
    571           execute("t1="+s3+";");
     569          a1=int(a);
     570          b1=int(b);
     571          t1=int(t);
    572572          g1=gcd(gcd(a1,b1),t1 div a1);
    573573          setring C;
  • Singular/LIB/chern.lib

    r6645dd r85c8c60  
    201201  }
    202202  ring r@ = create_ring(ringlist(basering)[1], l2, "dp", "no_minpoly");
    203   execute( "map F= br@,"+varstr(br@)+";" ); // define the corresponding inclusion of rings
    204   ideal V=F(V);
    205   poly f=F(f);
     203  ideal V = fetch(br@,V);// define the corresponding inclusion of rings
     204  poly f = fetch(br@,f);
    206205  int i;
    207206  for(i=1; i<=nV; i++)
     
    705704    }
    706705    ring r@ = create_ring(ringlist(basering)[1], l3, "lp", "no_minpoly");
    707     execute( "map F= br@,"+varstr(br@)+";" ); // define the corresponding inclusion of rings
    708     list c=F(c); // embed c in the bigger ring
     706   // execute( "map F= br@,"+varstr(br@)+";" ); // define the corresponding inclusion of rings
     707   // list c=F(c); // embed c in the bigger ring
     708    list c= fetch(br@,c);
    709709    poly rez;
    710710    list A=a@(1..n);
     
    721721    rez = -subst(I[1], c@, 0);
    722722    setring br@; // come back to the initial base ring
    723     execute( "map FF= r@,0,"+varstr(br@)+";" ); // define the specialization homomorphism
    724     poly rez=FF(rez); // bring the result to the base ring
     723    //execute( "map FF= r@,0,"+varstr(br@)+";" ); // define the specialization homomorphism
     724    poly rez=fetch(r@,rez); // bring the result to the base ring
    725725    return( (1/factorial(n))*rez);
    726726  }
  • Singular/LIB/dmod.lib

    r6645dd r85c8c60  
    459459  tmp = 0;
    460460  // block ord (dp(N),dp);
    461   s = "iv=";
    462   for (i=1; i<=Nnew-1; i++)
    463   {
    464     s = s+"1,";
    465   }
    466   s[size(s)]=";";
    467   execute(s);
     461  iv = 1:(Nnew-1);
    468462  tmp[1] = "dp";  // string
    469463  tmp[2] = iv;   // intvec
     
    471465  // continue with dp 1,1,1,1...
    472466  tmp[1] = "dp";  // string
    473   s[size(s)] = ",";
    474   s = s+"1;";
    475   execute(s);
     467  iv = 1:Nnew;
    476468  kill s;
    477469  kill NName;
     
    601593  // continue with dp 1,1,1,1...
    602594  tmp[1]  = "dp"; // string
    603   s       = "iv=";
    604   for (i=1; i<=Nnew; i++)
    605   {
    606     s = s+"1,";
    607   }
    608   s[size(s)]= ";";
    609   execute(s);
     595  iv=1:Nnew;
    610596  kill s;
    611597  tmp[2]    = iv;
     
    670656  L[2] = NName;
    671657  // dp ordering;
    672   string s = "iv=";
    673   for (i=1; i<=Nnew; i++)
    674   {
    675     s = s+"1,";
    676   }
    677   s[size(s)] = ";";
    678   execute(s);
     658  iv=1:Nnew;
    679659  kill s;
    680660  tmp     = 0;
     
    711691  export LD;
    712692  return(@R2);
     693  string  strr1 = "Printing from ALTannfsBM";     // comment by sachin
     694  strr1;    // comment by sachin 
    713695}
    714696example
     
    791773  // continue with dp 1,1,1,1...
    792774  tmp[1]  = "dp"; // string
    793   s       = "iv=";
    794   for (i=1; i<=Nnew; i++)
    795   {
    796     s = s+"1,";
    797   }
    798   s[size(s)]= ";";
    799   execute(s);
     775  iv=1:(Nnew-1);
    800776  kill s;
    801777  tmp[2]    = iv;
     
    859835  tmp = 0;
    860836  // block ord (dp(N),dp);
    861   string s = "iv=";
    862   for (i=1; i<=Nnew-1; i++)
    863   {
    864     s = s+"1,";
    865   }
    866   s[size(s)]=";";
    867   execute(s);
     837  iv=1:(Nnew-1);
    868838  tmp[1] = "dp";  //string
    869839  tmp[2] = iv;    //intvec
     
    871841  // continue with dp 1,1,1,1...
    872842  tmp[1] = "dp";  //string
    873   s[size(s)] = ",";
    874   s = s+"1;";
    875   execute(s);
    876   kill s;
     843  iv=1:Nnew;
    877844  kill NName;
    878845  tmp[2]  = iv;
     
    1021988  // continue with dp 1,1,1,1...
    1022989  tmp[1]  = "dp"; // string
    1023   s       = "iv=";
    1024   for (i=1; i<=Nnew; i++)
    1025   {
    1026     s = s+"1,";
    1027   }
    1028   s[size(s)]= ";";
    1029   execute(s);
     990  iv=1:Nnew;
    1030991  kill s;
    1031992  tmp[2]    = iv;
     
    10891050  tmp = 0;
    10901051  // block ord (dp(N),dp);
    1091   string s = "iv=";
    1092   for (i=1; i<=Nnew-1; i++)
    1093   {
    1094     s = s+"1,";
    1095   }
    1096   s[size(s)]=";";
    1097   execute(s);
     1052  iv=1:(Nnew-1);
    10981053  tmp[1] = "dp";  //string
    10991054  tmp[2] = iv;    //intvec
     
    11011056  // continue with dp 1,1,1,1...
    11021057  tmp[1] = "dp";  //string
    1103   s[size(s)] = ",";
    1104   s = s+"1;";
    1105   execute(s);
    1106   kill s;
     1058  iv=1:Nnew;
    11071059  kill NName;
    11081060  tmp[2]  = iv;
     
    11851137  L[2] = NName;
    11861138  // dp ordering;
    1187   string s = "iv=";
    1188   for (i=1; i<=Nnew; i++)
    1189   {
    1190     s = s+"1,";
    1191   }
    1192   s[size(s)] = ";";
    1193   execute(s);
    1194   kill s;
     1139  iv=1:Nnew;
    11951140  tmp     = 0;
    11961141  tmp[1]  = "dp";  //string
     
    13441289  L[2] = NName;
    13451290  // dp ordering;
    1346   string s = "iv=";
    1347   for (i=1; i<=Nnew; i++)
    1348   {
    1349     s = s+"1,";
    1350   }
    1351   s[size(s)] = ";";
    1352   execute(s);
     1291  iv=1:Nnew;
    13531292  tmp     = 0;
    13541293  tmp[1]  = "dp";  // string
    13551294  tmp[2]  = iv;  // intvec
    13561295  Lord[1] = tmp;
    1357   kill s;
    13581296  tmp[1]  = "C";
    13591297  iv = 0;
     
    15391477  L[2] = NName;
    15401478  // dp ordering;
    1541   string s = "iv=";
    1542   for (i=1; i<=Nnew; i++)
    1543   {
    1544     s = s+"1,";
    1545   }
    1546   s[size(s)] = ";";
    1547   execute(s);
     1479  iv=1:Nnew;
    15481480  tmp     = 0;
    15491481  tmp[1]  = "dp";  // string
    15501482  tmp[2]  = iv;  // intvec
    15511483  Lord[1] = tmp;
    1552   kill s;
    15531484  tmp[1]  = "C";
    15541485  iv = 0;
     
    16751606  // continue with dp 1,1,1,1...
    16761607  tmp[1]  = "dp"; // string
    1677   s       = "iv=";
    1678   for (i=1; i<=Nnew; i++)
    1679   {
    1680     s = s+"1,";
    1681   }
    1682   s[size(s)]= ";";
    1683   execute(s);
     1608  iv=1:Nnew;
    16841609  kill s;
    16851610  tmp[2]    = iv;
     
    17431668  tmp = 0;
    17441669  // block ord (dp(N),dp);
    1745   string s = "iv=";
    1746   for (i=1; i<=Nnew-1; i++)
    1747   {
    1748     s = s+"1,";
    1749   }
    1750   s[size(s)]=";";
    1751   execute(s);
     1670  iv=1:(Nnew-1);
    17521671  tmp[1] = "dp";  //string
    17531672  tmp[2] = iv;    //intvec
     
    17551674  // continue with dp 1,1,1,1...
    17561675  tmp[1] = "dp";  //string
    1757   s[size(s)] = ",";
    1758   s = s+"1;";
    1759   execute(s);
    1760   kill s;
     1676  iv=1:Nnew;
    17611677  kill NName;
    17621678  tmp[2]  = iv;
     
    18401756  L[2] = NName;
    18411757  // dp ordering;
    1842   string s = "iv=";
    1843   for (i=1; i<=Nnew; i++)
    1844   {
    1845     s = s+"1,";
    1846   }
    1847   s[size(s)] = ";";
    1848   execute(s);
    1849   kill s;
     1758  iv=1:Nnew;
    18501759  tmp     = 0;
    18511760  tmp[1]  = "dp";  //string
     
    21292038  // continue with dp 1,1,1,1...
    21302039  tmp[1]  = "dp"; // string
    2131   s       = "iv=";
    2132   for (i=1; i<=Nnew; i++)
    2133   {
    2134     s = s+"1,";
    2135   }
    2136   s[size(s)]= ";";
    2137   execute(s);
     2040  iv=1:Nnew;
    21382041  kill s;
    21392042  tmp[2]    = iv;
     
    22232126  L[2] = NName;
    22242127  // dp ordering;
    2225   string s = "iv=";
    2226   for (i=1; i<=Nnew; i++)
    2227   {
    2228     s = s+"1,";
    2229   }
    2230   s[size(s)] = ";";
    2231   execute(s);
    2232   kill s;
     2128  iv=1:Nnew;
    22332129  tmp     = 0;
    22342130  tmp[1]  = "dp";  //string
     
    24002296  // block ord (lp(P),dp);
    24012297  tmp[1] = "lp";  //string
    2402   s      = "iv=";
    2403   for (i=1; i<=2*P; i++)
    2404   {
    2405     s = s+"1,";
    2406   }
    2407   s[size(s)]= ";";
    2408   execute(s);
     2298  iv=1:(2*P);
    24092299  tmp[2] = iv;  //intvec
    24102300  Lord[1] = tmp;
    24112301  // continue with dp 1,1,1,1...
    24122302  tmp[1] = "dp";  //string
    2413   s      = "iv=";
    2414   for (i=1; i<=Nnew; i++)  //actually i<=2*N
    2415   {
    2416     s = s+"1,";
    2417   }
    2418   s[size(s)]= ";";
    2419   execute(s);
     2303  iv=1:Nnew;
    24202304  kill s;
    24212305  tmp[2]  = iv;
     
    25292413  tmp = 0;
    25302414  // block ord (dp(N),dp);
    2531   string s = "iv=";
    2532   for (i=1; i<=Nnew-P; i++)
    2533   {
    2534     s = s+"1,";
    2535   }
    2536   s[size(s)]=";";
    2537   execute(s);
     2415  iv=1:(Nnew-P);
    25382416  tmp[1] = "dp";  //string
    25392417  tmp[2] = iv;    //intvec
     
    25412419  // continue with dp 1,1,1,1...
    25422420  tmp[1] = "dp";  //string
    2543   s[size(s)] = ",";
    2544   for (j=1; j<=P; j++)
    2545   {
    2546     s = s+"1,";
    2547   }
    2548   s[size(s)]=";";
    2549   execute(s);
    2550   kill s;
     2421  iv=1:P;
    25512422  kill NName;
    25522423  tmp[2]  = iv;
     
    26602531  tmp = 0;
    26612532  // dp ordering;
    2662   string s = "iv=";
    2663   for (i=1; i<=Nnew; i++)
    2664   {
    2665     s = s+"1,";
    2666   }
    2667   s[size(s)]=";";
    2668   execute(s);
    2669   kill s;
     2533  iv=1:Nnew;
    26702534  kill NName;
    26712535  tmp[1] = "dp";  //string
     
    28152679  // continue with dp 1,1,1,1...
    28162680  tmp[1]  = "dp"; // string
    2817   s       = "iv=";
    2818   for(i=1;i<=Nnew;i++)
    2819   {
    2820     s = s+"1,";
    2821   }
    2822   s[size(s)]= ";";
    2823   execute(s);
     2681  iv=1:Nnew;
    28242682  tmp[2]    = iv;
    28252683  Lord[2]   = tmp;
     
    28922750  // continue with a(1,1,1,1)...
    28932751  tmp[1]  = "dp"; s  = "iv=";
    2894   for(i=1; i<= Nnew; i++)
    2895   {
    2896     s = s+"1,";
    2897   }
    2898   s[size(s)]= ";";  execute(s);
     2752  iv=1:Nnew;
    28992753  kill NName;
    29002754  tmp[2]    = iv;
     
    29552809  tmp    = 0;
    29562810  // block ord (a(1,1...),dp);
    2957   string  s = "iv=";
    2958   for(i=1; i<=Nnew-1; i++)
    2959   {
    2960     s = s+"1,";
    2961   }
    2962   s[size(s)]= ";";
    2963   execute(s);
     2811  iv=1:(Nnew-1);
    29642812  tmp[1]  = "a"; // string
    29652813  tmp[2]  = iv; //intvec
     
    29672815  // continue with dp 1,1,1,1...
    29682816  tmp[1]  = "dp"; // string
    2969   s[size(s)]=","; s= s+"1;";
    2970   execute(s);
    2971   kill s;
    29722817  kill NName;
    29732818  tmp[2]    = iv;
     
    30512896  L[2]   = NName;
    30522897  // dp ordering;
    3053   string   s       = "iv=";
    3054   for(i=1;i<=Nnew;i++)
    3055   {
    3056     s = s+"1,";
    3057   }
    3058   s[size(s)]= ";";
    3059   execute(s);
    30602898  tmp     = 0;
    30612899  tmp[1]  = "dp"; // string
    30622900  tmp[2]  = iv; //intvec
    30632901  Lord[1] = tmp;
    3064   kill s;
    30652902  tmp[1]    = "C";
    30662903  iv        = 0;
     
    31933030  // continue with dp 1,1,1,1...
    31943031  tmp[1]  = "dp"; // string
    3195   s       = "iv=";
    3196   for(i=1;i<=Nnew;i++)
    3197   {
    3198     s = s+"1,";
    3199   }
    3200   s[size(s)]= ";";
    3201   execute(s);
     3032  iv=1:Nnew;
    32023033  tmp[2]    = iv;
    32033034  Lord[2]   = tmp;
     
    32703101  // continue with a(1,1,1,1)...
    32713102  tmp[1]  = "dp"; s  = "iv=";
    3272   for(i=1; i<= Nnew; i++)
    3273   {
    3274     s = s+"1,";
    3275   }
    3276   s[size(s)]= ";";  execute(s);
     3103  iv=1:Nnew;
    32773104  kill NName;
    32783105  tmp[2]    = iv;
     
    33293156  tmp = 0;
    33303157  // block ord (dp(N),dp);
    3331   string s = "iv=";
    3332   for (i=1; i<=Nnew-1; i++)
    3333   {
    3334     s = s+"1,";
    3335   }
    3336   s[size(s)]=";";
    3337   execute(s);
     3158  iv=1:(Nnew-1);
    33383159  tmp[1] = "dp";  // string
    33393160  tmp[2] = iv;   // intvec
     
    33413162  // continue with dp 1,1,1,1...
    33423163  tmp[1] = "dp";  // string
    3343   s[size(s)] = ",";
    3344   s = s+"1;";
    3345   execute(s);
    3346   kill s;
     3164  iv=1:Nnew;
    33473165  kill NName;
    33483166  tmp[2]      = iv;
     
    33653183  dbprint(ppl-1, @R3);
    33663184  ideal MM = maxideal(1);
     3185  string  strr1 = "printing s";     // comment by sachin
     3186strr1;    // comment by sachin
     3187  s;
    33673188  MM = 0,s,MM;
    33683189  map R01 = @R2, MM;
     
    34653286  // continue with dp 1,1,1,1...
    34663287  tmp[1]  = "dp"; // string
    3467   s       = "iv=";
    3468   for(i=1;i<=Nnew;i++)
    3469   {
    3470     s = s+"1,";
    3471   }
    3472   s[size(s)]= ";";
    3473   execute(s);
     3288  iv=1:Nnew;
    34743289  kill s;
    34753290  tmp[2]    = iv;
     
    35373352  tmp = 0;
    35383353  // block ord (dp(N),dp);
    3539   string s = "iv=";
    3540   for (i=1; i<=Nnew-1; i++)
    3541   {
    3542     s = s+"1,";
    3543   }
    3544   s[size(s)]=";";
    3545   execute(s);
     3354  iv=1:(Nnew-1);
    35463355  tmp[1] = "dp";  // string
    35473356  tmp[2] = iv;   // intvec
     
    35493358  // continue with dp 1,1,1,1...
    35503359  tmp[1] = "dp";  // string
    3551   s[size(s)] = ",";
    3552   s = s+"1;";
    3553   execute(s);
    3554   kill s;
     3360  iv=1:Nnew;
    35553361  kill NName;
    35563362  tmp[2]      = iv;
     
    39813787  // continue with dp 1,1,1,1...
    39823788  tmp[1]  = "dp"; // string
    3983   s       = "iv=";
    3984   for(i=1;i<=Nnew;i++)
    3985   {
    3986     s = s+"1,";
    3987   }
    3988   s[size(s)]= ";";
    3989   execute(s);
     3789  iv=1:Nnew;
    39903790  kill s;
    39913791  tmp[2]    = iv;
     
    40543854  // just dp
    40553855  // block ord (dp(N),dp);
    4056   string s = "iv=";
    4057   for (i=1; i<=Nnew; i++)
    4058   {
    4059     s = s+"1,";
    4060   }
    4061   s[size(s)]=";";
    4062   execute(s);
     3856  iv=1:Nnew;
    40633857  tmp[1] = "dp";  // string
    40643858  tmp[2] = iv;   // intvec
     
    42184012  // continue with dp 1,1,1,1...
    42194013  tmp[1]  = "dp"; // string
    4220   s       = "iv=";
    4221   for(i=1;i<=Nnew;i++)
    4222   {
    4223     s = s+"1,";
    4224   }
    4225   s[size(s)]= ";";
    4226   execute(s);
     4014  iv=1:Nnew;
    42274015  tmp[2]    = iv;
    42284016  Lord[3]   = tmp;
     
    42974085  // block ord (dp(N),dp);
    42984086  // string s is already defined
    4299   s = "iv=";
    4300   for (i=1; i<=Nnew-1; i++)
    4301   {
    4302     s = s+"1,";
    4303   }
    4304   s[size(s)]=";";
    4305   execute(s);
     4087  iv=1:(Nnew-1);
    43064088  tmp[1] = "dp";  // string
    43074089  tmp[2] = iv;   // intvec
     
    43094091  // continue with dp 1,1,1,1...
    43104092  tmp[1] = "dp";  // string
    4311   s[size(s)] = ",";
    4312   s = s+"1;";
    4313   execute(s);
     4093  iv=1:Nnew;
    43144094  kill s;
    43154095  kill NName;
     
    44434223  // continue with dp 1,1,1,1...
    44444224  tmp[1]  = "dp"; // string
    4445   s       = "iv=";
    4446   for(i=1;i<=Nnew;i++)
    4447   {
    4448     s = s+"1,";
    4449   }
    4450   s[size(s)]= ";";
    4451   execute(s);
     4225  iv=1:Nnew;
    44524226  tmp[2]    = iv;
    44534227  Lord[2]   = tmp;
     
    45224296  // block ord (dp(N),dp);
    45234297  // string s is already defined
    4524   s = "iv=";
    4525   for (i=1; i<=Nnew-1; i++)
    4526   {
    4527     s = s+"1,";
    4528   }
    4529   s[size(s)]=";";
    4530   execute(s);
     4298  iv=1:(Nnew-1);
    45314299  tmp[1] = "dp";  // string
    45324300  tmp[2] = iv;   // intvec
     
    45344302  // continue with dp 1,1,1,1...
    45354303  tmp[1] = "dp";  // string
    4536   s[size(s)] = ",";
    4537   s = s+"1;";
    4538   execute(s);
    4539   kill s;
     4304  iv=1:Nnew;
     4305  kill s;
    45404306  kill NName;
    45414307  tmp[2]      = iv;
     
    47274493  L[2] = NName;
    47284494  // dp ordering;
    4729   string s = "iv=";
    4730   for (i=1; i<=Nnew; i++)
    4731   {
    4732     s = s+"1,";
    4733   }
    4734   s[size(s)] = ";";
    4735   execute(s);
     4495  iv=1:Nnew;
    47364496  tmp     = 0;
    47374497  tmp[1]  = "dp";  // string
    47384498  tmp[2]  = iv;  // intvec
    47394499  Lord[1] = tmp;
    4740   kill s;
    47414500  tmp[1]  = "C";
    47424501  iv = 0;
     
    53675126  // block ord (dp);
    53685127  tmp[1] = "dp"; // string
    5369   s = "iv=";
    5370   for (i=1; i<=Nnew; i++)
    5371   {
    5372     s = s+"1,";
    5373   }
    5374   s[size(s)]=";";
    5375   execute(s);
     5128  iv=1:Nnew;
    53765129  kill s;
    53775130  tmp[2]  = iv;
     
    55025255  // block ord (dp);
    55035256  tmp[1] = "dp"; // string
    5504   s = "iv=";
    5505   for (i=1; i<=N; i++)
    5506   {
    5507     s = s+"1,";
    5508   }
    5509   s[size(s)]=";";
    5510   execute(s);
     5257  iv=1:N;
    55115258  kill s;
    55125259  tmp[2]  = iv;
     
    57995546  tmp=0;
    58005547  // dp ordering;
    5801   string s = "iv=";
    5802   for (i=1; i<=Nnew; i++)
    5803   {
    5804     s = s+"1,";
    5805   }
    5806   s[size(s)] = ";";
    5807   execute(s);
     5548  iv=1:Nnew;
    58085549  tmp     = 0;
    58095550  tmp[1]  = "dp";  // string
    58105551  tmp[2]  = iv;  // intvec
    58115552  Lord[2] = tmp;
    5812   kill s;
    58135553  tmp     = 0;
    58145554  L[3]    = Lord;
  • Singular/LIB/finitediff.lib

    r6645dd r85c8c60  
    583583  svars=svars+","+strvar;            ////variables
    584584  string sord="(c,lp)";         ////ordering
    585   string sring="ring Q="+scf+","+svars+","+sord+";";
    586   execute(sring);
     585  ring Q = create_ring(scf, svars, sord);
    587586  ideal Id=i**2+1;
    588587  int j;
  • Singular/LIB/finvar.lib

    r6645dd r85c8c60  
    719719        execute
    720720        ("ring newring=(0"+chst+"),("+varstr(br)+"),("+ordstr(br)+")");
    721         string  strr2 = "First else";     // comment by sachin
    722         strr2;    // comment by sachin
    723         ringlist(newring);
    724721       
    725722      }
  • Singular/LIB/lejeune.lib

    r6645dd r85c8c60  
    422422     ERROR("n must range between 1 and 26!");
    423423   }
    424   string s="ring r=0,("+A_Z("a",n)+"),ds;";
    425   execute(s);
     424  ring r = create_ring(0, A_Z_L("a",n), "ds");
    426425  return (string(var(n)));
    427426}
  • Singular/LIB/ncrat.lib

    r6645dd r85c8c60  
    716716  string s;
    717717
    718   s = "ring NCRING = (0, I), (";
     718  s = "(";
    719719  for (i = 1; i <= size(NCVARIABLES); i++)
    720720  {
     
    728728    }
    729729  }
    730   s = s + "), dp; minpoly = I^2+1;";
    731   execute(s);
     730  s = s + ")";
     731  ring NCRING = create_ring("(0, I)", s, "dp");
     732  minpoly = I^2+1;
    732733  short = 0;
    733734  exportto(Top, NCRING);
  • Singular/LIB/polymake.lib

    r6645dd r85c8c60  
    617617  // create a ring whose variables are indexed by the points in
    618618  // boundary resp. by interior
    619  string rst="ring cyclering=0,(u"+string(interior[1])+string(interior[2]);
     619 string rst="(u"+string(interior[1])+string(interior[2]);
    620620 for (j=1;j<=size(boundary);j++)
    621621 {
    622622   rst=rst+",u"+string(boundary[j][1])+string(boundary[j][2]);
    623623 }
    624  rst=rst+"),lp;";
    625  execute(rst);
    626 /*string rst="(u"+string(interior[1])+string(interior[2]);
    627   for (j=1;j<=size(boundary);j++)
    628   {
    629    rst=rst+",u"+string(boundary[j][1])+string(boundary[j][2]);
    630   }
    631   rst=rst+")";
    632   ring cyclering = create_ring(0, rst, "lp"); */
     624 rst=rst+")";
     625 ring cyclering = create_ring(0, rst, "lp");
    633626  // add the first and second point at the end of boundary
    634627  boundary[size(boundary)+1]=boundary[1];
Note: See TracChangeset for help on using the changeset viewer.