Changeset ca51f3 in git
- Timestamp:
- Mar 12, 2019, 4:48:12 PM (5 years ago)
- Branches:
- (u'spielwiese', '8e0ad00ce244dfd0756200662572aef8402f13d5')
- Children:
- 78231d7139b2194b384eff921c267d7838ebfba8
- Parents:
- 0561fc58af61d0d8b4a15e3d4ea172191a712c49
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/VecField.lib
r0561fc rca51f3 469 469 470 470 list W; 471 VecField temp; 471 472 for(i=1;i<=size(l);i++) 472 473 { 473 VecFieldtemp = changeCoordinates(l[i],coordChange);474 temp = changeCoordinates(l[i],coordChange); 474 475 W = insert(W,temp,size(W)); 475 476 } … … 821 822 822 823 // find out how to write this in the basis again by brute force: 824 poly mon,lcoef; 823 825 while(p != 0) 824 826 { 825 polymon = leadmonom(p);826 polylcoef = leadcoef(p);827 mon = leadmonom(p); 828 lcoef = leadcoef(p); 827 829 p = p-mon*lcoef; 828 830 for(int h=1;h<=ncols(W);h++) … … 1069 1071 int pos=1; 1070 1072 int i,j,k; 1073 module M; 1074 matrix m; 1071 1075 for(i=1;i<=size(spaces);i++) 1072 1076 { 1073 moduleM = spaces[i];1074 m atrix m= M;1077 M = spaces[i]; 1078 m = M; 1075 1079 for(j=1;j<=ncols(m);j++) 1076 1080 { … … 1220 1224 // we compute this by brute force: for each monomial in H, 1221 1225 // check whether it is wh of lambda,weights. 1226 poly mon; 1227 matrix res; 1222 1228 for(int j=1;j<=ncols(H);j++) 1223 1229 { 1224 polymon = H[j];1230 mon = H[j]; 1225 1231 matrix expos[nvars(basering)][1] = leadexp(mon); 1226 1232 expos = transpose(expos); 1227 matrix res = expos*w; 1233 res = expos*w; 1234 kill expos; 1228 1235 if(res[1,1] == lambda) 1229 1236 { … … 1259 1266 ideal W; 1260 1267 matrix w[nvars(basering)][1] = weights; 1268 poly mon; 1269 matrix res; 1261 1270 for(int j=1;j<=ncols(H);j++) 1262 1271 { 1263 polymon = H[j];1272 mon = H[j]; 1264 1273 matrix expos[nvars(basering)][1] = leadexp(mon); 1265 1274 expos = transpose(expos); 1266 matrix res = expos*w; 1275 res = expos*w; 1276 kill expos; 1267 1277 if(res[1,1] != lambda) // this line is the only difference to weaklyHomogeneousSpace() 1268 1278 { … … 1424 1434 poly c = 0; 1425 1435 int origSize = size(g); 1436 poly mon; 1426 1437 for(int j=1;j <= origSize;j++) 1427 1438 { 1428 polymon = leadmonom(g);1439 mon = leadmonom(g); 1429 1440 mon = mon*leadcoef(g); 1430 1441 g = g-mon;
Note: See TracChangeset
for help on using the changeset viewer.