Changeset 921f32 in git
- Timestamp:
- Dec 18, 2018, 4:55:28 PM (5 years ago)
- Branches:
- (u'spielwiese', 'd1b01e9d51ade4b46b745d3bada5c5f3696be3a8')
- Children:
- 76b17c9f4370c2a5339d376eb75d3d06e5da4a72
- Parents:
- b8498566aa4a3ab88d5be05cde0f7e4a5652c890
- Location:
- Singular/LIB
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/fpadim.lib
rb849856 r921f32 136 136 int i; 137 137 for (i = 1; i <= size(L); i++) 138 {if (entryViolation(L[i], lp BlockSize(basering)))138 {if (entryViolation(L[i], lpVarBlockSize(basering))) 139 139 {ERROR("Not allowed monomial/intvec found!");} 140 140 } … … 907 907 // construct the Graph of normal words [Studzinski page 78] 908 908 // construct set of vertices 909 int v = lp BlockSize(basering); int d = lpDegBound(basering);909 int v = lpVarBlockSize(basering); int d = lpDegBound(basering); 910 910 ideal V; poly p,q,w; 911 911 ideal LG = lead(G); … … 1691 1691 EXAMPLE: example lpDHilbert; shows examples 1692 1692 " 1693 {int degbound = lpDegBound(basering);int n = lp BlockSize(basering);1693 {int degbound = lpDegBound(basering);int n = lpVarBlockSize(basering); 1694 1694 if (size(#) > 0){if (typeof(#[1])=="int"){if (#[1] >= 0){degbound = #[1];}}} 1695 1695 if (size(#) > 1){if (typeof(#[1])=="int"){if (#[2] > 0){n = #[2];}}} … … 1733 1733 EXAMPLE: example lpDHilbertSickle; shows examples 1734 1734 " 1735 {int degbound = lpDegBound(basering);int n = lp BlockSize(basering);1735 {int degbound = lpDegBound(basering);int n = lpVarBlockSize(basering); 1736 1736 if (size(#) > 0){if (typeof(#[1])=="int"){if (#[1] >= 0){degbound = #[1];}}} 1737 1737 if (size(#) > 1){if (typeof(#[1])=="int"){if (#[2] > 0){n = #[2];}}} … … 1778 1778 SEE ALSO: ncHilb_lib 1779 1779 " 1780 {int degbound = lpDegBound(basering);int n = lp BlockSize(basering);1780 {int degbound = lpDegBound(basering);int n = lpVarBlockSize(basering); 1781 1781 if (size(#) > 0){if (typeof(#[1])=="int"){if (#[1] >= 0){degbound = #[1];}}} 1782 1782 if (size(#) > 1){if (typeof(#[1])=="int"){if (#[2] > 0){n = #[2];}}} … … 1809 1809 EXAMPLE: example lpKDimCheck; shows examples 1810 1810 " 1811 {int n = lp BlockSize(basering);1811 {int n = lpVarBlockSize(basering); 1812 1812 list L; 1813 1813 ideal R; … … 1844 1844 EXAMPLE: example lpKDim; shows examples 1845 1845 " 1846 {int degbound = lpDegBound(basering);int n = lp BlockSize(basering);1846 {int degbound = lpDegBound(basering);int n = lpVarBlockSize(basering); 1847 1847 if (size(#) > 0){if (typeof(#[1])=="int"){if (#[1] >= 0){degbound = #[1];}}} 1848 1848 if (size(#) > 1){if (typeof(#[1])=="int"){if (#[2] > 0){n = #[2];}}} … … 1948 1948 EXAMPLE: example lpSickle; shows examples 1949 1949 " 1950 {int degbound = lpDegBound(basering); int n = lp BlockSize(basering);1950 {int degbound = lpDegBound(basering); int n = lpVarBlockSize(basering); 1951 1951 if (size(#) > 0){if (typeof(#[1])=="int"){if (#[1] >= 0){degbound = #[1];}}} 1952 1952 if (size(#) > 1){if (typeof(#[1])=="int"){if (#[2] > 0){n = #[2];}}} … … 1988 1988 EXAMPLE: example lpSickleDim; shows examples 1989 1989 " 1990 {int degbound = lpDegBound(basering);int n = lp BlockSize(basering);1990 {int degbound = lpDegBound(basering);int n = lpVarBlockSize(basering); 1991 1991 if (size(#) > 0){if (typeof(#[1])=="int"){if (#[1] >= 0){degbound = #[1];}}} 1992 1992 if (size(#) > 1){if (typeof(#[1])=="int"){if (#[2] > 0){n = #[2];}}} … … 2028 2028 EXAMPLE: example lpSickleHil; shows examples 2029 2029 " 2030 {int degbound = lpDegBound(basering);int n = lp BlockSize(basering);2030 {int degbound = lpDegBound(basering);int n = lpVarBlockSize(basering); 2031 2031 if (size(#) > 0){if (typeof(#[1])=="int"){if (#[1] >= 0){degbound = #[1];}}} 2032 2032 if (size(#) > 1){if (typeof(#[1])=="int"){if (#[2] > 0){n = #[2];}}} … … 2082 2082 if (m == 1) 2083 2083 {if (d == 0) 2084 {if (h == 0) {return(lpSickle(G,degbound,lp BlockSize(basering)));}2085 else {return(lpSickleHil(G,degbound,lp BlockSize(basering)));}2084 {if (h == 0) {return(lpSickle(G,degbound,lpVarBlockSize(basering)));} 2085 else {return(lpSickleHil(G,degbound,lpVarBlockSize(basering)));} 2086 2086 } 2087 2087 else 2088 {if (h == 0) {return(lpSickleDim(G,degbound,lp BlockSize(basering)));}2089 else {return(lpDHilbertSickle(G,degbound,lp BlockSize(basering)));}2088 {if (h == 0) {return(lpSickleDim(G,degbound,lpVarBlockSize(basering)));} 2089 else {return(lpDHilbertSickle(G,degbound,lpVarBlockSize(basering)));} 2090 2090 } 2091 2091 } … … 2093 2093 {if (d == 0) 2094 2094 {if (h == 0) {ERROR("You request to do nothing, so relax and do so");} 2095 else {return(lpHilbert(G,degbound,lp BlockSize(basering)));}2095 else {return(lpHilbert(G,degbound,lpVarBlockSize(basering)));} 2096 2096 } 2097 2097 else 2098 {if (h == 0) {return(lpKDim(G,degbound,lp BlockSize(basering)));}2099 else {return(lpDHilbert(G,degbound,lp BlockSize(basering)));}2098 {if (h == 0) {return(lpKDim(G,degbound,lpVarBlockSize(basering)));} 2099 else {return(lpDHilbert(G,degbound,lpVarBlockSize(basering)));} 2100 2100 } 2101 2101 } -
Singular/LIB/fpaprops.lib
rb849856 r921f32 91 91 // if longest word has length 1 we handle it as a special case 92 92 if (l == 1) { 93 int n = lp BlockSize(basering); // variable count93 int n = lpVarBlockSize(basering); // variable count 94 94 int k = size(G); 95 95 if (k == n) { // only the field left … … 746 746 ERROR("Ufnarovskij graph not implemented for l = 1"); 747 747 } 748 int lV = lp BlockSize(basering);748 int lV = lpVarBlockSize(basering); 749 749 // TODO: what if l <= 0? 750 750 dbprint("computing standard words"); … … 836 836 return (words); // no standard words 837 837 } 838 int lV = lp BlockSize(basering); // variable count838 int lV = lpVarBlockSize(basering); // variable count 839 839 list prevWords = ivStandardWords(G, length - 1); 840 840 list words; … … 907 907 // if longest word has length 1, or G is the zero ideal, we handle it as a special case 908 908 if (l == 1 || size(G) == 0) { 909 int n = lp BlockSize(basering); // variable count909 int n = lpVarBlockSize(basering); // variable count 910 910 int k = size(G); 911 911 if (k == n) { // V = {1} no edges … … 1004 1004 if (size(G1) > 0) { 1005 1005 while (size(G1) > 0) { 1006 if (lp BlockSize(R) > 1) {1006 if (lpVarBlockSize(R) > 1) { 1007 1007 def @R = R - string(G1[1]); 1008 1008 R = @R; … … 1082 1082 { 1083 1083 list LG = lpId2ivLi(lead(G)); 1084 int n = lp BlockSize(basering);1084 int n = lpVarBlockSize(basering); 1085 1085 int degbound = lpDegBound(basering); 1086 1086 … … 1384 1384 /* // removes a variable from a letterplace ring (a bit of a hack) */ 1385 1385 /* static proc lpDelVar(int index) { */ 1386 /* int lV = lp BlockSize(basering); // number of variables in the main block */1386 /* int lV = lpVarBlockSize(basering); // number of variables in the main block */ 1387 1387 /* int d = lpDegBound(basering); // degree bround */ 1388 1388 /* list LR = ringlist(basering); */ -
Singular/LIB/freegb.lib
rb849856 r921f32 26 26 isFreeAlgebra(r); check whether r is a letterplace ring (free algebra) 27 27 lpDegBound(R); returns the degree bound of a letterplace ring 28 lp BlockSize(R);returns the size of the letterplace blocks28 lpVarBlockSize(R); returns the size of the letterplace blocks 29 29 30 30 letplaceGBasis(I); (deprecated, use twostd) two-sided Groebner basis of a letterplace ideal I … … 126 126 ring r = 0,(x(1),y(1),x(2),y(2),x(3),y(3),x(4),y(4)),dp; 127 127 def R = setLetterplaceAttributes(r, 4, 2); setring R; 128 lp BlockSize(R);128 lpVarBlockSize(R); 129 129 lieBracket(x(1),y(1),2); 130 130 } … … 1018 1018 } 1019 1019 1020 proc lp BlockSize(def R)1021 "USAGE: lp BlockSize(R); R a letterplace ring1020 proc lpVarBlockSize(def R) 1021 "USAGE: lpVarBlockSize(R); R a letterplace ring 1022 1022 RETURN: int 1023 PURPOSE: returns the block size of the letterplace ring, that is the number of variables of the original ring.1024 EXAMPLE: example lp BlockSize; shows examples1023 PURPOSE: returns the variable block size of the letterplace ring, that is the number of variables of the original ring. 1024 EXAMPLE: example lpVarBlockSize; shows examples 1025 1025 " 1026 1026 { … … 1036 1036 ring r = 0,(x,y,z),dp; 1037 1037 def R = freeAlgebra(r, 7); 1038 lp BlockSize(R);1038 lpVarBlockSize(R); 1039 1039 } 1040 1040 … … 1128 1128 def A = makeLetterplaceRing(2); // same as makeLetterplaceRing(2,0) 1129 1129 setring A; A; 1130 lp BlockSize(A);1130 lpVarBlockSize(A); 1131 1131 lpDegBound(A); // degree bound 1132 1132 setring r; def B = makeLetterplaceRing(2,1); // to compare: 1133 1133 setring B; B; 1134 lp BlockSize(B);1134 lpVarBlockSize(B); 1135 1135 lpDegBound(B); // degree bound 1136 1136 setring r; def C = makeLetterplaceRing(2,2); // to compare: … … 1219 1219 setring A; 1220 1220 A; 1221 lp BlockSize(A);// number of variables in the main block1221 lpVarBlockSize(A);// number of variables in the main block 1222 1222 lpDegBound(A); // degree bound 1223 1223 } … … 1324 1324 setring A; 1325 1325 A; 1326 lp BlockSize(A); // number of variables in the main block1326 lpVarBlockSize(A); // number of variables in the main block 1327 1327 lpDegBound(A); // degree bound 1328 1328 } … … 1388 1388 setring A; 1389 1389 A; 1390 lp BlockSize(A); // number of variables in the main block1390 lpVarBlockSize(A); // number of variables in the main block 1391 1391 lpDegBound(A); // degree bound 1392 1392 } … … 1494 1494 setring A; 1495 1495 A; 1496 lp BlockSize(A); // number of variables in the main block1496 lpVarBlockSize(A); // number of variables in the main block 1497 1497 lpDegBound(A); // degree bound 1498 1498 } … … 2941 2941 // make new ring 2942 2942 def save = basering; 2943 int norigvars = lp BlockSize(save);2943 int norigvars = lpVarBlockSize(save); 2944 2944 def Rtagged; def temp = save; 2945 2945 for (int i = 1; i <= size(I); i++) { … … 2948 2948 } kill i; 2949 2949 // currently R + "var" doesn't preserve uptodeg 2950 Rtagged = setLetterplaceAttributes(Rtagged, lp BlockSize(Rtagged), lpDegBound(save));2950 Rtagged = setLetterplaceAttributes(Rtagged, lpVarBlockSize(Rtagged), lpDegBound(save)); 2951 2951 setring Rtagged; 2952 2952 … … 3060 3060 {if (I==intvec(0)) {return(intvec(0));} 3061 3061 int j,k,l; 3062 int n = lp BlockSize(basering); int d = lpDegBound(basering);3062 int n = lpVarBlockSize(basering); int d = lpDegBound(basering); 3063 3063 intvec w; j = 1; 3064 3064 while (j <= d) … … 3096 3096 " 3097 3097 {int i,j,k,r1,r2; intvec D; 3098 int n = lp BlockSize(basering);3098 int n = lpVarBlockSize(basering); 3099 3099 k = size(V) div n; r1 = 0; r2 = 0; 3100 3100 for (i=1; i<= k; i++) … … 3141 3141 3142 3142 int i,j,r; intvec T; list R; 3143 int n = lp BlockSize(basering);3143 int n = lpVarBlockSize(basering); 3144 3144 int k = size(V) - size(W) + 1; 3145 3145 if (intvec(V[1..size(W)])-W == 0){R[1]=0;} … … 3309 3309 } 3310 3310 3311 int lV = lp BlockSize(basering); // variable count3311 int lV = lpVarBlockSize(basering); // variable count 3312 3312 ideal prevMonomials = lpMonomialsWithHoles(d - 1); 3313 3313 … … 3335 3335 static proc getlpCoeffs(poly q, poly p) 3336 3336 {list R; intvec cq,t,lv,rv,bla; 3337 int n = lp BlockSize(basering); int d = lpDegBound(basering);3337 int n = lpVarBlockSize(basering); int d = lpDegBound(basering); 3338 3338 int i; 3339 3339 cq = leadexp(p)-leadexp(q); /* p/q */ … … 3383 3383 int i; 3384 3384 for (i = 1; i <= size(L); i++) 3385 {if (entryViolation(L[i], lp BlockSize(basering)))3385 {if (entryViolation(L[i], lpVarBlockSize(basering))) 3386 3386 {ERROR("Not allowed monomial/intvec found!");} 3387 3387 } … … 3427 3427 {int r = 0; intvec w; 3428 3428 intvec l = leadexp(p); 3429 int n = lp BlockSize(basering); int d = lpDegBound(basering);3429 int n = lpVarBlockSize(basering); int d = lpDegBound(basering); 3430 3430 int i,j,c,c1; 3431 3431 while (1 <= d) … … 3822 3822 if (p == 0) {ERROR("Monomial is not allowed to equal zero");} 3823 3823 intvec lep = leadexp(p); 3824 for ( i = 1; i <= lp BlockSize(basering); i++) {if (lep[i] == 1) {I = i; break;}}3825 for (i = (lp BlockSize(basering)+1); i <= size(lep); i++)3824 for ( i = 1; i <= lpVarBlockSize(basering); i++) {if (lep[i] == 1) {I = i; break;}} 3825 for (i = (lpVarBlockSize(basering)+1); i <= size(lep); i++) 3826 3826 {if (lep[i] == 1) 3827 { j = (i mod lp BlockSize(basering));3828 if (j == 0) {I = I,lp BlockSize(basering);}3827 { j = (i mod lpVarBlockSize(basering)); 3828 if (j == 0) {I = I,lpVarBlockSize(basering);} 3829 3829 else {I = I,j;} 3830 3830 } -
Singular/LIB/ncHilb.lib
rb849856 r921f32 75 75 def save = basering; 76 76 int sz=size(#); 77 int lV=lp BlockSize(save); // nvars(orig comm ring)77 int lV=lpVarBlockSize(save); // nvars(orig comm ring) 78 78 int ig=0; 79 79 int mgrad=0; … … 403 403 EXAMPLE: example rcolon; shows an example" 404 404 { 405 int lV = lp BlockSize(R); //nvars(save);405 int lV = lpVarBlockSize(R); //nvars(save); 406 406 if (lV == 0) 407 407 {
Note: See TracChangeset
for help on using the changeset viewer.