Changeset 360741 in git for Singular/LIB/zeroset.lib
- Timestamp:
- Jan 21, 2020, 4:38:40 PM (4 years ago)
- Branches:
- (u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', 'c5facdfddea2addfd91babd8b9019161dea4b695')
- Children:
- f47408b1fbaa5550e599a333fbef837cbffa286d
- Parents:
- b262f1218ae1112aa4dde91dcd16243c0b6acfeecbf01b413214dc520717e4af698cf40a2a2954e4
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/zeroset.lib
rb262f12 r360741 421 421 int n = nvars(basering); 422 422 def RINB = basering; 423 string ringSTR = "ring RINR = 0, " + string(var(n)) + ", dp;"; 424 execute(ringSTR); // new ring = Q[a] 423 ring RINR = create_ring(0, string(var(n)), "dp"); // new ring = Q[a] 425 424 426 425 list gcdList; … … 637 636 if(npars(basering) == 0) { @parName = "0";} 638 637 else { @parName = "(0, " + parstr(basering) + ")"; } 639 @str = "ring @RGCD = " + @parName + ", " + string(var(varIndex)) + ", dp;"; 640 execute(@str); 638 ring @RGCD = create_ring(@parName, string(var(varIndex)), "dp"); 641 639 if(@mp != "0") { execute ("minpoly = " + @mp + ";"); } 642 640 execute("poly @f = " + @sf + ";"); … … 717 715 g = f; 718 716 719 ringSTR = "ring SNRM1 = 0, " + string(var(1)) + ", dp;"; // univariate ring 720 execute(ringSTR); 717 ring SNRM1 = create_ring(0, string(var(1)), "dp"); // univariate ring 721 718 poly N, N1, h; // N, N1 do not contain 'a', use built-in gcd 722 719 h = gcd(imap(SNRMB, N), imap(SNRMB, N1)); … … 1100 1097 // create univariate ring 1101 1098 1102 ringSTR = "ring RIS1 = 0, (" + string(var(vars[1])) + ", " + string(var(n+1)) + "), lp;"; 1103 execute(ringSTR); 1099 ring RIS1 = create_ring(0, "(" + string(var(vars[1])) + ", " + string(var(n+1)) + ")", "lp"); 1104 1100 ideal mpoly = std(imap(NLZR, mpoly)); 1105 1101 list roots; … … 1405 1401 minPoly = string(minpoly); 1406 1402 } 1407 ringSTR = "ring TR = 0, (" + varstr(basering) + "," + parName + "), lp;"; 1408 1409 execute(ringSTR); 1403 ring TR = create_ring(0, "(" + varstr(basering) + "," + parName + ")", "lp"); 1410 1404 execute("ideal mpoly = std(" + minPoly + ");"); 1411 1405 export(mpoly); … … 1444 1438 varString = varString + string(var(n-1)); 1445 1439 1446 ringSTR = "ring TR = " + parName + ", (" + varString + "), lp;"; 1447 execute(ringSTR); 1440 ring TR = create_ring(parName, "(" + varString + ")", "lp"); 1448 1441 if(MP) { minpoly = number(imap(BR, mpoly)[1]); } 1449 1442 setring BR;
Note: See TracChangeset
for help on using the changeset viewer.