Changeset 52d301 in git
- Timestamp:
- Nov 4, 2022, 1:41:50 PM (5 months ago)
- Branches:
- (u'spielwiese', 'f875bbaccd0831e36aaed09ff6adeb3eb45aeb94')
- Children:
- 4969bd84679f00e0090770de53eae064004f1330
- Parents:
- 4c83c9834db9cbdfc54c7e750ea85bc9df089d8e
- Location:
- Singular/LIB
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/normal.lib
r4c83c98 r52d301 277 277 { 278 278 // We change to dp ordering. 279 list rl = ring list(origR);279 list rl = ring_list(origR); 280 280 list origOrd = rl[3]; 281 281 list newOrd = list("dp", intvec(1:nvars(origR))), list("C", 0); … … 298 298 // We define anyway a new ring, equal to R, to be able to return it. 299 299 setring origR; 300 list lR = ring list(origR);300 list lR = ring_list(origR); 301 301 def ROut = ring(lR); 302 302 setring ROut; … … 422 422 savebasering=basering; 423 423 setring newR; // must be in a compatible ring to newR 424 // as ring list may produce ring-dep. stuff424 // as ring_list may produce ring-dep. stuff 425 425 if(!defined(newRList)) { list newRList;} 426 newRList = ring list(newR);426 newRList = ring_list(newR); 427 427 setring savebasering; 428 428 U = norComp[j][1]; … … 908 908 else 909 909 { 910 list RL = ring list(newRing);910 list RL = ring_list(newRing); 911 911 def lastRing = ring(RL); 912 912 setring lastRing; … … 1176 1176 int ti=timer; 1177 1177 1178 gnirlist = ring list(BAS);1178 gnirlist = ring_list(BAS); 1179 1179 1180 1180 //----------- the trivial case of a zero ideal as input, RETURN ------------ … … 2324 2324 { 2325 2325 list newvar=N[1..ncols(N)]; 2326 list @rl0=ring list(R0);2326 list @rl0=ring_list(R0); 2327 2327 @rl0[2]=list(N[1..ncols(N)]); 2328 2328 @rl0[3]=list(list("dp",1:ncols(N)),list("C",0)); … … 2330 2330 setring R; 2331 2331 ideal I =imap(R0,I); 2332 list @rl0=ring list(R0);2332 list @rl0=ring_list(R0); 2333 2333 if((@rl0[3][1][1]=="dp")&&(size(@rl0[3])==2)) 2334 2334 { … … 2354 2354 else 2355 2355 { 2356 list @rl0=ring list(R0);2356 list @rl0=ring_list(R0); 2357 2357 if(dim(I)==0) 2358 2358 { … … 2695 2695 option(redSB); 2696 2696 def R=basering; 2697 list @rl=ring list(R);2697 list @rl=ring_list(R); 2698 2698 @rl[2]=list("x","y"); 2699 2699 @rl[3]=list(list("lp",1:2),list("C",0)); … … 2987 2987 2988 2988 int dblvl = printlevel-voice+2; 2989 list gnirlist = ring list(basering);2989 list gnirlist = ring_list(basering); 2990 2990 2991 2991 //---- Some auxiliary variables: … … 3012 3012 { 3013 3013 def R = L[1]; 3014 list Rlist = ring list(R);3014 list Rlist = ring_list(R); 3015 3015 def BAS = basering; 3016 3016 setring R; … … 3022 3022 ker=simplify(interred(ker),15); 3023 3023 // Rlist may be not defined in this new ring, so we define it again. 3024 list Rlist2 = ring list(R);3024 list Rlist2 = ring_list(R); 3025 3025 def R0 = ring(Rlist2); 3026 3026 setring R0; … … 3323 3323 { 3324 3324 //### noch abfragen ob Z(i) definiert ist 3325 list gnirlist = ring list(R(k));3325 list gnirlist = ring_list(R(k)); 3326 3326 n2 = size(gnirlist[2]); 3327 3327 n3 = size(gnirlist[3]); … … 3500 3500 { 3501 3501 //### noch abfragen ob Z(i) definiert ist 3502 list gnirlist = ring list(R(k));3502 list gnirlist = ring_list(R(k)); 3503 3503 n2 = size(gnirlist[2]); 3504 3504 n3 = size(gnirlist[3]); … … 3733 3733 result[2]=result[1]; 3734 3734 def BR = basering; 3735 list rl = ring list(basering);3735 list rl = ring_list(basering); 3736 3736 def rng = ring(rl); 3737 3737 setring rng; … … 4045 4045 } 4046 4046 4047 list gnirlist = ring list(basering);4047 list gnirlist = ring_list(basering); 4048 4048 endid = Le[1]; 4049 4049 //endphi;""; … … 4194 4194 def R = basering; 4195 4195 poly c = J[1]; // the denominator 4196 list gnirlist = ring list(basering);4196 list gnirlist = ring_list(basering); 4197 4197 string svars = varstr(basering); 4198 4198 int nva = nvars(basering); … … 4538 4538 // We define anyway a new ring, equal to R, to be able to return it. 4539 4539 setring R; 4540 list lR = ring list(R);4540 list lR = ring_list(R); 4541 4541 def ROut = ring(lR); 4542 4542 setring ROut; … … 4669 4669 { 4670 4670 // We change to global dp ordering. 4671 list rl = ring list(R);4671 list rl = ring_list(R); 4672 4672 list origOrd = rl[3]; 4673 4673 list newOrd = list("dp", intvec(1:nvars(R))), list("C", 0); … … 4871 4871 // We define anyway a new ring, equal to R, to be able to return it. 4872 4872 setring R; 4873 list lR = ring list(R);4873 list lR = ring_list(R); 4874 4874 def ROut = ring(lR); 4875 4875 setring ROut; … … 5065 5065 5066 5066 def R = basering; 5067 list RL = ring list(R);5067 list RL = ring_list(R); 5068 5068 // We keep everything from R but we change the ordering to lp, and we 5069 5069 // order the variables as needed. … … 5219 5219 if(isGlobal != 1) 5220 5220 { 5221 list rl = ring list(origEre);5221 list rl = ring_list(origEre); 5222 5222 list origOrd = rl[3]; 5223 5223 list newOrd = list("dp", intvec(1:nvars(origEre))), list("C", 0); … … 6783 6783 6784 6784 //--------- create new ring with one dp block keeping weights ------------ 6785 list LR = ring list(R);6785 list LR = ring_list(R); 6786 6786 list g3 = LR[3]; 6787 6787 int n3 = size(g3); … … 7164 7164 //### minor verbessern 7165 7165 //### preimage verbessern (Ist imm Kern map oder imap verwendet?) 7166 //### Gleich in Ordnung dp wechseln, ring list verwenden7166 //### Gleich in Ordnung dp wechseln, ring_list verwenden 7167 7167 //### interred ev nur zum Schluss 7168 7168 // (z.B. wenn nacher std; wenn nacher minor: testen ) -
Singular/LIB/paraplanecurves.lib
r4c83c98 r52d301 1 1 ////////////////////////////////////////////////////////////////////////////// 2 version="version paraplanecurves.lib 4. 1.2.0 Feb_2019"; // $Id$2 version="version paraplanecurves.lib 4.3.2.2 Nov_2022 "; // $Id$ 3 3 category="Algebraic Geometry"; 4 4 info=" … … 130 130 int m = size(phi); 131 131 /*phi: P^(n-1) --> P^(m-1)*/ 132 list rl = ring list(Roriginal);132 list rl = ring_list(Roriginal); 133 133 int k; 134 134 for(k = 1; k <= n; k++) … … 174 174 } 175 175 relpsi = compress(relpsi); 176 list rl = ring list(Rbig);176 list rl = ring_list(Rbig); 177 177 list rl2 = rl[2]; 178 178 rl[2] = list(rl2[n+1..n+m]); … … 300 300 /*checking assumptions and handling the conic case*/ 301 301 checkAssumptions(f); 302 list rl = ring list(Roriginal);302 list rl = ring_list(Roriginal); 303 303 rl[2] = list("x","y","z"); 304 304 rl[3] = list(list("dp",1:3),list("C",0)); … … 349 349 def rp1 = invertBirMap(PROJ,If); // ring containing ideal psi 350 350 setring rp1; 351 list rl1 = ring list(rp1);351 list rl1 = ring_list(rp1); 352 352 rl1[2] = list("s","t"); 353 353 def RP1 = ring(rl1); … … 391 391 def rp1 = invertBirMap(PROJ,If); // (ring, (J,psi)) 392 392 setring rp1; 393 list rl1 = ring list(rp1);393 list rl1 = ring_list(rp1); 394 394 rl1[2] = list("s","t"); 395 395 def RP1 = ring(rl1); … … 630 630 { 631 631 def Roriginal = basering; 632 list rl = ring list(Roriginal);632 list rl = ring_list(Roriginal); 633 633 rl[3] = list(list("dp",1:nvars(Roriginal)),list("C",0)); 634 634 def RP2 = ring(rl); … … 663 663 } 664 664 int j,k; 665 list rl = ring list(RP2);665 list rl = ring_list(RP2); 666 666 rl[2] = list(var(1), var(2)); 667 667 rl[3] = list(list("dp",1:2),list("C",0)); … … 794 794 poly dhf = subst(f,var(3),1); 795 795 def Roriginal = basering; 796 list rl = ring list(Roriginal);796 list rl = ring_list(Roriginal); 797 797 rl[2] = list(var(1), var(2)); 798 798 rl[3] = list(list("dp",1:2),list("C",0)); … … 803 803 list DATA = imap(Roriginal,choices); 804 804 /* Creating rings for later use */ 805 list rl = ring list(Rdummy);805 list rl = ring_list(Rdummy); 806 806 rl[2] = list(var(2), var(1)); 807 807 rl[3] = list(list("lp",1:2),list("C",0)); 808 808 def Rred = ring(rl); // make var(2) > var(1) 809 rl = ring list(Rdummy);809 rl = ring_list(Rdummy); 810 810 rl[1] = list(0,list(var(1)),list(list("dp",1)),ideal(0)); 811 811 rl[2] = list(var(2)); … … 901 901 poly dhf = subst(f,var(3),1); 902 902 def Roriginal = basering; 903 list rl = ring list(Roriginal);903 list rl = ring_list(Roriginal); 904 904 rl[2] = list(var(1), var(2)); 905 905 rl[3] = list(list("dp",1:2),list("C",0)); … … 955 955 def Roriginal = basering; 956 956 ideal IC = f; 957 list rl = ring list(Roriginal);957 list rl = ring_list(Roriginal); 958 958 /* begin workaround elimination*/ 959 959 for(k = 1; k <= 3; k++) … … 979 979 intvec HJ = hilb(SJ,1); 980 980 ideal RNC = eliminate(J,x(1)*x(2)*x(3),HJ); 981 list rl = ring list(Relim);981 list rl = ring_list(Relim); 982 982 list rl2 = rl[2]; 983 983 rl[2] = list(rl2[4..n+3]); … … 1023 1023 { 1024 1024 def Roriginal = basering; 1025 list rl = ring list(Roriginal);1025 list rl = ring_list(Roriginal); 1026 1026 rl[3] = list(list("dp",1:nvars(Roriginal)),list("C",0)); 1027 1027 def RoriginalDP = ring(rl); … … 1069 1069 if((n mod 2) == 1){ERROR("Pn has even dimension");} 1070 1070 def Roriginal = basering; 1071 list rlo = ring list(Roriginal);1071 list rlo = ring_list(Roriginal); 1072 1072 rlo[3]= list(list("dp",1:n),list("C",0)); 1073 1073 int k; … … 1094 1094 { 1095 1095 ideal PR = rncAntiCanonicalMap(J); 1096 list rl = ring list(Rold);1096 list rl = ring_list(Rold); 1097 1097 Rbig = Rold + RoriginalCopy; 1098 1098 setring Rbig; … … 1133 1133 HJJ = hilb(SJJ,1); 1134 1134 J = eliminate(JJ,pvar,HJJ); 1135 list rl = ring list(Relim);1135 list rl = ring_list(Relim); 1136 1136 rl2 = rl[2]; 1137 1137 rl[2] = list(rl2[n+1..2*n-2]); … … 1185 1185 if((n mod 2) == 0){ERROR("Pn has odd dimension");} 1186 1186 def Roriginal = basering; 1187 list rlo = ring list(Roriginal);1187 list rlo = ring_list(Roriginal); 1188 1188 rlo[3]= list(list("dp",1:n),list("C",0)); 1189 1189 int k; … … 1210 1210 { 1211 1211 ideal PR = rncAntiCanonicalMap(J); 1212 list rl = ring list(Rold);1212 list rl = ring_list(Rold); 1213 1213 Rbig = Rold + RoriginalCopy; 1214 1214 setring Rbig; … … 1249 1249 HJJ = hilb(SJJ,1); 1250 1250 J = eliminate(JJ,pvar,HJJ); 1251 list rl = ring list(Relim);1251 list rl = ring_list(Relim); 1252 1252 rl2 = rl[2]; 1253 1253 rl[2] = list(rl2[n+1..2*n-2]); … … 1308 1308 poly dhf = subst(f,var(3),1); 1309 1309 def Roriginal = basering; 1310 list rl = ring list(Roriginal);1310 list rl = ring_list(Roriginal); 1311 1311 rl[2] = list(var(1), var(2)); 1312 1312 rl[3] = list(list("dp",1:2),list("C",0)); … … 1800 1800 def rp1 = invertBirMap(Ipoint, ideal(fetch(Roriginal,q))); 1801 1801 setring rp1; 1802 list rl = ring list(rp1);1802 list rl = ring_list(rp1); 1803 1803 rl[2] = list("s","t"); 1804 1804 def RP1 = ring(rl); … … 1848 1848 { 1849 1849 def Roriginal = basering; 1850 list rl = ring list(Roriginal);1850 list rl = ring_list(Roriginal); 1851 1851 rl[3] = list(list("dp",1:3),list("C",0)); 1852 1852 def RP20 = ring(rl); … … 1858 1858 // field 1859 1859 setring RP21; 1860 list rl1 = ring list(RP21);1860 list rl1 = ring_list(RP21); 1861 1861 rl1[2] = list("u","v","w"); 1862 1862 rl1[3] = list(list("dp",1:3),list("C",0)); … … 2381 2381 /* X^2 + br*bs^2*Y^2 + cr*cs^2*Z^2 = 0 */ 2382 2382 def bRing = basering; 2383 list L = ring list(bRing);2383 list L = ring_list(bRing); 2384 2384 2385 2385 if (b != 0) … … 2689 2689 /* begin workaround elimination*/ 2690 2690 int k; 2691 list rl = ring list(rTT);2691 list rl = ring_list(rTT); 2692 2692 rl[2] = list("s","t","x","y","z"); 2693 2693 rl[3]= list(list("dp",1:5),list("C",0)); -
Singular/LIB/primdecint.lib
r4c83c98 r52d301 1 1 ////////////////////////////////////////////////////////////////////////////// 2 version="version primdecint.lib 4. 1.2.0 Feb_2019"; // $Id$2 version="version primdecint.lib 4.3.2.2 Nov_2022 "; // $Id$ 3 3 category = "Commutative Algebra"; 4 4 info=" … … 83 83 def R=basering; 84 84 number s; 85 list rl=ring list(R);85 list rl=ring_list(R); 86 86 int i,j,p,m,ex,nu,k_link; 87 87 list P,B,IS; … … 415 415 ideal K; 416 416 def R=basering; 417 list rl=ring list(R);417 list rl=ring_list(R); 418 418 int i,j,p,m; 419 419 list P,B; … … 552 552 ideal K=1; 553 553 def R=basering; 554 list rl=ring list(R);554 list rl=ring_list(R); 555 555 int i,j,p,m; 556 556 list P; … … 645 645 ideal K=1; 646 646 def R=basering; 647 list rl=ring list(R);647 list rl=ring_list(R); 648 648 int i,j,p,m; 649 649 list P; … … 750 750 def R=basering; 751 751 bigint s; 752 list rl=ring list(R);752 list rl=ring_list(R); 753 753 int i,j,p,m,ex; 754 754 list P,IS,B; … … 963 963 //=== decomposition (else) of J in Z/p and maps the result back to the basering 964 964 def R = basering; 965 list rp = ring list(R);965 list rp = ring_list(R); 966 966 rp[1] = p; 967 967 def Rp = ring(rp); … … 1115 1115 ideal I=imap(R,I); 1116 1116 I=stdZ(I); 1117 list rl=ring list(S);1117 list rl=ring_list(S); 1118 1118 rl[1]=0; 1119 1119 def Shelp =ring(rl); … … 1165 1165 C[i+1]=I[i]; 1166 1166 } 1167 list rl=ring list(R);1167 list rl=ring_list(R); 1168 1168 rl[1]=0; 1169 1169 def Shelp =ring(rl); … … 1399 1399 if(size(B)==1){return(M);} 1400 1400 I=std(I); 1401 list rl=ring list(R);1401 list rl=ring_list(R); 1402 1402 if(deg(I[1])==0) 1403 1403 { … … 1438 1438 module I=imap(R,I); 1439 1439 I=std(I); 1440 list rl=ring list(S);1440 list rl=ring_list(S); 1441 1441 rl[1]=0; 1442 1442 def Shelp =ring(rl); … … 1470 1470 C[i]=I[i]; 1471 1471 } 1472 list rl=ring list(R);1472 list rl=ring_list(R); 1473 1473 rl[1]=0; 1474 1474 def Shelp =ring(rl); -
Singular/LIB/ring.lib
r4c83c98 r52d301 1 1 ///////////////////////////////////////////////////////////////////////////// 2 version="version ring.lib 4. 1.2.0 Feb_2019"; // $Id$2 version="version ring.lib 4.3.1.2 Nov_2022 "; // $Id$ 3 3 category="General purpose"; 4 4 info=" … … 917 917 setring OldRing; 918 918 919 list rl = ring list( OldRing);919 list rl = ring_list( OldRing); 920 920 while ( j<=@sized ) 921 921 { -
Singular/LIB/sagbi.lib
r4c83c98 r52d301 1 1 /////////////////////////////////////////////////////////////////////////// 2 version="version sagbi.lib 4. 1.2.0 Feb_2019"; // $Id$2 version="version sagbi.lib 4.3.1.2 Nov_2022 "; // $Id$ 3 3 category="Commutative Algebra"; 4 4 info=" … … 76 76 // one for each new element in ideal algebra ------------- 77 77 setring r; 78 list l = ring list(r);78 list l = ring_list(r); 79 79 for (i=nvars(r)-nvars(br)+1; i<=numTotalAdditionalVars;i++) 80 80 { … … 234 234 //Create the preimage ring K[@y(1),...,@y(m)], where m=ncols(algebra). 235 235 string variableName=uniqueVariableName("@y"); 236 list l = ring list(basering);236 list l = ring_list(basering); 237 237 for (i=1; i<=m;i++) 238 238 { … … 918 918 919 919 string extendVarName=uniqueVariableName("@c"); 920 list l2 = ring list(basering);920 list l2 = ring_list(basering); 921 921 for (i=1; i<=ncols(I);i++) 922 922 { -
Singular/LIB/scheme.lib
r4c83c98 r52d301 1 1 // 2 version="version scheme.lib 1.0.0.0 Apr_2021"; // $Id$2 version="version scheme.lib 4.3.1.2 Nov_2022 "; // $Id$ 3 3 category="Algebraic geometry"; 4 4 info=" … … 753 753 for(i=1;i<=n;i=i+1) 754 754 { 755 list rl=ring list(R);755 list rl=ring_list(R); 756 756 rl[2]=delete(rl[2],i); 757 757 rl[3]=list(list("dp",1:(n-1)),list("C",0)); … … 1702 1702 { ERROR("Wrong number of arguments!!!");} 1703 1703 1704 list rl=ring list(S);1704 list rl=ring_list(S); 1705 1705 list newvars=list(); 1706 1706 for(int k=1; k<=n; k=k+1)
Note: See TracChangeset
for help on using the changeset viewer.