Changeset b45eec in git for Singular/LIB/graal.lib
- Timestamp:
- Apr 19, 2020, 9:57:14 PM (4 years ago)
- Branches:
- (u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
- Children:
- ffcf715a570941627bf020648e93955f1596f2ef
- Parents:
- d49a473d44aa25c5ecc3126360ab3179ddb8af7e636fa5455e699fa8c0a820dd60e00a059f96f8ef
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/graal.lib
rd49a47 rb45eec 1 ////////////////////////////////////////////////////////////////////////////1 graalMixed//////////////////////////////////////////////////////////////////////////// 2 2 version="version graal.lib 4.1.2.0 Feb_2019 "; // $Id$ 3 3 category="Commutative Algebra"; … … 332 332 **/ 333 333 ideal H = ringlist(Gr.A)[4]; 334 execute("ring Q = "+string(Gr.A)+";"); 334 list l3 = ringlist(Gr.A); 335 def Q = ring(l3); 336 setring Q; 335 337 ideal H = fetch(Gr.A,H); H = sort(std(H))[1]; 336 338 ideal J = fetch(Gr.A,L) + H; J = sort(std(J))[1]; … … 358 360 else 359 361 { string pars; string vars = varstr(basering); } 360 execute("ring Q0 = ("+charstr(basering)+pars+"),("+vars+"),dp;");362 ring Q0 = create_ring("("+charstr(basering)+pars+")", "("+vars+")", "dp"); 361 363 ideal J0 = imap(Q,J); 362 364 Gr.Q0 = Q0; … … 404 406 string ostring = ordstr(basering); 405 407 ostring = ostring[1..size(ostring)-2]; 406 list l2; 407 for (int ii = 1; ii <= s; ii++) 408 { 409 l2[ii] = "Y("+string(ii)+")"; 410 } 411 l2[s+1] = ringlist(basering)[2]; 412 execute("ring Q0y = ("+charstr(basering)+"),(Y(1..s),"+varstr(basering)+"),(ds(s),c,"+ostring+")"); 408 list l4; 409 for (int zz = 1; zz <= s; zz++) 410 { 411 l4[zz] = "Y("+string(zz)+")"; 412 } 413 ring Q0y = create_ring(ringlist(basering)[1], "("+string(l4)+","+varstr(basering)+")", "(ds("+string(s)+"),c,"+ostring+")", "no_minpoly"); 413 414 setring Q0y; 414 415 ideal H = imap(Gr.A,H); … … 428 429 * construct Ky and sigmainI 429 430 **/ 430 execute("ring KKy = ("+charstr(Q0)+"),(Y(1..s),"+varstr(Q0,nvars(Q0))+"),(c,dp(s),dp(1));");431 ring KKy = create_ring(ringlist(Q0)[1], "("+string(l4)+","+varstr(Q0,nvars(Q0))+")", "(c,dp("+string(s)+"),dp(1))", "no_minpoly"); 431 432 poly minpolyOfK = imap(Q0,g(n)); 432 433 qring Ky = std(minpolyOfK);
Note: See TracChangeset
for help on using the changeset viewer.