1 | // File: finvar2_l.tst |
---|
2 | // Long tests for finvar.lib |
---|
3 | |
---|
4 | LIB "tst.lib"; |
---|
5 | tst_init(); |
---|
6 | LIB "finvar.lib"; |
---|
7 | option(redSB); |
---|
8 | |
---|
9 | |
---|
10 | ring R2=29,(x,y,z),(L(255),lp); |
---|
11 | number e=7; |
---|
12 | number i=12; |
---|
13 | number s=6; |
---|
14 | number a=(e^5-e^2)/(i*s); |
---|
15 | number b=(e^3-e^4)/(i*s); |
---|
16 | number c=(e^6-e)/(i*s); |
---|
17 | matrix A[3][3]=a,b,c,b,c,a,c,a,b; |
---|
18 | matrix B[3][3]=e,0,0,0,e^4,0,0,0,e^2; |
---|
19 | list L=primary_invariants(A,B,intvec(1,0,1)); |
---|
20 | // the following is not unique! |
---|
21 | L; |
---|
22 | tst_status(); |
---|
23 | matrix S, IS1 = secondary_and_irreducibles_no_molien(L[1..size(L)],1); |
---|
24 | // the following is not unique! |
---|
25 | S; |
---|
26 | IS1; |
---|
27 | tst_status(); |
---|
28 | kill S; |
---|
29 | matrix IS2=irred_secondary_no_molien(L[1..size(L)],1); |
---|
30 | // the following is not unique! |
---|
31 | IS2; |
---|
32 | tst_status(); |
---|
33 | list L2 = invariant_ring(A,B,intvec(0,0,1)); |
---|
34 | // the following is not unique! |
---|
35 | L2; |
---|
36 | tst_status(); |
---|
37 | list L3 = invariant_ring(A,B,intvec(1,0,1)); |
---|
38 | // the following is not unique! |
---|
39 | L3; |
---|
40 | tst_status(); |
---|
41 | matrix TST2a = groebner(ideal(L[1])+ideal(IS1)); |
---|
42 | matrix TST2b = groebner(ideal(L[1])+ideal(IS2)); |
---|
43 | kill L,IS1,IS2; |
---|
44 | matrix TST2c = groebner(ideal(L2[1])+ideal(L2[3])); |
---|
45 | matrix TST2d = groebner(ideal(L3[1])+ideal(L3[3])); |
---|
46 | // the following is unique |
---|
47 | TST2a; |
---|
48 | TST2a==TST2b; |
---|
49 | TST2a==TST2c; |
---|
50 | TST2a==TST2d; |
---|
51 | tst_status(); |
---|
52 | kill R2; |
---|
53 | |
---|
54 | tst_status(1);$ |
---|