1 | LIB "tst.lib"; |
---|
2 | tst_init(); |
---|
3 | |
---|
4 | ///////////////////////////////////////////////////////////////////////////// |
---|
5 | /// Further examples for testing the main procedures |
---|
6 | /// Timings on wawa Sept 29 |
---|
7 | ///////////////////////////////////////////////////////////////////////////// |
---|
8 | LIB"paraplanecurves.lib"; |
---|
9 | // ------------------------------------------------------- |
---|
10 | // Example 1 |
---|
11 | // ------------------------------------------------------- |
---|
12 | ring RR = 0, (x,y,z), dp; |
---|
13 | poly f = 7*z2+11*y2+13*z*y+17*x2+19*x*y; // conic |
---|
14 | def RP1 = paraConic(f); |
---|
15 | setring RP1; PARACONIC; |
---|
16 | setring RR; |
---|
17 | RP1 = paraPlaneCurve(f); |
---|
18 | testParametrization(f,RP1); |
---|
19 | setring RP1; PARA; |
---|
20 | kill RR;kill RP1; |
---|
21 | // ------------------------------------------------------- |
---|
22 | // Example 21 |
---|
23 | // ------------------------------------------------------- |
---|
24 | ring RR = 0, (x,y,z), dp; |
---|
25 | poly f = -1965466244509920x5y+34871245546721380061760x4y2; |
---|
26 | f = f+104613747941595046117320x3y3+113331564241941002407560x2y4; |
---|
27 | f = f+52306876673313609259800xy5+8717812860780028397880y6; |
---|
28 | f = f+1040297748510024x5z+4468147845634872x4yz; |
---|
29 | f = f-22398508728211453743258x3y2z-33223996581074443306854x2y3z; |
---|
30 | f = f-10638598235041298082366xy4z+186886189971594356382y5z; |
---|
31 | f = f-1385078844909312x4z2-34893092731637052532683x3yz2; |
---|
32 | f = f-98591463214095439056609x2y2z2-92339459334829609336485xy3z2; |
---|
33 | f = f-24923289542522905755711y4z2+472440640471377x3z3; |
---|
34 | f = f+33821511925664516716011x2yz3+49745237303968344397437xy2z3; |
---|
35 | f = f+11040465960074786720475y3z3+8728735735878837099404x2z4; |
---|
36 | f = f+17676785754519678518537xyz4+17935885079051421934609y2z4; |
---|
37 | f = f-11314701999743172607075xz5-16164284825803158969425yz5; |
---|
38 | f = f+3666695988537425618750z6; |
---|
39 | // 4 nodes, 1 OMP of mult 4 |
---|
40 | adjointIdeal(f,2); |
---|
41 | kill RR; |
---|
42 | |
---|
43 | tst_status(1);$ |
---|
44 | |
---|