source: git/dyn_modules/syzextra/p_content.tst @ f115489

fieker-DuValspielwiese
Last change on this file since f115489 was f115489, checked in by Oleksandr Motsak <motsak@…>, 13 years ago
ADD: bugtest for ntIntDiv & p_Content using syzextra From: Oleksandr Motsak <motsak@mathematik.uni-kl.de> git-svn-id: file:///usr/local/Singular/svn/trunk@14091 2c84dea3-7e68-4137-9b89-c4e89433aadc
  • Property mode set to 100644
File size: 2.2 KB
Line 
1echo = 2;
2
3"ndebug?: ", system("with", "ndebug");
4"om_ndebug?: ", system("with", "om_ndebug");
5
6if( system("with", "ndebug") )
7{
8  "Loading the Release version!";
9  LIB("./syzextra.so");
10}
11else
12{
13  "Loading the Debug version!";
14  LIB("./syzextra_g.sog");
15}
16
17noop();
18
19
20option(redSB);option(redTail);
21
22ring r=(0,q),(U,V,S,T),lp;
23
24poly o = (-q4-5q3-2q2+10q-4)/(q)*ST2+(-q4-6q3-5q2+12q-4)*ST+(q5+5q4+2q3-10q2+4q)*T2+(q5+6q4+5q3-12q2+4q)*T;
25
26poly O = o;
27
28DetailedPrint(o);
29
30p_Content( O );
31
32O;
33
34simplify(O, 1 + 2 + 4 + 8);
35
36leadcoef(o) / leadcoef(O);
37
38o = o - lead(o);
39O = O - lead(O);
40
41leadcoef(o) / leadcoef(O);
42
43
44o = o - lead(o);
45O = O - lead(O);
46
47leadcoef(o) / leadcoef(O);
48
49o = o - lead(o);
50O = O - lead(O);
51
52leadcoef(o) / leadcoef(O);
53
54o;
55O;
56
57$$$
58
59
60/*
61matrix C[4][4];
62C[1,2]=1/q;
63C[1,3]=1;
64C[1,4]=1;
65C[2,3]=1;
66C[2,4]=1;
67C[3,4]=1/q;
68def R=nc_algebra(C,0);
69setring R;
70R;
71
72DetailedPrint(R);
73
74*/
75
76
77// characteristic : 0
78// 1 parameter : q
79// minpoly : 0
80// number of vars : 4
81// block 1 : ordering lp
82// : names U V S T
83// block 2 : ordering C
84// noncommutative relations:
85// VU=1/(q)*UV
86// TS=1/(q)*ST
87ideal I=U*V+S-V-1,U-V-V*T,U+V-V*S;
88
89I;
90// I[1]=UV-V+S-1
91// I[2]=U-VT-V
92// I[3]=U-VS+V
93
94slimgb(I);// correct
95
96
97option(mem);
98option(prot);
99option(teach);
100
101
102
103$$
104
105
106LIB "teachstd.lib";
107printlevel = 0;
108standard(I);
109
110simplify(_, 1+2+4);
111 
112ideal II = _;
113slimgb(II);
114
115
116interred(II);
117
118
119
120$
121
122
123
124// std(ideal(II[7], II[1]));
125
126std(ideal(II[7], II[3], II[9]));
127
128$$$
129
130
131std(I);
132// G[1]=T
133// G[2]=S-1
134// G[3]=V
135// G[4]=U
136
137
138
139matrix M;
140ideal H=liftstd(I,M);
141
142// _[1]=T
143// _[2]=S-1
144// _[3]=V
145// _[4]=U
146H;
147// H[1]=(q7-q6-3q5+5q4-2q3)*S+(-q7+q6+3q5-5q4+2q3)
148// H[2]=(-q11+2q10+2q9-8q8+7q7-2q6)*V
149// H[3]=(q6-3q4+2q3)*U
150std(H);
151// _[1]=S-1
152// _[2]=V
153// _[3]=U
154
155M;
156print(M);
157
158
159
160
161kill r, R;
162
163$$$$
164
165
166ring r=(0,Q),(U,V,S,T),lp;
167matrix C[4][4];
168C[1,2]=Q;
169C[1,3]=1;
170C[1,4]=1;
171C[2,3]=1;
172C[2,4]=1;
173C[3,4]=Q;
174print(C);
175def R=nc_algebra(C,0);
176setring R;
177R;
178ideal I=U*V+S-V-1,U-V-V*T,U+V-V*S;
179option(redSB);
180option(redTail);
181slimgb(I);// correct
182std(I); // is correct as well!!!
183
184
185// U, V, S-1
186
187
188matrix Z;
189liftstd(I,Z);
190
191
192// _[1]=(2Q11-5Q10+3Q9+Q8-Q7)*S+(-2Q11+5Q10-3Q9-Q8+Q7)
193// _[2]=(-2Q16+7Q15-8Q14+2Q13+2Q12-Q11)*V
194// _[3]=(2Q7-3Q6+Q4)*U
195
196
197simplify(_,1+2+4 );
198
199$$
200
Note: See TracBrowser for help on using the repository browser.