1% wu.tst
2
3% Russell Bradford, 8 June 90.
4
5% Some tests for the Wu algorithm
6% The order directives are not necessary for general use: they just
7% help tie things down for testing purposes.
8
9% run after loading Wu code: in "wu.red"$
10
11% test 1
12
13order x,y,a,b;
14
15wu({x^2+y^2-a,x*y-b}, {x,y});
16
17% test 2
18
19order x,y,a,b;
20
21wu({x^2+y^2-a,x*y-b},{x,y,a,b});
22
23% test 3
24
25order x,y,z,r;
26
27wu({x^2+y^2+z^2-r^2, x*y+z^2-1, x*y*z-x^2-y^2-z+1}, {x,y,z});
28
29% test 4
30
31order x,y,z,r;
32
33wu({x^2+y^2+z^2-r^2, x*y+z^2-1, x*y*z-x^2-y^2-z+1}, {x,y,z,r});
34
35% test 5
36
37order x,y,z;
38
39wu({(x-1)*(y-1)*(z-1), (x-2)*(y-2)*(z-2), (x-3)*(y-3)*(z-3)}, {x,y,z});
40
41% test 6
42
43order x,y,z;
44
45wu({(x-1)*(y-1)*(z-1), (x-2)*(y-2)*(z-2), (x-3)*(y-3)*(z-3)});
46
47% test 7
48
49order x1,x2,x3,x4;
50
51p1 := x1+x2+x3+x4;
52p2 := x1*x2+x2*x3+x3*x4+x4*x1;
53p3 := x1*x2*x3+x2*x3*x4+x3*x4*x1+x4*x1*x2;
54p4 := x1*x2*x3*x4 - 1;
55
56wu({p1,p2,p3,p4}, {x1,x2,x3,x4});
57
58% test 8
59
60order x,y,z;
61
62wu({z*z,y*z-1,x*z-1}, {x,y,z});
63
64end;
65