Changeset 3d5b277 in git for Singular/LIB/numerDecom.lib
- Timestamp:
- Jul 25, 2019, 1:23:11 PM (5 years ago)
- Branches:
- (u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
- Children:
- 35a153cf7acb82519d52d0b74581e9a193d15765
- Parents:
- c9824e54c705dcaedb3b991c6bf81f9f974e758d3f94cbdb24129f0002a999ab5fb63a34f0fc933d
- git-author:
- Hans Schoenemann <hannes@mathematik.uni-kl.de>2019-07-25 13:23:11+02:00
- git-committer:
- GitHub <noreply@github.com>2019-07-25 13:23:11+02:00
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/numerDecom.lib
r3f94cb r3d5b277 349 349 if((n-rn)!=dd) 350 350 { 351 execute("ring R=0,("+varstr(S)+",z(1..dd)),dp;"); 351 list l2 = ringlist(S)[2]; 352 for (int ii = 1; ii <= dd; ii++) 353 { 354 l2[size(l2)+1] = "z("+string(ii)+")"; 355 } 356 ring R = create_ring(0, l2, "dp"); 352 357 ideal I=imap(rs,I); 353 358 ideal H(0..n),L,LL,L(1..dd),LL(1..dd),h(1..dd),N(0..dd); … … 483 488 { 484 489 int w(q-1)=0; 485 execute("ring D(q)=(0,s,gamma),("+varstr(S)+",z(1..q)),dp;"); 490 list l2 = ringlist(S)[2]; 491 for (int ii = 1; ii <= q; ii++) 492 { 493 l2[size(l2)+1] = "z("+string(ii)+")"; 494 } 495 ring D(q) = create_ring("(0,s,gamma)", l2, "dp"); 486 496 string nonsin(q),stnonsin(q); 487 497 ideal H(1..q); … … 579 589 for(qq=q-1;qq>=1;qq--) 580 590 { 581 execute("ring T(qq)=(complex,16,I),("+varstr(S)+",z(1..qq)),dp;"); 591 list l2 = ringlist(S)[2]; 592 for (int ii = 1; ii <= qq; ii++) 593 { 594 l2[size(l2)+1] = "z("+string(ii)+")"; 595 } 596 ring T(qq) = create_ring("(complex,16,I)", l2, "dp"); 582 597 list W(qq-1)=var(1); 583 598 } … … 591 606 { 592 607 int w(qq-1); 593 execute("ring T(qq)=(complex,16,I),("+varstr(S)+",z(1..qq)),dp;"); 608 list l2 = ringlist(S)[2]; 609 for (int ii = 1; ii <= qq; ii++) 610 { 611 l2[size(l2)+1] = "z("+string(ii)+")"; 612 } 613 ring T(qq) = create_ring("(complex,16,I)", l2, "dp"); 594 614 list W(qq-1)=var(1); 595 615 } … … 600 620 for(qq=q;qq>=1;qq--) 601 621 { 602 execute("ring T(qq)=(complex,16,I),("+varstr(S)+",z(1..qq)),dp;"); 622 list l2 = ringlist(S)[2]; 623 for (int ii = 1; ii <= qq; ii++) 624 { 625 l2[size(l2)+1] = "z("+string(ii)+")"; 626 } 627 ring T(qq) = create_ring("(complex,16,I)", l2, "dp"); 603 628 list W(qq-1)=var(1); 604 629 } … … 1951 1976 export(J(i)); 1952 1977 } 1953 execute("ring RR=0,(x(1..n),"+varstr(S)+"),dp;"); 1978 list l2; 1979 for (int ii = 1; ii <= n; ii++) 1980 { 1981 l2[ii] = "x("+string(ii)+")"; 1982 } 1983 l2 = l2+ringlist(S)[2]; 1984 ring RR = create_ring(0, l2, "dp"); 1954 1985 for(i=1;i<=d;i++) 1955 1986 { … … 1995 2026 int zc=size(D(d)); 1996 2027 export(zc); 1997 execute("ring DR=0,("+varstr(S)+",x(1..zc)),dp;"); 2028 list l2 = ringlist(S)[2]; 2029 for (int ii = 1; ii <= zc; ii++) 2030 { 2031 l2[size(l2)+1] = "x("+string(ii)+")"; 2032 } 2033 ring DR = create_ring(0, l2, "dp"); 1998 2034 matrix TT(d)=imap(R,TT(d)); 1999 2035 ideal I=imap(S,I);
Note: See TracChangeset
for help on using the changeset viewer.