fieker-DuValspielwiese
|
Last change
on this file since fab637 was
0abc87,
checked in by Martin Lamm <lamm@…>, 25 years ago
|
|
* added check for proc primitive_extra
* deleted check for proc randomLast (no longer part of primitiv.lib)
git-svn-id: file:///usr/local/Singular/svn/trunk@3428 2c84dea3-7e68-4137-9b89-c4e89433aadc
|
-
Property mode set to
100644
|
|
File size:
1.2 KB
|
| Line | |
|---|
| 1 | // File: primitiv_s.tst |
|---|
| 2 | // Short tests for primitiv lib |
|---|
| 3 | LIB "tst.lib"; |
|---|
| 4 | tst_init(); |
|---|
| 5 | LIB "primitiv.lib"; |
|---|
| 6 | // ------------ test of primitive: ------------------- |
|---|
| 7 | printlevel=3; |
|---|
| 8 | ring r=0,x,dp; |
|---|
| 9 | primitive(ideal(x2+1)); |
|---|
| 10 | kill r; |
|---|
| 11 | ring r=0,(x,y),dp; |
|---|
| 12 | ideal i=x2-2,y2-3; |
|---|
| 13 | ideal j=primitive(i); |
|---|
| 14 | factorize(j[1]); |
|---|
| 15 | ring r1=(0,a),x,lp; |
|---|
| 16 | map phi=r,x,a; |
|---|
| 17 | string s=string(phi(j)[1]); |
|---|
| 18 | execute("minpoly="+s+";"); |
|---|
| 19 | ideal v=phi(j)[1],phi(j)[2]^2-2,phi(j)[3]^2-3; |
|---|
| 20 | v; |
|---|
| 21 | if (size(v)>0) |
|---|
| 22 | { "ERROR!!! Result of primitive(i) is not a primitive element!"; } |
|---|
| 23 | // else minpoly(a)==0, g_1(a) is squareroot of 2, g_2(a) is squareroot of 3 |
|---|
| 24 | // as it should be |
|---|
| 25 | setring r; |
|---|
| 26 | i=x2-3,y2-x; |
|---|
| 27 | primitive(i); |
|---|
| 28 | kill r; |
|---|
| 29 | ring r=3,(x,y,z),dp; |
|---|
| 30 | ideal i=x2+1,y3-y-1,z2+yz-1; |
|---|
| 31 | primitive(i); |
|---|
| 32 | // ------------ test of primitive_extra: ------------- |
|---|
| 33 | primitive_extra(i); |
|---|
| 34 | primitive_extra(ideal(i[1],i[2])); |
|---|
| 35 | kill r; |
|---|
| 36 | // ------------ test of splitring: ------------------- |
|---|
| 37 | ring r=0,(x,y),ds; |
|---|
| 38 | splitring(x2-5,"R1"); |
|---|
| 39 | R1; |
|---|
| 40 | list l=splitring(x2-7,"R2",a); |
|---|
| 41 | l[1]^2; |
|---|
| 42 | R2; |
|---|
| 43 | splitring(y2+1,"",a); |
|---|
| 44 | R2; |
|---|
| 45 | kill r,R1,R2; |
|---|
| 46 | ring r=(2003,j),(a,b,c,d),dp; |
|---|
| 47 | minpoly=j2-2; |
|---|
| 48 | list l=splitring(b3+b2+b-2,"R1",list(a,b,c,d,j)); |
|---|
| 49 | l; |
|---|
| 50 | number(l[5])^2; |
|---|
| 51 | factorize(b3+b2+b-2); |
|---|
| 52 | list L=splitring(c2-2j,"",list(a2+ja,j,l[5])); |
|---|
| 53 | L; |
|---|
| 54 | L[3]^2; |
|---|
| 55 | kill R1,r; |
|---|
| 56 | tst_status(1);$ |
|---|
| 57 | |
|---|
Note: See
TracBrowser
for help on using the repository browser.
