spielwiese
Last change
on this file since 1427f62 was
75f460,
checked in by Hans Schoenemann <hannes@…>, 9 years ago
|
format
|
-
Property mode set to
100644
|
File size:
704 bytes
|
Line | |
---|
1 | LIB "tst.lib"; |
---|
2 | tst_init(); |
---|
3 | |
---|
4 | // incorrect radical entry in a primdecSY |
---|
5 | LIB"primdec.lib"; |
---|
6 | ring rng = (0),(m,p,w,o,j,d,n,z),dp; |
---|
7 | |
---|
8 | ideal I = -2*j+19*o^2, |
---|
9 | 18*w*z+19*d, |
---|
10 | -13*p*d*z-16*o, |
---|
11 | -2*j*n^2+9*n, |
---|
12 | 23*m-8*j, |
---|
13 | -17*m*n+17*m*o, |
---|
14 | -18*n*z+3*j, |
---|
15 | 23*m*p*n-30*w*n, |
---|
16 | -6*m*w, |
---|
17 | 23*j*d*z+4*d*z, |
---|
18 | 16*m*w+14*j, |
---|
19 | -20*p*j*n+27*z^2, |
---|
20 | 23*p*w+14*m*d, |
---|
21 | 7*m*j+29*o*d*n; |
---|
22 | |
---|
23 | list L = primdecSY(I); |
---|
24 | ideal pc2 = L[2][1]; |
---|
25 | pc2; |
---|
26 | ideal rpc2 = L[2][2]; // incorrect!! |
---|
27 | rpc2; |
---|
28 | |
---|
29 | radical(L[2][1]); |
---|
30 | |
---|
31 | idealsEqual( radical(pc2), rpc2 ); // = 0! |
---|
32 | /////////////////////////////////////////// |
---|
33 | // same with groebner(I): |
---|
34 | L = primdecSY (groebner(I) ); |
---|
35 | pc2 = L[2][1]; |
---|
36 | rpc2 = L[2][2]; |
---|
37 | rpc2; |
---|
38 | radical(L[2][1]); |
---|
39 | idealsEqual( radical(pc2), rpc2 ); //ok! |
---|
40 | |
---|
41 | tst_status(1);$ |
---|
Note: See
TracBrowser
for help on using the repository browser.